wrf修改土壤湿度

Wrf修改土壤湿度思路与方法（初始场替换数据）思路：不同下垫面的土壤湿度不同，但NCEP资料初始化的时候因为模拟区域太小，湿度值基本不变，用不同下垫面类型修改对应的不同土壤湿度，用MODIS_landuse资料20种土地类型判别替换裸地，耕地，草地的土壤湿度。Read_wrf_nc可以查看，修改wrf模式的初始场，结果等netcdf格式的文件1.下载read_wrf_nc.f（wrf主页）2.pgf90read_wrf_nc.f-L/usr/local/netcdf/lib-lnetcdf-lm-I/usr/local/netcdf/include-Mfree-oread_wrf_nc3.命令与用法参考usersGuide，例如：-h/helpPrinthelpinformation.列出帮助文件-mPrintlistoffieldsavailableforeachtime,plustheminandmaxvaluesforeachfield.列出变量信息（相当于ncdump–h）-tslllatlonVARVAR…..willgeneratetimeseriesoutputforallVAR’satx/ylocationnearesttolat/lon列出时间变量序列值-VVAR查看当前时间变量值4.初始场文件wrfinput_d01拷贝到read_wrf_nc文件夹下5.打开read_wrf_nc.f到1600行，添加变量data_real1real,dimension(dim1,dim2,dim3)::data_real,data_real1到1620行添加语句：open(9999,file='/root/download/readwrfnc/1.txt')if(var=='LU_INDEX')thendata_real1=data_real!写入到文本文件1.txtwrite(9999,*)data_real1重新编译一下read_wrf_nc.f文件，修改变量的命令为：./read_wrf_ncwrfinput_d03–EditDataLU_INDEX选yes6.在文本编辑器里按照wrf的LANDUSE.TBL文件查看MODIS的相应的土壤湿度值，替换数据。例如：1,12.,.30,.95,50.,4.,3.33,29.2e5,'EvergreenNeedleleafForest'常绿针叶林2,12.,.50,.95,50.,5.,1.67,29.2e5,'EvergreenBroadleafForest'常绿阔叶林3,14.,.30,.94,50.,4.,2.86,25.0e5,'DeciduousNeedleleafForest'落叶针叶林4,16.,.30,.93,50.,4.,2.63,25.0e5,'DeciduousBroadleafForest'落叶阔叶林5,13.,.30,.97,50.,4.,2.11,41.8e5,'MixedForests'混交林6,22.,.10,.93,5.,3.,1.56,20.8e5,'ClosedShrublands'封闭灌丛7,20.,.15,.95,6.,3.,2.14,20.8e5,'OpenShrublands'开放灌丛8,22.,.10,.93,5.,3.,1.56,20.8e5,'WoodySavannas'热带稀疏草原9,20.,.15,.92,15.,3.,2.00,25.0e5,'Savannas'热带、亚热带稀疏大草原10,19.,.15,.96,12.,3.,2.37,20.8e5,'Grasslands'草场、牧场11,14.,.42,.95,30.,5.5,1.3235.5e5'Permanentwetlands'永久湿地12,17.,.30,.985,15.,4.,2.71,25.0e5,'Croplands'农田、耕地13,15.,.10,.88,80.,3.,1.67,18.9e5,'UrbanandBuilt-Up'城市和建筑1418.,.25,.98,14.,4.,2.56,25.0e5,'cropland/naturalvegetationmosaic'耕地、自然植被镶嵌15,55.,.95,.95,0.1,5.,0.,9.0e25,'SnowandIce'冰雪16,25.,.02,.90,1.,2.,0.81,12.0e5,'BarrenorSparselyVegetated'稀疏荒芜植被17,8.,1.0,.98,0.01,6.,0.,9.0e25,'Water'水体18,15.,.50,.93,30.,5.,2.67,9.0e25,'WoodedTundra'林地苔藓19,15.,.50,.92,15.,5.,2.67,9.0e25,'MixedTundra'混合苔藓20,25.,.02,.90,10.,2.,1.60,12.0e5,'BarrenTundra'贫瘠苔藓7.wrfinput_d03的SMOIS量有4层，LU_INDEX只有一层，把修改过的文本内容往下复制3次，这样把每层土壤湿度都用同一值了，有时间在乘以系数修订。8.替换：open(9999,file='/root/download/readwrfnc/1.txt')read(9999,*)data_real1if(var=='SMOIS')thendata_real=data_real1!替换再重新编译一下read_wrf_nc.f文件，修改变量的命令为：./read_wrf_ncwrfinput_d03–EditDataSMOIS选yes9.把修改过的wrfinput_d03替换原始文件，直接运行./wrf.exe10.昀保险的做法是把wps前处理的初始场文件met_em做修改，但是嵌套和时次太多，麻烦。