[ncl-talk] calcMonAnomTLL

Setareh Rahimi setareh.rahimi at gmail.com
Thu May 27 01:23:59 MDT 2021


Dear all NCL users,

I visited this link:
https://www.ncl.ucar.edu/Document/Functions/Contributed/calcMonAnomTLL.shtml ,
which is about monthly anomalies. But I could not understand the way it
calculates the anomalies.

_____________________________________________________________________________________________

Example 1


  fa     = *addfile*
<https://www.ncl.ucar.edu/Document/Functions/Built-in/addfile.shtml>("air.mean.nc","r")
  fb     = *addfile*
<https://www.ncl.ucar.edu/Document/Functions/Built-in/addfile.shtml>("air.nobs.nc","r")

  air  = *short2flt*
<https://www.ncl.ucar.edu/Document/Functions/Contributed/short2flt.shtml>(
fa->air)                           ; air temp
  nobs = *short2flt*
<https://www.ncl.ucar.edu/Document/Functions/Contributed/short2flt.shtml>(
fb->air)                           ; # obs for air
; filter out geographic locations that have < nMin observations
  nMin   = 2                                           ; min # obs
  air    = *mask*
<https://www.ncl.ucar.edu/Document/Functions/Built-in/mask.shtml>(air,nobs.ge.nMin,True)

;---Compute the climatology using a function in contributed.ncl
  yrStrt = 1950
  yrLast = 1979
  moStrt = (yrStrt-1800)*12              ; start subscript
  moLast = (yrLast-1800)*12 + 11         ; last  subscript
  clm = *clmMonTLL*
<https://www.ncl.ucar.edu/Document/Functions/Contributed/clmMonTLL.shtml>
(air(moStrt:moLast,:,:)) ; monthly climatology

;---Compute the anomalies from the climatology just computed
  xAnom = *calcMonAnomTLL* (air,clm)

_________________________________________________________________

May I ask you please explain subscript method? Why using 1800 and then
multiply 12?

what do moStrt and moLast show?

I appriciate any advse

All the best,


-- 
S.Rahimi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20210527/e199f70b/attachment.html>


More information about the ncl-talk mailing list