读取相同经纬度的站点资料和wrfout格点资料二者变量的代码

felin

经过一周的编写，调试，迷茫，惆怅后，终于实现的自己的一小步！
现在将程序贴出来，供有需要的同学使用！
压缩包中有五个文件，使用过程如下：
1.先确定你需要的站点的站点名（不知道的百度，搜：“气象台站名”），将站点名存放在station.txt中，保存后关闭！
2.将你的站点资料存放进data文件夹中（我用的是micaps的domain1资料，请解压后放入，文件名为“***.000”）!
3.启动“read_sta.f90”，编译、执行，如有报错，请仔细调试！（错误一般多是没有.txt文件啦啥的，自己新建个就行）
4.将输出的结果按“站点_经度_纬度”存放至“sta_lon_lat.txt”中，方便“read_wrf.f90”读取。（不要问为啥，自己笨，没把两个.f90合起来，有好心人下来合成下呗！）
5.我模拟后的结果“wrfout_d***.nc”用ARWPots处理的，所以放上读取".ctl"和“.dat”后处理文件的grads脚本“fwrite.gs”,修改下输出变量就好。（这里当然还可以选择NCL，我也编写了NCL脚本，那个更简单，只是我没想好，暂时没发）（这里涉及的模拟变量都是高度为1的，即“zlev 1”，4D变量我今后在尝试读取，写出来发给大家。）
6.启动“read_wrf.f90”，编译、执行，如有报错，请仔细调试！（错误一般和上一个类似）
7.最后，你生成了站点的变量和模拟资料的变量，想干啥就干啥（站点资料和模拟资料的时间不是一一对应哈！！！！还需各位需要处理下，我想应该不难，就当练手吧！）
感谢下面帖子的作者和清风（忘记清风的帖子存在什么地方，大约是12年的帖子！）
http://bbs.06climate.com/forum.p ... =6532&extra=&page=1

mezhe

辛苦！学习下！

放逐流年

没有ncl程序的吗？楼主

放逐流年

能否加个好友，向你学习学习，楼主，谢谢了
我的qq是51067585

felin

放逐流年 发表于 2017-4-9 22:00
没有ncl程序的吗？楼主

当时没想好怎么写NCL脚本，你可以试试，有个命令是wrf_user_ll_to_ij好像可以一步到位

泡泡妹儿

先看看吧，说不定以后用得上呢

饶平仔

楼主你好，看了你这个，感觉受益颇多。想多问问，就是这些都是Fortran的程序。
nf_get_var_real这个函数可以读取u/v/xlat/xlong等变量，因为都是real型。但Times是字符型，nf_get_var_char和nf_get_var_character都不能读取，有些不明白，网上一直找不到答案，求助！

felin

饶平仔 发表于 2017-9-21 11:05
楼主你好，看了你这个，感觉受益颇多。想多问问，就是这些都是Fortran的程序。
nf_get_var_real这个函数可 ...

nf_get_var_char是啥命令？

字俣丞

感谢楼主，非常实用

两米阳光

饶平仔 发表于 2017-9-21 11:05
楼主你好，看了你这个，感觉受益颇多。想多问问，就是这些都是Fortran的程序。
nf_get_var_real这个函数可 ...

你好，请问问题解决了吗？nf_get_var_char为什么不能读取，搞明白了吗？