Changeset 219 for trunk/SRC/ToBeReviewed/WIDGET

Timestamp:

03/14/07 16:44:26 (17 years ago)

Author:

smasson

Message:

get back to changeset:217

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 {default='T'}{type=scalar string}
-; Used to specify on which grid type are located the data
+; @keyword GRID
+;
 ;
 ; @keyword _EXTRA
-; Used to pass your keywords to isafile and ncdf_getaxis
+; Used to pass your keywords
 ; 
 ; @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
@@ -87 +88 @@
 ; What contains the file?
 ;------------------------------------------------------------
-  inside = ncdf_inquire(cdfid)  ;
+  infile = ncdf_inquire(cdfid)  ;
 ;------------------------------------------------------------
 ; name of all dimensions
 ;------------------------------------------------------------
-  namedim = strarr(inside.ndims)
-  for dimiq = 0, inside.ndims-1 do begin
+  namedim = strarr(infile.ndims)
+  for dimiq = 0, infile.ndims-1 do begin
     ncdf_diminq, cdfid, dimiq, tmpname, value 
     namedim[dimiq] = strlowcase(tmpname)
   ENDFOR
-;------------------------------------------------------------
-; x/y dimensions id
-;------------------------------------------------------------
-  ncdf_getaxis, cdfid, dimidx, dimidy, _extra = ex
+; 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
 ;------------------------------------------------------------
 ; name of all variables
 ;------------------------------------------------------------
 ; we keep only the variables containing at least x, y and time dimension (if existing...)
-  namevar = strarr(inside.nvars)
-  for varid = 0, inside.nvars-1 do begin
+  namevar = strarr(infile.nvars)
+  for varid = 0, infile.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 inside.recdim))[0] NE -1 OR inside.recdim EQ -1) $
+       ((where(invar.dim EQ infile.recdim))[0] NE -1 OR infile.recdim EQ -1) $
     THEN namevar[varid] = invar.name 
   ENDFOR
@@ -159 +210 @@
 ;------------------------------------------------------------
   date0fk = date2jul(19000101)
-  IF inside.recdim EQ -1 THEN BEGIN 
+  IF infile.recdim EQ -1 THEN BEGIN 
     jpt = 1
     time = date0fk
     fakecal = 1
   ENDIF ELSE BEGIN
-    ncdf_diminq, cdfid, inside.recdim, timedimname, jpt
+    ncdf_diminq, cdfid, infile.recdim, timedimname, jpt
 ; we look for the variable containing the time axis
-; we look for the first variable having for only dimension inside.recdim
+; we look for the first variable having for only dimension infile.recdim
     varid = 0
     repeat BEGIN
-      IF varid LT inside.nvars THEN BEGIN
+      IF varid LT infile.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 inside.recdim) OR (varid EQ 0)
+    endrep until (n_elements(invar.dim) EQ 1 AND invar.dim[0] EQ infile.recdim) OR (varid EQ 0)
     varid = varid-1
 ;