湿位涡剖面图

*湿位涡剖面图
'reinit'
'open D:\66\ctatcrhuv\662.ctl'
'enable print D:\66\ctatcrhuv\14pv1pm.gmf'

'set ylevs 995 900 800 700 600 500 400 300'
*选lev和z都会出现cannot contour grid-all undefined values的错误
*'set lev 995 300'
'set z 2 20'
'set zlog on'
'set lon 108 125'
'set lat 20 37'

'set t 37'

*假相当位温
'define prs=lev'
'define es=(6.112*exp((17.67*tc)/(tc+243.5)))'
'define qs=0.62197*es/(prs-0.378*es)'
'define q=rh*qs/100'
'define e=prs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(tc+273.16)-log(e)-4.805)'
'define theta=(tc+273.16)*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'


*湿位涡3
'define vo=hcurl(u,v)'
'define f=2*7.292*sin(lat*3.14159/180.0)*0.00001'
'define g=9.8'
'define dp=100*(lev(z-1)-lev(z+1))'
'define deqt=eqt(z-1)-eqt(z+1)'
'define du=u(z-1)-u(z+1)'
'define dv=u(z-1)-u(z+1)'
'define dx=2.0*6370949.0*cos(lat*3.14159/180.0)*3.14159/180.0'
'define dy=2.0*6370949.0*3.14159/180.0'
'define dtx=cdiff(eqt,x)'
'define dty=cdiff(eqt,y)'
'define pv1=-g*(vo+f)*deqt/dp' 
'define pv2=g*((dv/dp)*(dtx/dx)-(du/dp)*(dty/dy))'
'define pv=pv1+pv2'

'set lon 113'
'set lat 21 37'
'set gxout contour'
'set csmooth on'
'set cthick 7'
'set clopts -1 -1 0.10'
'd pv1*10e5'
*'d pv2*10e6'
*'d pv*10e5'
'print'
'disable print'
;


ctl文件
dset D:\66\ctatcrhuv\662.dat
options  byteswapped
undef 1.e30
title  OUTPUT FROM WRF V3.1.1 MODEL
pdef  200 200 lcc  28.500  116.500  100.500  100.500  28.50000  28.50000  116.50000  10000.000  10000.000
xdef  512 linear  104.95424   0.04504504
ydef  423 linear   18.79075   0.04504504
zdef   24 levels 
1000.00000
995.00000
990.00000
980.00000
970.00000
960.00000
950.00000
900.00000
850.00000
800.00000
750.00000
700.00000
650.00000
600.00000
550.00000
500.00000
450.00000
400.00000
350.00000
300.00000
250.00000
200.00000
150.00000
100.00000
tdef   48 linear 12Z12JUN2004      60MN     
VARS    6
U             24  0  x-wind component (m s-1)
V             24  0  y-wind component (m s-1)
W             24  0  z-wind component (m s-1)
tc            24  0  Temperature (C)
theta         24  0  Potential Temperature (K)
rh            24  0  Relative Humidity (%)
ENDVARS
@ global String comment TITLE =  OUTPUT FROM WRF V3.1.1 MODEL
@ global String comment SIMULATION_START_DATE = 2004-06-10_12:00:00
@ global String comment GRIDTYPE = C
@ global String comment MMINLU = USGS
@ global String comment WEST-EAST_GRID_DIMENSION =   201
@ global String comment SOUTH-NORTH_GRID_DIMENSION =   201
@ global String comment BOTTOM-TOP_GRID_DIMENSION =    35
@ global String comment MAP_PROJ =     1
@ global String comment DX =     10000.00
@ global String comment DY =     10000.00
@ global String comment CEN_LAT =        28.50
@ global String comment CEN_LON =       116.50
@ global String comment TRUELAT1 =        28.50
@ global String comment TRUELAT2 =        28.50
@ global String comment MOAD_CEN_LAT =        28.50
@ global String comment STAND_LON =       116.50
@ global String comment DIFF_OPT =     19
@ global String comment KM_OPT =     48
@ global String comment DAMP_OPT =     0
@ global String comment KHDIF =         0.00
@ global String comment KVDIF =         0.00
@ global String comment MP_PHYSICS =     2
@ global String comment RA_LW_PHYSICS =     1
@ global String comment RA_SW_PHYSICS =     1
@ global String comment SF_SFCLAY_PHYSICS =     2
@ global String comment SF_SURFACE_PHYSICS =     2
@ global String comment BL_PBL_PHYSICS =     2
@ global String comment CU_PHYSICS =     3
@ global String comment SURFACE_INPUT_SOURCE =     1
@ global String comment SST_UPDATE =     0
@ global String comment GRID_FDDA =     0
@ global String comment FEEDBACK =     1
@ global String comment SMOOTH_OPTION =     0
@ global String comment W_DAMPING =     0
@ global String comment OBS_NUDGE_OPT =     0
@ global String comment GRID_ID =     1
@ global String comment PARENT_ID =     0
@ global String comment I_PARENT_START =     1
@ global String comment J_PARENT_START =     1
@ global String comment PARENT_GRID_RATIO =     1
@ global String comment DT =        30.00
@ global String comment ISWATER =    16
@ global String comment ISICE =    24
@ global String comment ISURBAN =     1
@ global String comment ISOILWATER =    14

（注：以上内容转自MNMUC）