爱气象，爱气象家园！ 
·设为首页 开启辅助访问 充值积分新帖子新回复热门贴精华帖

气象家园

 找回密码
 立即注册

QQ登录

只需一步，快速开始

搜索
气象家园»论坛 › 编程作图 › GrADS › JAR-25 读取程序
12下一页
返回列表 发新帖
查看: 6952|回复: 10

[分享资料] JAR-25 读取程序

[复制链接]
kongfeng0824
发表于 2013-11-18 09:18:50 | 显示全部楼层 |阅读模式

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

您需要 登录 才可以下载或查看，没有帐号？立即注册

x
  1. ctl说明文件：

  3. dset  G:\jra25\%y4\%m2\fcst_phy2m.%y4%m2%d2%h2

  5. index G:\jra25\fcst_phy2m.1979_2004.idx

  7. undef 9.999E+20

  9. title fcst_phy2m

  11. dtype grib 255

  13. options yrev template

  15. xdef 320 linear 0.0 1.125

  17. ydef 160 levels

  19. -89.142 -88.029 -86.911 -85.791 -84.670 -83.549 -82.428 -81.307 -80.185 -79.064

  21. -77.943 -76.821 -75.700 -74.578 -73.457 -72.336 -71.214 -70.093 -68.971 -67.850

  23. -66.728 -65.607 -64.485 -63.364 -62.242 -61.121 -60.000 -58.878 -57.757 -56.635

  25. -55.514 -54.392 -53.271 -52.149 -51.028 -49.906 -48.785 -47.663 -46.542 -45.420

  27. -44.299 -43.177 -42.056 -40.934 -39.813 -38.691 -37.570 -36.448 -35.327 -34.205

  29. -33.084 -31.962 -30.841 -29.719 -28.598 -27.476 -26.355 -25.234 -24.112 -22.991

  31. -21.869 -20.748 -19.626 -18.505 -17.383 -16.262 -15.140 -14.019 -12.897 -11.776

  33. -10.654  -9.533  -8.411  -7.290  -6.168  -5.047  -3.925  -2.804  -1.682  -0.561

  35.   0.561   1.682   2.804   3.925   5.047   6.168   7.290   8.411   9.533  10.654

  37. 11.776  12.897  14.019  15.140  16.262  17.383  18.505  19.626  20.748  21.869

  39. 22.991  24.112  25.234  26.355  27.476  28.598  29.719  30.841  31.962  33.084

  41. 34.205  35.327  36.448  37.570  38.691  39.813  40.934  42.056  43.177  44.299

  43. 45.420  46.542  47.663  48.785  49.906  51.028  52.149  53.271  54.392  55.514

  45. 56.635  57.757  58.878  60.000  61.121  62.242  63.364  64.485  65.607  66.728

  47. 67.850  68.971  70.093  71.214  72.336  73.457  74.578  75.700  76.821  77.943

  49. 79.064  80.185  81.307  82.428  83.549  84.670  85.791  86.911  88.029  89.142

  51. tdef 37988 linear 00Z01jan1979  6hr

  53. zdef 1 linear 1 1

  55. * Level suffix abbreviation: sfc=Surface; col=Total Column; toa=Top of atmosphere

  57. * Last Column Abbreviation: AC=Accumulation; MN=Time-Mean; SNP=Snap-shot (instantaneous)

  59. vars 60

  61. CPRATsfc 0 63,1,0   ** 6-hourly convective precipitation rate       [mm/day] AC

  63. CSDLFsfc 0 163,1,0  ** Clear sky downward long wave flux at surface  [W/m^2] MN

  65. CSDSFsfc 0 161,1,0  ** Clear sky downward solar flux at surface      [W/m^2] MN

  67. CSULFtoa 0 162,8,0  ** Clear sky upward long wave flux at top        [W/m^2] MN

  69. CSUSFsfc 0 160,1,0  ** Clear sky upward solar flux at surface        [W/m^2] MN

  71. CSUSFtoa 0 160,8,0  ** Clear sky upward solar flux at top            [W/m^2] MN

  73. CWATcol  0 76,1,0   ** Total cloud water                            [kg/m^2] AC

  75. DISTsfc  0 8,1,0    ** Model surface height                              [m] SNP

  77. DLWRFsfc 0 205,1,0  ** Downward long wave radiation flux at surface  [W/m^2] MN

  79. DSWRFsfc 0 204,1,0  ** Downward solar radiation flux at surface      [W/m^2] MN

  81. DSWRFtoa 0 204,8,0  ** Downward solar radiation flux at top          [W/m^2] MN

  83. FGLUsfc  0 147,1,0  ** Zonal momemtum flux by long gravity wave      [N/m^2] MN

  85. FGLVsfc  0 148,1,0  ** Meridional momemtum flux by long gravity wave [N/m^2] MN

  87. FGSUsfc  0 159,1,0  ** Zonal momemtum flux by short gravity wave     [N/m^2] MN

  89. FGSVsfc  0 154,1,0  ** Meridional momemtum flux by short gravity wave [N/m^2] MN

  91. FRQCsfc  0 169,1,0  ** Frequency of cumulus precipitation                [%]

  93. FRQPsfc  0 168,1,0  ** Frequency of precipitation                        [%]

  95. HCDC     0 75,1,0   ** High level cloud cover                            [%] MN

  97. ICECsfc  0 91,1,0   ** Ice concentration (ice=1; no ice=0)             [1/0] MN

  99. LCDC     0 73,1,0   ** Low level cloud cover                             [%] MN

  101. LHTFLsfc 0 121,1,0  ** Latent heat flux                              [W/m^2] MN

  103. MCDC     0 74,1,0   ** Mid level cloud cover                             [%] MN

  105. LPRATsfc 0 62,1,0   ** 6-hourly large-scale precipitation rate      [mm/day] AC

  107. PRESsfc  0 1,1,0    ** Surface pressure                                 [Pa] SNP

  109. PRMSLmsl 0 2,102,0  ** Mean pressure reduced to MSL                     [Pa] MN

  111. PRMXsfc  0 220,1,0  ** Maximum hourly precipitaion                 [mm/hour]

  113. PWATcol  0 54,1,0   ** Precipitable water                           [kg/m^2] SNP

  115. SPFHMsfc 0 139,1,0  ** Surface specific humidity                     [kg/kg] MN

  117. RHsfc    0 52,1,0   ** Surface relative humidity                         [%] SNP

  119. SFCRsfc  0 83,1,0   ** Surface roughness                                 [m] SNP

  121. SHTFLsfc 0 122,1,0  ** Sensible heat flux                            [W/m^2] MN

  123. SPFHsfc  0 51,1,0   ** Specific humidity                             [kg/kg] SNP

  125. SPFHhbl  0 51,109,1 ** Specific humidity at the lowest model level   [kg/kg] SNP

  127. SRWEQsfc 0 64,1,0   ** 6-hourly snowfall rate water equivalent    [kg/m^2/s] AC

  129. TCDCcol  0 71,1,0   ** Total cloud cover                                 [%] MN

  131. TMAXsfc  0 15,1,0   ** Surface 2m air maximum temperature                [K]

  133. TMINsfc  0 16,1,0   ** Surface 2m air minimum temperature                [K]

  135. TMPsfc   0 11,1,0   ** Surface 2m air temperature                        [K] SNP

  137. TMPhbl   0 11,109,1 ** Temperature at the lowest model level             [K] SNP

  139. TMPMsfc  0 138,1,0  ** Surface 2m air temperature                        [K] MN

  141. TTHMcol  0 218,1,0  ** Column total of diabatic heat                 [W/m^2] MN

  143. UFLXsfc  0 124,1,0  ** Zonal momentum flux by friction               [N/m^2] MN

  145. UGRDsfc  0 33,1,0   ** Surface 10m zonal wind                          [m/s] SNP

  147. UGRDhbl  0 33,109,1 ** Zonal wind  at lowest model level               [m/s] SNP

  149. ULWRFsfc 0 212,1,0  ** Upward long wave radiation flux at surface    [W/m^2] MN

  151. ULWRFtoa 0 212,8,0  ** Upward long wave radiation flux at top        [W/m^2] MN

  153. UGRDMsfc 0 136,1,0  ** Surface 10m zonal wind                          [m/s] MN

  155. USWRFsfc 0 211,1,0  ** Upward solar radiation at surface             [W/m^2] MN

  157. USWRFtoa 0 211,8,0  ** Upward solar radiation at top                 [W/m^2] MN

  159. UWVcol   0 157,1,0  ** Column total of zonal water vapor flux       [kg/m/s] MN

  161. VFLXsfc  0 125,1,0  ** Meridional momentum flux by friction          [N/m^2] MN

  163. VGRDsfc  0 34,1,0   ** Surface 10m meridional wind                     [m/s] SNP

  165. VGRDhbl  0 34,109,1 ** Meridional wind at lowest model level           [m/s] SNP

  167. VGRDMsfc 0 137,1,0  ** Surface 10m meridional wind                     [m/s] MN

  169. VTUcol   0 200,1,0  ** Column total of TU                         [K*Pa*m/s] MN

  171. VTVcol   0 201,1,0  ** Column total of TV                         [K*Pa*m/s] MN

  173. VWVcol   0 152,1,0  ** Column total of meridional water vapor flux  [Kg/m/s] MN

  175. WSMXsfc  0 219,1,0  ** Surface 10m maximum surface wind                [m/s]

  177. WSMXhbl  0 219,109,1 ** Maximum surface wind at lowest model level     [m/s]

  179. WTMPsfc  0 80,1,0   ** SST and sea ice surface temperature               [K] SNP

  181. ENDVARS

  183. * WTMPsfc includes reversed weighted-mean ground and canopy tempereture

  185. * over land area. Please use GrADS commands like;

  187. * d abs(WTMPsfc)             ---> Effective skin temperature

  189. * d maskout(WTMPsfc,WTMPsfc) ---> SST and sea ice surface temperature


  192. 批量读取出图程序:

  194. 'reinit'

  196. 'open f:\fcst_phy2m.1979_2004.ctl'

  198. year=1979

  200. **设定开始的年

  202. month=1

  204. day=1

  206. hour=00

  208. t=1

  210. while(year<2011)

  212. while(month<13)

  214. while(day<32)

  216. while(hour<24)

  218. if(month<10)

  220.     if(day<10)

  222.         if(hour<10)

  224. ***  设定输出的字符，输出的格式是4位年+2位月

  226. aaa=year%0%month%0%day%0%hour

  228. ***如果月小于10,则在月前加0，保证月是两位

  230.         else

  232.           aaa=year%0%month%0%day%hour

  234.         endif

  236.     else

  238.      aaa=year%0%month%day%hour

  240.     endif

  242. else

  244. aaa=year%month%day%hour

  246. endif




  251. 'set t 't   




  256. 'set gxout shaded'

  258. 'set mpdset cnworld'

  260. 'set grads off'




  265. **30-36:  red (rr%: 0, 600%)

  267. *'set rgb    30 255 215 225'

  269. *'set rgb    31 255 190 190'

  271. *'set rgb    32 255 160 150'

  273. *'set rgb    33 255 130 115'

  275. *'set rgb     34 255 100 85'

  277. *'set rgb      35 255 70 50'

  279. *'set rgb       36 255 40 0'

  281. *'set rgb      37 161 1   1'

  283. *'set rgb      38 121  1  1'

  285. *'set rgb      39 80   1  1'

  287. **40-46:green (rr%: 0, 600%)

  289. *'set rgb 40 225 255 225'

  291. *'set rgb 41 190 255 190'

  293. *'set rgb 42 150 255 150'

  295. *'set rgb 43 115 255 115'

  297. *'set rgb 44 75  255 75 '

  299. *'set rgb 45 45  255 45 '

  301. *'set rgb 46 0   255   0'

  303. **50-56:  blue (rr%: 0,600%)

  305. *'set rgb   50  200 254 254'

  307. *'set rgb   51  178 251 254'

  309. *'set rgb   52  155 242 254'

  311. *'set rgb   53  130 230 254'

  313. *'set rgb   54  108 211 254'

  315. *'set rgb    55  90 191 254'

  317. *'set rgb   56   70 166 254'

  319. *'set rgb   57   42 136 254'

  321. *'set rgb   58   24 105 254'

  323. *'set rgb     59   0 64 254'

  325. **40-46: yello (rr%: 0,600%)

  327. *'set rgb   60  254 254 230'

  329. *'set rgb   61  254 254 205'

  331. *'set rgb   62  254 254 180'

  333. *'set rgb   63  254 253 155'

  335. *'set rgb   64  254 253 130'

  337. *'set rgb   65  254 250 105'

  339. *'set rgb   66  254 240  80'

  341. *'set rgb   67  254 238  50'

  343. **50-56:purple (rr%: 0, 600%)  

  345. *'set rgb 70 225 225 255'

  347. *'set rgb 71 190 190 255'

  349. *'set rgb 72 150 150 255'

  351. *'set rgb 73 115 115 255'

  353. *'set rgb 74 75 75 255'

  355. *'set rgb 75 45 45 255'

  357. *'set rgb 76 0 0 255'




  362. *'set clevs 0.001 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0'

  364. *'set ccols 0 71 51 54 55 42 45 62 64 66 33 35'

  366. 'set csmooth on'

  368. 'd TMAXsfc'

  370. 'draw title JAR_'%aaa%''

  372. 'cbarn.gs '

  374. 'enable print f:\1\JAR_'%aaa%'.gmf'

  376. 'printim  f:\1\JAR_'%aaa%'.png x1600 y1200 white'

  378. 'print'

  380. 'disable print'

  382. 'c'




  387. t=t+6

  389. hour=hour+6

  391. endwhile

  393. hour=0

  395. day=day+1

  397. endwhile




  402. day=1

  404. month=month+1

  406. endwhile




  411. month=1

  413. year=year+1

  415. endwhile



  419. 'reinit'
复制代码
以上程序仅供抛砖引玉。以前用到的时候写的，水平不是很高，请大家指正
程序

评分

参与人数 2金钱 +18 贡献 +4 收起 理由
river + 6 + 1 赞一个!
mofangbao + 12 + 3

查看全部评分

相关帖子

密码修改失败请联系微信：mofangbao
回复

使用道具 举报

linikiku
发表于 2013-11-18 10:01:49 | 显示全部楼层
先收藏，谢谢！
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

qxtlyf
发表于 2013-11-18 10:32:32 | 显示全部楼层
谢谢版主，学习了
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

北京小脚
发表于 2013-11-21 08:54:19 | 显示全部楼层
是JRA-25吧？谢谢分享
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

rainyday1218
发表于 2014-2-21 15:29:55 | 显示全部楼层
请问楼主fcst_phy2m的意思是预报场吗，那还是再分析资料吗？
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

kongfeng0824
 楼主| 发表于 2014-2-21 18:49:26 | 显示全部楼层
rainyday1218 发表于 2014-2-21 15:29

                               
登录/注册后可看大图

请问楼主fcst_phy2m的意思是预报场吗，那还是再分析资料吗？

是日本的再分系资料
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

dannyjoyride
发表于 2014-4-16 11:20:29 | 显示全部楼层
同求anl和fcst的区别
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

hxyj
发表于 2015-10-30 09:36:40 | 显示全部楼层
楼主能否说说下载JRA的流程？谢谢了
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

藕茭白空心菜
发表于 2016-11-28 16:33:07 | 显示全部楼层
dannyjoyride 发表于 2014-4-16 11:20
同求anl和fcst的区别

Analysis & forecast
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

Haryu
发表于 2017-6-21 18:50:19 | 显示全部楼层
请问这个资料是.grib还是.grd呀
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

下一页 »
12下一页
返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

小黑屋|帮助中心|清除痕迹

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

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

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