Changeset 49 for trunk/ToBeReviewed/WIDGET/AUTOUR_de_XXX/scanfile.pro

Timestamp:

05/09/06 11:47:44 (18 years ago)

Author:

pinsard

Message:

upgrade of AUTOUR_de_XXX according to cerbere.lodyc.jussieu.fr: /usr/home/smasson/SAXO_RD/ : files

File:

• trunk/ToBeReviewed/WIDGET/AUTOUR_de_XXX/scanfile.pro (copied) (copied from trunk/WIDGET/AUTOUR_de_XXX/scanfile.pro) (5 diffs)

trunk/ToBeReviewed/WIDGET/AUTOUR_de_XXX/scanfile.pro

@@ -28 +28 @@
 ;
 ;------------------------------------------------------------
-FUNCTION scanfile, nomficher, IODIRECTORY = iodirectory
+FUNCTION scanfile, nomficher, _extra = ex
 @common
 ;------------------------------------------------------------
@@ -34 +34 @@
 ; declarees ds un common)
 ;------------------------------------------------------------
-   res = -1
+  res = -1
 ;------------------------------------------------------------
 ; choix du nom du fichier
 ;------------------------------------------------------------
-   if NOT keyword_set(iodirectory) then iodirectory = iodir
-   nom = isafile(filename = nomficher, IODIRECTORY = iodirectory)
+  nom = isafile(filename = nomficher, IODIRECTORY = iodir, _extra = ex)
 ;------------------------------------------------------------
 ; ouverture du fichier nom
 ;------------------------------------------------------------
-   cdfid=ncdf_open(nom)
+  cdfid = ncdf_open(nom)
 ;------------------------------------------------------------
 ; que contient le fichier??
 ;------------------------------------------------------------
-   contient=ncdf_inquire(cdfid)
-   vargrid = strupcase(strmid(nom, strlen(nom)-9))
-   if vargrid EQ 'GRID.T.NC' OR vargrid EQ 'GRID.U.NC' $
-    OR vargrid EQ 'GRID.V.NC' OR vargrid EQ 'GRID.F.NC' $
-    OR vargrid eq 'GRID.W.NC' $
-    OR vargrid EQ 'GRID_T.NC' OR vargrid EQ 'GRID_U.NC' $
-    OR vargrid EQ 'GRID_V.NC' OR vargrid EQ 'GRID_F.NC' $
-    OR vargrid eq 'GRID_W.NC' then $
-    vargrid = strupcase(strmid(nom, strlen(nom)-4, 1)) ELSE vargrid='T'
+  contient = ncdf_inquire(cdfid);
+; find vargrid ...
+  vargrid = 'T'                 ; default definition
+  pattern = ['GRID.', 'GRID_', 'GRID', 'UPID_', '30ID_']
+  gdtype = ['T', 'U', 'V', 'W', 'F']
+  fnametest = strupcase(nom)
+  FOR i = 0, n_elements(pattern)-1 DO BEGIN
+    FOR j = 0, n_elements(gdtype)-1 DO BEGIN
+      substr = pattern[i]+gdtype[j]
+      pos = strpos(fnametest, substr)
+      IF pos NE -1 THEN $
+        vargrid = strmid(fnametest, pos+strlen(substr)-1, 1)
+    ENDFOR
+  ENDFOR
 ;------------------------------------------------------------
 ; nom des differentes variables
@@ -62 +66 @@
 ; contient au moins les dimensions qui doivent etre appelles x et y
 ; et la dimension infinie
-   nomdim   =strarr(contient.ndims)
-   for dimiq=0,contient.ndims-1 do begin
-      ncdf_diminq,cdfid,dimiq,name,value ; nom et valeur de la dimension
-      nomdim[dimiq]=strlowcase(name)
-   ENDFOR
-   indexdimx = where(nomdim EQ 'x' OR nomdim EQ 'lon' OR nomdim EQ 'longitude')
-   indexdimx = indexdimx[0]
-   if indexdimx EQ -1 then begin
-      print, 'one of the dimensions must have the name: ''x'' or ''lon'' or ''longitude'''
-      stop
-   endif
-   indexdimy = where(nomdim EQ 'y' OR nomdim EQ 'lat' OR nomdim EQ 'latitude')
-   indexdimy = indexdimy[0]
-   if indexdimy EQ -1 then begin
-      print, 'one of the dimensions must have the name: ''y'' or ''lat'' or ''latitude'''
-      stop
-   endif
+  nomdim   = strarr(contient.ndims)
+  for dimiq = 0, contient.ndims-1 do begin
+    ncdf_diminq, cdfid, dimiq, name, value ; nom et valeur de la dimension
+    nomdim[dimiq] = strlowcase(name)
+  ENDFOR
+  indexdimx = where(nomdim EQ 'x' OR nomdim EQ 'lon' OR nomdim EQ 'longitude' OR nomdim EQ 'xt_i7_156' OR nomdim EQ 'xi_rho' OR nomdim EQ 'xi_u' OR nomdim EQ 'xi_v' OR nomdim EQ 'xi_psi')
+  indexdimx = indexdimx[0]
+  if indexdimx EQ -1 then begin
+    print, 'one of the dimensions must have the name: ''x'' or ''lon'' or ''longitude'' or ''xt_i7_156'''
+    stop
+  endif
+  indexdimy = where(nomdim EQ 'y' OR nomdim EQ 'lat' OR nomdim EQ 'latitude' OR nomdim EQ 'yt_j6_75' OR nomdim EQ 'eta_rho' OR nomdim EQ 'eta_u' OR nomdim EQ 'eta_v' OR nomdim EQ 'eta_psi')
+  indexdimy = indexdimy[0]
+  if indexdimy EQ -1 then begin
+    print, 'one of the dimensions must have the name: ''y'' or ''lat'' or ''latitude'' or ''yt_j6_75'''
+    stop
+  endif
 ;
-   namevar = strarr(contient.nvars)
-   for varid=0,contient.nvars-1 do begin
-      varcontient=ncdf_varinq(cdfid,varid) ; que contient la variable
-      if (where(varcontient.dim EQ indexdimx))[0] NE -1 AND $
-       (where(varcontient.dim EQ indexdimy))[0] NE -1 AND $
-       (where(varcontient.dim EQ contient.recdim))[0] NE -1 THEN namevar[varid]=varcontient.name 
-   ENDFOR
-   namevar = namevar[where(namevar NE '')]
-   listgrid = replicate(vargrid,  n_elements(namevar))
+  namevar = strarr(contient.nvars)
+  for varid = 0, contient.nvars-1 do begin
+    varcontient = ncdf_varinq(cdfid, varid) ; que contient la variable
+    if (where(varcontient.dim EQ indexdimx))[0] NE -1 AND $
+      (where(varcontient.dim EQ indexdimy))[0] NE -1 AND $
+      (where(varcontient.dim EQ contient.recdim))[0] NE -1 THEN namevar[varid] = varcontient.name 
+  ENDFOR
+  namevar = namevar[where(namevar NE '')]
+  listgrid = replicate(vargrid,  n_elements(namevar))
 ;------------------------------------------------------------
 ; on recupere l''axe des temps
 ;------------------------------------------------------------
 ; on recupere le nom de la variable contenant l''axe des temps
-   varid = 0
-   repeat BEGIN                 ; tant que l''on a pas trouve une variable qui n''a qu''
+  varid = 0
+  repeat BEGIN ; tant que l''on a pas trouve une variable qui n''a qu''
                                 ; une dimension: la dimension infinie
-      varcontient=ncdf_varinq(cdfid,varid) ; que contient la variable
-      varid = varid+1
-   endrep until n_elements(varcontient.dim) EQ 1 AND varcontient.dim[0] EQ contient.recdim
+    varcontient = ncdf_varinq(cdfid, varid) ; que contient la variable
+    varid = varid+1
+  endrep until n_elements(varcontient.dim) EQ 1 AND varcontient.dim[0] EQ contient.recdim ;
+  varid = varid-1
+;
+; we want to know which attrributes are attached to the time variable...
+;
+  if varcontient.natts EQ 0 then BEGIN
+    ncdf_close, cdfid
+    return, report('the variable '+varcontient.name+' has no attribut.!C we need attribut ''units'' for the calendar ...')
+  endif
+  attnames = strarr(varcontient.natts)
+  for attiq = 0, varcontient.natts-1 do attnames[attiq] = ncdf_attname(cdfid, varid, attiq)
+  if (where(attnames EQ 'units'))[0] EQ -1 then BEGIN
+    ncdf_close, cdfid
+    return, report('Attribut ''units'' not found for the variable '+varid.name+'we need it for the calendar...')
+  endif
 ; on lit l''axe des temps
-   varid = varid-1
-   ncdf_varget, cdfid, varid, time
-   time = long(time)
-   ncdf_attget,cdfid,varid,'units',value
+  ncdf_varget, cdfid, varid, time
+  time = double(time)
+  ncdf_attget, cdfid, varid, 'units', value
 ; time_counter:units = "seconds since 0001-01-01 00:00:00" ;
 ; time_counter:units = "hours since 0001-01-01 00:00:00" ;
@@ -109 +126 @@
 ; time_counter:units = "months since 1979-01-01 00:00:00" ;
 ; time_counter:units = "years since 1979-01-01 00:00:00" ;
-   value = strtrim(strcompress(string(value)), 2)
-   mots = str_sep(value, ' ')
-   unite = mots[0]
-   debut = str_sep(mots[2], '-')
+  value = strtrim(strcompress(string(value)), 2)
+  mots = str_sep(value, ' ')
+  unite = mots[0]
+  debut = str_sep(mots[2], '-')
 ;
+; now we try to find the attribut called calendar...
+; the the attribute "calendar" exists?
+; If no, we suppose that the calendar is gregorian calendar
 ;
-; ATTENTION il faut recuperer l''atribut calendar et ajuster time en
+  if (where(attnames EQ 'calendar'))[0] NE -1 then BEGIN
+      ncdf_attget, cdfid, varid, 'calendar', value
+      value = string(value)
+        CASE value OF
+          '360d':key_caltype = '360d'
+          'greg':IF n_elements(key_caltype) EQ 0 THEN key_caltype = 'greg'
+          ELSE:BEGIN
+;            notused = report('Unknown calendar: '+value+', we use greg calendar.') 
+            key_caltype = 'greg'
+          END
+        ENDCASE
+      ENDIF ELSE BEGIN
+;        notused = report('Unknown calendar, we use '+key_caltype+' calendar.') 
+        IF n_elements(key_caltype) EQ 0 THEN key_caltype = 'greg'
+      ENDELSE
+;
+; ATTENTION il faut recuperer l''attribut calendar et ajuster time en
 ; consequense...
 ;
@@ -122 +158 @@
 ; on suppose qu''on utilise le vrai calendrier.
 ;
-   case strlowcase(unite) of
-      'seconds':time = julday(debut[1], debut[2], debut[0])+time/(long(24)*3600)
-      'hours':time = julday(debut[1], debut[2], debut[0])+time/(long(24))
-      'days':time = julday(debut[1], debut[2], debut[0])+time
-      'months':BEGIN 
-         for t = 0, n_elements(time)-1  do begin
-            time[t] = julday(debut[1]+time[t], debut[2], debut[0])
-         endfor
-      END
-      'years':BEGIN
-         for t = 0, n_elements(time)-1  do begin
-            time[t] = julday(debut[1], debut[2], debut[0]+time[t])
-         endfor
-      END
-   ENDCASE
+  unite = strlowcase(unite)
+  IF strpos(unite, 's', strlen(unite)-1) NE -1 THEN unite = strmid(unite, 0, strlen(unite)-1)
+  IF strpos(unite, 'julian_') NE -1 THEN unite = strmid(unite, 7)
+  case unite of
+    'second':time = julday(debut[1], debut[2], debut[0])+time/(long(24)*3600)
+    'hour':time = julday(debut[1], debut[2], debut[0])+time/(long(24))
+    'day':time = julday(debut[1], debut[2], debut[0])+time
+    'month':BEGIN 
+      if total(fix(time) NE time) NE 0 then $ ; we switch to days with 30d/m
+        time = julday(debut[1], debut[2], debut[0])+round(time*30) $
+      ELSE for t = 0, n_elements(time)-1 DO $
+        time[t] = julday(debut[1]+time[t], debut[2], debut[0])
+    END
+    'year':BEGIN
+      if total(fix(time) NE time) NE 0 then $ ; we switch to days with 365d/y
+        time = julday(debut[1], debut[2], debut[0])+round(time*365) $
+      ELSE for t = 0, n_elements(time)-1 do $
+        time[t] = julday(debut[1], debut[2], debut[0]+time[t])
+    END
+  ENDCASE
+  time = long(time)
 ;------------------------------------------------------------
-   return, {filename:nom, time_counter:time, listvar:namevar, listgrid:strupcase(listgrid)}
+  return, {filename:nom, time_counter:time, listvar:namevar, listgrid:strupcase(listgrid), calendartype:key_caltype}
 end