替换下垫面数据步骤
1.裁切MODIS数据
我们只需要研究区域的数据,因此需要对MODIS裁切,减少需要处理的数据量。以下所有的操作都是在linux下进行。
Tips:裁切的大小应该略大于WRF模拟中选中的区域,这样可以完全的覆盖住我们的研究区域。
步骤:“matlab”命令——左上角当前文件夹双击“MCD12C1.A2018001.006.2019200161458.hdf”进入“HDF Importort Tool”
上图总共4个位置需要修改,按照从上往下:
1.选择第一个
2.选择“geographic box”
3.Corner1是网格左下角的的经纬度(可略大),corner2是网格右上角的经纬度(可略大)
经纬度可通过geogrid步骤生成的geo_em.d01.nc文件,用ncl查看
4.修改为“Band3”
5.点击右下角“import”
在matlab右上角工作空间这里会出现Band3的大小802(y)*761(x)
2.转换成nc格式
vi landuse.m
上图总共3个位置需要修改,按照从上往下:
1.根据上一步骤标黄的x,y修改
2.同上
3.按照第一步修改的经纬度修改,注意各自对应
然后run lanuse.m
Tip:这里报错说明你目录里面存在上次替换留下的文件,这里每替换一次需要删除一次
3.执行fortran文件
修改一些位置的fname(共两处)以及nx与ny
执行f90文件
ifort -o test.exe Replacelanduse.f90 -L/software/netcdf/lib -lnetcdf -I/software/netcdf/include
./test.exe
然后就在目录的原位置覆盖掉原来的数据了,注意检查是否成功
依次修改d02和d03的数据。在下次修改之前,先把刚刚生成的红框文件删除掉。