Changeset 495 for trunk

Timestamp:

03/12/12 17:59:38 (12 years ago)

Author:

lelod

Message:

correction bugs dimensions n_scan nzon

File:

• trunk/src/extract_amsua.pro (modified) (11 diffs)

trunk/src/extract_amsua.pro

@@ -81 +81 @@
 ;   IDL> dd=13
 ;   IDL> resol=1
-;   IDL> lon_min=-60.
-;   IDL> lon_max=50.
-;   IDL> lat_min=-30.
-;   IDL> lat_max=45.
+;   IDL> lon_min=-20.
+;   IDL> lon_max=40.
+;   IDL> lat_min=-5.
+;   IDL> lat_max=40.
 ;   IDL> extract_amsua, numch, yyyy, mm, dd, resol, lon_min, lon_max, lat_min, lat_max
 ;
@@ -251 +251 @@
 tbmin=100
 tbmax=350
-PRINT, 'iii, debut programme extract_amsu',SYSTIME() 
+PRINT, 'www, debut programme extract_amsu',SYSTIME() 
 ; lecture fichier land - sea (S. Masson)
 file=project_id_env+'/MASK/ETOPO1_Ice_g_gmt4.nc'
@@ -286 +286 @@
 a = STRARR(nb_file)
    ; filea = STRARR(nb_filea)
-print, 'iii ouverture pour lecture de ', list_file
-print, 'iii traitement du jour ',yyyy,mm,dd
+;print, 'www ouverture pour lecture de ', list_file
+print, 'www traitement du jour ',yyyy,mm,dd
 
 openr, lun_a, list_file, /get_lun
@@ -340 +340 @@
    amzeni=REFORM(ama_scan.angles[0,*,*]/100.)
    ttt=REFORM(ama_scan.scnlintime/3600000.)
-
    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 dans le domaine geographique ",nzon, SYSTIME()
+   print,"www nb de points du fichier dans le domaine geographique +15deg ",nzon, SYSTIME()
    if nzon ne 0 then begin
       amalat=amalati[*,jnd]
@@ -366 +365 @@
  
 ; correction nadir des donnees
-         ch_nadir=fltarr(nbpix,n_scan)
-      ; ch_nadir=amch[ifo,isc]
-         landseamask=intarr(nbpix,n_scan)+2 ; valeur hors zone selectionnee
+         ch_nadir=fltarr(nbpix,nzon)
+         landseamask=intarr(nbpix,nzon)+2 ; valeur hors zone selectionnee
         ; PRINT,'boucle sur les points du fichier, correction nadir', SYSTIME() 
 
-         for isc=0L,n_scan-1L do begin
+         for isc=0L,nzon-1L do begin
 ;recherche de la zone xxe,yye englobant la fauchee
             for ifo=0,nbpix-1 do begin
@@ -395 +393 @@
         ; PRINT, 'fin boucle',SYSTIME() 
          moych=fltarr(nbpix)
-      
-;for ifo=0,nbpix-1 do begin
-;   jnd=where(ch_nadir[ifo,*] gt tbmin and ch_nadir[ifo,*] lt tbmax and amalon[ifo,*] ge lon_min and amalon[ifo,*] le lon_max $
-;         and amalat[ifo,*] ge lat_min and amalat[ifo,*] le lat_max,npts)
-;  if npts ne 0 then begin
-;     moych[ifo]= mean(ch_nadir[ifo,jnd])
-;  endif
-;endfor
-;print,'correction',moych
-;print,'cor_s',cor_s
-;print,'cor_l',cor_l
-;oplot,amafov,moych,psym=1,color=numwindow
-;print,'ecarts max a la valeur moyenne au nadir',min(moych)-mean(moych[12:17]),max(moych)-mean(moych[12:17])
-
 
 ; appel a interpolswath pour ajuster les pixels amsua sur une grille
@@ -419 +403 @@
          cont=0L
          chint=fltarr(1)
-         for i=0L,n_scan-1L do begin
+         for i=0L,nzon-1L do begin
             tb=ch_nadir[*,i]
             lon=amalon[*,i]
@@ -433 +417 @@
                   fovgrid=indgen(nbgrid)+1
                endif
-         ;longit=longrid
-         ;latit=latgrid
-         ;oplot,longit,latit,psym=3,color=0
 
 ; selection des taches au sol situees dans la zone d'interet
@@ -453 +434 @@
                      lonint=[lonint,longrid[zone]]
                      fovint=[fovint,fovgrid[zone]]
-                     maskint=[maskint,maskgrid[zone]]
+                     maskint=[maskint,maskgrid[zone]]               
                      timeint=[timeint,replicate(tt[i],npt)]
                   endelse
                endif else begin
-            print, 'www : aucun point dans la zone'
+            ;print, 'www : aucun point dans la zone'
                endelse
             endif
@@ -465 +446 @@
          nn=n_elements(chint)
 
-        ; print,'nb points dans la zone',nn
+         print,'www : nb points dans la zone en fin de traitement de l'orbite',nn
 
          for i=0L, nn-1L do begin
@@ -491 +472 @@
 ; fin boucle sur les fichiers lus
    endif
-   PRINT,'passage au fichier suivant', SYSTIME() 
+   PRINT,'www : passage au fichier suivant ', SYSTIME() 
 endwhile
 labfile: