Changeset 221 for trunk/SRC/ToBeReviewed/WIDGET

Timestamp:

03/14/07 17:00:19 (17 years ago)

Author:

smasson

Message:

add x/ydimname and cleaning to solve ticket:61

File:

• trunk/SRC/ToBeReviewed/WIDGET/AUTOUR_de_XXX/scanfile.pro (modified) (4 diffs)

trunk/SRC/ToBeReviewed/WIDGET/AUTOUR_de_XXX/scanfile.pro

@@ -9 +9 @@
 ;
 ; 
-; @keyword GRID
-;
+; @keyword GRID {default='T'}{type=scalar string}
+; Used to specify on which grid type are located the data
 ;
 ; @keyword _EXTRA
-; Used to pass your keywords
+; Used to pass your keywords to isafile and ncdf_getaxis
 ; 
 ; @returns
@@ -76 +76 @@
 @common
 ;------------------------------------------------------------
-  res = -1
-;------------------------------------------------------------
 ; filename
 ;------------------------------------------------------------
   fullname = isafile(filename = namefile, IODIRECTORY = iodir, _extra = ex)
+  IF size(fullname, /type) NE 7 THEN return, -1
 ;------------------------------------------------------------
 ; open file
@@ -88 +87 @@
 ; What contains the file?
 ;------------------------------------------------------------
-  infile = ncdf_inquire(cdfid)  ;
+  inside = ncdf_inquire(cdfid)  ;
 ;------------------------------------------------------------
 ; name of all dimensions
 ;------------------------------------------------------------
-  namedim = strarr(infile.ndims)
-  for dimiq = 0, infile.ndims-1 do begin
+  namedim = strarr(inside.ndims)
+  for dimiq = 0, inside.ndims-1 do begin
     ncdf_diminq, cdfid, dimiq, tmpname, value 
     namedim[dimiq] = strlowcase(tmpname)
   ENDFOR
-; roms file?
-  dimidx = where(namedim EQ 'xi_rho' OR namedim EQ 'xi_u' OR namedim EQ 'xi_v' OR namedim EQ 'xi_psi')
-  IF dimidx[0] EQ -1 THEN BEGIN
-; we are looking for a x dimension with a name matching one of the following regular expression:
-    testname = ['longitude', 'lon', 'x', 'longitude*', 'lon*', 'x*', '*longitude*', '*lon*', '*x*']
-    cnt = -1
-    ii = 0
-    WHILE cnt NE 1 AND ii LT n_elements(testname) DO BEGIN
-      dimidx = where(strmatch(namedim, testname[ii]) EQ 1, cnt)
-      ii = ii+1
-    ENDWHILE
-    CASE cnt OF
-      0:begin
-        dummy = report(['none of the dimensions name matches one of the following regular expression:' $
-                        , '''longitude'', ''lon'', ''x'', ''longitude*'', ''lon*'', ''x*'', ''*longitude*'', ''*lon*'', ''*x*''' $
-                        , ' => we cannot find the x dimension'])
-        stop
-      END
-      1:dimidx = dimidx[0]
-      ELSE:begin
-        dummy = report(['several (and not one unique) dimensions name matches the following regular expression:' $
-                        , '''longitude'', ''lon'', ''x'', ''longitude*'', ''lon*'', ''x*'', ''*longitude*'', ''*lon*'', ''*x*''' $
-                        , ' => we cannot find the x dimension'])
-        stop
-      ENDELSE
-    ENDCASE
-  ENDIF
-; roms file?
-  dimidy = where(namedim EQ 'eta_rho' OR namedim EQ 'eta_u' OR namedim EQ 'eta_v' OR namedim EQ 'eta_psi')
-  IF dimidy[0] EQ -1 THEN BEGIN
-; we are looking for a y dimension with a name matching one of the following regular expression:
-    testname = ['latitude', 'lat', 'y', 'latitude*', 'lat*', 'y*', 'eta_*', '*latitude*', '*lat*', '*y*']
-    cnt = -1
-    ii = 0
-    WHILE cnt NE 1 AND ii LT n_elements(testname) DO BEGIN
-      dimidy = where(strmatch(namedim, testname[ii]) EQ 1, cnt)
-      ii = ii+1
-    ENDWHILE
-    CASE cnt OF
-      0:begin
-        dummy = report(['none of the dimensions name matches one of the following regular expression:' $
-                        , '''latitude'', ''lat'', ''y'', ''latitude*'', ''lat*'', ''y*'', ''eta_*'', ''*latitude*'', ''*lat*'', ''*y*''' $
-                        , ' => we cannot find the y dimension'])
-        stop
-      END
-      1:dimidy = dimidy[0]
-      ELSE:begin
-        dummy = report(['several (and not one unique) dimensions name matches the following regular expression:' $
-                        , '''latitude'', ''lat'', ''y'', ''latitude*'', ''lat*'', ''y*'', ''eta_*'', ''*latitude*'', ''*lat*'', ''*y*''' $
-                        , ' => we cannot find the x dimension'])
-        stop
-      ENDELSE
-    ENDCASE
-  ENDIF
+;------------------------------------------------------------
+; x/y dimensions id
+;------------------------------------------------------------
+  ncdf_getaxis, cdfid, dimidx, dimidy, _extra = ex
 ;------------------------------------------------------------
 ; name of all variables
 ;------------------------------------------------------------
 ; we keep only the variables containing at least x, y and time dimension (if existing...)
-  namevar = strarr(infile.nvars)
-  for varid = 0, infile.nvars-1 do begin
+  namevar = strarr(inside.nvars)
+  for varid = 0, inside.nvars-1 do begin
     invar = ncdf_varinq(cdfid, varid) ; what contains the variable?
     if (inter(invar.dim, dimidx))[0] NE -1 AND $
        (inter(invar.dim, dimidy))[0] NE -1 AND $
-       ((where(invar.dim EQ infile.recdim))[0] NE -1 OR infile.recdim EQ -1) $
+       ((where(invar.dim EQ inside.recdim))[0] NE -1 OR inside.recdim EQ -1) $
     THEN namevar[varid] = invar.name 
   ENDFOR
@@ -210 +159 @@
 ;------------------------------------------------------------
   date0fk = date2jul(19000101)
-  IF infile.recdim EQ -1 THEN BEGIN 
+  IF inside.recdim EQ -1 THEN BEGIN 
     jpt = 1
     time = date0fk
     fakecal = 1
   ENDIF ELSE BEGIN
-    ncdf_diminq, cdfid, infile.recdim, timedimname, jpt
+    ncdf_diminq, cdfid, inside.recdim, timedimname, jpt
 ; we look for the variable containing the time axis
-; we look for the first variable having for only dimension infile.recdim
+; we look for the first variable having for only dimension inside.recdim
     varid = 0
     repeat BEGIN
-      IF varid LT infile.nvars THEN BEGIN
+      IF varid LT inside.nvars THEN BEGIN
         invar = ncdf_varinq(cdfid, varid) 
         varid = varid+1        
       ENDIF ELSE varid = 0
-    endrep until (n_elements(invar.dim) EQ 1 AND invar.dim[0] EQ infile.recdim) OR (varid EQ 0)
+    endrep until (n_elements(invar.dim) EQ 1 AND invar.dim[0] EQ inside.recdim) OR (varid EQ 0)
     varid = varid-1
 ;