Changeset 3

Timestamp:

09/26/07 10:53:22 (17 years ago)

Author:

kolasinski

Message:

use of ncdf_meshread routine instead of ncdf_meshlec obsolete routine

File:

• trunk/procs/meshes/mesh_orca.pro (modified) (1 diff)

trunk/procs/meshes/mesh_orca.pro

@@ -115 +115 @@
                   'V2': mesmsk = 'meshmask_ORCA_R2.nc'
                   'V3': mesmsk = 'meshmask_ORCA_R2_V2.nc'
-               ENDCASE 
-               print, '                         mask config '+orca_mask_version
+               ENDCASE
                IF whole_arrays EQ 0 THEN BEGIN 
-                  ncdf_meshlec, hom_idl+'grids/'+mesmsk, glam = glamboundary_box, onearth = onearth
+                   cmd_grid = 'ncdf_meshread, sm_file, GLAMBOUNDARY = glamboundary_box, onearth = onearth'
                ENDIF ELSE BEGIN 
-                  ncdf_meshlec, hom_idl+'grids/'+mesmsk, onearth = onearth
+                   cmd_grid = 'ncdf_meshread, sm_file, GLAMBOUNDARY = glamboundary_box'
                ENDELSE 
-                  
             END 
             'L300': BEGIN
-               ncdf_meshlec, hom_idl+'grids/meshmask.orca.2d.L300.nc', glam = glamboundary_box
+               mesmsk = 'meshmask.orca.2d.L300.nc'
+               cmd_grid = 'ncdf_meshread, sm_file, GLAMBOUNDARY = glamboundary_box'
             END 
          ENDCASE 
       END 
-      'ORCA_R4': BEGIN 
-         ncdf_meshlec, hom_idl+'grids/meshmask_orca4.nc', glam = glamboundary_box 
-      END 
-   ENDCASE 
-
-   IF keyword_set(WHOLE_ARRAYS) THEN print, '                         [whole array]'               
+      'ORCA_R4': BEGIN
+         mesmsk = 'meshmask_orca4.nc'
+         cmd_grid = 'ncdf_meshread, sm_file, GLAMBOUNDARY = glamboundary_box' 
+      END 
+   ENDCASE
+   
+   print, '                         mask config '+orca_mask_version
+   IF keyword_set(WHOLE_ARRAYS) THEN print, '                         [whole array]'          
+        
+   sm_file = hom_idl+'grids/'+mesmsk
+   res = find(sm_file)
+
+   IF res NE 'NOT FOUND' THEN BEGIN
+       res_grid =  execute(cmd_grid)
+   ENDIF ELSE BEGIN
+       stop, 'No meshmask file found for ORCA mesh config'
+   ENDELSE
 
 ; reduce grid to zonal mean