Changeset 69 for trunk/Grid/ncdf_meshread.pro

Timestamp:

05/11/06 12:35:53 (18 years ago)

Author:

smasson

Message:

debug + new xxx

File:

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

trunk/Grid/ncdf_meshread.pro

@@ -43 +43 @@
 ;    glamboundary (when defined) by using the first line of glamt. if
 ;    key_periodic=0 then in any case key_shift = 0. 
+;
+;    STRCALLING: a string containing the calling command used to call
+;    computegrid (this is used by xxx.pro)
 ;
 ;    STRIDE = : a 3 elements vector to specify the stride in x, y, z
@@ -78 +81 @@
                   , ONEARTH = onearth, GETDIMENSIONS = getdimensions $
                   , PERIODIC = periodic, SHIFT = shift, STRIDE = stride $
-                  , _EXTRA = ex
+                  , STRCALLING = strcalling, _EXTRA = ex
 ;
 ;---------------------------------------------------------
@@ -536 +539 @@
   key_gridtype = 'c'
 ;
-
+;====================================================
+; grid parameters used by xxx
+;====================================================
+;
+  IF NOT keyword_set(strcalling) THEN BEGIN 
+    IF n_elements(ccmeshparameters) EQ 0 THEN strcalling = 'ncdf_meshread' $
+    ELSE strcalling = ccmeshparameters.filename
+  ENDIF 
+  glaminfo = moment(glamt)
+  IF finite(glaminfo[2]) EQ 0 THEN glaminfo = glaminfo[0:1]
+  gphiinfo = moment(gphit)
+  IF finite(gphiinfo[2]) EQ 0 THEN gphiinfo = gphiinfo[0:1]
+  ccmeshparameters = {filename:strcalling  $
+          , glaminfo:glaminfo $
+          , gphiinfo:gphiinfo $
+          , jpiglo:jpiglo, jpjglo:jpjglo, jpkglo:jpkglo $
+          , jpi:jpi, jpj:jpj, jpk:jpk $
+          , ixminmesh:ixminmesh, ixmaxmesh:ixmaxmesh $
+          , iyminmesh:iyminmesh, iymaxmesh:iymaxmesh $
+          , izminmesh:izminmesh, izmaxmesh:izmaxmesh $
+          , key_shift:key_shift, key_periodic:key_periodic $
+          , key_stride:key_stride, key_gridtype:key_gridtype $
+          , key_yreverse:key_yreverse, key_zreverse:key_zreverse $
+          , key_partialstep:key_partialstep, key_onearth:key_onearth}
+;
   if keyword_set(key_performance) THEN $
     print, 'time ncdf_meshread', systime(1)-tempsun