[ncl-talk] rcm2rgrid issue
Dr. Suman Maity
suman.buie at gmail.com
Fri Jun 26 22:58:07 MDT 2020
Hi Dennis
Extremely sorry for late reply. I would like to thank you for your replies.
It worked fine as you suggested.
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 5:00 AM Dennis Shea <shea at ucar.edu> wrote:
> *https://www.ncl.ucar.edu/Applications/ESMF.shtml*
> <https://www.ncl.ucar.edu/Applications/ESMF.shtml>
>
> ESMF regridding works with 1D, 2D (lat,lon), 3D (time,lat,lon), 4D
> (time,level,lat,lon)
>
> b = addfile("temp_seltime.nc","r") ; source curvilinear grid
> ;;t = b->TSOI ; 4D: (time, levgrnd, y, x)
> ;;t = b->TSOI(:,0,:,:) ; 3D: (time, y, x)
> t = b->TSOI(0,0,:,:) ; 2D: (y, x)
> printVarSummary(t)
> printMinMax(t,0)
> print("=====")
>
> The error message you mention tells us nothing. There is no context.
>
> fatal:Eq: Dimension size, for dimension number 0, of operands does not
> match, can't continue
> Will you please comment?
>
> On Wed, Jun 24, 2020 at 3:34 PM Dr. Suman Maity <suman.buie at gmail.com>
> wrote:
>
>> 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/20200627/11768021/attachment.html>
More information about the ncl-talk
mailing list