Changeset 583 for trunk/src/extract_amsu.pro

Timestamp:

07/09/12 16:22:12 (12 years ago)

Author:

lelod

Message:

adaptation chaine aux autres canaux

File:

• trunk/src/extract_amsu.pro (modified) (6 diffs)

trunk/src/extract_amsu.pro

@@ -305 +305 @@
 ; appel au ssprgm qui interpole les fonctions de correction (inutile
 ; de l'appeler a chaque fichier....)
+; ou appeler pixelsize et reduire interpolcorrec a la partie
+; correction nadir
 interpol_correc,numch,nbpixel,cor_l,cor_s,swath,track
 ; ouverture des fichiers liste (annee, mois, jour, tous satellites) pour
@@ -314 +316 @@
 ;PRINT,'demarrage boucle sur lichiers', SYSTIME()
 FOR ifile = 0, nb_file - 1 do begin
-   filea = files_list[ifile]
    if nomcanal eq 'a' then begin
+      filea = files_list[ifile]
       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
+      ;openr,lu1,filea,Error=erra,/get_lun
       read_amsua1c,filea, flag1
-      close,lu1
-      free_lun,lu1
+     ;close,lu1
+     ;free_lun,lu1
       if (flag1 eq 0) then goto, labfile
       na=SIZE(reform(ama_scan.btemps[0,*,*]))
@@ -337 +339 @@
       amalong=REFORM(ama_scan.latlon[1,*,*]/1.E4)
       amalati=REFORM(ama_scan.latlon[0,*,*]/1.E4)
-      amcha=REFORM(ama_scan.btemps[nocanal,*,*]/100.)
+      amcha=REFORM(ama_scan.btemps[nocanal-1,*,*]/100.)
       amzeni=REFORM(ama_scan.angles[0,*,*]/100.)
       ttt=REFORM(ama_scan.scnlintime/3600000.)
@@ -345 +347 @@
       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
+     ;openr,lu1,fileb,Error=errb,/get_lun
       read_amsub1c,fileb, flag2
-      close,lu1
-      free_lun,lu1
+      ;close,lu1
+     ;free_lun,lu1
       if (flag2 eq 0) then goto, labfile
       nb=SIZE(reform(amb_scan.btemps[0,*,*]))
@@ -359 +362 @@
       amalong=REFORM(amb_scan.latlon[1,*,*]/1.E4)
       amalati=REFORM(amb_scan.latlon[0,*,*]/1.E4)
-      amcha=REFORM(amb_scan.btemps[nocanal,*,*]/100.)
+      amcha=REFORM(amb_scan.btemps[nocanal-1,*,*]/100.)
       amzeni=REFORM(amb_scan.angles[0,*,*]/100.)
       ttt=REFORM(amb_scan.scnlintime/3600000.)
@@ -365 +368 @@
    endif
 
-   fov=indgen(nbpix)            ; nb pixels dans la fauchee AMSU
+   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)