登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
在编译WPS成功后,会生成三个可执行文件-geogrid.exe,ungrib.exe,metgrid.exe。第一个是用于静态地形数据进行水平插值,第二个用于气象数据(FNL等数据)的水平插值,第三个则是把前二者结合。
./geogrid.exe后,会生成n个(n=你的嵌套层数)geo_em.d0*.nc.
./ungrib.exe后,会生成n个(n=你的嵌套层数)FILE*, GRIB*.
./metgrid.exe后,会生成n个(n=你的嵌套层数)met_em.d0*.nc
同样,在./real.exe后,会有wrfinput,wrfbdy文件,最后的输出结果有wrfout文件,对于这些NC文件,我们是否可以随时查看呢?.
当然可以,但是前提是你已成功安装NETCDF,
命令:ncdump *.nc>*.txt(第一个*表示你要读的NC文件,第二个表示你要存的txt文件名),这样,你查看txt看看,就可以看到所有的变量名以及其数值了。
但是,有时候我们不需要查看那么多变量,因为这些NC文件中都有几十种变量,我们不一定都关心,而且数据多了,查看也慢。再来一个看到某个变量的命令:
ncdump –v T2 *.nc>*.txt,这就是查看T2的数值了。如果要查看两个及其以上的呢:
ncdump –v U;V *.nc>*.txt
好了,看看效果:
我提取的是T2,
1
前面N行的表头,写明了该NC文件中所有的变量名,以及变量名的意义,存放顺序,单位等信息。
2
这就是T2的描述,三维的,时间,南北,东西。注意循环次序。2m的温度,单位K。
3
从上面的数据,我们可以清晰的看出,第一大块是在第一个纬度值上,各个经度点上的数值。也就是先固定南北,从东到西循环,然后再从南到北循环。最后,时次的循环。
同样的,我们同样也可以把数据写进去,替换原来的数据。这个就需要附件中的程序了。
4
重点看一下这一段,这个是自己加的,这是要改变其中的某个变量是用的。
5
这个是自己把文件读进去,然后完全替换原来的数据,当然这个就需要你先把数据处理好了。
改后怎么运行,程序前几行有说明,一般为
pgf90 read_wrf_nc.f -L/usr/local/netcdf/lib -lnetcdf -lm \ -I/usr/local/netcdf/include -Mfree -o read_wrf_nc 然后改变数据;
./read_wrf_nc geo_em.d04.nc -EditData LU_INDEX
有提示,输入yes即可。注意,要改什么就把命令改成相应的变量和文件。
当然,这个才是第一步,你会了这个,还有很多东西要做,慢慢来吧!
read_wrf_nc.f
(65.99 KB, 下载次数: 350)
|