[ncl-talk] Contour map in ncl 6.6. 2
Zilore Mumba
zmumba at gmail.com
Mon Nov 2 16:08:54 MST 2020
I wonder what happened to climo_0.ncl. The link on the page gives the
message " The requested URL was not found on this server." I would
appreciate if anyone has a copy to share.
On Sun, Nov 1, 2020 at 10:48 AM Zilore Mumba <zmumba at gmail.com> wrote:
> Thank you very much indeed for this very useful information.
> I appreciate your time and effort.
>
> On Sun, Nov 1, 2020 at 10:39 AM Gus Correa via ncl-talk <
> ncl-talk at mailman.ucar.edu> wrote:
>
>> Your input file name suggests that it is daily data.
>> However, clmMonTLL expects an input file with monthly data,
>> and a time dimension (number of months) that is a multiple of 12:
>> https://www.ncl.ucar.edu/Document/Functions/Contributed/clmMonTLL.shtml
>>
>> You can try first to calculate monthly values from the daily input using
>> calculate_monthly_values:
>>
>> https://www.ncl.ucar.edu/Document/Functions/Contributed/calculate_monthly_values.shtml
>> But you need to make sure that your input (daily data) has a time
>> dimension with a recognized
>> calendar (I guess it would be gregorian).
>> A "ncdump -h " of your input file should tell.
>>
>>
>> On Sun, Nov 1, 2020 at 3:19 AM Zilore Mumba via ncl-talk <
>> ncl-talk at mailman.ucar.edu> wrote:
>>
>>> I hope I can get assistance from Rick or Dennis or indeed anyone who can
>>> help. I am working with climate simulations using cordex data (20+ regional
>>> models). I have files of the 30 year period (2021-2050), though I tried
>>> climo_0.ncl <https://www.ncl.ucar.edu/Applications/Scripts/climo_0.ncl> with
>>> only one model simulation.
>>> When I run the script, I get "*(0) contributed.ncl: clmMonTLL:
>>> dimension must be a multiple of 12*"
>>>
>>> *The following are modifications I made to the original script:*
>>> a =
>>> addfile("../../Data/rcp45/precip/pr_AFR-44_NOAA-GFDL-GFDL-ESM2M_rcp45_r1i1p1_SMHI-RCA4_v1_day_20210101_20501231.nc","r")
>>> prc = a->pr ; (time,lat,lon)
>>> time = a->time ; (time) ==> YYMM
>>> ntim = dimsizes(time)
>>>
>>> ;---Read lat/lon arrays and attach as special "lat2d" / "lon2d"
>>> attributes
>>> lat2d = a->rlat ; (y,x)
>>> lon2d = a->rlon ; (y,x)
>>>
>>> prc&rlat at units = "degrees_north"
>>> prc&rlon at units = "degrees_east"
>>>
>>> months = (/"January", "February", "March", "April" \
>>> ,"May", "June", "July", "August" \
>>> ,"September", "October", "November" \
>>> ,"December" /)
>>>
>>>
>>> printVarSummary(prc) ; note "coordinates" attribute
>>> printVarSummary(lat2d)
>>> printVarSummary(lon2d)
>>> ;*************************************************************
>>>
>>> *I also used (see last line below):*
>>>
>>> ; Compute the climatology using a function in contributed.ncl
>>> ;************************************************
>>> prcClm = clmMonTLL (prc) ; monthly climatology
>>> prcStd = stdMonTLL (prc) ; monthly interannual variability
>>> ;************************************************
>>>
>>> ;To get the climatology for Jan 2021 through Dec 2030 for this dataset:
>>> ;prcClm = clmMonTLL (prc(12:131,:,:)), using conventional subscripts.
>>> prcClm = clmMonTLL (prc({202101:203012},:,:)) ;using coordinate
>>> subscripting.
>>>
>>> ;*******************************************************************************
>>>
>>> I will appreciate any help.
>>> Zilore
>>>
>>>
>>> On Sat, Oct 31, 2020 at 2:50 AM Md. Jalal Uddin via ncl-talk <
>>> ncl-talk at mailman.ucar.edu> wrote:
>>>
>>>> Hi all,
>>>>
>>>> Why am I not getting the contour colour in the map in ncl 6.6.2?
>>>>
>>>> Any suggestions would be appreciated?
>>>>
>>>> Best wishes,
>>>> Jalal
>>>> --
>>>> *Md. Jalal Uddin*
>>>> Founder and Director of Research Society (
>>>> https://researchsociety20.org/founder-and-director/)
>>>> PhD candidate - Atmospheric Physics and Atmospheric Environment (NUIST)
>>>> M.Sc. in Applied Meteorology (NUIST)
>>>> B.Sc. in Disaster Management (PSTU)
>>>> Email: dmjalal90 at nuist.edu.cn
>>>> Personal website: https://dmjalal90.weebly.com/
>>>> Research gate profile: https://www.researchgate.net/profile/Md_Uddin125
>>>>
>>>> YouTube channel: http://www.youtube.com/c/MdJalalUddindmjalal90
>>>> Cell: +8613260859092, +8801792052662
>>>> _______________________________________________
>>>> 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
>>
>> _______________________________________________
>> 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/20201103/93ad3ef2/attachment.html>
More information about the ncl-talk
mailing list