[ncl-talk] a bug about geopt in WRFUserARW.ncl

[Bill Ladwig]

Thank you for the bug report.  This issue has already been addressed in the
upcoming 6.4.0 release of NCL.

Regards,

Bill

On Tue, Jun 7, 2016 at 9:50 AM, Jiang Ziqiang <zqjiang86 at hotmail.com> wrote:

> 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"
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160607/15de39a3/attachment.html