[ncl-talk] rcm2rgrid issue
Rashed Mahmood
rashidcomsis at gmail.com
Wed Jun 24 14:32:36 MDT 2020
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/20200624/e6acd3d0/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Suman_test.ncl
Type: application/octet-stream
Size: 1866 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200624/e6acd3d0/attachment-0001.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: plot.pdf
Type: application/pdf
Size: 557914 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200624/e6acd3d0/attachment-0001.pdf>
More information about the ncl-talk
mailing list