[ncl-talk] rcm2rgrid issue

Rashed Mahmood rashidcomsis at gmail.com
Tue Jun 23 15:29:57 MDT 2020 

did you check what the min/max of your data values are?

printMinMax(xgrid,0)

What does printVarSummary(xgrid) tell?

Without data file(s) there is not much people can do to help. You could
also try using ESMF regridding:
http://www.ncl.ucar.edu/Applications/ESMF.shtml




On Tue, Jun 23, 2020 at 9:18 PM Dr. Suman Maity via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> Hi Dennis
> Thank you for your suggestions. As you mentioned, I modified my script as
> beow:
> begin
>   a = addfile("EIN5.soilT.0-7cm.JJASmean.1982-2016.remap.nc","r")
>   soil = a->stl1
>   soil = soil(:,::-1,:)
>   b = addfile("temp_seltime.nc","r")
>   t = rm_single_dims(b->TSOI(:,0,:,:))
>   lat2d = b->lat
>   lon2d = b->lon
>   xgrd  = rcm2rgrid_Wrap(lat2d,lon2d,t,soil&latitude,soil&longitude,0)
>
>   wks  = gsn_open_wks ("x11", "plot")
>   sres                      = True
>   sres at gsnAddCyclic         = False
>   sres at mpMinLatF            = 5
>   sres at mpMaxLatF            = 40
>   sres at mpMinLonF            = 65
>   sres at mpMaxLonF            = 100
>   sres at cnFillOn = True
>   sres at mpFillOn = False
>   sres at cnLevelSelectionMode = "ManualLevels"
>   sres at cnMinLevelValF = 280.
>   sres at cnMaxLevelValF = 320.
>   sres at cnLevelSpacingF = 10
>   plot = gsn_csm_contour_map(wks,xgrd,sres)
> end
> But unfortunately, all the efforts goes to vain and the plot remains same
> as observed earlier. Is there any other issue? Please comment.
> Thanking you,
>
> Best
> Suman
> ++++++++++++++++++++++++++++++++++++++++++
> Dr. Suman Maity
> Research Associate
> School of Atmospheric Science
> Sun Yat-sen University, Zhuhai campus,Tangjiawan Town
> Zhuhai, Guangdong, Postcode-519082,China.
> Phone:+86-13543861045 (China), +91-9732636778 (India)
> skype: reach2suman
> https://www.researchgate.net/profile/Suman_Maity
> ++++++++++++++++++++++++++++++++++++++++++
>
>
> On Tue, Jun 23, 2020 at 11:02 PM Dennis Shea <shea at ucar.edu> wrote:
>
>>
>> *http://www.ncl.ucar.edu/Document/Functions/Contributed/rcm2rgrid_Wrap.shtml*
>> <http://www.ncl.ucar.edu/Document/Functions/Contributed/rcm2rgrid_Wrap.shtml>
>>
>> and *rcm2rgrid*
>> <http://www.ncl.ucar.edu/Document/Functions/Built-in/rcm2rgrid.shtml>
>> require that the destination (rectilinear) grid have the coordinates
>> (lat[*] and lon[*])  in *ascending order*.
>>
>> You rectilinear grid looks like it is in descending order:  *latitude:
>> [49.5..-19.5]*
>>
>> *f = addfile("EIN5.soilT.0-7cm.JJASmean.1982-2016.remap.nc
>> <http://EIN5.soilT.0-7cm.JJASmean.1982-2016.remap.nc>","r")*
>> *x = f->SOIL*
>> *printVarSummary(x)*
>>
>> *print("========")*
>>
>>
>> *x = x(::-1,:)    *
>> *printVarSummary(x)*
>> *print("========")*
>>
>>   xgrd  = *rcm2rgrid*(lat2d,lon2d,*x,x&lat,x&lon*,0)
>>
>>
>> On Tue, Jun 23, 2020 at 9:32 AM Dr. Suman Maity via ncl-talk <
>> ncl-talk at mailman.ucar.edu> wrote:
>>
>>> Hi Barry
>>> Thanks for your prompt reply. As you mentioned, here are the details:
>>>  See http://www.ncl.ucar.edu/ for more details.
>>> (0) ============Before Interpolation===========
>>> Variable: t
>>> Type: float
>>> Total Size: 265960 bytes
>>>             66490 values
>>> Number of Dimensions: 2
>>> Dimensions and sizes: [218] x [305]
>>> Coordinates:
>>> Number Of Attributes: 6
>>>   _FillValue : 1e+36
>>>   time : 287808
>>>   levgrnd : 0.007100635417193535
>>>   long_name : soil temperature (vegetated landunits only)
>>>   units : K
>>>   cell_method : time: mean
>>> (0) ============After Interpolation===========
>>> Variable: xgrd
>>> Type: float
>>> Total Size: 52080 bytes
>>>             13020 values
>>> Number of Dimensions: 2
>>> Dimensions and sizes: [latitude | 93] x [longitude | 140]
>>> Coordinates:
>>>             latitude: [49.5..-19.5]
>>>             longitude: [20.25..124.5]
>>> Number Of Attributes: 7
>>>   _FillValue : 1e+36
>>>   time : 287808
>>>   levgrnd : 0.007100635417193535
>>>   long_name : soil temperature (vegetated landunits only)
>>>   units : K
>>>   cell_method : time: mean
>>>   ncl : rcm2rgrid used for interpolation
>>> [suman at localhost ght]$
>>> Please reply if you need further information. I am unable to share the
>>> data as ncl community allow only 2 MB.
>>> Thanking you once again.
>>>
>>> Best
>>> Suman
>>> ++++++++++++++++++++++++++++++++++++++++++
>>> Dr. Suman Maity
>>> Research Associate
>>> School of Atmospheric Science
>>> Sun Yat-sen University, Zhuhai campus,Tangjiawan Town
>>> Zhuhai, Guangdong, Postcode-519082,China.
>>> Phone:+86-13543861045 (China), +91-9732636778 (India)
>>> skype: reach2suman
>>> https://www.researchgate.net/profile/Suman_Maity
>>> ++++++++++++++++++++++++++++++++++++++++++
>>>
>>>
>>> On Tue, Jun 23, 2020 at 8:48 PM Barry Lynn <barry.h.lynn at gmail.com>
>>> wrote:
>>>
>>>> Hello Dr. Suman Maity:
>>>>
>>>> It would be helpful if you could attach printouts from your
>>>> interpolation, as well as printVarSummary output of variables you are using
>>>> and/or calculating during the interpolation.
>>>>
>>>> Barry
>>>>
>>>> On Tue, Jun 23, 2020 at 5:58 PM Dr. Suman Maity via ncl-talk <
>>>> ncl-talk at mailman.ucar.edu> wrote:
>>>>
>>>>> Dear All
>>>>> I am facing an issue from "rcm2rgrid" which is mentioned below:
>>>>> I have two netcdf file "temp_seltime.nc" and "
>>>>> EIN5.soilT.0-7cm.JJASmean.1982-2016.remap.nc" and I want to
>>>>> interpolate "temp_seltime.nc" to the lon/lat of "
>>>>> EIN5.soilT.0-7cm.JJASmean.1982-2016.remap.nc" for evaluation. Now "
>>>>> temp_seltime.nc" is an curviliear data with lat and lon information
>>>>> in 2D and the lat2D and lon2D contains some missing points. On the
>>>>> otherhand, "EIN5.soilT.0-7cm.JJASmean.1982-2016.remap.nc" is a normal
>>>>> netdf file (rectilinear grid).  Although interpolation doesn't through any
>>>>> error but the resultant plot is blank. I couldn't understand where I did
>>>>> the mistake? Please scripts and plots as attached.
>>>>> Any sort of suggestion/comments are highly appreciable.
>>>>> Best
>>>>> Suman
>>>>> ++++++++++++++++++++++++++++++++++++++++++
>>>>> Dr. Suman Maity
>>>>> Research Associate
>>>>> School of Atmospheric Science
>>>>> Sun Yat-sen University, Zhuhai campus,Tangjiawan Town
>>>>> Zhuhai, Guangdong, Postcode-519082,China.
>>>>> Phone:+86-13543861045 (China), +91-9732636778 (India)
>>>>> skype: reach2suman
>>>>> https://www.researchgate.net/profile/Suman_Maity
>>>>> ++++++++++++++++++++++++++++++++++++++++++
>>>>> _______________________________________________
>>>>> ncl-talk mailing list
>>>>> ncl-talk at mailman.ucar.edu
>>>>> List instructions, subscriber options, unsubscribe:
>>>>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>>>
>>>>
>>>> --
>>>> Barry H. Lynn, Ph.D
>>>> Senior Associate Scientist, Lecturer,
>>>> The Institute of the Earth Science,
>>>> The Hebrew University of Jerusalem,
>>>> Givat Ram, Jerusalem 91904, Israel
>>>> Tel: 972 547 231 170
>>>> Fax: (972)-25662581
>>>>
>>>> C.E.O, Weather It Is, LTD
>>>> Weather and Climate Focus
>>>> http://weather-it-is.com
>>>> Jerusalem, Israel
>>>> Local: 02 930 9525
>>>> Cell: 054 7 231 170
>>>> Int-IS: x972 2 930 9525
>>>>
>>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at mailman.ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200623/26f89ae3/attachment.html>

More information about the ncl-talk mailing list