[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