Changeset 19 for trunk/initorca2_bab.pro

Timestamp:

05/15/07 15:07:51 (17 years ago)

Author:

pinsard

Message:

some improvements on validation of step 1 and begining of parametrisation to orcares

File:

• trunk/initorca2_bab.pro (modified) (3 diffs)

trunk/initorca2_bab.pro

@@ -10 +10 @@
 ; @todo
 ; parametrisation for ORCA025
+; can I use CASE statement in a file called by @...++
+; can I use RETURN or STOP statement in a file called by @...++
 ;
 ; zminmesh and izmaxmesh may be set to 0,0 because we only work with 
@@ -31 +33 @@
 @cm_4data
 ;
-;---------------------------------------------------------
+
+; test if ${GEOMAG_ID} defined
+geomag_id_env=GETENV('GEOMAG_ID')
+; ++ pb syntaxe case dans @...
+
+iodirin = isadirectory(geomag_id_env)
+;
+; existence and protection of ${GEOMAG_ID} ++
+
+filename_oce='meshmask_bab.nc'
+; check if this file exists
+fullfilename_oce = isafile(iodirin + filename_oce, NEW=0, /MUST_EXIST, $
+                     RECURSIVE=0)
+print,fullfilename_oce[0]
+
+;--------------------------------------------------------
 ; boundaries of the read grid regarding to the original grid
 ;----------------------------------------------------------
@@ -49 +66 @@
 ; read the grid
 ;------------------------------------------------------
-iodir = '/usr/work/sur/fvi/OPA/ORCA2/'
 
 key_stride = [1, 1, 1]
-ncdf_meshlec, 'meshmask_bab.nc', glamboundary = [20, 380]
+ncdf_meshlec, fullfilename_oce[0], glamboundary = [20, 380]
 ;-------------------------------------------------------------
 domdef