[ncl-talk] Contour map in ncl 6.6. 2

[Dennis Shea]

[0] The *climo_0.ncl* script is attached

[2] Gus: THX for answering

[3] Zilore:  (a) Please DO NOT put a new question under a previous ncl-talk
question.
                 (b) Create a new 'thread' with a Subject

On Mon, Nov 2, 2020 at 4:09 PM Zilore Mumba via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:

> 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
>>
>> _______________________________________________
> 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/20201102/fef7a3ca/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: climo_0.ncl
Type: application/octet-stream
Size: 2602 bytes
Desc: not available
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20201102/fef7a3ca/attachment.obj>