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

Timestamp:

11/27/07 17:13:57 (16 years ago)

Author:

kolasinski

Message:

Add var_read_grd_file in def_grid and mesh_from_file - Define valmask if needed when the grid is read in mesh_from_file

File:

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

trunk/procs/meshes/mesh_from_file.pro

@@ -6 +6 @@
 ; masks made in data_read
 ;
+   IF debug_w THEN print, '  ENTER mesh_from_file...'
 
    IF strpos(ncdf_db, ':') GE 1 THEN directory = (str_sep(ncdf_db, ':'))[1] $
@@ -12 +13 @@
    s_file = directory+file_name
 
-   print,' Model inits for ', model, ' from file', s_file
+   print,'   Model inits for ', model, ' from file: ', s_file
 
    sm_file = hom_idl+'grids/grids_'+model+'.nc' 
@@ -21 +22 @@
       initncdf, s_file, GLAMBOUNDARY = glamboundary_box
 
-      IF debug_w THEN print, ' Found mask from ',sm_file 
+      print, '    Found mask from ',sm_file 
 
       tmask = byte(read_ncdf('sftlf', 0, 0, /timestep, file = sm_file, /nostruct))
@@ -35 +36 @@
 
    ENDIF ELSE BEGIN 
-      initncdf, s_file, USEASMASK = var_name, GLAMBOUNDARY = glamboundary_box, ZAXISNAME = 'depth'
+      CASE var_name OF
+         '@@voenergy': var_local = 'so'
+         ELSE: var_local = var_read_grd_file
+      ENDCASE 
+
+      ; find valmask
+      cdfidl=ncdf_open(s_file) 
+      contient=ncdf_inquire(cdfidl)
+      varidl = ncdf_varid(cdfidl, var_local)
+      varcontient=ncdf_varinq(cdfidl, var_local)
+      valmask = 1.e20
+      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 
+      ENDFOR 
+      ncdf_close, cdfidl
+      ; read grid
+      initncdf, s_file, USEASMASK = var_local, missing_value = valmask, GLAMBOUNDARY = glamboundary_box, ZAXISNAME = 'depth'
+
    ENDELSE 
-
-   IF debug_w THEN print, 'vargrid = ',vargrid 
 
    key_offset = [0, 0, 0]
@@ -45 +63 @@
 ;
    diaznl_idx = 1
+   IF debug_w THEN print, '  ...EXIT mesh_from_file'
 
   return