求助 read_wrf_nc 修改风场资料

hdswy

要将西安站（108.91°E，34.30°N）周围50km范围内850-300hPa的风速资料（nc格式）增大10%，是使用 read_wrf_nc 吗？如果用这个，用-EditData VAR选项，编辑 read_wrf_nc .f 中的subroutine USER_CODE，在if块中如何改写？

ps：我在网上查到 已知地球上两地的地理坐标，求两地间的实际距离-----设两点A、B的经、纬度分别为(jA,wA)(jB,wB)，则半径为R的球面上两点间的距离(大圆弧)为： 弧AB=R*arccos[sin(wA)sin(wB)+cos(wA)cos(wB)*cos(jA-jB)]
不知是不是要用到这个公式来确定50km内的范围。

天目神眉

850-300hPa的风速资料，也就是说你要修改垂直高度上不止一层的变量，所以：
1.当然就是你说的区域问题，你得首先确定你要修改的区域，找到对应的（x,y）对应的值，可以用你所说的公式；
2.找到850-300hPa这个高度对应的是垂直层哪几层；
3.相关变量应该就是U;V吧，自己再查看确认一下。

godenflame135

首先要明白：
1.wrfout的文件中高度层不是1000hPa，850hPa等是p坐标，是sigma坐标，我们通常所说的850hPa、1000hPa高度是经过ARWpost后处理插值得到的。其次：
2.不知道为什么要将850-300hPa高度的风增加，但结果应该还要进入wrf数值模式？
   如果全部增加倒是可以办到，不过是sigma高度的。
  如果不进入wrf数值模式，用ARWpost后处理得到的dat文件即可。

hdswy

godenflame135 发表于 2012-3-25 12:02

                               
登录/注册后可看大图

首先要明白：
1.wrfout的文件中高度层不是1000hPa，850hPa等是p坐标，是sigma坐标，我们通常所说的850hPa、 ...

我修改的是进wrf前的资料，改完之后再进模式运行看结果有什么不同。850-300hPa在资料中对应的是垂直第7-17层。

godenflame135

那就更好整了，思路：打开旧文件，创建新文件。
把旧文件中需要修改的风场数组读出来并修改写入新文件（nc格式），其余的copy到新文件。WPS在没有grib1格式之前是可以读取nc格式的文件的吧！
工作量稍微有点大！
嘿嘿！可参考http://bbs.06climate.com/forum.php?mod=viewthread&tid=4904

cuixindong

hdswy 发表于 2012-3-25 12:12

                               
登录/注册后可看大图

我修改的是进wrf前的资料，改完之后再进模式运行看结果有什么不同。850-300hPa在资料中对应的是垂直第7-1 ...

怎么找WPS处理完的数据中气压场对应的是数据中的第几层

cuixindong

天目神眉 发表于 2012-3-25 10:24

                               
登录/注册后可看大图

850-300hPa的风速资料，也就是说你要修改垂直高度上不止一层的变量，所以：
1.当然就是你说的区域问题，你 ...

天目神眉大侠，您上次发的读取WPS处理完的数据的程序中，怎么找经纬度对应的XY？谢谢帮助！

hdswy

cuixindong 发表于 2012-4-3 09:13

                               
登录/注册后可看大图

怎么找WPS处理完的数据中气压场对应的是数据中的第几层

我的是met_开头的文件，一共27层，第1至27层分别是： 地面 1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 30 20 10

ztftom

你可以用你的方法算好longitude和latitude的范围，然后算好model level的范围。之后，你看一下你的grids中有哪些在这个范围内， 对这些格点上的数值修改。

天目有一个帖子讲model level转换的。你可以翻翻看。

ztftom

hdswy 发表于 2012-4-3 10:25

                               
登录/注册后可看大图

我的是met_开头的文件，一共27层，第1至27层分别是： 地面 1000 975 950 925 900 850 800 750 700 650 60 ...

哦，你修改met_文件了。那就直接改这些压力层上的数据就行了。