[ncl-talk] Deriving monthly mean from files of ten years daily data
Rick Brownrigg
brownrig at ucar.edu
Fri Jun 16 11:48:41 MDT 2017
Hi,
There must be more going on to this script than is shown. The printed
output does not align with print statements in the script, and I don't
believe an addfiles command works with ".txt" files. The error message is
coming from asciiread(), which does not appear anywhere in the script
you've shown. In any case, the error message is telling you that there is
no variable names AOT_550 in the files. Are there other error messages from
the script?
Rick
On Fri, Jun 16, 2017 at 8:41 AM, <najib.yusuf at carnasrda.com> wrote:
>
>
> -------- 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
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170616/961ba2cc/attachment.html
More information about the ncl-talk
mailing list