Changeset 501 for trunk/SRC/ToBeReviewed/LECTURE/read_ncdf.pro

Timestamp:

09/20/16 20:15:00 (8 years ago)

Author:

smasson

Message:

set of bugfixes...

File:

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

trunk/SRC/ToBeReviewed/LECTURE/read_ncdf.pro

@@ -120 +120 @@
                     , GRID = grid, CALLITSELF = callitself, DIREC = direc $
                     , ZETAFILENAME = zetafilename, ZETAZERO = zetazero $
-                    , ZINVAR = zinvar, _EXTRA = ex
+                    , ZINVAR = zinvar, ISROMS = isroms, _EXTRA = ex
 ;
   compile_opt idl2, strictarrsubs
@@ -164 +164 @@
     dimnames[i] = tmp
   ENDFOR
-;------------------------------------------------------------
-; we check if the variable has a vertical dimension and/or
-; a record dimension. This is useful to define boxzoom
-; keyword in domdef
-;------------------------------------------------------------
-  dummy = where(varinq.dim EQ inq.recdim, tinvar)
-; check the presence of a vertical dimension according to the
-; number of dimensions and the presence of a record dimension
-  zinvar = (varinq.ndims EQ 3 AND tinvar NE 1) OR varinq.ndims EQ 4
-;------------------------------------------------------------
-; shall we redefine the grid parameters
-;------------------------------------------------------------
-  IF keyword_set(init) THEN initncdf, filename, _extra = ex
 ;------------------------------------------------------------
 ; check the time axis and the debut and ending dates
@@ -248 +235 @@
     END
   ENDCASE
+  IF inq.recdim EQ -1 AND jpt NE 1 THEN inq.recdim = varinq.dim[varinq.ndims-1]
+;------------------------------------------------------------
+; we check if the variable has a vertical dimension and/or
+; a record dimension. This is useful to define boxzoom
+; keyword in domdef
+;------------------------------------------------------------
+  dummy = where(varinq.dim EQ inq.recdim, tinvar)
+; check the presence of a vertical dimension according to the
+; number of dimensions and the presence of a record dimension
+  zinvar = (varinq.ndims EQ 3 AND tinvar NE 1) OR varinq.ndims EQ 4
+;------------------------------------------------------------
+; shall we redefine the grid parameters
+;------------------------------------------------------------
+  IF keyword_set(init) THEN initncdf, filename, _extra = ex
 ;------------------------------------------------------------
 ; Name of the grid on which the field refer to.
@@ -405 +406 @@
 ; if it is roms outputs, we need to get additional infos...
   IF NOT keyword_set(callitself) THEN BEGIN
-    IF strmid(dimnames[0], 0, 3) EQ 'xi_' AND strmid(dimnames[1], 0, 4) EQ 'eta_' THEN BEGIN
+    IF (strmid(dimnames[0], 0, 3) EQ 'xi_' AND strmid(dimnames[1], 0, 4) EQ 'eta_') OR keyword_set(isroms) THEN BEGIN
       ncdf_attget, cdfid, 'theta_s', theta_s, /global
       ncdf_attget, cdfid, 'theta_b', theta_b, /global