Changeset 302 for trunk/SRC/Grid/ncdf_meshread.pro

Timestamp:

10/10/07 14:29:38 (17 years ago)

Author:

smasson

Message:

compatibility with new meshmask using gdept_O, ...

File:

• trunk/SRC/Grid/ncdf_meshread.pro (modified) (3 diffs)

trunk/SRC/Grid/ncdf_meshread.pro

@@ -12 +12 @@
 ; @param filename {in}{optional}{default='meshmask.nc'}{type=scalar string}
 ;    Name of the meshmask file to read. If this name does not contain any "/"
-;    and if iodirectory keyword is not specify, then the common variable
+;    and if iodirectory keyword is not specified, then the common variable
 ;    iodir will be used to define the mesh file path.
 ;
@@ -479 +479 @@
 ; 1d arrays
 ;-------------------------------------------------------
+  IF (where(allvarname EQ 'e3t_O'))[0] NE -1 THEN fnamevar = ['e3t_O', 'e3w_O', 'gdept_O', 'gdepw_O'] $
+  ELSE                                            fnamevar = ['e3t', 'e3w', 'gdept', 'gdepw'] 
   namevar = ['e3t', 'e3w', 'gdept', 'gdepw']
   for i = 0, n_elements(namevar)-1 do begin
@@ -484 +486 @@
     CASE n_elements(varinq.dim) OF
       4:BEGIN
-        command = 'ncdf_varget,cdfid,namevar[i],'+namevar[i] $
+        command = 'ncdf_varget,cdfid,fnamevar[i],'+namevar[i] $
                    +',offset = [0,0,izminmesh,0], count = [1,1,jpk,1]'
         if key_stride[2] NE 1 then command = command+', stride=[1,1,key_stride[2],1]'
       END
       2:BEGIN
-        command = 'ncdf_varget,cdfid,namevar[i],'+namevar[i] $
+        command = 'ncdf_varget,cdfid,fnamevar[i],'+namevar[i] $
                    +',offset = [izminmesh,0], count = [jpk,1]'
         if key_stride[2] NE 1 then command = command+', stride=key_stride[2]'
       END
       1:BEGIN
-        command = 'ncdf_varget,cdfid,namevar[i],'+namevar[i] $
+        command = 'ncdf_varget,cdfid,fnamevar[i],'+namevar[i] $
                    +',offset = [izminmesh], count = [jpk]'
         if key_stride[2] NE 1 then command = command+', stride=key_stride[2]'