用grads从WRf模拟结果中提取数据，写出对应的ctl文件，但是画不出图。

shirly

从WRF模拟结果中提取变量，写出对应的ctl文件，但是总出错。将原始数据的ctl文件、提取变量的gs文件、提取出的数据的ctl文件贴出：

WRf模拟结果的ctl文件

（注：原始文件有101个变量，我截图时候没有截全。
另外endvars之后有很多行类似这样的语句：@ global String comment TITLE =  OUTPUT FROM WRF V3.5.1 MODEL）

提取变量的gs文件（只提取一个点的数据：lon=104.12，lat=35.85）：
'reinit'
   'open E:\data\winter\wmai.ctl'
   'set fwrite E:\data\winter\tiqu\u.v.q.t.dat'
   'set gxout fwrite'
   'set lon 104.12'
   'set lat 35.85'
    tt=1
   while(tt<=193)
   'set t 'tt
   zz=1
   while(zz<=29)
   'set z 'zz
  'd u'   'd v'  
  'd t'
  'd qvapor'
   zz=zz+1
   endwhile
   tt=tt+1
   endwhile
  'disable fwrite'
  'reinit'


提取出的数据ctl文件：
dset E:/data/winter/tiqu/u.v.q.t.dat
options  byteswapped
undef 1.e30
title  OUTPUT FROM WRF V3.5.1 MODEL
pdef   73  58 lcc  35.850  104.120   37.000   29.500 60.00000  30.00000  104.00000  30000.000  30000.000
Xdef  1 linear   104.12   0.13513514        
Ydef  1 linear   35.85   0.13513514
Zdef   29 linear 1 1  
Tdef  193 linear 00Z09JAN2009      15MN      
VARS  4
U             29  0  x-wind component (m s-1)
V             29  0  y-wind component (m s-1)
T             29  0  perturbation potential temperature (theta-t0) (K)
QVAPOR        29  0  Water vapor mixing ratio (kg kg-1)
ENDVARS

（注：**xdef和ydef的初始经纬度若改成和原始数据的一样，画图时不出错，显示说都是缺省值，没有图；改成要提取数据的经纬度时，画图显示错误如下）

运行出错结果

第一次转化WRF模拟结果，没有说清楚的地方还请谅解，希望有前辈可以帮忙看一下~不胜感激！

andrewsoong

不清楚。帮你顶。为什么不用NCL ？

shirly

@godenflame135找到了你发的贴子，提取的时候选取的水平区域是所有的点，然后用fortran提取出来某一个点的。我是grads提取变量的时候选取了一个点，这样是不是不行呢？必须得提取出一个区域的才行？

river

把你编的ctl 文件里的options项和pdef项去掉试试。只有一个点的数据还有什么字节顺序与本机相反啊，还有什么特殊投影啊。

shirly

river 发表于 2013-11-30 21:38

                               
登录/注册后可看大图

把你编的ctl 文件里的options项和pdef项去掉试试。只有一个点的数据还有什么字节顺序与本机相反啊，还有什么 ...

谢谢你的指导，我今晚按照godenflame135的办法提取出来了，我明天按你的方法再试试看，谢谢river啦~！

river

shirly 发表于 2013-11-30 22:12

                               
登录/注册后可看大图

谢谢你的指导，我今晚按照godenflame135的办法提取出来了，我明天按你的方法再试试看，谢谢river啦~！{:e ...

先试试看行不行吧，我没怎么弄过wrf的资料

做个霸气的木头

不是只提取了u的数据吗，为嘛ctl里还有v什么的

shirly

做个霸气的木头 发表于 2013-12-1 23:26

                               
登录/注册后可看大图

不是只提取了u的数据吗，为嘛ctl里还有v什么的

首先感谢你的回复，那个当时是想同时提取四个时候写的，不出图就先提取了一个，ctl粘成之前写的了，我在这改一下gs就对应上了。

奋豆儿

楼主你好我现在也在做同样的尝试，把wrf后处理出来的结果用grads同时提取多个需要的变量，但是提取出来感觉有些问题，仅仅提取了9个变量生成的.grd数据比wrf结果本身都大，而且在某个变量的某一层画图的时候发现值明显不对。之前用同样的方法提取过单个变量没有问题。想请教您，您最后是怎么处理的？谢谢！

shirly

奋豆儿 发表于 2014-1-4 10:55

                               
登录/注册后可看大图

楼主你好我现在也在做同样的尝试，把wrf后处理出来的结果用grads同时提取多个需要的变量，但是提取出来感觉 ...

不好意思啊，我才看到，那我估计是你的循环写错了，你具体的gs我也不知道，所以不是很明白你的问题所在。