ncl做敏感性试验

用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

欢迎关注榴莲气象