Changeset 470 for trunk

Timestamp:

12/15/11 19:09:21 (12 years ago)

Author:

lelod

Message:

divers

File:

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

trunk/src/extract_amsua.pro

@@ -277 +277 @@
 nlist=0
 ilist=0
-loadct,39
-numwindow=1
+;loadct,39
+tek_color
+numwindow=2
+window,numwindow
+plot,indgen(30),cor_l+200,psym=2,color=0,yrange=[150,250]
+oplot,indgen(30),cor_s+200,psym=2,color=100
+
 ;while ilist le nlist do begin
 while (not eof(lun_a)) do begin
@@ -360 +365 @@
       endfor
    endfor
-;moych=fltarr(nbpix)
-;for ifo=0,nbpix-1 do moych[ifo]= mean(ch_nadir[ifo,*])
+
+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
-;plot,amafov,moych,psym=1,color=0
-;oplot,amafov,cor_l+200,psym=2,color=200
+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
@@ -380 +393 @@
 ;endfor
 ;plots,lolo,lala,psym=5,symsize=1,color=colscale
-;numwindow=numwindow+1
+numwindow=numwindow+1
 
 ; appel a interpolswath pour ajuster les pixels amsua sur une grille
@@ -429 +442 @@
             print, 'www : aucun point dans la zone'
          endelse
-      endif else begin
-         print,'fauchee incomplete',nbon
-      endelse
+      endif ;else begin
+        ; print,'fauchee incomplete',nbon
+     ; endelse
 
    endfor
@@ -491 +504 @@
 
 endif else begin
-    print, 'www : no data to write'
+  ;  print, 'www : no data to write'
 endelse
 loadct,39