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

[Gus Correa]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20201101/45d106d2/attachment.html>