Changeset 399

Timestamp:

07/02/09 10:38:54 (15 years ago)

Author:

smasson

Message:

wrong bugfix in previous changeset:398

Location:

trunk/SRC

Files:

• Grid/computegrid.pro (modified) (7 diffs)
• ReadWrite/ncdf_getaxis.pro (modified) (7 diffs)
• ReadWrite/ncdf_getmask.pro (modified) (3 diffs)
• ToBeReviewed/LECTURE/read_ncdf.pro (modified) (1 diff)

trunk/SRC/Grid/computegrid.pro

@@ -353 +353 @@
 ;====================================================
 ;
-  jpiglo = long(nx)
-  jpjglo = long(ny)
-  jpkglo = long(nz)
-  IF keyword_set(romsh) THEN BEGIN
-    jpiglo = jpiglo - 1
-    jpjglo = jpjglo - 1
-    fullcgrid = 1
-  ENDIF
+  IF keyword_set(romsh) THEN fullcgrid = 1
+  CASE 1 OF
+    keyword_set(fbase2tbase):key_gridtype = 'c_f'
+    keyword_set(ubase2tbase):key_gridtype = 'c_u'
+    keyword_set(vbase2tbase):key_gridtype = 'c_v'
+    else:key_gridtype = 'c'
+  ENDCASE
+  IF strlen(key_gridtype) EQ 3 THEN fullcgrid = 1
+;
+  IF n_elements(xminmesh) NE 0 AND n_elements(xmaxmesh) NE 0 THEN BEGIN
+    IF nx EQ jpi AND xminmesh EQ ixminmesh AND xmaxmesh EQ ixmaxmesh THEN xalreadycut = 1
+  ENDIF
+  IF keyword_set(xalreadycut) THEN BEGIN
+    xmin = 0 
+    xmax = jpi - 1 
+    nxx = jpi
+  ENDIF ELSE BEGIN
+    jpiglo = long(nx)
+    IF keyword_set(romsh) THEN jpiglo = jpiglo - 1    
+    IF n_elements(xminmesh) NE 0 THEN ixminmesh = long(xminmesh[0]) ELSE ixminmesh  = 0l
+    IF n_elements(xmaxmesh) NE 0 THEN ixmaxmesh = long(xmaxmesh[0]) ELSE ixmaxmesh  = jpiglo-1
+    IF ixmaxmesh LT 0 THEN ixmaxmesh = jpiglo -1 + ixmaxmesh
+    ixmaxmesh = 0 > ixmaxmesh < (jpiglo-1)
+    ixminmesh = 0 > ixminmesh < ixmaxmesh
+    jpi = ixmaxmesh-ixminmesh+1
+    xmin = ixminmesh 
+    xmax = ixmaxmesh
+    nxx = jpiglo
+  ENDELSE
+
+  IF n_elements(yminmesh) NE 0 AND n_elements(ymaxmesh) NE 0 THEN BEGIN
+    IF ny EQ jpj AND yminmesh EQ iyminmesh AND ymaxmesh EQ iymaxmesh THEN yalreadycut = 1
+  ENDIF
+  IF keyword_set(yalreadycut) THEN BEGIN
+    ymin = 0 
+    ymax = jpj - 1 
+    nyy = jpj
+  ENDIF ELSE BEGIN 
+    jpjglo = long(ny)
+    IF keyword_set(romsh) THEN jpjglo = jpjglo - 1
+    IF n_elements(yminmesh) NE 0 THEN iyminmesh = long(yminmesh[0]) ELSE iyminmesh  = 0l
+    IF n_elements(ymaxmesh) NE 0 THEN iymaxmesh = long(ymaxmesh[0]) ELSE iymaxmesh  = jpjglo-1
+    IF key_gridtype EQ 'c_v' OR key_gridtype EQ 'c_f' THEN iymaxmesh = iymaxmesh-1
+    IF iymaxmesh LT 0 THEN iymaxmesh = jpjglo -1 + iymaxmesh
+    iymaxmesh = 0 > iymaxmesh < (jpjglo-1)
+    iyminmesh = 0 > iyminmesh < iymaxmesh
+    jpj = iymaxmesh-iyminmesh+1
+    ymin = iyminmesh 
+    ymax = iymaxmesh
+    nyy = jpjglo
+  ENDELSE
+
+  IF n_elements(zminmesh) NE 0 AND n_elements(zmaxmesh) NE 0 THEN BEGIN
+    IF nz EQ jpk AND zminmesh EQ izminmesh AND zmaxmesh EQ izmaxmesh THEN zalreadycut = 1
+  ENDIF
+  IF keyword_set(zalreadycut) THEN BEGIN
+    zmin = 0 
+    zmax = jpk - 1 
+    nzz = jpk
+  ENDIF ELSE BEGIN 
+   jpkglo = long(nz)
+    IF n_elements(zminmesh) NE 0 THEN izminmesh = long(zminmesh[0]) ELSE izminmesh  = 0l
+    IF n_elements(zmaxmesh) NE 0 THEN izmaxmesh = long(zmaxmesh[0]) ELSE izmaxmesh  = jpkglo-1
+    IF izmaxmesh LT 0 THEN izmaxmesh = jpkglo -1 + izmaxmesh
+    izmaxmesh = 0 > izmaxmesh < (jpkglo-1)
+    izminmesh = 0 > izminmesh < izmaxmesh
+    jpk = izmaxmesh-izminmesh+1
+    zmin = izminmesh 
+    zmax = izmaxmesh
+    nzz = jpkglo
+  ENDELSE
 ;
 ; impact of plain keyword:
@@ -372 +435 @@
   ENDIF
 ;
-  IF n_elements(xminmesh) NE 0 THEN ixminmesh = long(xminmesh[0]) ELSE ixminmesh  = 0l
-  IF n_elements(xmaxmesh) NE 0 THEN ixmaxmesh = long(xmaxmesh[0]) ELSE ixmaxmesh  = jpiglo-1
-  IF n_elements(yminmesh) NE 0 THEN iyminmesh = long(yminmesh[0]) ELSE iyminmesh  = 0l
-  IF n_elements(ymaxmesh) NE 0 THEN iymaxmesh = long(ymaxmesh[0]) ELSE iymaxmesh  = jpjglo-1
-  IF n_elements(zminmesh) NE 0 THEN izminmesh = long(zminmesh[0]) ELSE izminmesh  = 0l
-  IF n_elements(zmaxmesh) NE 0 THEN izmaxmesh = long(zmaxmesh[0]) ELSE izmaxmesh  = jpkglo-1
-;
-  CASE 1 OF
-    keyword_set(fbase2tbase):key_gridtype = 'c_f'
-    keyword_set(ubase2tbase):key_gridtype = 'c_u'
-    keyword_set(vbase2tbase):key_gridtype = 'c_v'
-    else:key_gridtype = 'c'
-  ENDCASE
-  IF key_gridtype EQ 'c_v' OR key_gridtype EQ 'c_f' THEN BEGIN
-    iymaxmesh = iymaxmesh-1
-  ENDIF
-  IF strlen(key_gridtype) EQ 3 THEN fullcgrid = 1
-;
-  IF ixmaxmesh LT 0 THEN ixmaxmesh = jpiglo -1 + ixmaxmesh
-  IF iymaxmesh LT 0 THEN iymaxmesh = jpjglo -1 + iymaxmesh
-  IF izmaxmesh LT 0 THEN izmaxmesh = jpkglo -1 + izmaxmesh
 ; avoid basics errors...
-  ixmaxmesh = 0 > ixmaxmesh < (jpiglo-1)
-  ixminmesh = 0 > ixminmesh < ixmaxmesh
-  iymaxmesh = 0 > iymaxmesh < (jpjglo-1)
-  iyminmesh = 0 > iyminmesh < iymaxmesh
-  izmaxmesh = 0 > izmaxmesh < (jpkglo-1)
-  izminmesh = 0 > izminmesh < izmaxmesh
-;
-  jpi = ixmaxmesh-ixminmesh+1
-  jpj = iymaxmesh-iyminmesh+1
-  jpk = izmaxmesh-izminmesh+1
 ;
   jpidta = jpiglo
@@ -468 +500 @@
 ; force glamt to have 2 dimensions
 ;
-  CASE size(reform(glamt), /n_dimensions) OF
-    0:glamt = replicate(glamt, jpi, jpj)
-    1:glamt = glamt[ixminmesh:ixmaxmesh]#replicate(1, jpj)
-    2:glamt = glamt[ixminmesh:ixmaxmesh, iyminmesh:iymaxmesh]
+  IF n_elements(glamt) EQ nxx*nyy THEN glamt = reform(glamt, nxx, nyy, /over) $
+  ELSE glamt = reform(glamt, nxx, /over)
+  CASE size(glamt, /n_dimensions) OF
+    1:BEGIN 
+      IF n_elements(glamt) EQ 1 THEN glamt = replicate(glamt[0], jpi, jpj) $
+      ELSE glamt = glamt[xmin:xmax]#replicate(1, jpj)
+    END
+    2:glamt = glamt[xmin:xmax, ymin:ymax]
   ENDCASE
 ; keep 2d array even with degenerated dimension
@@ -498 +534 @@
 ; force gphit to have 2 dimensions
 ;
-  CASE size(reform(gphit), /n_dimensions) OF
-    0:gphit = replicate(gphit, jpi, jpj)
-    1:gphit = replicate(1, jpi)#gphit[iyminmesh:iymaxmesh]
-    2:gphit = gphit[ixminmesh:ixmaxmesh, iyminmesh:iymaxmesh]
+  IF n_elements(gphit) EQ nxx*nyy THEN gphit = reform(gphit, nxx, nyy, /over) $
+  ELSE gphit = reform(gphit, nyy, /over)
+  CASE size(gphit, /n_dimensions) OF
+    1:BEGIN 
+      IF n_elements(gphit) EQ 1 THEN gphit = replicate(gphit[0], jpi, jpj) $
+      ELSE gphit = replicate(1, jpi)#gphit[ymin:ymax]
+    END
+    2:gphit = gphit[xmin:xmax, ymin:ymax]
   ENDCASE
 ; keep 2d array even with degenerated dimension
@@ -990 +1030 @@
       END
       ELSE:BEGIN
-        gdept = zaxis[izminmesh:izmaxmesh]
+        gdept = zaxis[zmin:zmax]
         IF n_elements(zreverse) EQ 0 THEN BEGIN
           IF jpk GT 1 THEN BEGIN
@@ -1028 +1068 @@
 ;
   if tmask[0] NE -1 then BEGIN
+    tmask = byte(temporary(tmask))
     IF keyword_set(romsh) THEN tmask = tmask[0:jpiglo-1, 0:jpjglo-1]
-    IF n_elements(mask) EQ jpiglo*jpjglo AND jpkglo GT 1 THEN BEGIN
-      tmask = tmask[*]#replicate(1, jpkglo)
-      tmask = reform(tmask, jpiglo, jpjglo, jpkglo, /overwrite)
+    IF n_elements(mask) EQ nxx*nyy AND nzz GT 1 THEN BEGIN
+      tmask = tmask[*]#replicate(1b, nzz)
+      tmask = reform(tmask, nxx, nyy, nzz, /overwrite)
     ENDIF
-    IF jpiglo EQ 1 OR jpjglo EQ 1 THEN tmask = reform(tmask, jpiglo, jpjglo, jpkglo, /overwrite)
-    tmask = byte(tmask[ixminmesh:ixmaxmesh, iyminmesh:iymaxmesh, izminmesh:izmaxmesh])
-    tmask = reform(tmask, jpi, jpj, jpk, /over)
+    IF nxx EQ 1 OR nyy EQ 1 OR nzz EQ 1 THEN tmask = reform(tmask, nxx, nyy, nzz, /overwrite)
+    tmask = tmask[xmin:xmax, ymin:ymax, zmin:zmax]
+    IF jpi EQ 1 OR jpj EQ 1 OR jpk EQ 1 THEN tmask = reform(tmask, jpi, jpj, jpk, /over)
     if key_shift NE 0 then tmask = shift(tmask, key_shift, 0, 0)
 ; because tmask = reverse(tmask, 2) is not working if the 3rd
@@ -1242 +1283 @@
   ENDELSE
   IF keyword_set(romsh) THEN $
-     romszinfos = {h:romsh[ixminmesh:ixmaxmesh, iyminmesh:iymaxmesh], zeta:replicate(0., jpi, jpj), theta_s:-1, theta_b:-1, hc:-1} $
+     romszinfos = {h:romsh[xmin:xmax, ymin:ymax], zeta:replicate(0., jpi, jpj), theta_s:-1, theta_b:-1, hc:-1} $
   ELSE romszinfos = {h:-1, zeta:-1, theta_s:-1, theta_b:-1, hc:-1}
 

trunk/SRC/ReadWrite/ncdf_getaxis.pro

@@ -51 +51 @@
 ; x/yaxis = keyword_set(start1) + findgen(jpi/jpj)
 ;
+; @keyword XMINMESH {default=0L}{type=scalar}
+; Define common (cm_4mesh) variables ixminmesh used to define the localization
+; of the first point of the grid along the x direction in a zoom of the original grid
+;
+; @keyword YMINMESH {default=0L}{type=scalar}
+; Define common (cm_4mesh) variables iyminmesh used to define the localization
+; of the first point of the grid along the y direction in a zoom of the original grid
+;
+; @keyword XMAXMESH {default=jpiglo-1}{type=scalar}
+; Define common (cm_4mesh) variables ixmaxmesh used to define the localization
+; of the last point of the grid along the x direction in a zoom of the original grid
+; Note that if XMAXMESH < 0 then ixmaxmesh is defined as ixmaxmesh = jpiglo -1 + xmaxmesh
+;
+; @keyword YMAXMESH {default=jpjglo-1}{type=scalar}
+; Define common (cm_4mesh) variables iymaxmesh used to define the localization
+; of the last point of the grid along the y direction in a zoom of the original grid
+; Note that if YMAXMESH < 0 then iymaxmesh is defined as iymaxmesh = jpjglo -1 + ymaxmesh
+;
 ; @keyword _EXTRA
 ; Used to be able to call ncdf_getaxis with _extra
@@ -62 +80 @@
 ;-
 PRO ncdf_getaxis, fileid, dimidx, dimidy, xaxis, yaxis $
-                  , XAXISNAME=xaxisname, YAXISNAME=yaxisname $
-                  , XDIMNAME=xdimname, YDIMNAME=ydimname $
-                  , XYINDEX=xyindex, START1=start1 $
-                  , ROMSGRID=romsgrid, _EXTRA=ex
+                  , XAXISNAME = xaxisname, YAXISNAME = yaxisname $
+                  , XDIMNAME = xdimname, YDIMNAME = ydimname $
+                  , XYINDEX = xyindex, START1 = start1 $
+                  , XMINMESH = xminmesh, XMAXMESH = xmaxmesh $
+                  , YMINMESH = yminmesh, YMAXMESH = ymaxmesh $
+                  , ROMSGRID = romsgrid, _EXTRA = ex
 ;
   compile_opt idl2, strictarrsubs
 ;
-  @cm_general                      ;   needed for iodir
-;
+@cm_general                   ;   needed for iodir
+@cm_4mesh
+;
+  IF n_elements(xminmesh) NE 0 THEN ixminmesh = 0L > long(xminmesh[0]) ELSE ixminmesh  = 0l
+  IF n_elements(yminmesh) NE 0 THEN iyminmesh = 0L > long(yminmesh[0]) ELSE iyminmesh  = 0l
+  IF n_elements(zminmesh) NE 0 THEN izminmesh = 0L > long(zminmesh[0]) ELSE izminmesh  = 0l
 ; should we open the file?
   IF size(fileid, /type) EQ 7 THEN $
@@ -140 +164 @@
     xaxisname = xinq.name
     dimidx = xinq.dim[0]
-    IF xinq.ndims GE 2 THEN ncdf_diminq, cdfid, xinq.dim[1], blabla, jpjfromx
+    xoffset = lonarr(xinq.ndims)
+    xcount = lonarr(xinq.ndims)
+    FOR i = 0, xinq.ndims -1 DO BEGIN
+      ncdf_diminq, cdfid, xinq.dim[i], blabla, tmpsz
+      xcount[i] = tmpsz
+    ENDFOR
+    jpiglo = xcount[0]
+    IF n_elements(xmaxmesh) NE 0 THEN BEGIN 
+      IF xmaxmesh GE 0 THEN ixmaxmesh = long(xmaxmesh[0]) ELSE ixmaxmesh = xcount[0] -1 + long(xmaxmesh[0])
+    ENDIF ELSE ixmaxmesh = xcount[0] - 1L
+    ixmaxmesh  = 0 > ixmaxmesh < (xcount[0] - 1L) ; make sure ixmaxmesh is not too big
+    xoffset[0] = ixminmesh
+    xcount[0] = ixmaxmesh - ixminmesh + 1L
+    jpifromx = xcount[0]
+    jpi = jpifromx
+    IF xinq.ndims GE 2 THEN BEGIN 
+      jpjglo = xcount[1]
+      IF n_elements(ymaxmesh) NE 0 THEN BEGIN 
+        IF ymaxmesh GE 0 THEN iymaxmesh = long(ymaxmesh[0]) ELSE iymaxmesh = xcount[1] - 1L + long(ymaxmesh[0])
+      ENDIF ELSE iymaxmesh = xcount[1] - 1L
+      iymaxmesh  = 0 > iymaxmesh < (xcount[1] - 1L) ; make sure ixmaxmesh is not too big
+      xoffset[1] = iyminmesh
+      xcount[1] = iymaxmesh - iyminmesh + 1L
+      jpjfromx = xcount[1]
+      jpj = jpjfromx
+    ENDIF
   ENDELSE
-  IF arg_present(xdimname) THEN ncdf_diminq, cdfid, dimidx,  xdimname, jpifromx
+  IF arg_present(xdimname) THEN ncdf_diminq, cdfid, dimidx,  xdimname, tmp
 ;
   IF arg_present(xaxis) THEN BEGIN
-    ncdf_diminq, cdfid, dimidx, blabla, jpifromx
 ; should we read or compute the xaxis?
     IF keyword_set(xyindex) OR xvarid EQ -1 THEN BEGIN
@@ -152 +200 @@
     ENDIF ELSE BEGIN
 ; read the xaxis
-      ncdf_varget, cdfid, xvarid, xaxis
+      ncdf_varget, cdfid, xvarid, xaxis, offset = xoffset, count = xcount
       ncdf_getatt, cdfid, xvarid, ADD_OFFSET = add_offset, SCALE_FACTOR = scale_factor
       IF scale_factor NE 1. THEN xaxis = temporary(xaxis) * scale_factor
@@ -208 +256 @@
     yinq = ncdf_varinq(cdfid, yvarid)
     yaxisname = yinq.name
-    IF yinq.ndims GE 2 THEN BEGIN
-      ncdf_diminq, cdfid, yinq.dim[0], blabla, jpifromy
-      dimidy = yinq.dim[1]
-    ENDIF ELSE dimidy = yinq.dim[0]
+    IF yinq.ndims GE 2 THEN dimidy = yinq.dim[1] ELSE dimidy = yinq.dim[0]
+    yoffset = lonarr(yinq.ndims)
+    ycount = lonarr(yinq.ndims)
+    FOR i = 0, yinq.ndims -1 DO BEGIN
+      ncdf_diminq, cdfid, yinq.dim[i], blabla, tmpsz
+      ycount[i] = tmpsz
+    ENDFOR
+    idy = yinq.ndims GE 2
+    jpjglo = ycount[idy]
+    IF n_elements(ymaxmesh) NE 0 THEN BEGIN 
+      IF ymaxmesh GE 0 THEN iymaxmesh = long(ymaxmesh[0]) ELSE iymaxmesh = ycount[idy] - 1L + long(ymaxmesh[0])
+    ENDIF ELSE iymaxmesh = ycount[idy] - 1L
+    iymaxmesh  = iymaxmesh < (ycount[idy] - 1L) ; make sure ixmaxmesh is not too big
+    yoffset[idy] = iyminmesh
+    ycount[idy] = iymaxmesh - iyminmesh + 1L
+    jpjfromy = ycount[idy]
+    jpj = jpjfromy
+    IF yinq.ndims GE 2 THEN BEGIN 
+      jpiglo = ycount[0]
+      IF n_elements(xmaxmesh) NE 0 THEN BEGIN 
+        IF xmaxmesh GE 0 THEN ixmaxmesh = long(xmaxmesh[0]) ELSE ixmaxmesh = ycount[0] - 1L + long(xmaxmesh[0])
+      ENDIF ELSE ixmaxmesh = ycount[0] - 1L
+      ixmaxmesh  = ixmaxmesh < (ycount[0] - 1L) ; make sure ixmaxmesh is not too big
+      yoffset[0] = ixminmesh
+      ycount[0] = ixmaxmesh - ixminmesh + 1L
+      jpifromy = xcount[0]
+      jpi = jpifromy
+    ENDIF
   ENDELSE
-  IF arg_present(ydimname) THEN ncdf_diminq, cdfid, dimidy,  ydimname, jpjfromy
+  IF arg_present(ydimname) THEN ncdf_diminq, cdfid, dimidy,  ydimname, tmp
 ;
   IF arg_present(yaxis) THEN BEGIN
     IF n_elements(jpifromy) NE 0 THEN BEGIN
+      IF n_elements(jpifromx) EQ 0 THEN ncdf_diminq, cdfid, dimidx, blabla, jpifromx
       IF jpifromy NE jpifromx THEN BEGIN
         dummy = report('x/y axes do not have the same x dimension...')
@@ -222 +295 @@
       ENDIF
     ENDIF
-    ncdf_diminq, cdfid, dimidy, blabla, jpjfromy
     IF n_elements(jpjfromx) NE 0 THEN BEGIN
+      IF n_elements(jpjfromy) EQ 0 THEN ncdf_diminq, cdfid, dimidy, blabla, jpjfromy
       IF jpjfromy NE jpjfromx THEN BEGIN
         dummy = report(' x/y axes do not have the same y dimension...')
@@ -234 +307 @@
     ENDIF ELSE BEGIN
 ; read the yaxis
-      ncdf_varget, cdfid, yvarid, yaxis
+      ncdf_varget, cdfid, yvarid, yaxis, offset = yoffset, count = ycount
       ncdf_getatt, cdfid, yvarid, ADD_OFFSET = add_offset, SCALE_FACTOR = scale_factor
       IF scale_factor NE 1. THEN yaxis = temporary(yaxis) * scale_factor

trunk/SRC/ReadWrite/ncdf_getmask.pro

@@ -56 +56 @@
 ;   3) !Values.f_nan (can be used only with NE and EQ operators)
 ;
+; @keyword XMINMESH {default=0L}{type=scalar}
+;       Define common (cm_4mesh) variables ixminmesh used to define the localization
+;       of the first point of the grid along the x direction in a zoom of the original grid
+;
+; @keyword YMINMESH {default=0L}{type=scalar}
+;       Define common (cm_4mesh) variables iyminmesh used to define the localization
+;       of the first point of the grid along the y direction in a zoom of the original grid
+;
+; @keyword ZMINMESH {default=0L}{type=scalar}
+;       Define common (cm_4mesh) variables izminmesh used to define the localization
+;       of the first point of the grid along the z direction in a zoom of the original grid
+;
+; @keyword XMAXMESH {default=jpiglo-1}{type=scalar}
+;       Define common (cm_4mesh) variables ixmaxmesh used to define the localization
+;       of the last point of the grid along the x direction in a zoom of the original grid
+;       Note that if XMAXMESH < 0 then ixmaxmesh is defined as ixmaxmesh = jpiglo -1 + xmaxmesh
+;
+; @keyword YMAXMESH {default=jpjglo-1}{type=scalar}
+;       Define common (cm_4mesh) variables iymaxmesh used to define the localization
+;       of the last point of the grid along the y direction in a zoom of the original grid
+;       Note that if YMAXMESH < 0 then iymaxmesh is defined as iymaxmesh = jpjglo -1 + ymaxmesh
+;
+; @keyword ZMAXMESH {default=jpkglo-1}{type=scalar}
+;       Define common (cm_4mesh) variables izmaxmesh used to define the localization
+;       of the last point of the grid along the z direction in a zoom of the original grid
+;       Note that if ZMAXMESH < 0 then izmaxmesh is defined as izmaxmesh = jpkglo -1 + maxmesh
+;
 ; @keyword 
 ; _EXTRA to be able to call ncdf_getmask with _extra keyword
@@ -80 +107 @@
                        , MASKNAME = maskname, USEASMASK = useasmask $
                        , MISSING_VALUE = missing_value, INVMASK = invmask $
+                       , XMINMESH = xminmesh, XMAXMESH = xmaxmesh $
+                       , YMINMESH = yminmesh, YMAXMESH = ymaxmesh $
+                       , ZMINMESH = zminmesh, ZMAXMESH = zmaxmesh $
                        , TESTOP = testop, _EXTRA = ex
 ;
   compile_opt idl2, strictarrsubs
 ;
-  @cm_general                      ;   needed for iodir
+  @cm_general                   ;   needed for iodir
+;
+  IF n_elements(xminmesh) NE 0 THEN ixminmesh = 0L > long(xminmesh[0]) ELSE ixminmesh  = 0l
+  IF n_elements(yminmesh) NE 0 THEN iyminmesh = 0L > long(yminmesh[0]) ELSE iyminmesh  = 0l
+  IF n_elements(zminmesh) NE 0 THEN izminmesh = 0L > long(zminmesh[0]) ELSE izminmesh  = 0l
 ;
   IF NOT (keyword_set(maskname) OR keyword_set(useasmask)) AND keyword_set(romsgrid) THEN maskname = 'mask_rho'
@@ -111 +145 @@
 ; is the mask variable containing the record dimension?
     withrcd = (where(mskinq.dim EQ inq.recdim))[0]
-    IF withrcd NE -1 THEN BEGIN
 ; in order to read only the first record
 ; we need to get the size of each dimension
-      count = replicate(1L, mskinq.ndims)
-      FOR d = 0, mskinq.ndims -1 DO BEGIN
-        IF d NE withrcd THEN BEGIN
-          ncdf_diminq, cdfid, mskinq.dim[d], name, size
-          count[d] = size
-        ENDIF
-      ENDFOR
+    offset = lonarr(mskinq.ndims)
+    count = replicate(1L, mskinq.ndims)
+    FOR d = 0, mskinq.ndims -1 DO BEGIN
+      IF d NE withrcd THEN BEGIN
+        ncdf_diminq, cdfid, mskinq.dim[d], name, size
+        count[d] = size
+      ENDIF
+    ENDFOR
+    
+    IF n_elements(xmaxmesh) NE 0 THEN BEGIN 
+      IF xmaxmesh GE 0 THEN ixmaxmesh = long(xmaxmesh[0]) ELSE ixmaxmesh = count[0] - 1L + long(xmaxmesh[0])
+    ENDIF ELSE ixmaxmesh = count[0] - 1L
+    ixmaxmesh  = 0 > ixmaxmesh < (count[0] - 1L) ; make sure ixmaxmesh is not too big
+    jpiglo = count[0]
+    offset[0] = ixminmesh
+    count[0] = ixmaxmesh - ixminmesh + 1L
+    jpi = count[0]
+
+    IF n_elements(ymaxmesh) NE 0 THEN  BEGIN 
+      IF ymaxmesh GE 0 THEN iymaxmesh = long(ymaxmesh[0]) ELSE iymaxmesh = count[1] - 1L + long(ymaxmesh[0])
+    ENDIF ELSE iymaxmesh = count[1] - 1L 
+    iymaxmesh  = 0 > iymaxmesh < (count[1] - 1L) ; make sure ixmaxmesh is not too big
+    jpjglo = count[1]
+    offset[1] = iyminmesh
+    count[1] = iymaxmesh - iyminmesh + 1L
+    jpj = count[1]
+
+    IF (mskinq.ndims EQ 3 AND withrcd EQ -1) OR (mskinq.ndims EQ 4) THEN BEGIN
+      IF n_elements(zmaxmesh) NE 0 THEN  BEGIN 
+        IF zmaxmesh GE 0 THEN izmaxmesh = long(zmaxmesh[0]) ELSE izmaxmesh = count[2] - 1L + long(zmaxmesh[0])
+      ENDIF ELSE izmaxmesh = count[2] - 1L
+      izmaxmesh  = 0 > izmaxmesh < (count[2] - 1L) ; make sure ixmaxmesh is not too big
+      jpkglo = count[2]
+      offset[2] = izminmesh
+      count[2] = izmaxmesh - izminmesh + 1L
+      jpk = count[2]
+    ENDIF
+
 ; read the variable for the first record
-      ncdf_varget, cdfid, mskid, mask, count = count
-    ENDIF ELSE ncdf_varget, cdfid, mskid, mask
-
+    ncdf_varget, cdfid, mskid, mask, count = count,  offset = offset
+    
 ; get add_offset, scale factor and missing value attributes
     ncdf_getatt, cdfid, mskid, add_offset = add, scale_factor = scl, missing_value = miss

trunk/SRC/ToBeReviewed/LECTURE/read_ncdf.pro

@@ -340 +340 @@
 ; vardate
 ; We make a legible date
-  caldat, time[0], year, month, day
+  caldat, time[0], month, day, year
   vardate = string(format = '(C(CMoA))', 31*(month-1))+' '+strtrim(day, 1)+', '+strtrim(year, 1)
   varexp = file_basename(filename)