[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