Changeset 588

Timestamp:

07/18/12 14:19:02 (12 years ago)

Author:

lelod

Message:

chaine prete pour amsub et mhs???

Location:

trunk/src

Files:

• extract_amsu.pro (modified) (24 diffs)
• interpolswath.pro (modified) (2 diffs)
• search_amsufiles.pro (modified) (3 diffs)
• traite_amsuab.sh (modified) (3 diffs)

trunk/src/extract_amsu.pro

@@ -13 +13 @@
 ; prgm de lecture des fichiers AMSU (A et B)
 ;
-; decode les noms des fichiers donnés en argument dans **files_list**
+; decode les noms des fichiers donnés en argument dans **files_list**
 ; dans la date choisie, puis appelle le prgm de lecture
 ;
@@ -29 +29 @@
 ; no satellite, le no de fov (1 a 40 ou 1 a 90), long, lat, tb
 ;
-; Les valeurs manquantes sont codées par NaN.
+; Les valeurs manquantes sont codées par NaN.
 ;
 ; .. only:: man
@@ -96 +96 @@
 ;   lat_max=20.
 ;   files_list = search_amsufiles(numch, yyyy, mm, dd)
+;   extract_amsu, numch, files_list, yyyy, mm, dd, resol, lon_min, lon_max, lat_min, lat_max
+;   result = file_bathy_to_mem(lon_min, lon_max, lat_min, lat_max)
 ;   extract_amsua, numch, files_list, yyyy, mm, dd, resol, lon_min, lon_max, lat_min, lat_max
 ;
@@ -125 +127 @@
 ; interpolswath ne gere pas cette situation
 ;
-; lever le doute sur le contenu du fichier écrit par le printf (je (fplod)
+; lever le doute sur le contenu du fichier écrit par le printf (je (fplod)
 ; crains des zeros louches!!).louche
 ;
@@ -133 +135 @@
 ; eviter la creation d'un fichier vide (si pas de points dans la zone)
 ;
-; tester la cohérence entre les fichiers présents dans files_list
-; en terme de date avec yyyy, mm et dd donnés en paramètre
+; tester la cohérence entre les fichiers présents dans files_list
+; en terme de date avec yyyy, mm et dd donnés en paramÚtre
 ;
 ; improve log
@@ -154 +156 @@
 ; - lelod 20120420 
 ;
-;   * extension de extract_amsua à amsub
+;   * extension de extract_amsua à amsub
 ;
 ; - lelod 20120418
@@ -272 +274 @@
 ;-
 PRO extract_amsu, numch, files_list, yyyy, mm, dd, resol $
-                 , lon_min, lon_max, lat_min, lat_max
+                 , lonmin, lonmax, latmin, latmax
 ;
 compile_opt idl2, strictarrsubs
@@ -278 +280 @@
 @cm_project
 @common
-
+s = size(SCOPE_TRACEBACK(/STRUCTURE),/DIMENSION)
+routine = (SCOPE_TRACEBACK(/STRUCTURE))[s-1].Routine
+key_performance=1
+
+if key_performance EQ 1 THEN BEGIN
+    time1 = SYSTIME(1)
+ENDIF
+;
 ; test if some files to read
 CASE size(files_list,/DIMENSION) OF
@@ -312 +321 @@
 ; boucle sur les elements de la liste
 desc=0
-print, 'iii : traitement du jour ',yyyy,mm,dd
+msg = report(['iii : ' + routine + ' : traitement du jour' $
+      + string(yyyy, format='(1x,I4.4)') $
+      + string(mm, format='(I2.2)') $
+      + string(dd, format='(I2.2)')])
 a = STRARR(nb_file)
-;PRINT,'demarrage boucle sur lichiers', SYSTIME()
+
+IF key_performance EQ 1 THEN BEGIN
+    msg = report(['ppp : ' + routine + ' :  duree intermediaire demarrage boucle sur fichier ' $
+                   + string(SYSTIME(1)-time1,format='(F12.6)')])
+ENDIF
+
 FOR ifile = 0, nb_file - 1 do begin
    if nomcanal eq 'a' then begin
@@ -320 +337 @@
       COMMON amsua_header,ama_head
       COMMON amsua_data  ,ama_scan
-      print, 'ouverture et lecture du fichier ', filea, SYSTIME()
-      ;openr,lu1,filea,Error=erra,/get_lun
+   IF key_performance EQ 1 THEN BEGIN
+       msg = report(['ppp : ' + routine + ' : duree intermediaire ouverture et lecture du fichier' $
+                   + string(SYSTIME(1)-time1,format='(F12.6)')])
+   ENDIF
       read_amsua1c,filea, flag1
-     ;close,lu1
-     ;free_lun,lu1
       if (flag1 eq 0) then goto, labfile
       na=SIZE(reform(ama_scan.btemps[0,*,*]))
@@ -333 +350 @@
       nbpix= na_fov
       n_fov=na_fov
-      nosat=ama_head.h_satid                  
+      nosat=ama_head.h_satid    ; numero du satellite NOAA (15 a 19) ou METOP 2
+      print,'NUMERO DU SATELLITE ',nosat
+      print,'INSTRUMENT',ama_head.h_instrument
 ; on lit la structure et on extrait les infos: temps, longitude,
 ; latitude, et les Tbs des differents canaux (15 pour AMSUA)
-; extraction du canal traite dans un domaine enveloppant la zone choisie
       amalong=REFORM(ama_scan.latlon[1,*,*]/1.E4)
       amalati=REFORM(ama_scan.latlon[0,*,*]/1.E4)
@@ -345 +363 @@
    endif
    if nomcanal eq 'b' then begin
+      fileb = files_list[ifile]
       COMMON amsub_header,amb_head
       COMMON amsub_data  ,amb_scan
-      fileb = files_list[ifile]
-      print, 'ouverture et lecture du fichier ', fileb, SYSTIME()
-     ;openr,lu1,fileb,Error=errb,/get_lun
-      read_amsub1c,fileb, flag2
-      ;close,lu1
-     ;free_lun,lu1
-      if (flag2 eq 0) then goto, labfile
+   IF key_performance EQ 1 THEN BEGIN
+       msg = report(['ppp : ' + routine + ' : duree intermediaire ouverture et lecture du fichier' $
+                   + string(SYSTIME(1)-time1,format='(F12.6)')])
+   ENDIF
+    read_amsub1c,fileb, flag2
+     if (flag2 eq 0) then goto, labfile
       nb=SIZE(reform(amb_scan.btemps[0,*,*]))
       nb_scan=nb[2]             ; nb lignes (ou fauchees)
@@ -359 +377 @@
       nb_fov=nb[1]              ; nb pixels dans la fauchee
       nbpix = nb_fov          ; nb pixels dans la fauchee AMSUB
-      nosat=amb_head.h_satid  ; on le lit dans le fichier, donc pas besoin de l'info dans la liste
+      nosat=amb_head.h_satid    ; numero du satellite NOAA (15 a 19) ou METOP 2
+      print,'NUMERO DU SATELLITE ',nosat
+      print,'INSTRUMENT',amb_head.h_instrument
+; on lit la structure et on extrait les infos: temps, longitude,
+; latitude, et les Tbs des differents canaux (5 pour AMSUB)
       amalong=REFORM(amb_scan.latlon[1,*,*]/1.E4)
       amalati=REFORM(amb_scan.latlon[0,*,*]/1.E4)
@@ -368 +390 @@
    endif
 
-   fov=indgen(nbpix)            ; nb pixels dans la fauchee AMSU ne sert a rien???
    ntime_1=fltarr(nbpix,n_scan)
    fovy_1=fltarr(nbpix,n_scan)
@@ -374 +395 @@
       fovy_1[j,*]=j+1           ; on remplit un tableau de nscan lignes avec les valeurs de no de pixel
    endfor
+; extraction du canal traite dans un domaine enveloppant la zone choisie
    midpix=fix(nbpix/2)
-   jnd=where(amalong[midpix,*] gt lon_min-15 and amalong[midpix,*] lt lon_max+15 and amalati[midpix,*] gt lat_min-15 and amalati[midpix,*] lt lat_max+15 ,nzon)
-   print,"iii : nb de points du fichier dans le domaine geographique +15deg ",nzon, SYSTIME()
+   jnd=where(amalong[midpix,*] gt lonmin-15 and amalong[midpix,*] lt lonmax+15 and amalati[midpix,*] gt latmin-15 and amalati[midpix,*] lt latmax+15 ,nzon)
+   msg = report(['iii : ' + routine + ' : nb de points du fichier dans le domaine geographique +15deg ' + string(nzon,format='(I8.8)')])
+   IF key_performance EQ 1 THEN BEGIN
+       msg = report(['ppp : ' + routine + ' : durï¿œe intermediaire avant test nzon' $
+                   + string(SYSTIME(1)-time1,format='(F12.6)')])
+   ENDIF
    if nzon gt 1 then begin
       amalat=amalati[*,jnd]
@@ -389 +415 @@
       amafov= fovy_1
       jour_ama=amaday[0]        ; jour du passage
-
-;
-; correction nadir des donnees
+      desc=desc+1
+;
+; correction nadir des donnees.En principe le 17 est elimine des
+; listes
+         isat=0
+         if nosat eq 15 then begin
+             isat=0
+             print,'no satellite',nosat, isat
+         endif
+         if nosat eq 16 then isat=1
+         if nosat eq 18 then isat=2
+         if nosat eq 19 then isat=3
+         if nosat eq 2 then isat=4
+ 
       ch_nadir=fltarr(nbpix,nzon)
       landseamask=intarr(nbpix,nzon)+2 ; valeur hors zone selectionnee
-; PRINT,'boucle sur les points du fichier, correction nadir', SYSTIME()
+         IF key_performance EQ 1 THEN BEGIN
+              msg = report(['ppp : ' + routine + ' : durï¿œe intermediaire boucle sur les points du fichier, correction nadir ' $
+                   + string(SYSTIME(1)-time1,format='(F12.6)')])
+         ENDIF
+
       for isc=0L,nzon-1L do begin
 ;recherche de la zone xxe,yye englobant la fauchee
@@ -417 +458 @@
                      bate2=bate1[*,yind]
                      if (mean(bate2) le 0.5) then begin
-                        ch_nadir[ifo,isc]=amch[ifo,isc]-cor_s[ifo]
+                        ch_nadir[ifo,isc]=amch[ifo,isc]-cor_s[ifo,isat]
                         landseamask[ifo,isc]=0 ; mer
                      endif else begin
-                        ch_nadir[ifo,isc]=amch[ifo,isc]-cor_l[ifo]
+                        ch_nadir[ifo,isc]=amch[ifo,isc]-cor_l[ifo,isat]
                         landseamask[ifo,isc]=1 ; terre
                      endelse
@@ -429 +470 @@
 
       endfor
-        ; PRINT, 'fin boucle ',SYSTIME()
+        IF key_performance EQ 1 THEN BEGIN
+             msg = report(['ppp : ' + routine + ' : durï¿œe intermediaire fin boucle ' $
+                   + string(SYSTIME(1)-time1,format='(F12.6)')])
+        ENDIF
       moych=fltarr(nbpix)
-
+print,'valeurs extremes tb',min(amch),min(ch_nadir),max(amch),max(ch_nadir)
 ; appel a interpolswath pour ajuster les pixels amsua sur une grille
 ; reguliere dans la fauchee
 ; et selection de la zone conservee
-; attention test sur la longitude code en dur:
-; pbs avec les fauchees qui traversent le meridien 180deg!!!!
-; l'interpolation ajoute des points parasites dans la zone...
       nbgrid=0
       cont=0L
@@ -447 +488 @@
          mask=landseamask[*,i]
          ind=where(tb gt tbmin,nbon)
+         print,'nbon',nbon
 ; test pour eviter pb d'interpolation de longitude (meridien 180)
          if nbon eq nbpix then begin
@@ -455 +497 @@
 ;print,'wwwtest: nb points fov interpole ',n_elements(tbgrid)
 ; selection des taches au sol situees dans la zone d'interet
-            zone=where((longrid ge lon_min) and (longrid le lon_max) $
-                          and (latgrid ge lat_min) and (latgrid le lat_max) and (tbgrid gt tbmin) and (tbgrid lt tbmax), npt)
+            zone=where((longrid ge lonmin) and (longrid le lonmax) $
+                          and (latgrid ge latmin) and (latgrid le latmax) and (tbgrid gt tbmin) and (tbgrid lt tbmax), npt)
 ;print,'wwwtest: nb points zone d interet ',npt
             if npt ne 0 then begin
@@ -479 +521 @@
          endif
       endfor
- ; PRINT,'fin interpolation et selection des points de la zone ', SYSTIME()
+         IF key_performance EQ 1 THEN BEGIN
+             msg = report(['ppp : ' + routine + ' : durï¿œe intermediaire fin interpolation et selection des points de la zone ' $
+                   + string(SYSTIME(1)-time1,format='(F12.6)')])
+         ENDIF
 
       nn=n_elements(chint)
-
-      print,'www : nb points dans la zone en fin de traitement de l''orbite ',nn
+         msg = report(['www : ' + routine + ' : nb points dans la zone en fin de traitement de l''orbite ' + string(nn,format='(I10.10)')])
       if nn gt 1 then begin
          for i=0L, nn-1L do begin
@@ -511 +555 @@
 ; fin boucle sur les fichiers lus
    endif
-   PRINT,'www : passage au fichier suivant ', SYSTIME()
+   IF key_performance EQ 1 THEN BEGIN
+       msg = report(['ppp : ' + routine + ' : durï¿œe intermediaire au passage au fichier suivant ' $
+                   + string(SYSTIME(1)-time1,format='(F12.6)')])
+   ENDIF
 endfor
 labfile:
@@ -523 +570 @@
              ,  resol : resol $
              }
-header2= { lon_min: lon_min $
-                      , lon_max: lon_max $
-                      , lat_min: lat_min $
-                      , lat_max: lat_max $
+header2= { lonmin: lonmin $
+                      , lonmax: lonmax $
+                      , latmin: latmin $
+                      , latmax: latmax $
                      }
 
@@ -535 +582 @@
 
 endif else begin
-    print, 'www : no data to write'
+    msg = report(['www : ' + routine + ' : no data to write'])
     goto, realend
 endelse
 
 nofile:
-    print, 'www : no data to read'
+    msg = report(['www : ' + routine + ' : no data to read'])
     goto, realend
 
 realend:
-    PRINT, 'www : fin programme extract_amsua ',SYSTIME()
+    IF key_performance EQ 1 THEN BEGIN
+       msg = report(['ppp : ' + routine + ' : durï¿œe totale ' $
+                   + string(SYSTIME(1)-time1,format='(F12.6)')])
+       ENDIF
 end

trunk/src/interpolswath.pro

@@ -65 +65 @@
 ;
 pixelsize,pixatot,pixbtot,alongatot,alongbtot
-if type='a' then begin
+if type eq 'a' then begin
    na=max(size(pixatot))
    nnadir=na/2
    pixa=pixatot[nnadir:(na-1)]
-endif else begin
+endif
+if type eq 'b' then begin
    na=max(size(pixbtot))
    nnadir=na/2
    pixa=pixbtot[nnadir:(na-1)]
-endelse
+endif
 fova=fltarr(na)
 fova[nnadir]=pixa[0]/2.
@@ -158 +159 @@
    endif
 endelse
-
+print,mean(lonlu),mean(longrid)
 ;print,tbint
 IF key_performance EQ 1 THEN BEGIN

trunk/src/search_amsufiles.pro

@@ -51 +51 @@
 ; ==========
 ;
+; lelod 20120716 
+; ajout de MHS dans la liste
 ; $Id$
 ;
@@ -173 +175 @@
 nomcanal_upper = STRUPCASE(nomcanal)
 ;
-files_expr = project_id_env $
+files_expr1 = project_id_env $
            + 'AMSU/AMSU' + nomcanal_upper + '*/L1C/' $
            + string(yyyy,format='(i4.4)') + '/' $
@@ -179 +181 @@
            + string(mm,format='(i2.2)') + '_' $
            + string(dd,format='(i2.2)') +'/'+'*.L1C'
+files_expr2 = project_id_env $
+           + 'AMSU/MHS' + '*/L1C/' $
+           + string(yyyy,format='(i4.4)') + '/' $
+           + string(yyyy,format='(i4.4)') + '_' $
+           + string(mm,format='(i2.2)') + '_' $
+           + string(dd,format='(i2.2)') +'/'+'*.L1C'
+files_expr=[files_expr1,files_expr2]
 ;
 ;print, files_expr

trunk/src/traite_amsuab.sh

@@ -315 +315 @@
 # !! when change
 yyyymmddb_min=19991230
-yyyymmdde_max=20130101
+yyyymmdde_max=20150102
 #
 # default
@@ -579 +579 @@
             print, 'iii : nb de fichiers', size(files_list,/DIMENSION)
             resol=1
-            extract_amsua, '${numch}', files_list, ${yyyy}, ${mm}, ${dd}, resol, ${lonmin}, ${lonmax}, ${latmin}, ${latmax}
+            extract_amsu, '${numch}', files_list, ${yyyy}, ${mm}, ${dd}, resol, ${lonmin}, ${lonmax}, ${latmin}, ${latmax}
         END
 ENDCASE
@@ -601 +601 @@
 ; compile users procedures to be scan by Profiler
 .compile interpol_correc
-.compile extract_amsua
+.compile extract_amsu
 .compile read_amsua1c
+.compile read_amsub1c
 .compile define_amsua_header_struct
 .compile define_amsua_struct