[ncl-talk] rcm2rgrid issue
Dr. Suman Maity
suman.buie at gmail.com
Wed Jun 24 15:33:56 MDT 2020
Dear Rashed
I have another query. Will ESMF regrid functionworks only in 2D? I mean if
for 3D data like t(time,lat,lon) will it be applicable? If yes, what to do
especially?
Because it is mentioned in the website that the data to be of any
dimensionality with the rightmost dimension lat x lon, but I am getting
error like:
fatal:Eq: Dimension size, for dimension number 0, of operands does not
match, can't continue
Will you 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 Thu, Jun 25, 2020 at 2:30 AM Dr. Suman Maity <suman.buie at gmail.com>
wrote:
> Dear Rashed
> Yes, you are absolutely right.It works fine for me.
> 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 Thu, Jun 25, 2020 at 2:02 AM Rashed Mahmood <rashidcomsis at gmail.com>
> wrote:
>
>> ESMF regrid function is my personal choice when it comes to regriding:
>> See the attached script and the resulting plot. I turned off zoom
>> resources to see the full picture. Your "t" data seems to have missing
>> values which are masked before regriding. Please make sure that this looks
>> right to you.
>>
>> Cheers,
>> Rashed
>>
>> On Wed, Jun 24, 2020 at 6:28 PM Dr. Suman Maity <suman.buie at gmail.com>
>> wrote:
>>
>>> Hi Dennis
>>> As you directed I uploaded both the files "temp_seltime.nc" and "
>>> EIN5.soilT.0-7cm.JJASmean.1982-2016.remap.nc" at ftp.cgd.ucar.edu
>>> I would like to request you to have a look and comment where I am doing
>>> mistake.
>>> @Rashed
>>> As directed by Dennis, the data are now available at ftp.cgd.ucar.edu.
>>> Please have a look with the data and comment.
>>> Thanking you all.
>>>
>>> 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 Wed, Jun 24, 2020 at 8:27 PM Dennis Shea <shea at ucar.edu> wrote:
>>>
>>>> Perhaps, you could ftp the two netcdf file "temp_seltime.nc" and "
>>>> EIN5.soilT.0-7cm.JJASmean.1982-2016.remap.nc"
>>>>
>>>> ftp ftp.cgd.ucar.edu
>>>> anonymous
>>>> cd incoming
>>>> put temp_seltime.nc
>>>> put EIN5.soilT.0-7cm.JJASmean.1982-2016.remap.nc
>>>> quit
>>>>
>>>> If the EINS file is large, sInce NCL needs only rom the latitude and
>>>> longitude variables from the EINS file, you could use (say) the netCDF
>>>> 'ncks' operator
>>>>
>>>> %> ncks -v latitude,longitude
>>>> EIN5.soilT.0-7cm.JJASmean.1982-2016.remap.nc EINS_latlon.nc
>>>>
>>>> and ftp that file
>>>>
>>>>
>>>> On Tue, Jun 23, 2020 at 3:30 PM Rashed Mahmood <rashidcomsis at gmail.com>
>>>> wrote:
>>>>
>>>>> 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/20200625/7b3a0a7a/attachment.html>
More information about the ncl-talk
mailing list