[ncl-talk] Problem with linint2_Wrap
Ipsita Putatunda
ipsita.putatunda at gmail.com
Mon Aug 6 11:20:24 MDT 2018
Hi Toni,
Please find the attachment.
Thanks,
Ipsita
On Mon, Aug 6, 2018 at 10:43 PM, Toni Klemm <toni-klemm at tamu.edu> wrote:
> Ipsita,
>
> Could you run the command
>
> ncdump -h MERRA2_400.inst3_3d_asm_Np.20170101.SUB.nc
>
> and share the result of that? It should give you a summary of the header
> of your netCDF file, i.e. the metadata of the file and each variable, like
> latitude and longitude range, units, etc.
>
> Toni
>
>
>
> *Toni Klemm, Ph.D.*Postdoctoral Research Associate
> Department of Ecosystem Science and Management
> College of Agriculture and Life Sciences
> Texas A&M University, College Station, TX
> Contributor at the Early Career Climate Forum <http://www.eccforum.org>
> www.toni-klemm.de | @toniklemm <http://twitter.com/toniklemm>
>
>
>
>
>
>
>
> On Aug 6, 2018, at 12:05 PM, Ipsita Putatunda <ipsita.putatunda at gmail.com>
> wrote:
>
> Hi all,
> Attaching my script herewith, though can't send the data file through
> email(already sent in ftp) due to its huge volume.
>
> Thanks,
> Ipsita
>
> On Mon, Aug 6, 2018 at 10:33 PM, Ipsita Putatunda <
> ipsita.putatunda at gmail.com> wrote:
>
>> Hi Toni,
>> I have already sent my data and script through ftp, as the data file is
>> big I can't attach it in mail.
>> I first lonFlip my data from -180 to +180 to 0 to 360, till that point
>> the plot was coming fine.
>> After that I tried to regrid from 0.5 * 0.625 to 0.25*0.25, which is
>> giving this plot.
>> The point you mentioned for shifting 180 degree, I have noticed that
>> also, but somehow not able to fix it.
>>
>> Thanks,
>> Ipsita
>>
>>
>> On Mon, Aug 6, 2018 at 10:25 PM, Toni Klemm <toni-klemm at tamu.edu> wrote:
>>
>>> Ipsita,
>>>
>>> Would you mind sharing your NCL script and the metadata of your
>>> variables, i.e. an ncdump -h of your data files, with the listserv? I’m not
>>> sure many on the listserv have access to your data files and your script on
>>> the ftp server.
>>>
>>> From which to which resolution did you regrid? Maybe check if your
>>> latitude range before and after the regrinding is the same, for example
>>> from 0 to 359 degrees or from -180 to +179 degrees? Judging from your plot
>>> (the position of Africa, for example) there might be a 180 degree shift
>>> between your data and your map overlay.
>>>
>>> Toni
>>>
>>>
>>>
>>> *Toni Klemm, Ph.D.*Postdoctoral Research Associate
>>> Department of Ecosystem Science and Management
>>> College of Agriculture and Life Sciences
>>> Texas A&M University, College Station, TX
>>> Contributor at the Early Career Climate Forum
>>> <https://urldefense.proofpoint.com/v2/url?u=http-3A__www.eccforum.org&d=DwMFaQ&c=ODFT-G5SujMiGrKuoJJjVg&r=WsLJ5ywbW3XnTzqo_ETV62H6s1a5bcGgnrpl7xPRfic&m=XVNuJ4lF58XrFDJa3KuL2JhhOUYC9qhu5P0f2FbS4jg&s=pwUb1iORHQjLkcIoGsCduDCmKc1xCPHjvAGzpZk2S1E&e=>
>>> www.toni-klemm.de
>>> <https://urldefense.proofpoint.com/v2/url?u=http-3A__www.toni-2Dklemm.de&d=DwMFaQ&c=ODFT-G5SujMiGrKuoJJjVg&r=WsLJ5ywbW3XnTzqo_ETV62H6s1a5bcGgnrpl7xPRfic&m=XVNuJ4lF58XrFDJa3KuL2JhhOUYC9qhu5P0f2FbS4jg&s=AyO6EOxiMiJtuJjCG8E6gJxcyuY_DhvETISLN3Ik-Ns&e=>
>>> | @toniklemm
>>> <https://urldefense.proofpoint.com/v2/url?u=http-3A__twitter.com_toniklemm&d=DwMFaQ&c=ODFT-G5SujMiGrKuoJJjVg&r=WsLJ5ywbW3XnTzqo_ETV62H6s1a5bcGgnrpl7xPRfic&m=XVNuJ4lF58XrFDJa3KuL2JhhOUYC9qhu5P0f2FbS4jg&s=fCKqNGjGZ_eVYsnTBfGzcaC9dMplBdDod9UH7gzlW9Y&e=>
>>>
>>>
>>>
>>>
>>>
>>>
>>>
>>> On Aug 6, 2018, at 11:45 AM, Ipsita Putatunda <
>>> ipsita.putatunda at gmail.com> wrote:
>>>
>>> Hi Barry,
>>>
>>> There is actually no particular reason for not using other regridding
>>> options. As I have used this method several times in my previous scripts
>>> and am much familiar with this method, and also feel it is quicker process
>>> compare to other options, that is why I have chosen it. You are always
>>> welcome to suggest something different as well.
>>>
>>> Thanks,
>>> Ipsita
>>>
>>> On Mon, Aug 6, 2018 at 10:04 PM, Barry Lynn <barry.h.lynn at gmail.com>
>>> wrote:
>>>
>>>> Hi:
>>>>
>>>> Is there a particular reason you don't use
>>>>
>>>> https://www.ncl.ucar.edu/Applications/regrid.shtml
>>>> <https://urldefense.proofpoint.com/v2/url?u=https-3A__www.ncl.ucar.edu_Applications_regrid.shtml&d=DwMFaQ&c=ODFT-G5SujMiGrKuoJJjVg&r=WsLJ5ywbW3XnTzqo_ETV62H6s1a5bcGgnrpl7xPRfic&m=z_8XH5Nfko--z8kQcoLyR7TgUbSIqNJios95j_gBDhU&s=Rhqzff0XfY92jYqWoyOdJvMWFgdzGZP1L5Ngpclzknc&e=>
>>>> ?
>>>>
>>>> On Mon, Aug 6, 2018 at 7:20 PM Ipsita Putatunda <
>>>> ipsita.putatunda at gmail.com> wrote:
>>>>
>>>>> Dear all,
>>>>> I am not able to plot correctly after I regrid my data, using
>>>>> linint2_Wrap. I am attaching the plot herewith what I am getting.
>>>>> I have already sent my data file and script via ftp. Any help in this
>>>>> issue will be appreciated.
>>>>>
>>>>> Thanks in advance,
>>>>> Ipsita
>>>>> _______________________________________________
>>>>> ncl-talk mailing list
>>>>> ncl-talk at ucar.edu
>>>>> List instructions, subscriber options, unsubscribe:
>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>> <https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.ucar.edu_mailman_listinfo_ncl-2Dtalk&d=DwMFaQ&c=ODFT-G5SujMiGrKuoJJjVg&r=WsLJ5ywbW3XnTzqo_ETV62H6s1a5bcGgnrpl7xPRfic&m=z_8XH5Nfko--z8kQcoLyR7TgUbSIqNJios95j_gBDhU&s=0-9HbNHh-gfdMU8beZJdMGrGf7Mo8tSSWdBpUaEjMGo&e=>
>>>>>
>>>>
>>>>
>>>> --
>>>> 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
>>>> <https://urldefense.proofpoint.com/v2/url?u=http-3A__weather-2Dit-2Dis.com&d=DwMFaQ&c=ODFT-G5SujMiGrKuoJJjVg&r=WsLJ5ywbW3XnTzqo_ETV62H6s1a5bcGgnrpl7xPRfic&m=z_8XH5Nfko--z8kQcoLyR7TgUbSIqNJios95j_gBDhU&s=7H-awogmLMDx41uEHgNM5reEnkQstsgQ6tdGZAekd88&e=>
>>>> Jerusalem, Israel
>>>> Local: 02 930 9525
>>>> Cell: 054 7 231 170
>>>> Int-IS: x972 2 930 9525
>>>>
>>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>> <https://urldefense.proofpoint.com/v2/url?u=http-3A__mailman.ucar.edu_mailman_listinfo_ncl-2Dtalk&d=DwMFaQ&c=ODFT-G5SujMiGrKuoJJjVg&r=WsLJ5ywbW3XnTzqo_ETV62H6s1a5bcGgnrpl7xPRfic&m=XVNuJ4lF58XrFDJa3KuL2JhhOUYC9qhu5P0f2FbS4jg&s=S9k6flfF-_MuGpqVi3aMoYm1p6ZxWbqEFW99Mqmy_pM&e=>
>>>
>>>
>>>
>>
> <merra1.ncl>_______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180806/cb27495b/attachment.html>
-------------- next part --------------
netcdf MERRA2_400.inst3_3d_asm_Np.20170101.SUB {
dimensions:
time = UNLIMITED ; // (8 currently)
lev = 25 ;
lat = 121 ;
lon = 576 ;
variables:
float RH(time, lev, lat, lon) ;
RH:standard_name = "relative_humidity_after_moist" ;
RH:long_name = "relative_humidity_after_moist" ;
RH:units = "1" ;
RH:_FillValue = 1.e+15f ;
RH:missing_value = 1.e+15f ;
RH:fmissing_value = 1.e+15f ;
RH:vmax = 1.e+15f ;
RH:vmin = -1.e+15f ;
double lat(lat) ;
lat:standard_name = "latitude" ;
lat:long_name = "latitude" ;
lat:units = "degrees_north" ;
lat:axis = "Y" ;
double lev(lev) ;
lev:standard_name = "air_pressure" ;
lev:long_name = "vertical level" ;
lev:units = "hPa" ;
lev:positive = "down" ;
lev:axis = "Z" ;
double lon(lon) ;
lon:standard_name = "longitude" ;
lon:long_name = "longitude" ;
lon:units = "degrees_east" ;
lon:axis = "X" ;
double time(time) ;
time:standard_name = "time" ;
time:long_name = "time" ;
time:units = "minutes since 2017-01-01 00:00:00" ;
time:calendar = "standard" ;
// global attributes:
:CDI = "Climate Data Interface version 1.6.9 (http://mpimet.mpg.de/cdi)" ;
:history = "Thu Jul 26 05:43:51 2018: /usr/bin/ncks -O -L 1 --cnk_plc=g2d --cnk_dmn=lon,144 --cnk_dmn=lat,91 --cnk_dmn=lev,1 --cnk_dmn=time,1 /tmpdata/regridder/services_56125/cdoMERRA2_400.inst3_3d_asm_Np.20170101.SUB.nc4 /tmpdata/regridder/services_56125/chunkdeflatecdoMERRA2_400.inst3_3d_asm_Np.20170101.SUB.nc4\nThu Jul 26 05:43:10 2018: cdo -s -L -f nc4 -sellonlatbox,-180.0,179.375,-30.0,30.0 -select,name=RH,levidx=1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25 -seltimestep,1,2,3,4,5,6,7,8 /ftp/data/s4pa/MERRA2/M2I3NPASM.5.12.4/2017/01/MERRA2_400.inst3_3d_asm_Np.20170101.nc4 /tmpdata/regridder/services_56125/MERRA2_400.inst3_3d_asm_Np.20170101.SUB.nc4" ;
:Conventions = "CF-1" ;
:History = "Original file generated: Fri Jan 13 16:47:50 2017 GMT" ;
:Comment = "GMAO filename: d5124_m2_jan10.inst3_3d_asm_Np.20170101.nc4" ;
:Filename = "MERRA2_400.inst3_3d_asm_Np.20170101.nc4" ;
:Institution = "NASA Global Modeling and Assimilation Office" ;
:References = "http://gmao.gsfc.nasa.gov" ;
:Format = "NetCDF-4/HDF-5" ;
:SpatialCoverage = "global" ;
:VersionID = "5.12.4" ;
:TemporalRange = "1980-01-01 -> 2016-12-31" ;
:identifier_product_doi_authority = "http://dx.doi.org/" ;
:ShortName = "M2I3NPASM" ;
:GranuleID = "MERRA2_400.inst3_3d_asm_Np.20170101.nc4" ;
:ProductionDateTime = "Original file generated: Fri Jan 13 16:47:50 2017 GMT" ;
:LongName = "MERRA2 inst3_3d_asm_Np: 3d,3-Hourly,Instantaneous,Pressure-Level,Assimilation,Assimilated Meteorological Fields" ;
:Title = "MERRA2 inst3_3d_asm_Np: 3d,3-Hourly,Instantaneous,Pressure-Level,Assimilation,Assimilated Meteorological Fields" ;
:SouthernmostLatitude = "-90.0" ;
:NorthernmostLatitude = "90.0" ;
:WesternmostLongitude = "-180.0" ;
:EasternmostLongitude = "179.375" ;
:LatitudeResolution = "0.5" ;
:LongitudeResolution = "0.625" ;
:DataResolution = "0.5 x 0.625 (42 pressure levels)" ;
:Contact = "http://gmao.gsfc.nasa.gov" ;
:identifier_product_doi = "10.5067/QBZ6MG944HW0" ;
:RangeBeginningDate = "2017-01-01" ;
:RangeBeginningTime = "00:00:00.000000" ;
:RangeEndingDate = "2017-01-01" ;
:RangeEndingTime = "21:00:00.000000" ;
:CDO = "Climate Data Operators version 1.6.9 (http://mpimet.mpg.de/cdo)" ;
:NCO = "20180726" ;
}
More information about the ncl-talk
mailing list