NCL 计算整层水汽通量

hxyj

请问楼主是否已解决该问题了？能否将计算的完整脚本贴出共享？

奚子惠

hxyj 发表于 2015-6-21 09:51
请问楼主是否已解决该问题了？能否将计算的完整脚本贴出共享？

是的已经解决了
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
begin
;************************************************
; open file and read in variable
;***********************************************
;levels 1000 925 850 700 600 500 400 300  
; 层数    0   1   2   3   4   5   6   7      
;***********************************************
     nlat=73
     nlon=144
         nyear=35   
     nlev=8         
         
         p       = (/ 1000.,925.,850.,700.,600.,500., \
               400.,300.,250.,200.,150.,100., \
                70.,50.,30.,20.,10. /)
         linlog = 1
     pbot   = 1100.                                             
     ptop   = 300.
         u_file=addfile("uwnd_Dtrend.nc","r")
     v_file=addfile("vwnd_Dtrend.nc","r")
         p_file=addfile("pres_Dtrend.nc","r")
         shum_file=addfile("shum_Dtrend.nc","r")
         xx_file=addfile("uwnd.mon.mean.nc","r")
         
         uu=u_file->uwnd_Dtrend
         vv=v_file->vwnd_Dtrend
         pres=p_file->pres_Dtrend
         shum=shum_file->shum_Dtrend
       
         qu = new((/nyear,nlev,nlat,nlon/),float)
     qv = new((/nyear,nlev,nlat,nlon/),float)
       
         qu=shum*uu
         qv=shum*vv
         copy_VarCoords(uu,qu)
         copy_VarCoords(vv,qv)
         qu@_FillValue=-9.99E+33  
         qv@_FillValue=-9.99E+33  
       
         qu_new = new((/nyear,nlev+9,nlat,nlon/),float)
     qv_new = new((/nyear,nlev+9,nlat,nlon/),float)

         do i=0,7
            qu_new(:,i,:,:)= qu(:,i,:,:)
                qv_new(:,i,:,:)= qv(:,i,:,:)
         end do
         
         do j=8,nlev+9-1
           qu_new(:,j,:,:)= 0.0
           qv_new(:,j,:,:)= 0.0
         end do
         qu_new&level = xx_file->uwnd&level
         qv_new&level = xx_file->uwnd&level
        ; printVarSummary(qu_new)
         
         u = vibeta (p,qu_new(time|:,lat|:,lon|:,level|:),linlog,pres,pbot,ptop)/9.8  ; returns u(time,lat,lon)
         v = vibeta (p,qv_new(time|:,lat|:,lon|:,level|:),linlog,pres,pbot,ptop)/9.8  
       
         copy_VarCoords(uu(:,0,:,:),u)
         copy_VarCoords(vv(:,0,:,:),v)
         u@_FillValue=-9.99E+33  
         v@_FillValue=-9.99E+33  
       
        ; printVarSummary(v)
         asciiwrite ("allqu.txt" , u)
         delete(uu)
         delete(vv)
         delete(pres)
     delete(shum)

end

hxyj

奚子惠 发表于 2015-6-21 20:41
是的已经解决了
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/nc ...

感谢楼主的分享。再请教一个问题，不知你计算出的水汽通量即qu,qv的量级是多少？我计算出来的是10000，这是在夏季的。

奚子惠

hxyj 发表于 2015-6-21 20:48
感谢楼主的分享。再请教一个问题，不知你计算出的水汽通量即qu,qv的量级是多少？我计算出来的是10000，这 ...

qu/9.8大概在10的0次方，整层积分后的qu最大达到10的4次方的量级，我做的是春季，比湿g/kg

麦田_smile

奚子惠 发表于 2015-6-21 20:41
是的已经解决了
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/nc ...

楼主pres资料能否发一份？要到2014年12月的，现在noaa下不了了 谢谢啊

奚子惠

麦田_smile 发表于 2015-6-25 19:55
楼主pres资料能否发一份？要到2014年12月的，现在noaa下不了了 谢谢啊

sorry哦我只有2014年9月的的

麦田_smile

奚子惠 发表于 2015-6-26 09:21
sorry哦我只有2014年9月的的

不过还是多谢啦

麦田_smile

奚子惠 发表于 2015-6-21 20:41
是的已经解决了
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/nc ...

楼主，为什么用这个脚本还是提示你帖子里的错误啊

奚子惠

麦田_smile 发表于 2015-6-26 10:38
楼主，为什么用这个脚本还是提示你帖子里的错误啊

什么错误呢~

麦田_smile

奚子惠 发表于 2015-6-26 20:34
什么错误呢~

和你帖子里一开始一样的错误