[ncl-talk] Problem with linint2_Wrap

Toni Klemm toni-klemm at tamu.edu
Mon Aug 6 12:34:58 MDT 2018 

Ah, I see. You’re welcome! I’m glad you found it.

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 <http://www.toni-klemm.de/> | @toniklemm <http://twitter.com/toniklemm>







> On Aug 6, 2018, at 1:00 PM, Ipsita Putatunda <ipsita.putatunda at gmail.com> wrote:
> 
> Hi Toni,
>  I was using "r" while doing linint2_Wrap instead of "rnew", that is why it was plotting differently.
> 
> Thanks a lot.
> Ipsita
> 
> On Mon, Aug 6, 2018 at 11:22 PM, Ipsita Putatunda <ipsita.putatunda at gmail.com <mailto:ipsita.putatunda at gmail.com>> wrote:
> Hi Toni,
>  If you look in my script, in line 20 I have already lonFlip the data longitude. lonFlip is a function which flips the longitude from -180/+180 to 0/360 or vice versa. So though my initial data lon is in -180/+180 format, I have changed it to 0/360 via lonFlip and there after creating the variable for regridding.
> a1=fspan(0.,359.75,1441)
> Thanks,
> Ipsita
> 
> On Mon, Aug 6, 2018 at 11:06 PM, Toni Klemm <toni-klemm at tamu.edu <mailto:toni-klemm at tamu.edu>> wrote:
> Ipsita,
> 
> It looks indeed like your longitudes are shifted. In your script you created an array that spans 0 to 359
> 
> 	a1=fspan(0.,359.75,1441)
> 
> but the metadata of your file lists as the longitudes of your data ranging from -180 to 179.375 degrees.
> 
> 	:WesternmostLongitude = "-180.0" ;
> 	:EasternmostLongitude = "179.375" ;
> 
> 
> I’m not familiar with the function lining_Wrap, but try setting up the array a1 so it ranges from -180 to +179.375 and see what happens.
> 
> 	a1=fspan(-180,179.375,1441)
> 
> 
> Also, at first I wondered why the plot only covers the region around the equator (30 S to 30 N). But from your script I assume you only want to plot that region.
> 
> 	b1=fspan(-30.,30.,241)
> 
> If not, that try changing the fspan from -90 to 90 and the steps from 241 to 723 (3x 241).
> 
> 
> 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=hC7Jy6J4NpNHflaHzXEiWxJBsu50UiPokV6f3h8mg18&s=KWMDmP_MUoEz6emOPruGgN7E1YLxiGL7BGdQi6KNY7Q&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=hC7Jy6J4NpNHflaHzXEiWxJBsu50UiPokV6f3h8mg18&s=dFvrD9TVPAcV2MZ-qe8fhEIfl9tMagRhNwdbeJeDHcU&e=> | @toniklemm <https://urldefense.proofpoint.com/v2/url?u=http-3A__twitter.com_toniklemm&d=DwMFaQ&c=ODFT-G5SujMiGrKuoJJjVg&r=WsLJ5ywbW3XnTzqo_ETV62H6s1a5bcGgnrpl7xPRfic&m=hC7Jy6J4NpNHflaHzXEiWxJBsu50UiPokV6f3h8mg18&s=xmgrg4QQp1XCAcid9QdF38s4ifaVxad2DFNcWyQdNf0&e=>
> 
> 
> 
> 
> 
> 
> 
>> On Aug 6, 2018, at 12:20 PM, Ipsita Putatunda <ipsita.putatunda at gmail.com <mailto:ipsita.putatunda at gmail.com>> wrote:
>> 
>> Hi Toni,
>>  Please find the attachment.
>> 
>> Thanks,
>> Ipsita
>> 
>> On Mon, Aug 6, 2018 at 10:43 PM, Toni Klemm <toni-klemm at tamu.edu <mailto:toni-klemm at tamu.edu>> wrote:
>> Ipsita,
>> 
>> Could you run the command 
>> 
>> ncdump -h MERRA2_400.inst3_3d_asm_Np.20170101.SUB.nc <https://urldefense.proofpoint.com/v2/url?u=http-3A__MERRA2-5F400.inst3-5F3d-5Fasm-5FNp.20170101.SUB.nc&d=DwMFaQ&c=ODFT-G5SujMiGrKuoJJjVg&r=WsLJ5ywbW3XnTzqo_ETV62H6s1a5bcGgnrpl7xPRfic&m=SLg59Bo2un5Gc1gbYkbJx-YvAwqss_Dek2CXQ2STNHE&s=dxd4uHhDSIhFa7kqqt2VPyDxzvaTTJHrtiFXR-JfcmQ&e=>
>> 
>> 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 <https://urldefense.proofpoint.com/v2/url?u=http-3A__www.eccforum.org&d=DwMFaQ&c=ODFT-G5SujMiGrKuoJJjVg&r=WsLJ5ywbW3XnTzqo_ETV62H6s1a5bcGgnrpl7xPRfic&m=SLg59Bo2un5Gc1gbYkbJx-YvAwqss_Dek2CXQ2STNHE&s=484HpwSun6C0ZiMj1-0Kg5OgKfF5MPH9TD5mFbp-Neg&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=SLg59Bo2un5Gc1gbYkbJx-YvAwqss_Dek2CXQ2STNHE&s=HkRbQg4yW_BcdF8HHwfqByYj8e4_s3t6-apgDERmspQ&e=> | @toniklemm <https://urldefense.proofpoint.com/v2/url?u=http-3A__twitter.com_toniklemm&d=DwMFaQ&c=ODFT-G5SujMiGrKuoJJjVg&r=WsLJ5ywbW3XnTzqo_ETV62H6s1a5bcGgnrpl7xPRfic&m=SLg59Bo2un5Gc1gbYkbJx-YvAwqss_Dek2CXQ2STNHE&s=Lz6FlJH03nFr_L2Yqb2bVwWMICIpwxF41tbkRjASn5Y&e=>
>> 
>> 
>> 
>> 
>> 
>> 
>> 
>>> On Aug 6, 2018, at 12:05 PM, Ipsita Putatunda <ipsita.putatunda at gmail.com <mailto: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 <mailto: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 <mailto: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 <mailto: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 <mailto: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 <mailto: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 <mailto: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 <mailto: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 <mailto: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=SLg59Bo2un5Gc1gbYkbJx-YvAwqss_Dek2CXQ2STNHE&s=nNYd_0WKTngaTKj5RQIssYSpNFpOuKsISGVrzxYuUfk&e=>
>> 
>> 
>> <ncdump1.txt>
> 
> 
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180806/4724649e/attachment.html>

More information about the ncl-talk mailing list