Changeset 435 for trunk/SRC

Timestamp:

06/23/10 10:06:05 (14 years ago)

Author:

smasson

Message:

update computegrid for closea domain

File:

• trunk/SRC/Grid/computegrid.pro (modified) (1 diff)

trunk/SRC/Grid/computegrid.pro

@@ -1109 +1109 @@
 ;
   IF jpi GT 2 AND jpj GT 2 AND NOT keyword_set(plain) $
-     AND ixminmesh EQ 0l AND ixmaxmesh eq jpiglo-1 $
-     AND iyminmesh EQ 0l AND iymaxmesh eq jpjglo-1 $
+;;      AND ixminmesh EQ 0l AND ixmaxmesh eq jpiglo-1 $
+;;      AND iyminmesh EQ 0l AND iymaxmesh eq jpjglo-1 $
      AND total(tmask[*, 0, *]) EQ 0 AND total(tmask[*, jpj-1, *]) EQ 0 $
      AND total(tmask[0, *, *]) EQ 0 AND total(tmask[jpi-1, *, *]) EQ 0 THEN BEGIN
-        xminmesh = 1
-        xmaxmesh = -1
-        yminmesh = 1
-        ymaxmesh = -1
-        computegrid, XAXIS = glamt, YAXIS = gphit, ZAXIS = zaxis $
-                     , MASK = mask, GLAMBOUNDARY = glamboundary $
-                     , XMINMESH = xminmesh, XMAXMESH = xmaxmesh $
-                     , YMINMESH = yminmesh, YMAXMESH = ymaxmesh $
-                     , ZMINMESH = zminmesh, ZMAXMESH = zmaxmesh $
-                     , ONEARTH = onearth, PERIODIC = periodic $
-                     , PLAIN = plain, SHIFT = shift, STRIDE = stride $
-                     , FULLCGRID = fullcgrid, XYINDEX = xyindex $
-                     , UBASE2TBASE = ubase2tbase, VBASE2TBASE = vbase2tbase $
-                     , FBASE2TBASE = fbase2tbase, STRCALLING = strcalling $
-                     , ROMSH = romsh, _extra = ex
-        return
+    IF NOT keyword_set(key_periodic) THEN BEGIN 
+      xminmesh = ixminmesh + 1
+      xmaxmesh = ixmaxmesh - 1
+    ENDIF
+    yminmesh = iyminmesh + 1
+    ymaxmesh = iymaxmesh - 1
+; come back to the original grid before calling computegrid with the
+; new parameters...
+    IF keyword_set(yreverse) THEN BEGIN
+      gphit = reverse(gphit, 2)
+      glamt = reverse(glamt, 2)
+      tmask = reverse(tmask, 2)    
+      IF jpk EQ 1 THEN tmask = reform(tmask, jpi, jpj, jpk, /over)
+    ENDIF
+    IF keyword_set(zreverse) THEN BEGIN
+      zaxis = reverse(zaxis)
+      tmask = reverse(tmask, 3)
+      IF jpk EQ 1 THEN tmask = reform(tmask, jpi, jpj, jpk, /over)
+    ENDIF
+    IF key_shift NE 0 THEN BEGIN
+      glamt = shift(glamt, -key_shift, 0)
+      gphit = shift(gphit, -key_shift, 0)
+      tmask = shift(tmask, -key_shift, 0, 0)
+      IF jpk EQ 1 THEN tmask = reform(tmask, jpi, jpj, jpk, /over)
+    ENDIF
+    IF ixminmesh NE 0 THEN BEGIN
+      glamt = [fltarr(ixminmesh, jpj), glamt]
+      gphit = [fltarr(ixminmesh, jpj), gphit]
+      tmask = [fltarr(ixminmesh, jpj), tmask]
+      jpi = jpi+ixminmesh
+      IF jpk EQ 1 THEN tmask = reform(tmask, jpi, jpj, jpk, /over)
+    ENDIF
+    IF ixmaxmesh NE jpiglo-1 THEN BEGIN
+      glamt = [glamt, fltarr(jpiglo-1-ixmaxmesh, jpj)]
+      gphit = [gphit, fltarr(jpiglo-1-ixmaxmesh, jpj)]
+      tmask = [tmask, fltarr(jpiglo-1-ixmaxmesh, jpj)]
+      jpi = jpi+jpiglo-1-ixmaxmesh
+      IF jpk EQ 1 THEN tmask = reform(tmask, jpi, jpj, jpk, /over)
+    ENDIF
+    IF iyminmesh NE 0 THEN BEGIN
+      glamt = [[fltarr(jpi, iyminmesh)], [glamt]]
+      gphit = [[fltarr(jpi, iyminmesh)], [gphit]]
+      tmask = [[fltarr(jpi, iyminmesh)], [tmask]]
+      jpj = jpj+iyminmesh
+      IF jpk EQ 1 THEN tmask = reform(tmask, jpi, jpj, jpk, /over)
+    ENDIF
+    IF iymaxmesh NE jpjglo-1 THEN BEGIN
+      glamt = [[glamt], [fltarr(jpi, jpjglo-1-iymaxmesh)]]
+      gphit = [[gphit], [fltarr(jpi, jpjglo-1-iymaxmesh)]]
+      tmask = [[tmask], [fltarr(jpi, jpjglo-1-iymaxmesh)]]
+      jpj = jpj+jpjglo-1-iymaxmesh
+      IF jpk EQ 1 THEN tmask = reform(tmask, jpi, jpj, jpk, /over)
+    ENDIF
+    computegrid, XAXIS = glamt, YAXIS = gphit, ZAXIS = zaxis $
+                 , MASK = tmask, GLAMBOUNDARY = glamboundary $
+                 , XMINMESH = xminmesh, XMAXMESH = xmaxmesh $
+                 , YMINMESH = yminmesh, YMAXMESH = ymaxmesh $
+                 , ZMINMESH = zminmesh, ZMAXMESH = zmaxmesh $
+                 , ONEARTH = onearth, PERIODIC = key_periodic $
+                 , PLAIN = plain, SHIFT = key_shift, STRIDE = key_stride $
+                 , FULLCGRID = fullcgrid, XYINDEX = xyindex $
+                 , UBASE2TBASE = ubase2tbase, VBASE2TBASE = vbase2tbase $
+                 , FBASE2TBASE = fbase2tbase, STRCALLING = strcalling $
+                 , ROMSH = romsh, _extra = ex
+    return
   ENDIF
 ;