[ncl-talk] rcm2rgrid issue

Dennis Shea shea at ucar.edu
Tue Jun 23 11:32:41 MDT 2020


*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/20200623/ceff3c7d/attachment.html>


More information about the ncl-talk mailing list