 # Delivered to the U.S. Government in Accordance with Contract FA8650-10-C-7080.
 #
 # EXPORT CONTROLLED
 # ECCN#   EAR99  
 # Warning:  This document contains technical data which is controlled under the
 # Bureau of Industry & Security’s Export Administration Regulations and may not be
 # exported or re-exported without authorization by the Bureau of Industry & Security.
 # Diversion contrary to US law is prohibited.
 #
 # (c) 2011 United Technologies Corp.
 # Developed with the sponsorship of the Defense Advanced Research Projects Agency (DARPA).
 #
 # Permission is hereby granted, free of charge, to any person obtaining a copy of this data,
 # including any software or models in source or binary form, as well as any drawings, 
 # specifications, and documentation (collectively "the Data"), to deal in the Data without 
 # restriction, including without limitation the rights to use, copy, modify, merge, publish,
 # distribute, sublicense, and/or sell copies of the Data, and to permit persons to whom the Data
 # is furnished to do so, subject to the following conditions:
 #  
 # The above copyright notice and this permission notice shall be included in all copies or
 # substantial portions of the Data.
 #  
 # THE DATA IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING
 # BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
 # NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS, SPONSORS, DEVELOPERS, CONTRIBUTORS, OR
 # COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
 # OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE DATA OR THE
 # USE OR OTHER DEALINGS IN THE DATA.


#### Atmosphere1:
option ipopt options "halt_on_ampl_error=yes";
param Mac;
param Tac;
param Pac;
param R;
  
var M1>=0;
var Vac>=0;
var To1>=1;
var Po1>=1;
var cp1>=900,<=1400;
var gam1>=1.2,<=1.5;


subject to conM1: M1 = Mac;
subject to conTo1: To1 = Tac;
subject to conPo1: Po1 = Pac;
subject to conVac: Vac = Mac*sqrt(gam1*R*Tac);

subject to concp1:cp1 = (2.506*10^(-11)*To1^2+1.454*10^(-7)*To1^1.5
-4.246*10^(-7)*To1+3.162*10^(-5)*To1^(0.5)+1.3303-1.512*10^4*To1^(-1.5)
+3.063*10^5*To1^(-2)-2.212*10^7*To1^(-3))*10^3;
subject to congam1: gam1 = cp1/(cp1-R);
