Changeset 212 for trunk/SRC

Timestamp:

02/14/07 11:47:10 (17 years ago)

Author:

smasson

Message:

include [yz]reverse in read_ncdf_varget

Location:

trunk/SRC

Files:

• Grid/ncdf_meshread.pro (modified) (2 diffs)
• Grid/ncdf_meshroms.pro (modified) (2 diffs)
• ToBeReviewed/LECTURE/read_ncdf.pro (modified) (1 diff)
• ToBeReviewed/LECTURE/read_ncdf_varget.pro (modified) (1 diff)

trunk/SRC/Grid/ncdf_meshread.pro

@@ -269 +269 @@
   izminmeshsauve = izminmesh
   izminmesh = 0
+;
+  key_yreverse = 0
+  key_zreverse = 0
 ;-------------------------------------------------------
 ; 2d arrays:
@@ -569 +572 @@
   widget_control, noticebase, bad_id = nothing, /destroy
 ;
-  key_yreverse = 0
-  key_zreverse = 0
   key_gridtype = 'c'
 ;

trunk/SRC/Grid/ncdf_meshroms.pro

@@ -254 +254 @@
   izminmeshsauve = izminmesh
   izminmesh = 0
+;
+  key_yreverse = 0
+  key_zreverse = 1
 ;-------------------------------------------------------
 ; 2d arrays:
@@ -490 +493 @@
   widget_control, noticebase, bad_id = nothing, /destroy
 ;
-  key_yreverse = 0
-  key_zreverse = 1
   key_gridtype = 'c'
 ;

trunk/SRC/ToBeReviewed/LECTURE/read_ncdf.pro

@@ -491 +491 @@
   varexp = file_basename(filename)
 
-; we apply reverse
-  IF keyword_set(key_yreverse) AND ny NE 1 THEN BEGIN 
-    IF varcontient.ndims - ((where(varcontient.dim EQ contient.recdim))[0] NE -1) EQ 2 THEN $
-       res = reverse(reform(res, nx, ny, jpt, /overwrite),  2) $
-    ELSE res = reverse(reform(res, nx, ny, nz, jpt, /overwrite),  2)
-  ENDIF
-  if keyword_set(key_zreverse) AND nz NE 1 $
-     AND varcontient.ndims - ((where(varcontient.dim EQ contient.recdim))[0] NE -1) EQ 3 THEN $
-        res = reverse(reform(res, nx, ny, nz, jpt, /overwrite),  3)
 ; We apply the value valmask on land points.
   if NOT keyword_set(cont_nofill) then begin

trunk/SRC/ToBeReviewed/LECTURE/read_ncdf_varget.pro

@@ -446 +446 @@
       endcase
    END
-endcase
+ ENDCASE
+
+; we apply reverse
+  IF keyword_set(key_yreverse) AND ny NE 1 THEN BEGIN 
+    IF varcontient.ndims - ((where(varcontient.dim EQ contient.recdim))[0] NE -1) EQ 2 THEN $
+       res = reverse(reform(res, nx, ny, jpt, /overwrite),  2) $
+    ELSE res = reverse(reform(res, nx, ny, nz, jpt, /overwrite),  2)
+  ENDIF
+  if keyword_set(key_zreverse) AND nz NE 1 $
+     AND varcontient.ndims - ((where(varcontient.dim EQ contient.recdim))[0] NE -1) EQ 3 THEN $
+        res = reverse(reform(res, nx, ny, nz, jpt, /overwrite),  3)