[ncl-talk] Deriving monthly mean from files of ten years daily data

najib.yusuf at carnasrda.com najib.yusuf at carnasrda.com
Fri Jun 16 08:41:17 MDT 2017



-------- Original Message --------
Subject: Deriving monthly mean from files of ten years data
Date: 2017-06-16 10:00
 From: najib.yusuf at carnasrda.com
To: ncl-talk at ucar.edu
Cc: shea at ucar.edu

     Dear ncl experts,


     I am trying to derived monthly mean (for example all January, all 
February....all December) mean from 10 files (0-9) years of Julian date, 
so I can plot the mean on the bulk of the data. I tried but I got an 
error;




        ;  to derive the mean of the data
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
   wkdir = "./"
   fil_dat  = "AOT_550"
; list all AOT files for different years
   files = systemfunc("ls " + wkdir +"AOT_550*.txt")
   print(files)

; read years (files)
   files10 = files(0:9)
   print(files10)
   f = addfiles(files10, "r")
AOT = addfiles_GetVar(f,files10,"AOT_550")
time = addfiles_GetVar(f,files10,"time")

; get date
t1=calendar_decode2(time,0)
; get day of year
printVarSummary(t1)
ndays = 365
doy = ndays
doy(9) = 365
  print(doy)

; read only  monthly mean (not sure?)
AOT_550_all12 = AOT_550_all12(0:11)
print(AOT_550_all12)

f = addfiles(files12, "r")
AOT = addfiles_GetVar(f,AOT_550_all12,"AOT_550_all")
time = addfiles_GetVar(f,AOT_550_all12,"time")
; get date
t1=calendar_decode2(time,0)
; get day of year
printVarSummary(t1)
doy=day_of_year(floattoint(t1(:,0)),floattoint(t1(:,1)),floattoint(t1(:,2)))
; move to mid-months
doy=doy-15
doy(11:9) = 365 ;
print(doy)
AOT_550_all12 = AOT

printVarSummary(AOT)

Variable: variables
Type: string
Total Size: 56 bytes
             7 values
Number of Dimensions: 1
Dimensions and sizes:	[7]
Coordinates:
Number Of Attributes: 1
   _FillValue :	missing
(0)	Day
(1)	AOT_675
(2)	AOT_500
(3)	AOT_440
(4)	440-675Angstrom
(5)	╬▒
(6)	AOT_550


Variable: files (subsection)
Type: string
Total Size: 8 bytes
             1 values
Number of Dimensions: 1
Dimensions and sizes:	[1]
Coordinates:
(0)	/home/model-user/AOT_550_2014.txt
warning:asciiread: End of file reached and only (157) elements were read 
from the file, filling remaining elements with the default missing value 
for the requested type

Variable: z1
Type: string
Total Size: 2928 bytes
             366 values
Number of Dimensions: 1
Dimensions and sizes:	[366]
Coordinates:
Number Of Attributes: 1
   _FillValue :	missing

Variable: new_str
Type: string
Total Size: 20496 bytes
             2562 values
Number of Dimensions: 2
Dimensions and sizes:	[366] x [7]
Coordinates:
Number Of Attributes: 1
   _FillValue :	missing


Variable: variables
Type: string
Total Size: 56 bytes
             7 values
Number of Dimensions: 1
Dimensions and sizes:	[7]
Coordinates:
Number Of Attributes: 1
   _FillValue :	missing
(0)	Day
(1)	AOT_675
(2)	AOT_500
(3)	AOT_440
(4)	440-675Angstrom
(5)	╬▒
(6)	AOT_550

Variable: AOT_550_all
Type: float
Total Size: 14600 bytes
             3650 values
Number of Dimensions: 2
Dimensions and sizes:	[years | 10] x [days | 365]
Coordinates:
Number Of Attributes: 1
   _FillValue :	9.96921e+36


Variable: files
Type: string
Total Size: 80 bytes
             10 values
Number of Dimensions: 1
Dimensions and sizes:	[10]
Coordinates:
(0)	/home/model-user/AOT_550_2005.txt
(1)	/home/model-user/AOT_550_2006.txt
(2)	/home/model-user/AOT_550_2007.txt
(3)	/home/model-user/AOT_550_2008.txt
(4)	/home/model-user/AOT_550_2009.txt
(5)	/home/model-user/AOT_550_2010.txt
(6)	/home/model-user/AOT_550_2011.txt
(7)	/home/model-user/AOT_550_2012.txt
(8)	/home/model-user/AOT_550_2013.txt
(9)	/home/model-user/AOT_550_2014.txt


Variable: files10
Type: string
Total Size: 80 bytes
             10 values
Number of Dimensions: 1
Dimensions and sizes:	[10]
Coordinates:
(0)	/home/model-user/AOT_550_2005.txt
(1)	/home/model-user/AOT_550_2006.txt
(2)	/home/model-user/AOT_550_2007.txt
(3)	/home/model-user/AOT_550_2008.txt
(4)	/home/model-user/AOT_550_2009.txt
(5)	/home/model-user/AOT_550_2010.txt
(6)	/home/model-user/AOT_550_2011.txt
(7)	/home/model-user/AOT_550_2012.txt
(8)	/home/model-user/AOT_550_2013.txt
(9)	/home/model-user/AOT_550_2014.txt
fatal:No valid instance of variable AOT_550 found in file list
fatal:["Execute.c":8575]:Execute: Error occurred at or near line 6095 in 
file 
/home/model-user/NCL/ncl-6.3.0/lib/ncarg/nclscripts/csm/contributed.ncl


    Please assist to put me through.

   Thank you for your time.


   Najib


More information about the ncl-talk mailing list