Changeset 465 for trunk

Timestamp:

12/14/11 19:15:56 (12 years ago)

Author:

lelod

Message:

tests

File:

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

trunk/src/extract_amsua.pro

@@ -85 +85 @@
 ;   IDL> lat_min=-30.
 ;   IDL> lat_max=45.
-;   IDL> extract_amsua, numch, yyyy, mm, dd, resol $
-;                     , lon_min, lon_max, lat_min, lat_max
+;   IDL> extract_amsua, numch, yyyy, mm, dd, resol, lon_min, lon_max, lat_min, lat_max
 ;
 ; :file:`${PROJECT_ID}/AMSU/2006/08/a5_20060801_060w30s_050e45n.dat`
@@ -269 +268 @@
 ilist=0
    ; while (not eof(lun_a)) do begin
+loadct,39
+numwindow=1
 while ilist le nlist do begin
    ilist=ilist+1
@@ -349 +350 @@
       endfor
    endfor
-
+window,numwindow
+plot,amalon,amalat,xrange=[lon_min,lon_max],yrange=[lat_min,lat_max],xstyle=1,ystyle=1,psym=2,color=0
+jnd=where(ch_nadir gt tbmin and ch_nadir lt tbmax)
+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
+print,min(colscale),max(colscale)
+plots,lolo,lala,psym=5,symsize=2,color=colscale
+numwindow=numwindow+1
 
 ; appel a interpolswath pour ajuster les pixels amsua sur une grille
@@ -365 +378 @@
       interpolswath,tb,lat,lon,mask,resol,nbgrid,tbgrid,latgrid,longrid,maskgrid
       if i eq 0L then begin
-       ; plot,lon,tb,psym=1
-       ; oplot,longrid,tbgrid,psym=4;
          fovgrid=indgen(nbgrid)+1
       endif
-  ; if nbgrid ne nbgrid1 then begin
-  ;    oplot,lon,tb,psym=1,color=3
-   ;   oplot,longrid,tbgrid,psym=4,color=5
-   ;   endif
+     ; ind=where(tbgrid gt tbmin and tbgrid lt tbmax,nbon)
+     ; if (nbon ne 0) then begin
+     ;    longit=longrid(ind)
+     ;    latit=latgrid(ind)
+     ;    oplot,longit,latit,psym=4,color=0
+     ; endif else begin
+     ;    print,'pas de point valable dans la fauchee'
+     ; endelse
+
 ; selection des taches au sol situees dans la zone d'interet
       zone=where((longrid ge lon_min) and (longrid le lon_max) $
@@ -455 +471 @@
     print, 'www : no data to write'
 endelse
-
-window,1
-plot,lonint,latint,xrange=[lon_min,lon_max],yrange=[lat_min,lat_max],xstyle=1,ystyle=1,psym=2
+loadct,39
+window,0
+plot,lonint,latint,xrange=[lon_min,lon_max],yrange=[lat_min,lat_max],xstyle=1,ystyle=1,psym=3,color=0
+ind=where(maskint eq 1)
+oplot,lonint[ind],latint[ind],psym=1
+jnd=where(chint gt tbmin and chint lt tbmax)
+tb=chint[jnd]
+lolo=lonint[jnd]
+lala=latint[jnd]
+colscale=intarr(n_elements(tb))
+for j=0,n_elements(tb)-1 do begin
+   colscale[j]=tb[j]*255/(max(tb)-min(tb))
+endfor
+print,min(colscale),max(colscale)
+   plots,lolo,lala,psym=5,symsize=2,color=colscale
 end