请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 66111|回复: 82

[分享资料] 湿位涡剖面图gs示例

  [复制链接]

新浪微博达人勋

发表于 2011-8-19 22:59:53 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 Mouse 于 2011-8-30 20:47 编辑

*湿位涡剖面图
'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)

评分

参与人数 2威望 +5 金钱 +26 贡献 +13 收起 理由
mofangbao + 3 + 10 + 10 必须狠狠的给分,这个问题几乎是每个群里都.
传说中的谁 + 2 + 16 + 3 很好!感谢鼠标无私的分享!

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2011-8-19 23:02:34 | 显示全部楼层
收下了,好东西啊,还有木有私货啊,别藏着了。感谢鼠标!!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2011-8-19 23:04:30 | 显示全部楼层

点评

你的在线时间涨了哇 哈哈  发表于 2011-8-20 11:35
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2011-8-19 23:29:47 | 显示全部楼层
Mouse 发表于 2011-8-19 23:04
知道你最近在整一些物理量的gs,加油啊

我昨天刚编好,你今天就发了,诚心气我的,不跟你玩儿了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2011-8-20 19:13:41 | 显示全部楼层
谢谢鼠标,我也正在收集呢,好东西
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2011-8-21 09:54:55 | 显示全部楼层
感谢分享。。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2011-8-21 10:37:42 | 显示全部楼层
感谢分享,很好的指导教本,收藏了,找时间用下。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2011-8-22 08:48:59 | 显示全部楼层
非常感谢分享,收藏了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2011-8-22 10:58:16 | 显示全部楼层
好东西,收藏了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2011-8-30 18:15:03 | 显示全部楼层
这个算出来对吗,我也这样算的
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表