Changeset 481 for trunk

Timestamp:

01/04/12 14:37:35 (12 years ago)

Author:

lelod

Message:

nettoyage extract et interpolswath

Location:

trunk/src

Files:

• extract_amsua.pro (modified) (4 diffs)
• interpolswath.pro (modified) (1 diff)

trunk/src/extract_amsua.pro

@@ -278 +278 @@
 ilist=0
 ;loadct,39
-tek_color
+;tek_color
 numwindow=2
-window,numwindow
-plot,indgen(30)+1,cor_l+200,psym=2,color=0,yrange=[150,250]
-oplot,indgen(30)+1,cor_s+200,psym=2,color=100
+;window,numwindow
+;plot,indgen(30)+1,cor_l+200,psym=2,color=0,yrange=[150,250]
+;oplot,indgen(30)+1,cor_s+200,psym=2,color=100
 
 ;while ilist le nlist do begin
@@ -378 +378 @@
 ;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])
-
-;window,numwindow
-;plot,amalon,amalat,xrange=[lon_min,lon_max],yrange=[lat_min,lat_max],xstyle=1,ystyle=1,/nodata
-;jnd=where(ch_nadir gt tbmin and ch_nadir lt tbmax and amalon ge lon_min and amalon le lon_max $
-;         and amalat ge lat_min and amalat le lat_max)
-;tb=ch_nadir[jnd]
-;lolo=amalon[jnd]
-;lala=amalat[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
-;plots,lolo,lala,psym=5,symsize=1,color=colscale
+;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])
+
+
 numwindow=numwindow+1
 
@@ -408 +397 @@
       lat=amalat[*,i]
       mask=landseamask[*,i]
-      ind=where(tb gt 0 and lon gt -60 and lon lt 100,nbon)
+      ind=where(tb gt tbmin and lon gt lon_min-15 and lon lt lon_max+15,nbon) 
+; test pour eviter pb d'interpolation de longitude (meridien 180)
       if nbon eq nbpix then begin
          cont=cont+1L
@@ -440 +430 @@
             endelse
          endif else begin
-            print, 'www : aucun point dans la zone'
+            ;print, 'www : aucun point dans la zone'
          endelse
-      endif ;else begin
-        ; print,'fauchee incomplete',nbon
-     ; endelse
+      endif
 
    endfor

trunk/src/interpolswath.pro

@@ -117 +117 @@
   ;     
 ;endfor
+
+
 latgrid=interpol(latlu,fova,grid)
 longrid=interpol(lonlu,fova,grid)
 
-ind=where(longrid gt 35 and longrid lt 45 and latgrid gt 10 and latgrid lt 25,nbind)
-if nbind ne 0 then begin
-  ; print,'longitude intiniale',lonlu
-  ; print,'long. interpolee',longrid
-endif
 ;print,tbint
 end