Changeset 474 for trunk

Timestamp:

12/16/11 18:36:23 (12 years ago)

Author:

pinsard

Message:

try to fix usage of file reading

Location:

trunk/src

Files:

• cresamsu.pro (modified) (11 diffs)
• idl_amsu_netcdf.pro (modified) (1 diff)

trunk/src/cresamsu.pro

@@ -50 +50 @@
 ; ====
 ;
+; fix pb with fov_int2 and fov_int3
+;
 ; check last line for uncomment endfor (not sure of the use of this line)
+;
+; aligner lal/lon de l'exemple sur les autres modules amsu
 ;
 ; lelod 2011/12/15
 ; modifier la creation du ncdf pour ajouter la variable "mask"
-; 
+;
 ; lelod 2011/10/24
 ; actuellement codes en dur
@@ -64 +68 @@
 ; past (temporelle, en jour decimal)
 ; et enfin l'heure d'interpolation heurinterp (en heures)
-; 
+;
 ; coding rules
 ;
@@ -78 +82 @@
 ;  $URL$
 ;
+; - pinsard 20111216T172940Z
+;
+;   * start to fix usage of result of file_amsu_t2_to_mem
+;
 ; - fplod 20111216T164041Z cratos (Linux)
 ;
+;   * change the name of the procedure to avoid ambiguity
 ;   * fix file_amsu_t2_to_mem call
 ;   * start to forgot idldoc
 ;   * fix for loop syntax error
 ;
-; - lelod 20121215 
+; - lelod 20121215
 ;
 ;   * reprise du programme apres refonte de la chaine (extract_amsua):
@@ -101 +110 @@
 ;     pixels en bordure de fauchee de facon correcte
 ;
-;
 ; - fplod 20110505T123230Z cratos.locean-ipsl.upmc.fr (Linux)
 ;
@@ -115 +123 @@
 
 ; grille de sortie en long / lat
-; recoit en entree les jours de l'annee jdeb et jfin correspondant au debut et la fin du mois 
+; recoit en entree les jours de l'annee jdeb et jfin correspondant au debut et la fin du mois
 ;=============================================================================================
 
@@ -163 +171 @@
 endif
 hrdeb=jdeb*24
-hrfin=jfin*24  
+hrfin=jfin*24
 print,'debut et fin', jdeb,jfin
 
@@ -223 +231 @@
                                  , lon_min, lon_max, lat_min, lat_max $
                                  , testfilename)
-            desc=result.data[0]
-            hour=result.data[1]
-            fov=result.data[2]
-            lon=result.data[3]
-            lat=result.data[4]
-            mask=result.data[5]
-            tb1=result.data[6]
-; decodage du nb de points (nn)
+            desc=result.data.desc
+            hour=result.data.hour
+            fov=result.data.fov
+            lon=result.data.lon
+            lat=result.data.lat
+            mask=result.data.landseamask
+            tb1=result.data.tb
+            ; decodage du nb de points (nn)
             nn=n_elements(tb1)
-; egalement aussi des bornes en lon / lat (a retirer des parametres d'entree)
+; egalement aussi des bornes en lon / lat (a retirer des parametres d'entree)++
 
 ; boucle sur les points du fichier
@@ -238 +246 @@
             for i=0,nn do begin
 
-; recup de tb1,jour,hour,lon,lat
-               if (tb1 ge tbmin and tb1 le tbmax) then begin
+; recup de tb1,hour,lon,lat
+               if (tb1[i] ge tbmin and tb1[i] le tbmax) then begin
 ;if  ((tb1 lt tbmin) or (tb1 gt tbmax)) then goto,suite0 ;; on ne s interesse pas aux donnees aberrantes
 ; a remplacer par if tb1 ne Nan then begin
 ;if (cont lt 10) then print,lon,lat,fov,tb1,hour
 ; test sur le debut et la fin des donnes utiles tenant compte du rayon temporel
-                  jour=julday(mm,dd,yyyy,0)*24.+hour ; en heures decimales
+                  jour=julday(mm,dd,yyyy,0)*24.+hour[i] ; en heures decimales
                                 ;if (cont lt 10) then print,mm,dd,yyyy,'jour julien lu (en heure)',jour
-                  if ((jour ge hrdeb-rtemphor) and (jour le hrfin+rtemphor) and (hour le hfin) and (hour ge hdeb) ) then begin
+                  if ((jour ge hrdeb-rtemphor) and (jour le hrfin+rtemphor) and (hour[i] le hfin) and (hour[i] ge hdeb) ) then begin
                                 ; print,'on cherche les donnees dans la fenetre temporelle'
                                 ; utilisation seulement des orbites du
@@ -261 +269 @@
                                 ; calcul de la fonction de Cressman
                      if (fov gt fov_int2 and fov lt fov_int3) then begin
-                        yy=lat*coef 
+                        yy=lat*coef
                         xx=lon*coef
                         coyy=cos(yy)
@@ -271 +279 @@
                         for nj=jdp,jfn do begin
                            dij=(jour-t[nj])/rtemphor ; normalisee par le rayon de Cressmann
-                           disj=dij*dij 
+                           disj=dij*dij
                            for l=ldp,lfn do begin
                               coco=coyy*cos(las[l])
@@ -281 +289 @@
                                  ac=acos(vv)
                                  dist=rterre*rterre*ac*ac/r2
-                                 if (dist lt 1. and disj lt 1.) then begin 
+                                 if (dist lt 1. and disj lt 1.) then begin
                                     rnorm=2.
-                                    poids=(rnorm-dist-disj)/(rnorm+dist+disj)       
+                                    poids=(rnorm-dist-disj)/(rnorm+dist+disj)
                                     stb1[k,l,nj]=stb1[k,l,nj]+tb1*poids
                                     z1[k,l,nj]=z1[k,l,nj]+poids

trunk/src/idl_amsu_netcdf.pro

@@ -207 +207 @@
 ncdf_attput, id, varid[2], 'units', 'days since 2000-01-01 00:00:00'
 ncdf_attput, id, varid[2], 'long_name', 'Time axis'
+ncdf_attput, id, varid[3], 'time_origin', '2000-01-01 00:00:00'
 ncdf_attput, id, varid[3], 'short_name', 'amsu_tb'
 ncdf_attput, id, varid[3], 'missing_value', '-999'