用ncl修改wrfinput做敏感性实验（受兰溪之水启发）

东风急流

前几天看了兰溪之水利用ncl做敏感性试验的帖子http://bbs.06climate.com/forum.php?mod=viewthread&tid=4592&extra=page%3D1&page=1
觉得很受启发，就结合自己的例子进行了尝试，并且做了一些改进，也在论坛上进行了询问http://bbs.06climate.com/forum.php?mod=viewthread&tid=38276&page=1#pid526317可惜没人搭理我啊。。。。。。
现在自己终于调好了，分享给大家吧~~（第一次分享自己东西，不太正确的地方还请大家指正）

我的个例是用自己同化的土壤湿度资料（soilana.dat）去代替wrfinput中的第一层土壤湿度
;======================================================
; This script is used to read a binary file and write the var to the netCDF file.
;======================================================
begin
    a=addfile("/home/lw/soil/wrfinput.nc","w");自己可以用ncdump看一下里面的土壤湿度到底是什么样子存放的
    s=new((/1,4,300,400/),"float")
    s=a->SMOIS
  print(s(0,1,97,66))；检验是否读入

    soil=new((/300,400/),"float")
    soil!0="Time"
    soil!1="soil_layers_stag"；声明维数名称的命令
    soil=fbindirread("/home/lw/soil/soilana.dat",0,(/nx,ny/),"float")；根据写入时的大小读入   
  print(soil(66,97))；检验是否读入
;-----------------------
; change the field  
;-----------------------
   do i=1,300,1
     do j=1,400,1
       s(0,0,i-1,j-1)=soil(i-1,j-1)
     end do
   end do
  print(s(0,0,66,97))；检验是否和上面一样
;-----------------------
;write the field
;-----------------------
   a->SMOIS=s
end

还有一个不太懂的问题想请教大家，就是维数名称的命名，我是根据提示错误才这么做的。但是其实命名的应该是纬度和经度，不知道这样后面会不会出错啊？

andrewsoong

东风急流 发表于 2016-9-19 10:15
我是对wrfinput进行处理的
方法就是这个帖子

这样不会出现“漏风”现象么

东风急流

想@一下兰溪之水的，可惜和大神不是好友啊

兰溪之水

猫小哩

求问你的土壤资料是怎么同化的啊？

东风急流

不好意思，我不是做资料同化的啊。这个资料是别人给我，让我用wrf给他模拟结果。

八维伊布

维度名称叫什么都行，自己不要搞乱就好，只是个名称，除非后面有函数要寻找特定的维度名称……
本例里面，你2D变量soil的维度名称用在什么地方了？

东风急流

本来没有用的，但那个循环里面的赋值语句（s(0,0,i-1,j-1)=soil(i-1,-1)）编译不过去，提示我左边的变量有维数名称，右边的没有，让我把左边的去掉。但我不会去，就把右边的给写上了......

东风急流

s这个变量本来是有维数名称的，其中前两维就是我给soil的那两个

东风急流

八维伊布 发表于 2015-10-19 07:51
维度名称叫什么都行，自己不要搞乱就好，只是个名称，除非后面有函数要寻找特定的维度名称……
本例里面， ...

请问两个变量相互赋值，是不是要维数名称对应呢？

小预报员

好专业啊不懂