Changeset 483 for trunk

Timestamp:

01/06/12 16:47:54 (12 years ago)

Author:

pinsard

Message:

timing dans traitement amsutraite_amsuab.sh

Location:

trunk/src

Files:

• extract_amsua.pro (modified) (7 diffs)
• interpol_correc.pro (modified) (4 diffs)

trunk/src/extract_amsua.pro

@@ -14 +14 @@
 ; generaliser pour amsub aussi)
 ;
-; decode les noms des fichiers dans la liste de fichiers, puis
-; appelle le prgm de lecture L1C
+; decode les noms des fichiers dans la date choisie, puis
+; appelle le prgm de lecture
 ;
 ; appelle :ref:`interpol_correc`, qui fournit les fonctions de correction au
@@ -22 +22 @@
 ;
 ; applique une interpolation tenant compte de la dimension de la tache
-; au sol selon la position dans la fauchee : appel :ref:`interpolswath`
+; au sol selon la position dans la fauchee
 ;
 ; ecrit les donnees pour la zone (lat/long) choisie
-; la structure de sortie contient un en-tete (date, zone, resol)
-; puis les donnees du jour et du canal traite
-;
+;
+; ++ le fichier de sortie contient dans l'ordre mm, jj (jour de
+; l'annee), le temps (identique pour les points d'une meme fauchee), le
+; no satellite, le no de fov (1 a 40 ou 1 a 90), long, lat, tb
+;
+; Les valeurs manquantes sont codées par NaN.
 ;
 ; .. only:: man
@@ -68 +71 @@
 ;
 ;   $ mkdir -p ${PROJECT_ID}/AMSU/2006/08
-;   $ echo "${PROJECT_ID}/AMSU/AMSUAN15/L1C/2006/2006_08_01/NSS.AMAX.NK.D06213.S0110.E0255.B4271112.GC.L1C" > ${PROJECT_ID}/list_file
-;   $ echo "${PROJECT_ID}/AMSU/AMSUAN15/L1C/2006/2006_08_01/NSS.AMAX.NK.D06213.S0250.E0437.B4271213.GC.L1C" >> ${PROJECT_ID}/list_file
+;   $ echo "${PROJECT_ID}/AMSU/AMSUAN15/L1C/2006/2006_08_01/NSS.AMAX.NK.D06213.S0110.E0255.B4271112.GC.L1C" > ${PROJECT_ID}/list_filea
+;   $ echo "${PROJECT_ID}/AMSU/AMSUAN15/L1C/2006/2006_08_01/NSS.AMAX.NK.D06213.S0250.E0437.B4271213.GC.L1C" >> ${PROJECT_ID}/list_filea
 ;
 ; Using AMSU-A channel 5::
@@ -189 +192 @@
 ; called by :ref:`traite_amsuab.sh`
 ;
-; use :ref:`read_amsua1c`, :ref:`read_amsub1c.pro`,
+; use :ref:`read_amsua1c.pro`, :ref:`read_amsub1c.pro`,
 ; :func:`geolocation_to_string_idl`, :func:`mem_to_file_amsu_t2`
-; use :ref:`interpol_correc` 
-; use :ref:`interpolswath`
-;
-; next step :  au choix: amsu2ncdf ou cresamsu pour transformation en netcdf
+;
+; next step :  ++
 ;
 ; TODO
@@ -202 +203 @@
 ; interpolswath ne gere pas cette situation
 ;
+; comprendre les points bizarres autour de -40,30 visibles avec l'exemple
+; de 20060801
+;
+; lever le doute sur le contenu du fichier écrit par le printf (je (fplod)
+; crains des zesors louches!!).louche
+;
 ; vérifier que pas NaN dans les lignes écrites
 ;
 ; faire un test avec un zone géographique dans le pacifique pour voir comment
-; le système réagit
+; le système réagi (devrait dire : rien à faire !)
 ;
 ; changer le nom des fichiers de listes pour éviter des conflits lors de runs
@@ -214 +221 @@
 ; eviter la creation d'un fichier vide (si pas de points dans la zone)
 ;
-; tester la cohérence entre les fichiers présents dans list_file
+; tester la cohérence entre les fichiers présents dans list_filea
 ; en terme de date avec yyyy, mm et dd donnés en paramètre
 ;
@@ -481 +488 @@
   ;  print, 'www : no data to write'
 endelse
-loadct,39
-window,0
-plot,lonint,latint,xrange=[lon_min,lon_max],yrange=[lat_min,lat_max],xstyle=1,ystyle=1,/nodata
-ind=where(maskint eq 1)
-;oplot,lonint[ind],latint[ind],psym=1
-jnd=where(chint gt tbmin and chint lt tbmax)
-tb=chint[jnd]
-lolo=lonint[jnd]
-lala=latint[jnd]
-colscale=intarr(n_elements(tb))
-for j=0,n_elements(tb)-1 do begin
-   colscale[j]=(tb[j]-min(tb))*255/(max(tb)-min(tb))
-endfor
-print,min(colscale),max(colscale)
-   plots,lolo,lala,psym=5,symsize=1,color=colscale
+; loadct,39
+; window,0
+; plot,lonint,latint,xrange=[lon_min,lon_max],yrange=[lat_min,lat_max],xstyle=1,ystyle=1,/nodata
+; ind=where(maskint eq 1)
+; ;oplot,lonint[ind],latint[ind],psym=1
+; jnd=where(chint gt tbmin and chint lt tbmax)
+; tb=chint[jnd]
+; lolo=lonint[jnd]
+; lala=latint[jnd]
+; colscale=intarr(n_elements(tb))
+; for j=0,n_elements(tb)-1 do begin
+;    colscale[j]=(tb[j]-min(tb))*255/(max(tb)-min(tb))
+; endfor
+; print,min(colscale),max(colscale)
+;    plots,lolo,lala,psym=5,symsize=1,color=colscale
 end

trunk/src/interpol_correc.pro

@@ -10 +10 @@
 ; ===========
 ;
-; interpolation des tableaux de correction amsub (calcules sur 
+; interpolation des tableaux de correction amsub (calcules sur
 ; echantillonnage amsua) sur les 90 points de la fauchee amsub
 ;
-; read :file:`${PROJECT_ID}/CORR_[LAND|SEA]_AMSU[AB]_JUIL2006.DAT`
+; read :file:`${PROJECT}/src/dataref_amsu/CORR_[LAND|SEA]_AMSU[AB]_JUIL2006.DAT`
 ; containing correction for land and sea. ++ ref. algo fatima karbou
 ;
@@ -77 +77 @@
 ; ==========
 ;
-; lelod 2011/12/09 modif organisation avec separation des traitements
-; un seul canal traite
 ; $Id$
 ;
 ; $URL$
 ;
+; - fplod 20120106
+;
+;   * reprise des modifs lelod : changement de localisation des fichiers de correction
+;
+; - lelod 20111209
+;
+;   * modif organisation avec separation des traitements un seul canal traite
+;
 ; - lelod 20111019
-; 
+;
 ;   * fonction pixelsize fait le calcul des tailles des taches au sol a et
 ;     b, et ressort le tableau correspondant sur la fauchee
@@ -160 +166 @@
 
 
-path=project_env + 'src/dataref_amsu/'
+path = project_env + '/src/dataref_amsu/'
 ; lecture des fichiers d'ajustement au nadir pour AMSUA et B
 ; cas amsua
 if nomcanal eq 'a' then begin
-   cor_landa=path + '/CORR_LAND_AMSUA_JUIL2006.DAT'
-   cor_seaa=path + '/CORR_SEA_AMSUA_JUIL2006.DAT'
+   cor_landa = path + 'CORR_LAND_AMSUA_JUIL2006.DAT'
+   cor_seaa = path +'CORR_SEA_AMSUA_JUIL2006.DAT'
   ;test de l'autre fichier de correction: quel est le bon???
   ; cor_landa=project_id_env+'AMSU/CORR_LAND_AMSUA.DAT'
@@ -213 +219 @@
 ;cas amsub
 if nomcanal eq 'b' then begin
-   cor_landb=path+'CORR_LAND_AMSUB_JUIL2006.DAT'
-   cor_seab=path+'CORR_SEA_AMSUB_JUIL2006.DAT'
-   ;cor_landb=path+'AMSU/CORR_LAND_AMSUB.DAT'
-   ;cor_seab=path+'AMSU/CORR_SEA_AMSUB.DAT'
+   cor_landb = path + 'CORR_LAND_AMSUB_JUIL2006.DAT'
+   cor_seab = path + 'CORR_SEA_AMSUB_JUIL2006.DAT'
+   ;cor_landb=path+'CORR_LAND_AMSUB.DAT'
+   ;cor_seab=path+'CORR_SEA_AMSUB.DAT'
    cor_l=fltarr(nfovb)
    cor_s=fltarr(nfovb)