Changeset 332 for trunk

Timestamp:

02/11/08 17:48:24 (16 years ago)

Author:

smasson

Message:

adapt for very large array

File:

• trunk/SRC/ReadWrite/write_oasis.pro (modified) (2 diffs)

trunk/SRC/ReadWrite/write_oasis.pro

@@ -23 +23 @@
 ; data to be appended.
 ;
+; @keyword RECSIZE
+; define the size of the full data array to be written. Usefull when
+; you want to save memory and write the data in several write_oasis
+; instructions. see example
+;
+; @keyword TEMPORARY
+; activate undefine z2d when yo write it (to save memory) -> z2d will
+; be lost once write_oasis is returning.
+
+; @keyword HEADER
+; activate to write the header ("character*8" contained in varname)
+; before writting the data. Used when recsize is defined and /= 0, see example
+;
+; @keyword ENDING
+; activate when you write the last part of the data. Used when recsize
+; is defined and /= 0, see example
+;
+; @keyword POSITION
+; specify the position (in byte) at which you want to write the
+; data. Used when recsize is defined and /= 0, see example
+;
 ; @restrictions
 ; varname is automatically written as a "character*8"
 ; by default z2d is written as an R8 array
+;
+; @examples
+;
+;   write_oasis, fa2of, 'WEIGHTS5', weig
+;
+; or in several call so save memory
+;
+;   ysz = 100L
+;   recsz8 = 16L * jpio * jpjo * 8L
+;   FOR i = 0L, jpjo-1L, ysz DO BEGIN 
+;     ii = (i+ysz-1L) < (jpjo-1L)
+;     position = (4L + 8L + 4L + 4L)*(i NE 0) + 16L * jpio * i * 8L
+;     weig = .... 
+;     write_oasis, fa2ou, 'WEIGHTS3', temporary(weig), /temporary, append = i NE 0, header = i EQ 0 $
+;                  , ending = ii EQ jpjo-1, recsize = recsz8, position = position
+;   ENDFOR
+;
 ;
 ; @history
@@ -34 +72 @@
 ;
 ;-
-PRO write_oasis, filename, varname, z2d, I2=i2, I4=i4, I8=i8, R4=r4, APPEND=append
+PRO write_oasis, filename, varname, z2d, I2 = i2, I4 = i4, I8 = i8, R4 = r4, APPEND = append $
+                 , RECSIZE = recsize, TEMPORARY = temporary, HEADER = header, ENDING = ending, POSITION = position
 ;
   compile_opt idl2, strictarrsubs
 ;
-   openw, unit, filename , /F77_UNFORMATTED, /GET_LUN, /SWAP_IF_LITTLE_ENDIAN $
-    , error=err, APPEND = append
-   if err ne 0 then begin
-      ras = report(!err_string)
-      return
-   endif
+  openw, unit, filename, F77_UNFORMATTED = keyword_set(recsize) EQ 0, /GET_LUN, /SWAP_IF_LITTLE_ENDIAN $
+         , error = err, APPEND = append
+  IF err NE 0 THEN BEGIN
+    ras = report(!err_string)
+    return
+  ENDIF
+  IF n_elements(position) NE 0 THEN point_lun, unit, position
 
-   writeu, unit, string(varname, format='(a8)')
-   case 1 of
-      keyword_set(i2):writeu, unit, fix(z2d)
-      keyword_set(i4):writeu, unit, long(z2d)
+  IF keyword_set(recsize) THEN BEGIN
+    IF keyword_set(header) THEN BEGIN
+      writeu, unit, 8L
+      writeu, unit, string(varname, format = '(a8)')
+      writeu, unit, 8L
+      writeu, unit, long(recsize)
+    ENDIF
+  ENDIF ELSE BEGIN 
+    writeu, unit, string(varname, format = '(a8)')
+  ENDELSE 
+
+  IF keyword_set(temporary)  THEN BEGIN
+    CASE 1 OF
+      keyword_set(i2):writeu, unit, fix(   temporary(z2d))
+      keyword_set(i4):writeu, unit, long(  temporary(z2d))
+      keyword_set(i8):writeu, unit, long64(temporary(z2d))
+      keyword_set(r4):writeu, unit, float( temporary(z2d))
+      ELSE:           writeu, unit, double(temporary(z2d))
+    ENDCASE
+  ENDIF ELSE BEGIN 
+    CASE 1 OF
+      keyword_set(i2):writeu, unit, fix(   z2d)
+      keyword_set(i4):writeu, unit, long(  z2d)
       keyword_set(i8):writeu, unit, long64(z2d)
-      keyword_set(r4):writeu, unit, float(z2d)
-      ELSE:writeu, unit, double(z2d)
-   endcase
+      keyword_set(r4):writeu, unit, float( z2d)
+      ELSE:           writeu, unit, double(z2d)
+    ENDCASE
+  ENDELSE 
 
-   free_lun,unit
-   return
-end
+  IF keyword_set(recsize) AND keyword_set(ending) THEN writeu, unit, long(recsize)
+
+  free_lun, unit
+  return
+END