.MCAD 303010000 1 0 16547 0
.CMD PLOTFORMAT
0 0 1 1 0 0 1 
0 0 1 1 0 0 1 
0 1 0 0 NO-TRACE-STRING
0 2 1 0 NO-TRACE-STRING
0 3 2 0 NO-TRACE-STRING
0 4 3 0 NO-TRACE-STRING
0 1 4 0 NO-TRACE-STRING
0 2 5 0 NO-TRACE-STRING
0 3 6 0 NO-TRACE-STRING
0 4 0 0 NO-TRACE-STRING
0 1 1 0 NO-TRACE-STRING
0 2 2 0 NO-TRACE-STRING
0 3 3 0 NO-TRACE-STRING
0 4 4 0 NO-TRACE-STRING
0 1 5 0 NO-TRACE-STRING
0 2 6 0 NO-TRACE-STRING
0 3 0 0 NO-TRACE-STRING
0 4 1 0 NO-TRACE-STRING
0 1 21 15 0 3 
.CMD FORMAT  rd=d ct=10 im=i et=5 zt=15 pr=3 mass length time charge temperature tr=0 vm=0
.CMD SET ORIGIN 1
.CMD SET TOL 0.001000000000000
.CMD SET PRNCOLWIDTH 8
.CMD SET PRNPRECISION 4
.CMD PRINT_SETUP 1.200000 0.000000 0.000000 0.000000 0
.CMD HEADER_FOOTER 1 1 *empty* *empty* *empty* 0 1 *empty* *empty* *empty*
.CMD HEADER_FOOTER_FONT fontID=14 family=Arial points=10 bold=0 italic=0 underline=0
.CMD HEADER_FOOTER_FONT fontID=15 family=Arial points=10 bold=0 italic=0 underline=0
.CMD DEFINE_FONTSTYLE_NAME fontID=0 name=Variables
.CMD DEFINE_FONTSTYLE_NAME fontID=1 name=Constants
.CMD DEFINE_FONTSTYLE_NAME fontID=2 name=Text
.CMD DEFINE_FONTSTYLE_NAME fontID=4 name=User^1
.CMD DEFINE_FONTSTYLE_NAME fontID=5 name=User^2
.CMD DEFINE_FONTSTYLE_NAME fontID=6 name=User^3
.CMD DEFINE_FONTSTYLE_NAME fontID=7 name=User^4
.CMD DEFINE_FONTSTYLE_NAME fontID=8 name=User^5
.CMD DEFINE_FONTSTYLE_NAME fontID=9 name=User^6
.CMD DEFINE_FONTSTYLE_NAME fontID=10 name=User^7
.CMD DEFINE_FONTSTYLE fontID=0 family=Arial points=10 bold=0 italic=0 underline=0
.CMD DEFINE_FONTSTYLE fontID=1 family=Arial points=10 bold=0 italic=0 underline=0
.CMD DEFINE_FONTSTYLE fontID=2 family=Arial points=10 bold=0 italic=0 underline=0
.CMD DEFINE_FONTSTYLE fontID=4 family=Arial points=10 bold=0 italic=0 underline=0
.CMD DEFINE_FONTSTYLE fontID=5 family=Courier^New points=10 bold=0 italic=0 underline=0
.CMD DEFINE_FONTSTYLE fontID=6 family=System points=10 bold=0 italic=0 underline=0
.CMD DEFINE_FONTSTYLE fontID=7 family=Script points=10 bold=0 italic=0 underline=0
.CMD DEFINE_FONTSTYLE fontID=8 family=Roman points=10 bold=0 italic=0 underline=0
.CMD DEFINE_FONTSTYLE fontID=9 family=Modern points=10 bold=0 italic=0 underline=0
.CMD DEFINE_FONTSTYLE fontID=10 family=Times^New^Roman points=10 bold=0 italic=0 underline=0
.CMD UNITS U=1
.CMD DIMENSIONS_ANALYSIS 0 0
.TXT 4 1 16167 0
Cg a50.375000,85.375000,109
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain {ATMOSPHERIC THERMODYNAMICS FUNCTIONS
 with examples. \par Atmospheric thermodynamics}
 for the practical person.}
}
.TXT 7 0 16170 0
Cg b84.375000,84.375000,427
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain {This section contains the standard
 thermodynamic functions normally placed
 at the beginning \par of a program for 
use later in specific programs.\par The 
functions are not documented in specific
 programs and are used like standard mathematic
 functions.\par In this FUNCTION EXAMPLE
 file the functions are used at the right
 of their definition to show how \par they
 work.  The numbers in the examples can 
be changed to see how they work.}}
}
.TXT 16 0 1 0
Cg a41.875000,61.000000,284
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain {1.0 Define CONSTANTS\par \par TFC
 = Top of freezing band (C)\par FB = Freezing
 band (K)\par \par CPA = Specific heat of
 air at constant pressure (J/kg)\par RA 
= Gas constant of air \par KA = RA/CPA\par 
\par }{CPV = Specific heat of vapor at constant
 pressure (J/kg)\par RV = Gas constant of
 vapor\par KV = RV/CPV\par }}
}
.EQN 34 0 16478 0
{0:TFC}NAME~-10
.EQN 0 14 16479 0
{0:FB}NAME~20
.EQN 5 -14 3 0
{0:CPA}NAME~1004.675
.EQN 4 0 2 0
{0:KA}NAME~(2)/(7)
.EQN 0 41 16425 0
{0:KA}NAME={0}?_n_u_l_l_
.EQN 8 -41 4 0
{0:RA}NAME~{0:CPA}NAME*{0:KA}NAME
.EQN 0 41 16423 0
{0:RA}NAME={0}?_n_u_l_l_
.EQN 8 -41 6 0
{0:CPV}NAME~1846.04
.EQN 4 0 5 0
{0:KV}NAME~(1)/(4)
.EQN 5 0 7 0
{0:RV}NAME~{0:CPV}NAME*{0:KV}NAME
.EQN 0 42 16424 0
{0:RV}NAME={0}?_n_u_l_l_
.TXT 5 -42 16429 0
Cg a37.125000,83.000000,56
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain {Ratio of gas constant of air to
 gas constant of vapor}}
}
.EQN 4 42 16422 0
{0:\e}NAME={19027}?_n_u_l_l_
.EQN 1 -42 8 0
{0:\e}NAME~({0:RA}NAME)/({0:RV}NAME)
.TXT 5 0 16431 0
Cg a55.750000,83.000000,84
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain {Functions for specific heat, gas
 constant, and their ratio for air-vapor
 mixtures}}
}
.TXT 2 40 16480 0
Cg a40.375000,43.000000,107
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Specific at constant pressure, gas
 constant \par and their ratio for 1 kg 
air containing 10 g of water vapor.}
}
.EQN 6 -40 16330 0
{0:fCPA}NAME({0:M}NAME):{0:CPA}NAME+{0:CPV}NAME*{0:M}NAME
.EQN 0 41 16427 0
{0:fCPA}NAME(0.01)={0}?_n_u_l_l_
.EQN 5 -41 16331 0
{0:fRA}NAME({0:M}NAME):{0:RA}NAME+{0:RV}NAME*{0:M}NAME
.EQN 0 41 16428 0
{0:fRA}NAME(0.01)={0}?_n_u_l_l_
.EQN 5 -41 16334 0
{0:fKA}NAME({0:M}NAME):({0:fRA}NAME({0:M}NAME))/({0:fCPA}NAME({0:M}NAME))
.EQN 1 41 16426 0
{0:fKA}NAME(.01)={0}?_n_u_l_l_
.EQN 1 23 16460 0
{0:MX}NAME:0.055
.EQN 3 0 16461 0
{0:fKA}NAME({0:MX}NAME)={0}?_n_u_l_l_
.TXT 3 -64 16432 0
Cg b83.000000,83.000000,453
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain {Other Constants\par \par P0 = Base
 pressure for air = 100 kPa\par PV0 = Base
 pressure for vapor = 0.61068 kPa\par T0
 = Base temperature = 273.15 K\par TF = 
Freezing temperature in K\par G = Accelleration
 of gravity = 9.8 m/s2\par PG, TG, MG Starting
 point guess for solver blocks.\par LV0 
= Latent heat of vaporization of water at
 T0\par LF0 = Latent heat of solidification
 of water at T0\par CW = Specific heat of
 liquid water = 4190 J/kg\par CI = Specific
 heat of ice = 2090 J/kg\par }}
}
.EQN 31 0 9 0
{0:P0}NAME~100
.EQN 0 19 10 0
{0:PV0}NAME~0.61068
.EQN 3 -19 12 0
{0:T0}NAME~273.15
.EQN 0 19 14 0
{0:TF}NAME~{0:T0}NAME+{0:TFC}NAME
.EQN 0 20 18 0
{0:G}NAME~9.8
.EQN 4 -39 11 0
{0:PG}NAME~70
.EQN 0 19 15 0
{0:TG}NAME~280
.EQN 0 18 16 0
{0:MG}NAME~0.01
.EQN 4 -37 19 0
{0:LV0}NAME~2500840
.EQN 0 19 20 0
{0:LF0}NAME~333660
.EQN 4 -19 21 0
{0:CW}NAME~4190
.EQN 0 19 22 0
{0:CI}NAME~2090
.TXT 4 -19 16434 0
Cg b83.000000,83.000000,79
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain {Constants used in calculation of
 the vapor pressure of liquid water and 
ice.}}
}
.EQN 5 0 23 0
{0:\aL}NAME~52.91688
.EQN 0 19 24 0
{0:\aI}NAME~25.48597
.EQN 3 -19 25 0
{0:\bL}NAME~6806.171
.EQN 0 19 26 0
{0:\bI}NAME~6286.1912
.EQN 3 -19 27 0
{0:\GL}NAME~5.078893
.EQN 0 19 28 0
{0:\GI}NAME~0.528613
.TXT 10 -19 16180 0
C x1,1,0,0
.TXT 5 0 29 0
Cg a49.250000,65.250000,191
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain {2.0  FUNCTIONS: Vapor pressure,
 mixing ratio, virtual temperature\par fPVL(T)
 Vapor pressure of liquid water\par }{fPVI(T)
 Vapor pressure of ice\par }{fPV(T) Vapor
 pressure of liquid water}{ or ice}}
}
.TXT 10 45 16481 0
Cg a30.750000,38.000000,55
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Vapor pressure of water and ice 
at 270 K\par both in kPa.}
}
.EQN 7 -45 30 0
{0:fPVL}NAME({0:T}NAME)~{0:exp}NAME({0:\aL}NAME-({0:\bL}NAME)/({0:T}NAME)-{0:\GL}NAME*{0:ln}NAME({0:T}NAME))
.EQN 0 45 16417 0
{0:fPVL}NAME(270)={19029}?_n_u_l_l_
.EQN 7 -45 31 0
{0:fPVI}NAME({0:T}NAME)~{0:exp}NAME({0:\aI}NAME-({0:\bI}NAME)/({0:T}NAME)-{0:\GI}NAME*{0:ln}NAME({0:T}NAME))
.EQN 0 45 16418 0
{0:fPVI}NAME(270)={19029}?_n_u_l_l_
.EQN 7 -45 32 0
{0:fPV}NAME({0:T}NAME)~{0:if}NAME({0:T}NAME>{0:TF}NAME,{0:fPVL}NAME({0:T}NAME),{0:fPVI}NAME({0:T}NAME))
.EQN 0 45 16419 0
{0:fPV}NAME(270)={19029}?_n_u_l_l_
.TXT 5 0 16482 0
Cg a26.125000,38.000000,38
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Define conditons for use in functions}
}
.TXT 1 -45 16435 0
Cg b83.000000,83.000000,15
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain {Mixing Ratio}}
}
.EQN 2 46 16443 0
{0:P1}NAME:100
.EQN 0 12 16444 0
{0:T1}NAME:290
.EQN 0 10 16445 0
{0:M1}NAME:0.01
.EQN 2 -22 16487 0
{0:P2}NAME:50
.EQN 0 12 16488 0
{0:T2}NAME:260
.EQN 0 10 16489 0
{0:M2}NAME:0.02
.EQN 2 -22 16525 0
{0:Z2}NAME:1000
.EQN 0 12 16484 0
{0:P3}NAME:20
.EQN 5 -58 33 0
{0:fMVS}NAME({0:P}NAME,{0:T}NAME)~{0:\e}NAME*(({0:fPV}NAME({0:T}NAME))/({0:P}NAME-{0:fPV}NAME({0:T}NAME)))
.TXT 0 46 16485 0
Cg a27.000000,37.000000,47
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Saturation Mixing Ratio at 100 kPa
 \par and 290 K}
}
.EQN 6 0 16447 0
{0:MS1}NAME:{0:fMVS}NAME({0:P1}NAME,{0:T1}NAME)
.EQN 1 -46 16129 0
{0:fMVL}NAME({0:P}NAME,{0:T}NAME)~{0:\e}NAME*(({0:fPVL}NAME({0:T}NAME))/({0:P}NAME-{0:fPVL}NAME({0:T}NAME)))
.EQN 1 46 16448 0
{0:MS1}NAME*1000={0}?_n_u_l_l_
.TXT 3 1 16486 0
Cg a24.625000,36.000000,68
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Vapor pressure of air at 100 kPa\par 
with a misxing ratio of 10 g/kg. }
}
.EQN 5 -47 34 0
{0:fPV2}NAME({0:P}NAME,{0:M}NAME)~{0:M}NAME*({0:P}NAME)/({0:M}NAME+{0:\e}NAME)
.EQN 0 47 16439 0
{0:fPV2}NAME({0:P1}NAME,{0:M1}NAME)={0}?_n_u_l_l_
.EQN 6 -47 35 0
{0:fMS3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)~{0:fMVS}NAME({0:P}NAME,{0:T}NAME)
.TXT 1 48 16490 0
Cg a26.375000,35.000000,38
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Mixing ratio water in the vapor 
phase}
}
.EQN 2 -48 36 0
{0:fMV3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)~{0:if}NAME({0:fMVS}NAME({0:P}NAME,{0:T}NAME)>{0:M}NAME,{0:M}NAME,{0:fMVS}NAME({0:P}NAME,{0:T}NAME))
.EQN 1 48 16453 0
{0:fMV3}NAME({0:P1}NAME,{0:T2}NAME,{0:M1}NAME)*1000={0}?_n_u_l_l_
.TXT 3 1 16493 0
Cg a12.500000,34.000000,18
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Relative Humidity}
}
.EQN 3 -49 37 0
{0:fU3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME):100*(({0:fPV2}NAME({0:P}NAME,{0:M}NAME))/({0:fPVL}NAME({0:T}NAME)))
.EQN 0 49 16442 0
{0:fU3}NAME({0:P1}NAME,{0:T1}NAME,{0:M1}NAME)={0}?_n_u_l_l_
.TXT 6 0 16495 0
Cg a11.000000,34.000000,15
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Vapor pressure}
}
.EQN 2 -49 38 0
{0:fPV3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)~{0:fMV3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)*({0:P}NAME)/({0:\e}NAME+{0:fMV3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME))
.EQN 1 49 16465 0
{0:fPV3}NAME({0:P1}NAME,{0:T1}NAME,{0:M1}NAME)={0}?_n_u_l_l_
.TXT 4 1 16496 0
Cg a15.375000,33.000000,24
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Partial pressure of air}
}
.EQN 3 -50 39 0
{0:fPA3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)~{0:P}NAME-{0:fPV3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)
.EQN 0 50 16451 0
{0:fPA3}NAME({0:P1}NAME,{0:T1}NAME,{0:M1}NAME)={0}?_n_u_l_l_
.EQN 4 -50 40 0
{0:fMC3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)~{0:if}NAME({0:fMVS}NAME({0:P}NAME,{0:T}NAME)<{0:M}NAME,{0:M}NAME-{0:fMVS}NAME({0:P}NAME,{0:T}NAME),0)
.EQN 4 0 41 0
{0:fMLB}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)~{0:if}NAME({0:T}NAME>{0:TF}NAME-{0:FB}NAME,{0:fMC3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)*({0:T}NAME-{0:TF}NAME+{0:FB}NAME)/({0:FB}NAME),0)
.TXT 4 52 16497 0
Cg a23.625000,31.000000,35
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Mixing ratio water in liquid phase}
}
.EQN 4 -52 42 0
{0:fML3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)~{0:if}NAME({0:T}NAMEň{0:TF}NAME,{0:fMC3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME),{0:fMLB}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME))
.EQN 0 51 16452 0
{0:fML3}NAME({0:P1}NAME,{0:T2}NAME,{0:M1}NAME)*1000={0}?_n_u_l_l_
.EQN 6 -51 43 0
{0:fMIB}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)~{0:if}NAME({0:T}NAME<{0:TF}NAME,{0:fMC3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)*({0:TF}NAME-{0:T}NAME)/({0:FB}NAME),0)
.TXT 2 52 16498 0
Cg a22.000000,31.000000,32
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Mixing ratio water in ice phase}
}
.EQN 4 -52 44 0
{0:fMI3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)~{0:if}NAME({0:T}NAME<{0:TF}NAME-{0:FB}NAME,{0:fMC3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME),{0:fMIB}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME))
.EQN 0 53 16454 0
{0:fMI3}NAME({0:P1}NAME,{0:T2}NAME,{0:M1}NAME)*1000={0}?_n_u_l_l_
.TXT 2 -53 16179 0
C x1,1,0,0
.TXT 4 0 45 0
Cg a10.500000,11.250000,14
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain {3.0 ENTROPY}}
}
.TXT 0 53 16508 0
Cg a13.250000,30.000000,20
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Entropy of pure air}
}
.EQN 4 -53 46 0
{0:fSA}NAME({0:P}NAME,{0:T}NAME)~{0:CPA}NAME*{0:ln}NAME(({0:T}NAME)/({0:T0}NAME))-{0:RA}NAME*{0:ln}NAME(({0:P}NAME)/({0:P0}NAME))
.EQN 0 53 16502 0
{0:fSA}NAME({0:P0}NAME,{0:T0}NAME)={0}?_n_u_l_l_
.EQN 2 0 16507 0
{0:fSA}NAME(100,{0:T0}NAME+1)={0}?_n_u_l_l_
.TXT 4 0 16509 0
Cg a16.000000,30.000000,23
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Entropy of water vapor}
}
.EQN 2 -53 47 0
{0:fSV}NAME({0:P}NAME,{0:T}NAME)~{0:CPV}NAME*{0:ln}NAME(({0:T}NAME)/({0:T0}NAME))-{0:RV}NAME*{0:ln}NAME(({0:P}NAME)/({0:PV0}NAME))+({0:LV0}NAME)/({0:T0}NAME)
.EQN 1 53 16506 0
{0:fSV}NAME({0:PV0}NAME,{0:T0}NAME)={0}?_n_u_l_l_
.EQN 3 0 16503 0
{0:fSV}NAME(1,300)={0}?_n_u_l_l_
.EQN 4 -53 48 0
{0:fSL}NAME({0:T}NAME)~{0:CW}NAME*{0:ln}NAME(({0:T}NAME)/({0:T0}NAME))
.TXT 0 53 16510 0
Cg a16.000000,30.000000,24
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Entropy of liquid water}
}
.EQN 4 0 16504 0
{0:fSL}NAME({0:T0}NAME)={0}?_n_u_l_l_
.EQN 2 0 16511 0
{0:fSL}NAME({0:T0}NAME+1)={0}?_n_u_l_l_
.EQN 2 -53 49 0
{0:fSI}NAME({0:T}NAME)~{0:CI}NAME*{0:ln}NAME(({0:T}NAME)/({0:T0}NAME))-({0:LF0}NAME)/({0:T0}NAME)
.TXT 3 53 16512 0
Cg a10.000000,30.000000,15
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Entropy of ice}
}
.EQN 4 0 16505 0
{0:fSI}NAME({0:T2}NAME)={0}?_n_u_l_l_
.EQN 1 -53 50 0
{0:fSA3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)~{0:fSA}NAME({0:fPA3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME),{0:T}NAME)
.EQN 3 0 51 0
{0:fSV3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)~{0:fMV3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)*{0:fSV}NAME({0:fPV3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME),{0:T}NAME)
.EQN 3 0 52 0
{0:fSL3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)~{0:fML3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)*{0:fSL}NAME({0:T}NAME)
.EQN 3 0 53 0
{0:fSI3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)~{0:fMI3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)*{0:fSI}NAME({0:T}NAME)
.EQN 3 0 54 0
{0:fST3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)~{0:fSA3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)+{0:fSV3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)+{0:fSL3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)+{0:fSI3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)
.TXT 3 53 16499 0
Cg a23.125000,30.000000,35
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Entropy of moist air per kg of air}
}
.EQN 2 -53 16130 0
{0:fSM3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME):({0:fST3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME))/((1+{0:M}NAME))
.EQN 2 53 16466 0
{0:S1}NAME:{0:fST3}NAME({0:P1}NAME,{0:T1}NAME,{0:M1}NAME)
.EQN 2 -53 55 0
{0:fSC}NAME({0:T}NAME)~{0:if}NAME({0:T}NAME>{0:TF}NAME,{0:fSL}NAME({0:T}NAME),{0:fSI}NAME({0:T}NAME))
.EQN 0 53 16469 0
{0:S1}NAME={0}?_n_u_l_l_
.EQN 5 -53 56 0
{0:fAA3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)~{0:fSA}NAME({0:P}NAME,{0:T}NAME)+{0:M}NAME*{0:fSC}NAME({0:T}NAME)
.TXT 1 53 16500 0
Cg a27.125000,30.000000,38
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Entropy moist air per kg of substance}
}
.EQN 3 0 16456 0
{0:fSM3}NAME({0:P1}NAME,{0:T1}NAME,{0:M1}NAME)={0}?_n_u_l_l_
.TXT 3 -3 16501 0
Cg a29.125000,32.000000,69
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Potential temperature of air containing
 \par water in vapor phase only.}
}
.EQN 1 -50 57 0
{0:f\qM4}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME,{0:PC}NAME)~{0:T}NAME*((({0:PC}NAME)/({0:P}NAME)))^(({0:RA}NAME+{0:RV}NAME*{0:M}NAME)/({0:CPA}NAME+{0:CPV}NAME*{0:M}NAME))
.EQN 5 50 16457 0
{0:f\qM4}NAME(50,250,0.01,100)={0}?_n_u_l_l_
.TXT 4 -50 16173 0
C x1,1,0,0
.TXT 4 52 16517 0
Cg a23.875000,31.000000,58
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain enthalpy of pure air, pure vapor,\par 
liquid water, and ice.}
}
.TXT 1 -52 58 0
Cg a11.250000,12.250000,15
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain {4.0 ENTHALPY}}
}
.EQN 4 0 59 0
{0:fHA}NAME({0:T}NAME)~{0:CPA}NAME*({0:T}NAME-{0:T0}NAME)
.EQN 0 52 16513 0
{0:fHA}NAME({0:T1}NAME)={0}?_n_u_l_l_
.EQN 3 -52 60 0
{0:fHV}NAME({0:T}NAME)~({0:CPV}NAME*({0:T}NAME-{0:T0}NAME)+{0:LV0}NAME)
.EQN 2 52 16514 0
{0:fHV}NAME({0:T1}NAME)={0}?_n_u_l_l_
.EQN 1 -52 61 0
{0:fHL}NAME({0:T}NAME)~({0:CW}NAME*({0:T}NAME-{0:T0}NAME))
.EQN 3 0 62 0
{0:fHI}NAME({0:T}NAME)~({0:CI}NAME*({0:T}NAME-{0:T0}NAME)-{0:LF0}NAME)
.EQN 0 52 16515 0
{0:fHL}NAME({0:T1}NAME)={0}?_n_u_l_l_
.EQN 4 -52 63 0
{0:fHA3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)~{0:fHA}NAME({0:T}NAME)
.EQN 0 52 16516 0
{0:fHI}NAME({0:T2}NAME)={0}?_n_u_l_l_
.EQN 3 -52 64 0
{0:fHV3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)~{0:fMV3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)*{0:fHV}NAME({0:T}NAME)
.EQN 3 0 65 0
{0:fHL3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)~{0:fML3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)*{0:fHL}NAME({0:T}NAME)
.TXT 1 63 16519 0
Cg a16.125000,20.000000,38
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Enthalpy of moist air\par per kg
 of air.}
}
.EQN 2 -63 66 0
{0:fHI3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)~{0:fMI3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)*{0:fHI}NAME({0:T}NAME)
.EQN 4 0 67 0
{0:fHT3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)~{0:fHA3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)+{0:fHV3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)+{0:fHL3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)+{0:fHI3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)
.EQN 3 63 16471 0
{0:H1}NAME:{0:fHT3}NAME({0:P1}NAME,{0:T1}NAME,{0:M2}NAME)
.EQN 1 -63 68 0
{0:fHC}NAME({0:T}NAME)~{0:if}NAME({0:T}NAME>{0:TF}NAME,{0:fHL}NAME({0:T}NAME),{0:fHI}NAME({0:T}NAME))
.EQN 3 0 69 0
{0:fEE3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)~{0:fHA}NAME({0:T}NAME)+{0:M}NAME*{0:fHC}NAME({0:T}NAME)
.EQN 0 63 16474 0
{0:H1}NAME={19024}?_n_u_l_l_
.EQN 4 -63 70 0
{0:fHW4}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME,{0:W}NAME)~{0:fHA}NAME({0:T}NAME)+{0:M}NAME*{0:fHV}NAME({0:T}NAME)+{0:if}NAME({0:W}NAME>{0:TF}NAME,{0:fHL}NAME({0:W}NAME),{0:fHI}NAME({0:W}NAME))*({0:fMVS}NAME({0:P}NAME,{0:W}NAME)-{0:M}NAME)
.TXT 3 50 16523 0
Cg a32.625000,33.000000,47
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Static energy of air at 0 and 1000
 m elevation}
}
.EQN 1 -50 71 0
{0:f\s4}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME,{0:Z}NAME)~{0:fHT3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)+{0:G}NAME*(1+{0:M}NAME)*{0:Z}NAME
.EQN 3 51 16524 0
{0:f\s4}NAME({0:P1}NAME,{0:T1}NAME,{0:M1}NAME,{0}0)={0}?_n_u_l_l_
.TXT 3 -51 72 0
Cg a32.250000,34.250000,37
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain {5.0 VIRTUAL TEMPERATURE AND HEIGHT}}
}
.EQN 0 51 16526 0
{0:f\s4}NAME({0:P1}NAME,{0:T1}NAME,{0:M1}NAME,{0:Z2}NAME)={0}?_n_u_l_l_
.TXT 9 -1 16542 0
Cg a21.000000,32.000000,32
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Virtual temperature of dry air.}
}
.EQN 3 1 16541 0
{0:fTV}NAME({0:T1}NAME,{0:M1}NAME)={0}?_n_u_l_l_
.EQN 3 -51 73 0
{0:fTV}NAME({0:T}NAME,{0:M}NAME)~{0:T}NAME*((1+({0:M}NAME)/({0:\e}NAME))/(1+{0:M}NAME))
.EQN 9 0 74 0
{0:fTV3}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME)~{0:if}NAME({0:M}NAME<{0:fMVS}NAME({0:P}NAME,{0:T}NAME),{0:fTV}NAME({0:T}NAME,{0:M}NAME),{0:fTV}NAME({0:T}NAME,{0:fMVS}NAME({0:P}NAME,{0:T}NAME))*((1+{0:fMVS}NAME({0:P}NAME,{0:T}NAME))/(1+{0:M}NAME)))
.TXT 6 53 16543 0
Cg a27.375000,30.000000,60
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Virtual temperature of air containing
 \par water in any phase.}
}
.EQN 4 -53 75 0
{0:f\a}NAME({0:P1}NAME,{0:P2}NAME,{0:TV1}NAME,{0:TV2}NAME)~({0:G}NAME)/({0:RA}NAME)*(({0:ln}NAME(({0:TV2}NAME)/({0:TV1}NAME)))/({0:ln}NAME(({0:P2}NAME)/({0:P1}NAME))))
.EQN 3 52 16475 0
{0:fTV3}NAME({0:P1}NAME,{0:T1}NAME,{0:M2}NAME)={0}?_n_u_l_l_
.TXT 3 0 16545 0
Cg a31.000000,31.000000,67
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Expansion exponent and Vertical 
distance between sounding levels.}
}
.EQN 6 -52 76 0
{0:fZ}NAME({0:P1}NAME,{0:P2}NAME,{0:TV1}NAME,{0:TV2}NAME)~({0:TV2}NAME-{0:TV1}NAME)/({0:f\a}NAME({0:P1}NAME,{0:P2}NAME,{0:TV1}NAME,{0:TV2}NAME))
.EQN 2 52 16547 0
{0:f\a}NAME(100,80,300,280)={0}?_n_u_l_l_
.TXT 3 -52 16172 0
C x1,1,0,0
.EQN 0 52 16546 0
-{0:fZ}NAME(100,80,300,280)={0}?_n_u_l_l_
.TXT 8 -52 16163 0
Cg a13.750000,85.375000,16
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain {SOLVER BLOCKS}}
}
.TXT 3 0 16131 0
Cg a32.000000,32.250000,35
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain {ISENTROPIC EXPANSION TEMPERATURE}}
}
.EQN 4 0 16132 0
{0:GIVEN}NAME
.TXT 0 48 16527 0
Cg a34.000000,35.000000,158
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain True adiabatic expansion temperature
 in K.\par (reversible isentropic expansion
 without \par separation of the condensate
 and with freezing of the condensate)   
}
}
.EQN 8 -48 16133 0
{0:fST3}NAME({0:P}NAME,{0:TG}NAME,{0:M}NAME)÷{0:S}NAME
.EQN 3 0 16134 0
{0:fTSOL}NAME({0:S}NAME,{0:P}NAME,{0:M}NAME):{0:FIND}NAME({0:TG}NAME)
.EQN 0 48 16476 0
{0:fTSOL}NAME({0:S1}NAME,20,{0:M1}NAME)={0}?_n_u_l_l_
.TXT 4 -48 16135 0
Cg a28.875000,29.250000,32
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain {ISENTROPIC EXPANSION PRESSURE}}
}
.EQN 4 0 16136 0
{0:GIVEN}NAME
.TXT 2 48 16528 0
Cg a18.125000,35.000000,25
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Freezing pressure in kPa}
}
.EQN 1 -48 16137 0
{0:fST3}NAME({0:PG}NAME,{0:T}NAME,{0:M}NAME)÷{0:S}NAME
.EQN 3 0 16138 0
{0:fPSOL}NAME({0:S}NAME,{0:T}NAME,{0:M}NAME):{0:FIND}NAME({0:PG}NAME)
.EQN 0 48 16477 0
{0:fPSOL}NAME({0:S1}NAME,{0:T0}NAME,{0:M1}NAME)={0}?_n_u_l_l_
.TXT 4 -48 16139 0
Cg a33.500000,34.250000,37
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain {ISENTROPIC DESSICATION TEMPERATURE}}
}
.EQN 4 0 16140 0
{0:GIVEN}NAME
.EQN 3 0 16141 0
{0:fAA3}NAME({0:P}NAME,{0:TG}NAME,{0:M}NAME)÷{0:S}NAME
.EQN 3 0 16142 0
{0:fASOL}NAME({0:S}NAME,{0:P}NAME,{0:M}NAME):{0:FIND}NAME({0:TG}NAME)
.TXT 4 0 16143 0
Cg a58.000000,59.250000,62
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain {ISENTHALPIC DESSICATION TEMPERATURE,
 EQUIVALENT TEMPERATURE}}
}
.EQN 4 0 16144 0
{0:GIVEN}NAME
.EQN 3 0 16145 0
{0:fEE3}NAME({0:P}NAME,{0:TG}NAME,{0:M}NAME)÷{0:H}NAME
.EQN 3 0 16146 0
{0:fESOL}NAME({0:H}NAME,{0:P}NAME,{0:M}NAME):{0:FIND}NAME({0:TG}NAME)
.TXT 4 0 16147 0
Cg a24.750000,26.250000,29
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain {MIXING RATIO FROM WET BULB}}
}
.TXT 2 52 16530 0
Cg a24.875000,31.000000,58
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Mixing ratio of air with a wet bulb\par 
temperature of 20 C.}
}
.EQN 2 -52 16148 0
{0:GIVEN}NAME
.EQN 3 0 16149 0
{0:fHT3}NAME({0:P}NAME,{0:W}NAME,{0:fMVS}NAME({0:P}NAME,{0:W}NAME))÷{0:fHW4}NAME({0:P}NAME,{0:T}NAME,{0:MG}NAME,{0:W}NAME)
.EQN 3 0 16150 0
{0:fMSOL}NAME({0:P}NAME,{0:T}NAME,{0:W}NAME):{0:FIND}NAME({0:MG}NAME)
.EQN 1 52 16529 0
{0:fMSOL}NAME({0:P0}NAME,{0:T0}NAME+30,{0:T0}NAME+20)*1000={0}?_n_u_l_l_
.TXT 3 -52 16151 0
Cg a24.750000,26.250000,29
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain {WET BULB FROM MIXING RATIO}}
}
.EQN 4 0 16152 0
{0:GIVEN}NAME
.TXT 0 53 16535 0
Cg a27.125000,30.000000,50
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Wet bulb of air with a mixing ratio
 of\par 10.7 g/kg}
}
.EQN 3 -53 16153 0
{0:fHT3}NAME({0:P}NAME,{0:TG}NAME,{0:fMVS}NAME({0:P}NAME,{0:TG}NAME))÷{0:fHW4}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME,{0:TG}NAME)
.EQN 3 0 16154 0
{0:fWSOL}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME):{0:FIND}NAME({0:TG}NAME)
.EQN 1 53 16533 0
{0:fWSOL}NAME({0:P0}NAME,{0:T0}NAME+30,0.0107)-{0:T0}NAME={0}?_n_u_l_l_
.TXT 3 -53 16155 0
Cg a24.750000,26.250000,29
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain {LIFTING CONDENSATION LEVEL}}
}
.EQN 4 0 16156 0
{0:GIVEN}NAME
.TXT 1 54 16538 0
Cg a24.375000,29.000000,35
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Lifting Condensation Level in kPa.}
}
.EQN 2 -54 16157 0
{0:fU3}NAME({0:PG}NAME,{0:f\qM4}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME,{0:PG}NAME),{0:M}NAME)-100÷0
.EQN 3 0 16158 0
{0:fCSOL}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME):{0:FIND}NAME({0:PG}NAME)
.EQN 0 53 16537 0
{0:fCSOL}NAME({0:P1}NAME,{0:T1}NAME,{0:M1}NAME)={0}?_n_u_l_l_
.TXT 4 -53 16159 0
Cg a9.125000,9.250000,12
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain {DEW POINT}}
}
.TXT 2 54 16539 0
Cg a16.625000,30.000000,23
{\rtf1\ansi \deff0
{\fonttbl
{\f0\fnil Arial;}
}
{\plain Dew Point in degree C.}
}
.EQN 2 -54 16160 0
{0:GIVEN}NAME
.EQN 3 0 16161 0
{0:fPV2}NAME({0:P}NAME,{0:M}NAME)÷{0:fPV}NAME({0:TG}NAME)
.EQN 3 0 16162 0
{0:fDSOL}NAME({0:P}NAME,{0:T}NAME,{0:M}NAME):{0:FIND}NAME({0:TG}NAME)
.EQN 0 53 16540 0
{0:fDSOL}NAME({0:P1}NAME,{0:T1}NAME,{0:M1}NAME)-{0:T0}NAME={0}?_n_u_l_l_