<html><head><meta http-equiv="Content-Type" content="text/html; charset=UTF-8"><title></title><style type="text/css">.felamimail-body-blockquote {margin: 5px 10px 0 3px;padding-left: 10px;border-left: 2px solid #000088;} </style></head><body>THank you.<br><br>The cdo function calculates the monthly values.<br>&nbsp;yrStrt = 1970<br>&nbsp; yrLast = 2000<br><br>; ==============================================================<br>; Open the file: Read only the user specified period first observations then model<br>; ==============================================================<br>f = addfile("ERAIOBS_monthly.nc", "r") ;<br>TIME&nbsp;&nbsp; = f-&gt;time<br>YYYY&nbsp;&nbsp; = cd_calendar(TIME,-1)/100&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; entire file<br>iYYYY&nbsp; = ind(YYYY.ge.yrStrt .and. YYYY.le.yrLast)<br>e&nbsp;&nbsp; = f-&gt;e(iYYYY,:,:)<br>printVarSummary(e)&nbsp; <br><br>SO the e output is:<br><br>Type: short<br>Total Size: 2635776 bytes<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1317888 values<br>Number of Dimensions: 3<br>Dimensions and sizes:&nbsp;&nbsp;&nbsp; [time | 264] x [latitude | 64] x [longitude | 78]<br>Coordinates: <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; time: [693216..885336]<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; latitude: [64..32.5]<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; longitude: [30..68.5]<br>Number Of Attributes: 6<br>&nbsp; long_name :&nbsp;&nbsp;&nbsp; Evaporation<br>&nbsp; standard_name :&nbsp;&nbsp;&nbsp; lwe_thickness_of_water_evaporation_amount<br>&nbsp; units :&nbsp;&nbsp;&nbsp; m of water equivalent<br>&nbsp; add_offset :&nbsp;&nbsp;&nbsp; -0.01280028930916859<br>&nbsp; scale_factor :&nbsp;&nbsp;&nbsp; 4.884610965422988e-07<br>&nbsp; _FillValue :&nbsp;&nbsp;&nbsp; -32767<br><br><br>BUT: <br><br>;===================================================<br>; compute desired global annual or monthly climatology (12,nlat,mlon)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>; ==============================================================<br>&nbsp; Tann=clmMonTLL(e) <br>&nbsp;&nbsp; printVarSummary(Tann) <br><br>The output is:<br><br>fatal:Type Mismatch: The type of missing value could not be converted to type of variable (e)<br>fatal:Assignment type mismatch, right hand side can't be coerced to type of left hand side<br><br><br>Is it wrong to do clmMonTLL on the monthly data output from cdo?<br><br></body></html>