
<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>Hello<br><br>Im using daily data (for evaporation) and would like to convert/calculate monthly values to plot my annual cycle.<br>The evaporation data is from 1979-2010 in daily format.<br>My script is:<br><br>; ==============================================================<br>; Open the file: Read only the user specified period first observations then model<br>; ==============================================================<br>f = addfile("ERAIL_SF_SRO_SSRO_TP_E_1979_2010.nc", "r") ;<br>;***********************************************************<br>; Read time and create required yyyyddd&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<br>;***********************************************************<br>&nbsp;&nbsp; time&nbsp;&nbsp;&nbsp; = f-&gt;time&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; time:units = "hours since 1-1-1 00:00:0.0"&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<br>&nbsp;&nbsp; TIME&nbsp;&nbsp;&nbsp; = cd_calendar(time, 0)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; type float <br>&nbsp;&nbsp; year&nbsp;&nbsp;&nbsp; = toint( TIME(:,0) )<br>&nbsp;&nbsp; month&nbsp;&nbsp; = toint( TIME(:,1) )<br>&nbsp;&nbsp; day&nbsp;&nbsp;&nbsp;&nbsp; = toint( TIME(:,2) ) <br>&nbsp;&nbsp; ddd&nbsp;&nbsp;&nbsp;&nbsp; = day_of_year(year, month, day) <br>&nbsp;&nbsp; yyyyddd = year*1000 + ddd&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; needed for input<br>;***********************************************************<br>; Read data: short2flt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<br>;*********************************************************** <br>&nbsp;&nbsp; e0&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; short2flt( f-&gt;e(:,:,:) )&nbsp;&nbsp;&nbsp; ; convert to float <br>&nbsp;&nbsp; printVarSummary( e0 ) <br>;***********************************************************<br>; Compute daily climatology: raw and then 'smoothed' &nbsp;<br>;***********************************************************<br>&nbsp;&nbsp; EClmDay&nbsp;&nbsp;&nbsp; = clmDayTLL(e0, yyyyddd)&nbsp;&nbsp;&nbsp;&nbsp; ; daily climatology at each grid point&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<br>&nbsp;&nbsp; printVarSummary(EClmDay)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; printVarSummary(EClmDay)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; (time, lat,lon)&nbsp; <br><br><br>and the output is:<br><br>Variable: e0<br>Type: float<br>Total Size: 233385984 bytes<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 58346496 values<br>Number of Dimensions: 3<br>Dimensions and sizes:&nbsp;&nbsp;&nbsp; [time | 11688] x [latitude | 64] x [longitude | 78]<br>Coordinates: <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; time: [692520..973008]<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: 7<br>&nbsp; units :&nbsp;&nbsp;&nbsp; m of water equivalent<br>&nbsp; long_name :&nbsp;&nbsp;&nbsp; Evaporation<br>&nbsp; standard_name :&nbsp;&nbsp;&nbsp; lwe_thickness_of_water_evaporation_amount<br>&nbsp; _FillValue_original :&nbsp;&nbsp;&nbsp; -32767<br>&nbsp; _FillValue :&nbsp;&nbsp;&nbsp; -32767<br>&nbsp; missing_value_original :&nbsp;&nbsp;&nbsp; -32767<br>&nbsp; missing_value :&nbsp;&nbsp;&nbsp; -32767<br><br>Variable: EClmDay<br>Type: float<br>Total Size: 7308288 bytes<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1827072 values<br>Number of Dimensions: 3<br>Dimensions and sizes:&nbsp;&nbsp;&nbsp; [year_day | 366] x [latitude | 64] x [longitude | 78]<br>Coordinates: <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; year_day: [1..366]<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: 4<br>&nbsp; long_name :&nbsp;&nbsp;&nbsp; Daily Climatology: Evaporation<br>&nbsp; units :&nbsp;&nbsp;&nbsp; m of water equivalent<br>&nbsp; information :&nbsp;&nbsp;&nbsp; Raw daily averages across all years<br>&nbsp; smoothing :&nbsp;&nbsp;&nbsp; None<br><br>Would really appreciate any advice<br></body></html>