前几天看了兰溪之水利用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
欢迎关注榴莲气象
