用ncl修改wrfinput做敏感性實驗（受蘭溪之水啟發(fā)

前幾天看了蘭溪之水利用ncl做敏感性試驗的帖子http://bbs.06climate.com/forum.php?mod=viewthread&tid=4592&extra=page%3D1&page=1
覺得很受啟發(fā)，就結(jié)合自己的例子進行了嘗試，并且做了一些改進，也在論壇上進行了詢問http://bbs.06climate.com/forum.php?mod=viewthread&tid=38276&page=1#pid526317可惜沒人搭理我啊。。。。。。
現(xiàn)在自己終于調(diào)好了，分享給大家吧~~（第一次分享自己東西，不太正確的地方還請大家指正）

我的個例是用自己同化的土壤濕度資料（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"；聲明維數(shù)名稱的命令 
    soil=fbindirread("/home/lw/soil/soilana.dat",0,(/nx,ny/),"float")；根據(jù)寫入時的大小讀入   
  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

歡迎關(guān)注榴蓮氣象