Changeset 41 for trunk/procs/meshes/mesh_from_file.pro

Timestamp:

01/26/08 10:18:52 (16 years ago)

Author:

ericg

Message:

Misc modifs, incl. read_from_grid in yfx, debug_w prints, removal of mesh_lmdz and small bug corrections

File:

• trunk/procs/meshes/mesh_from_file.pro (modified) (3 diffs)

trunk/procs/meshes/mesh_from_file.pro

@@ -40 +40 @@
          ELSE: var_local = var_read_grd_file
       ENDCASE 
+      IF strpos(var_local, '=f(next)') NE -1 THEN BEGIN
+         idx = strpos(var_local, '=f(')
+         var_local = strmid(var_local, 0, idx)
+      ENDIF 
+
+      IF debug_w THEN print, '   var_local = ', var_local
 
       ; find valmask
@@ -49 +55 @@
       FOR i = 0, varcontient.natts-1 DO BEGIN
          att_txt = ncdf_attname(cdfidl, varidl, i)
-         IF att_txt EQ 'missing_value' OR att_txt EQ 'mask value' OR att_txt EQ '_FillValue' THEN ncdf_attget, cdfidl, varidl, att_txt, valmask
-         IF debug_w THEN print, 'valmask found = ',valmask 
+         IF att_txt EQ 'missing_value' OR att_txt EQ 'mask value' OR att_txt EQ '_FillValue' THEN BEGIN 
+            ncdf_attget, cdfidl, varidl, att_txt, valmask
+            IF debug_w THEN print, '     valmask found = ',valmask, ' in attribute ', att_txt
+         ENDIF 
       ENDFOR 
       ncdf_close, cdfidl
-      ; read grid
-      initncdf, s_file, USEASMASK = var_local, missing_value = valmask, GLAMBOUNDARY = glamboundary_box, ZAXISNAME = 'depth'
-
+      ; build grid
+      CASE mesh_type of
+         'atm': initncdf, s_file, GLAMBOUNDARY = glamboundary_box, ZAXISNAME = name_level
+         ELSE: initncdf, s_file, USEASMASK = var_local, missing_value = valmask, GLAMBOUNDARY = glamboundary_box, ZAXISNAME = 'depth'
+      ENDCASE 
+      IF debug_w THEN print,  '     Grid info after initncdf :'
+      IF debug_w THEN print,  '     glamt :', size(glamt)
+      IF debug_w THEN print,  '     gphit :', size(gphit)
+      IF debug_w THEN print,  '     gdept :', size(gdept)
+      
    ENDELSE 
 
@@ -65 +80 @@
 ;
    diaznl_idx = 1
-   IF debug_w THEN print, '     model, vargrid, varexp at exit = ',model, vargrid,  varexp
+   IF debug_w THEN print, '     model, vargrid, varexp at exit = ',model, ' ', vargrid,' ',  varexp
    IF debug_w THEN print, '  ...EXIT mesh_from_file'
+   IF debug_w THEN print, '  '
 
   return