[ncl-talk] Calculate Annual cycle from daily data

Sri Nandini snandini at marum.de
Tue Jun 13 07:35:08 MDT 2017


Hello

Im using daily data (for evaporation) and would like to convert/calculate monthly values to plot my annual cycle.
The evaporation data is from 1979-2010 in daily format.
My script is:

; ==============================================================
; Open the file: Read only the user specified period first observations then model
; ==============================================================
f = addfile("ERAIL_SF_SRO_SSRO_TP_E_1979_2010.nc", "r") ;
;***********************************************************
; Read time and create required yyyyddd                    
;***********************************************************
   time    = f->time                          ; time:units = "hours since 1-1-1 00:00:0.0"                               
   TIME    = cd_calendar(time, 0)          ; type float 
   year    = toint( TIME(:,0) )
   month   = toint( TIME(:,1) )
   day     = toint( TIME(:,2) ) 
   ddd     = day_of_year(year, month, day) 
   yyyyddd = year*1000 + ddd                                  ; needed for input
;***********************************************************
; Read data: short2flt                                     
;*********************************************************** 
   e0     =  short2flt( f->e(:,:,:) )    ; convert to float 
   printVarSummary( e0 ) 
;***********************************************************
; Compute daily climatology: raw and then 'smoothed'  
;***********************************************************
   EClmDay    = clmDayTLL(e0, yyyyddd)     ; daily climatology at each grid point                                      
   printVarSummary(EClmDay)                                                                  
         printVarSummary(EClmDay)                                ; (time, lat,lon)  


and the output is:

Variable: e0
Type: float
Total Size: 233385984 bytes
            58346496 values
Number of Dimensions: 3
Dimensions and sizes:    [time | 11688] x [latitude | 64] x [longitude | 78]
Coordinates: 
            time: [692520..973008]
            latitude: [64..32.5]
            longitude: [30..68.5]
Number Of Attributes: 7
  units :    m of water equivalent
  long_name :    Evaporation
  standard_name :    lwe_thickness_of_water_evaporation_amount
  _FillValue_original :    -32767
  _FillValue :    -32767
  missing_value_original :    -32767
  missing_value :    -32767

Variable: EClmDay
Type: float
Total Size: 7308288 bytes
            1827072 values
Number of Dimensions: 3
Dimensions and sizes:    [year_day | 366] x [latitude | 64] x [longitude | 78]
Coordinates: 
            year_day: [1..366]
            latitude: [64..32.5]
            longitude: [30..68.5]
Number Of Attributes: 4
  long_name :    Daily Climatology: Evaporation
  units :    m of water equivalent
  information :    Raw daily averages across all years
  smoothing :    None

Would really appreciate any advice

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170613/0b0d5285/attachment.html 


More information about the ncl-talk mailing list