Changeset 302 for trunk/SRC/Grid/micromeshmask.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/micromeshmask.pro (modified) (3 diffs)

trunk/SRC/Grid/micromeshmask.pro

@@ -102 +102 @@
   ENDELSE
 ; get the variables list related to the partial steps
-  varlist_ps = ncdf_listvars(cdfid)
-  varlist_ps = strtrim(strlowcase(varlist_ps), 2)
+  zgr_varlist = ncdf_listvars(cdfid)
+  zgr_varlist = strtrim(strlowcase(zgr_varlist), 2)
 ;------------------------------------------------------
 ;------------------------------------------------------
@@ -135 +135 @@
     varid[16+z] = ncdf_vardef(cdfidout, zgrlist[z], [dimidz], /float)
 ; variables related to the partial steps
-  IF (where(varlist_ps EQ 'hdept'))[0] NE -1 THEN $
+  IF (where(zgr_varlist EQ 'hdept'))[0] NE -1 THEN $
     varid = [varid, ncdf_vardef(cdfidout, 'hdept', [dimidx, dimidy], /float)]
-  IF (where(varlist_ps EQ 'hdepw'))[0] NE -1 THEN $
+  IF (where(zgr_varlist EQ 'hdepw'))[0] NE -1 THEN $
     varid = [varid, ncdf_vardef(cdfidout, 'hdepw', [dimidx, dimidy], /float)]
 ; old variable name. keep for compatibility with old run. Change e3tp to e3t_ps
-  IF (where(varlist_ps EQ 'e3tp'))[0] NE -1 THEN $
+  IF (where(zgr_varlist EQ 'e3tp'))[0] NE -1 THEN $
     varid = [varid, ncdf_vardef(cdfidout, 'e3t_ps', [dimidx, dimidy], /float)]
 ; old variable name. keep for compatibility with old run. Change e3wp to e3w_ps
-  IF (where(varlist_ps EQ 'e3wp'))[0] NE -1 THEN $
+  IF (where(zgr_varlist EQ 'e3wp'))[0] NE -1 THEN $
     varid = [varid, ncdf_vardef(cdfidout, 'e3w_ps', [dimidx, dimidy], /float)]
 ;
-  IF (where(varlist_ps EQ 'e3t_ps'))[0] NE -1 THEN $
+  IF (where(zgr_varlist EQ 'e3t_ps'))[0] NE -1 THEN $
     varid = [varid, ncdf_vardef(cdfidout, 'e3t_ps', [dimidx, dimidy], /float)]
-  IF (where(varlist_ps EQ 'e3w_ps'))[0] NE -1 THEN $
+  IF (where(zgr_varlist EQ 'e3w_ps'))[0] NE -1 THEN $
     varid = [varid, ncdf_vardef(cdfidout, 'e3w_ps', [dimidx, dimidy], /float)]
-;   IF (where(varlist_ps EQ 'e3u_ps'))[0] NE -1 THEN $
+;   IF (where(zgr_varlist EQ 'e3u_ps'))[0] NE -1 THEN $
 ;     varid = [varid, ncdf_vardef(cdfidout, 'e3u_ps', [dimidx, dimidy], /float)]
-;   IF (where(varlist_ps EQ 'e3v_ps'))[0] NE -1 THEN $
+;   IF (where(zgr_varlist EQ 'e3v_ps'))[0] NE -1 THEN $
 ;     varid = [varid, ncdf_vardef(cdfidout, 'e3v_ps', [dimidx, dimidy], /float)]
+  IF (where(zgr_varlist EQ 'mbathy'))[0] NE -1 THEN $
+    varid = [varid, ncdf_vardef(cdfidout, 'mbathy', [dimidx, dimidy], /short)]
 ; mask variable
   msklist = ['tmask', 'umask', 'vmask', 'fmask']
@@ -185 +187 @@
   ENDIF
 ;
+  inzgrlist =  zgrlist
+  IF (where(varlist EQ 'gdept_0'))[0] NE -1 THEN inzgrlist =  inzgrlist+'_0'
   FOR z = 0, n_elements(zgrlist)-1 DO $
-    ncdf_transfer, cdfid, cdfidout, zgrlist[z]
+    ncdf_transfer, cdfid, cdfidout, inzgrlist[z], zgrlist[z]
 ; partial step variables
-  IF (where(varlist_ps EQ 'hdept'))[0] NE -1 THEN $
+  IF (where(zgr_varlist EQ 'hdept'))[0] NE -1 THEN $
     ncdf_transfer, cdfid, cdfidout, 'hdept'
-  IF (where(varlist_ps EQ 'hdepw'))[0] NE -1 THEN $
+  IF (where(zgr_varlist EQ 'hdepw'))[0] NE -1 THEN $
     ncdf_transfer, cdfid, cdfidout, 'hdepw'
-  IF (where(varlist_ps EQ 'e3tp'))[0] NE -1 THEN $
+  IF (where(zgr_varlist EQ 'e3tp'))[0] NE -1 THEN $
     ncdf_transfer, cdfid, cdfidout, 'e3tp', 'e3t_ps'
-  IF (where(varlist_ps EQ 'e3wp'))[0] NE -1 THEN $
+  IF (where(zgr_varlist EQ 'e3wp'))[0] NE -1 THEN $
     ncdf_transfer, cdfid, cdfidout, 'e3wp', 'e3w_ps'
-  IF (where(varlist_ps EQ 'e3t_ps'))[0] NE -1 THEN $
+  IF (where(zgr_varlist EQ 'e3t_ps'))[0] NE -1 THEN $
     ncdf_transfer, cdfid, cdfidout, 'e3t_ps'
-  IF (where(varlist_ps EQ 'e3w_ps'))[0] NE -1 THEN $
+  IF (where(zgr_varlist EQ 'e3w_ps'))[0] NE -1 THEN $
     ncdf_transfer, cdfid, cdfidout, 'e3w_ps'
-;   IF (where(varlist_ps EQ 'e3u_ps'))[0] NE -1 THEN $
+;   IF (where(zgr_varlist EQ 'e3u_ps'))[0] NE -1 THEN $
 ;     ncdf_transfer, cdfid, cdfidout, 'e3u_ps'
-;   IF (where(varlist_ps EQ 'e3v_ps'))[0] NE -1 THEN $
+;   IF (where(zgr_varlist EQ 'e3v_ps'))[0] NE -1 THEN $
 ;     ncdf_transfer, cdfid, cdfidout, 'e3v_ps'
+  IF (where(zgr_varlist EQ 'mbathy'))[0] NE -1 THEN $
+    ncdf_transfer, cdfid, cdfidout, 'mbathy'
 ;
 ; mask