[ncl-talk] adding multiple files once to extract data at desired time period

[Rick Brownrigg]

Are you saying you are having trouble coming up with the right filename
expression?

What about something like:

ls MXLDEPTH.1440x720.[0-9][0-9][0-9][0-9]12*.nc
MXLDEPTH.1440x720.[0-9][0-9][0-9][0-9]0[1-2]*.nc    # DJF, all years
MXLDEPTH.1440x720.[0-9][0-9][0-9][0-9]0[3-5]*.nc    # MAM
MXLDEPTH.1440x720.[0-9][0-9][0-9][0-9]0[6-8]*.nc    # JJA
MXLDEPTH.1440x720.[0-9][0-9][0-9][0-9]09*.nc
 MXLDEPTH.1440x720.[0-9][0-9][0-9][0-9]1[0-1]*.nc    # SON

Its kind of crude and there are likely better ways, but I think that works.

Rick

On Sun, Mar 11, 2018 at 11:54 PM, Wickramage Chathurika <
wickramagechathurika at rocketmail.com> wrote:

> Hello,
>
>
> I have 7671 (21 years) Mixlayer depth files with dimensions, time, lat and
> lon. I want to extract mix layer depth data for "JJA". but I cannot read
> such amount of files at once and give errors;
>
> for that I have used;
>
> diri="/mnt/g/Ocean/ECCO2/MXLDEPTH/"
> fili=systemfunc("cd "+diri+" ;ls MXLDEPTH.1440x720.*.nc")
>
> "warning:_NclOpenFile: Can not open file <*.nc>; file format not supported
> or file is corrupted"
>
> if I add only one-month data, it runs properly and gives the output. for
> that I have used;
>
> diri="/mnt/g/Ocean/ECCO2/MXLDEPTH/"
> fili=systemfunc("cd "+diri+" ;ls MXLDEPTH.1440x720.201106*.nc")
>
> And one more thing, when I try to use only the 21 years June data, It
> gives not only June data but other months data also (the part of the
> printed output has been displayed below). for that I have used;
>
> diri="/mnt/g/Ocean/ECCO2/MXLDEPTH/"
> fili=systemfunc("cd "+diri+" ;ls MXLDEPTH.1440x720.*06*.nc")
>
> (1141)  MXLDEPTH.1440x720.20111006.nc
> (1142)  MXLDEPTH.1440x720.20111106.nc
> (1143)  MXLDEPTH.1440x720.20111206.nc
> (1144)  MXLDEPTH.1440x720.20120106.nc
> (1145)  MXLDEPTH.1440x720.20120206.nc
> (1146)  MXLDEPTH.1440x720.20120306.nc
> (1147)  MXLDEPTH.1440x720.20120406.nc
> (1148)  MXLDEPTH.1440x720.20120506.nc
> (1149)  MXLDEPTH.1440x720.20120601.nc
> (1150)  MXLDEPTH.1440x720.20120602.nc
> (1151)  MXLDEPTH.1440x720.20120603.nc
> (1152)  MXLDEPTH.1440x720.20120604.nc
> (1153)  MXLDEPTH.1440x720.20120605.nc
> (1154)  MXLDEPTH.1440x720.20120606.nc
> (1155)  MXLDEPTH.1440x720.20120607.nc
> (1156)  MXLDEPTH.1440x720.20120608.nc
> (1157)  MXLDEPTH.1440x720.20120609.nc
> (1158)  MXLDEPTH.1440x720.20120610.nc
> (1159)  MXLDEPTH.1440x720.20120611.nc
> (1160)  MXLDEPTH.1440x720.20120612.nc
> (1161)  MXLDEPTH.1440x720.20120613.nc
> (1162)  MXLDEPTH.1440x720.20120614.nc
> (1163)  MXLDEPTH.1440x720.20120615.nc
> (1164)  MXLDEPTH.1440x720.20120616.nc
> (1165)  MXLDEPTH.1440x720.20120617.nc
> (1166)  MXLDEPTH.1440x720.20120618.nc
> (1167)  MXLDEPTH.1440x720.20120619.nc
> (1168)  MXLDEPTH.1440x720.20120620.nc
> (1169)  MXLDEPTH.1440x720.20120621.nc
> (1170)  MXLDEPTH.1440x720.20120622.nc
> (1171)  MXLDEPTH.1440x720.20120623.nc
> (1172)  MXLDEPTH.1440x720.20120624.nc
> (1173)  MXLDEPTH.1440x720.20120625.nc
> (1174)  MXLDEPTH.1440x720.20120626.nc
> (1175)  MXLDEPTH.1440x720.20120627.nc
> (1176)  MXLDEPTH.1440x720.20120628.nc
> (1177)  MXLDEPTH.1440x720.20120629.nc
> (1178)  MXLDEPTH.1440x720.20120630.nc
> (1179)  MXLDEPTH.1440x720.20120706.nc
> (1180)  MXLDEPTH.1440x720.20120806.nc
> (1181)  MXLDEPTH.1440x720.20120906.nc
> (1182)  MXLDEPTH.1440x720.20121006.nc
> (1183)  MXLDEPTH.1440x720.20121106.nc
> (1184)  MXLDEPTH.1440x720.20121206.nc
>
> The variable summary for one file is,
>
> Variable: mix
> Type: float
> Total Size: 560 bytes
>             140 values
> Number of Dimensions: 2
> Dimensions and sizes:   [TIME | 1] x [LONGITUDE_T | 140]
> Coordinates:
>             TIME: [2507.5..2507.5]
>             LONGITUDE_T: [45.125..79.875]
> Number Of Attributes: 6
>   missing_value :       -9.9999998e+22
>   FillValue :   -9.9999998e+22
>   long_name :   Mixed layer depth, dT=.8degC density criterion
>   units :       m
>   LATITUDE_T :  -0.125
>   _FillValue :  1e+20
>
>
> I want to plot the mix layer depth climatology for seasons wise (
> DJF,MAM,JJA,SON). Please help me to solve this problem and I stuck here so
> long time. Any help appriciate in advance. I have tried so may methods up
> to now, but nothing work properly.
>
> Thank you so much
>
> Kind regards,
> Chathu
>
> ........................................
>
>
>
>
> _______________________________________________
> 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/20180312/e3bb89dc/attachment.html>