[ncl-talk] WRF interpolation (amendment-output)

[Ehsan Taghizadeh]

  td2  = wrf_user_getvar(fw,"td2",-1) "***works***"and td2=fw->td2(0,:,:) "***didn't work***"likely, "td2" is diagnostic and is not in the file. It should be computed via "wrf_user_getvar".My problem is why stations out of WRF domain, got interpolated values?!!
SincerelyEhsan    On Wednesday, September 11, 2019, 07:24:50 PM GMT+4:30, Dennis Shea <shea at ucar.edu> wrote:  
 
  ---
  td2  = wrf_user_getvar(fw,"td2",-1) ; td2=fw->td2(0,:,:) ******didn't work****  printVarSummary(td2)
  printMinMax(td2,0)---Likely: 

    [1] td2' is not  on the file 
    [2] one of the variables needed to derive 'td2' is not present on the file.    [3] you have an old version of 'wrf_user_getvar'
%> less WRFUserARW.ncl       <=== 6.5.0

  if( variable .eq. "td2" ) then
       ;; function wrf_td needs qv and pressure (Pa) on input and returns dewpoint temperature on return
       PSFC = _get_wrf_var(file_handle,"PSFC",time)
       Q2   = _get_wrf_var(file_handle,"Q2",time)
       td   = wrf_td( PSFC , Q2 )
       copy_VarAtts_except(Q2,td,(/"description","units"/))
       td at description = "2m Dewpoint Temperature"   ; Overwrite return description
       return(td)
  end if

On Wed, Sep 11, 2019 at 6:22 AM Ehsan Taghizadeh via ncl-talk <ncl-talk at ucar.edu> wrote:

Dears,I forgot to attach output file from interpolation.At the attachment, stations like "40899" or "40890" got "9969209968386869046778552952102584320.0" value. These stations don't locate in wrf domain 2. However, stations like "40745" has values, in spite of locating out of domain 2.
SincerelyEhsan_______________________________________________
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/20190911/49586a02/attachment.html>