search

MITgcm / gcmfaces

gcmfaces is a Matlab / Octave toolbox that handles gridded earth variables in generic fashion. Read more at:
http://gcmfaces.readthedocs.io/en/latest/
MIT License
24 stars 21 forks source link

Question about the use of read_nctiles #19

Open GC0225 opened 1 year ago

GC0225 commented 1 year ago

I have copy the update commit 'read_nctiles', but when I read the variables to calculate the volume budget like this:
oceFWflx = read_nctiles( [dir1,'oceFWflx/','oceFWflx'],'oceFWflx',ii ) It will show that '/ECCOv4r4/nctiles_monthly/oceFWflx/oceFWflx.0001.nc not found in the current directory'. Actually, I noticed that the file named by oceFWflx.0001.nc didn't exist in the ECCOv4r4, in which the files are stored in sub-directories of different year and named by the months, rather than divided into 13 tiles like v4r3. So how can I solve it and read the variables successfully to calculate the volume budget with v4r4?

owang01 commented 1 year ago

@Bitterground I can use the updated read_nctiles in PR #17 to read V4r4's files (organized by year and month) without issue (see below). Could it be in your case somehow the path is not specified correctly?

dir1='ECCO/Release4/nctiles_monthly/'; ii=3; oceFWflx = read_nctiles( [dir1,'oceFWflx/','oceFWflx'],'oceFWflx',ii) nFaces: 5

  f1: [90x270 double]
  f2: [90x270 double]
  f3: [90x90 double] 
  f4: [270x90 double]
  f5: [270x90 double]