Changeset 275

Timestamp:

09/03/07 15:59:09 (17 years ago)

Author:

smasson

Message:

reduce memory use

File:

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

trunk/SRC/ToBeReviewed/LECTURE/read_ncdf.pro

@@ -324 +324 @@
     valmask = 1.e20
     case 1 of
-      varinq.ndims eq 2:BEGIN ;xy array
-        earth = where(mask[*, *, 0] EQ 0)
-      END
-      varinq.ndims eq 3 AND (where(varinq.dim EQ inq.recdim))[0] EQ -1:BEGIN ;xyz array
-        earth = where(mask EQ 0)
-      END
-      varinq.ndims eq 3 AND (where(varinq.dim EQ inq.recdim))[0] NE -1:BEGIN ;xyt array    
-        earth = where(mask[*, *, 0] EQ 0)
-        if earth[0] NE -1 then BEGIN
-          earth = earth#replicate(1, jpt)+replicate(nx*ny, n_elements(earth))#lindgen(jpt)
-        END
-      END
-      varinq.ndims eq 4:BEGIN ;xyzt array
-        earth = where(mask EQ 0)
-        if earth[0] NE -1 then BEGIN
-          earth = earth#replicate(1, jpt)+replicate(nx*ny*nz, n_elements(earth))#lindgen(jpt)
-        END
-      END
+      varinq.ndims eq 2:                                               earth = where(mask[*, *, 0] EQ 0) ;xy   array
+      varinq.ndims eq 3 AND (where(varinq.dim EQ inq.recdim))[0] EQ -1:earth = where(mask EQ 0)          ;xyz  array
+      varinq.ndims eq 3 AND (where(varinq.dim EQ inq.recdim))[0] NE -1:earth = where(mask[*, *, 0] EQ 0) ;xyt  array    
+      varinq.ndims eq 4:                                               earth = where(mask EQ 0)          ;xyzt array
     ENDCASE
   ENDIF ELSE earth = -1
@@ -347 +333 @@
 ; we apply add_offset, scale_factor and missing_value
   IF keyword_set(addscl_before) THEN BEGIN
-    if scale_factor NE 1 then res = temporary(res)*scale_factor
-    if add_offset NE 0 then res = temporary(res)+add_offset
-  ENDIF
-  IF size(missing_value, /type) NE 7 then BEGIN
+    IF scale_factor NE 1 THEN res = temporary(res)*scale_factor
+    IF add_offset   NE 0 THEN res = temporary(res)+add_offset
+  ENDIF
+  IF size(missing_value, /type) NE 7 THEN BEGIN
     IF finite(missing_value) EQ 1 THEN BEGIN
         CASE 1 OF
@@ -364 +350 @@
     if add_offset NE 0 then res = temporary(res)+add_offset
   ENDIF
-  if missing[0] NE -1 then res[temporary(missing)] = !values.f_nan
-  if earth[0] NE -1 then res[temporary(earth)] = 1.e20
+  IF missing[0] NE -1 THEN res[temporary(missing)] = !values.f_nan
+  IF earth[0] NE -1 THEN BEGIN 
+    IF varinq.ndims eq 3 AND (where(varinq.dim EQ inq.recdim))[0] NE -1 THEN $   ;xyt array    
+       earth = earth#replicate(1, jpt)+replicate(nx*ny, n_elements(earth))#lindgen(jpt)
+    IF varinq.ndims eq 4 THEN earth = earth#replicate(1, jpt)+replicate(nx*ny*nz, n_elements(earth))#lindgen(jpt)
+    res[temporary(earth)] = 1.e20
+  ENDIF 
 ;---------------------------------------------------------------------
 ; if it is roms outputs, we need to get additionals infos...