[ncl-talk] a bug about geopt in WRFUserARW.ncl
Jiang Ziqiang
zqjiang86 at hotmail.com
Tue Jun 7 09:50:12 MDT 2016
Hi, ncl-talk
The time in four position should be replace by time_in.
==>> time_in <<==
if( any( variable .eq. (/"geopt","geopotential","z","height"/) ) ) then
; Height [=full geopotentail height / 9.81]
if(isfilevar(nc_file,"PH")) then
if ( time .eq. -1 ) then
if(ISFILE) then
var = nc_file->PH
PHB = nc_file->PHB
else
var = file_handle[:]->PH
PHB = file_handle[:]->PHB
end if
else
if(ISFILE) then
var = nc_file->PH(time,:,:,:) ==>> time_in <<==
PHB = nc_file->PHB(time,:,:,:) ==>> time_in <<==
else
var = file_handle[:]->PH(time,:,:,:) ==>> time_in <<==
PHB = file_handle[:]->PHB(time,:,:,:) ==>> time_in <<==
end if
end if
var = var + PHB
z = wrf_user_unstagger(var,var at stagger)
z at description = "Geopotential"
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160607/947c091e/attachment.html
More information about the ncl-talk
mailing list