[ncl-talk] rcm2rgrid issue

Dr. Suman Maity suman.buie at gmail.com
Tue Jun 23 13:18:03 MDT 2020 

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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200624/a30b448a/attachment.html>

More information about the ncl-talk mailing list