Changeset 276

Timestamp:

09/04/07 11:55:29 (17 years ago)

Author:

smasson

Message:

add direc keyword to read_ncdf

File:

• trunk/SRC/ToBeReviewed/LECTURE/read_ncdf.pro (modified) (3 diffs)

trunk/SRC/ToBeReviewed/LECTURE/read_ncdf.pro

@@ -41 +41 @@
 ; For ROMS outputs. Use by read_ncdf itself to access auxilliary data (h and zeta).
 ;
+; @keyword DIREC
+; a string used to specify the direction along which we want to make
+; spatial and/or temporal mean. It could be: 'x' 'y' 'z' 't' 'xy' 'xz'
+; 'yz' 'xyz' 'xt' 'yt' 'zt' 'xyt' 'xzt' 'yzt' or 'xyzt' 
+;
 ; @keyword FILENAME {required}{type=string}
 ; It contains he file's name.
@@ -96 +101 @@
                     , PARENTIN = parentin, TIMESTEP = timestep, ADDSCL_BEFORE = addscl_before $
                     , TOUT = tout, NOSTRUCT = nostruct, CONT_NOFILL = CONT_NOFILL, INIT = init $
-                    , GRID = grid, CALLITSELF = callitself $
+                    , GRID = grid, CALLITSELF = callitself, DIREC = direc $
                     , ZETAFILENAME = zetafilename, ZETAZERO = zetazero $
                     , _EXTRA = ex
@@ -386 +391 @@
   ENDIF
 ;---------------------------------------------------------------------
+  IF keyword_set(direc) THEN BEGIN
+    IF jpt EQ 1 THEN res = moyenne(temporary(res), direc) $
+    ELSE BEGIN
+      res = grossemoyenne(temporary(res), direc)
+      IF ( strpos(strlowcase(direc), 't') ge 0 ) THEN BEGIN
+        vardate = strtrim(jul2date(time[0]), 1)+' - '+strtrim(jul2date(time[jpt-1]), 1)
+        time = total(time)/float(jpt)
+        jpt = 1
+      ENDIF
+    ENDELSE 
+  ENDIF
+;---------------------------------------------------------------------
   ncdf_close, cdfid