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

[Zilore Mumba]

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