Changeset 449 for trunk

Timestamp:

12/07/11 23:25:20 (13 years ago)

Author:

lelod

Message:

tests extractamsua

Location:

trunk/src

Files:

• extract_amsua.pro (modified) (2 diffs)
• interpolswath.pro (modified) (4 diffs)

trunk/src/extract_amsua.pro

@@ -286 +286 @@
 ; reguliere dans la fauchee
 ; et selection de la zone conservee
+;tek_color
 resol=1
+;nbgrid=0
 chint=fltarr(100)
 latint=fltarr(100)
@@ -293 +295 @@
 timeint=fltarr(100)
 time=fltarr(100)
-for i=0L,n_scan-1L do begin
+;n_scan-1L
+for i=0L,10L do begin
    tb=amch[*,i]
    lon=amalon[*,i]
    lat=amalat[*,i]
+  ; nbgrid1=nbgrid
    interpolswath,tb,lat,lon,resol,nbgrid,tbgrid,latgrid,longrid
-   if i eq 0L then begin
-      fovgrid=indgen(nbgrid)+1
-      plot,lon,tb,psym=1
-      oplot,longrid,tbgrid,psym=4
-  endif
+     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
   zone=where((longrid ge lon_min) and (longrid le lon_max) $
              and (latgrid ge lat_min) and (latgrid le lat_max),npt)

trunk/src/interpolswath.pro

@@ -41 +41 @@
 ;
 ;-
-pro interpolswath, tb, latlu,lonlu,resol, nbgrid,tbint,latgrid,longrid
+pro interpolswath, tb, latlu, lonlu, resol, nbgrid, tbint, latgrid, longrid
+
 pixelsize,pixatot,pixbtot,alongatot,alongbtot
 na=max(size(pixatot))
@@ -54 +55 @@
 lon=lonlu*deg
 mnadir=nnadir-1
-
+ifov=indgen(na)
 for i=nnadir+1,na-1 do begin
    dist=rterre*acos(sin(lat(i-1))*sin(lat(i)) + cos(lat(i-1))*cos(lat(i))*cos(lon(i)-lon(i-1)))
@@ -64 +65 @@
    fova[i]=fova[i+1]-dist
 endfor
-
-print,fova
+;print,'fova'
+;print,fova
+;print,'pixatot'
+;print,pixatot
+window,1
+ploterr,ifov,fova,pixatot
 
 distmax=-min(fova)+max(fova)+max(pixa)
@@ -84 +89 @@
 ;for i=ndemigrid-1,0,-1 do grid1[i]=-(ndemigrid-i-1)*delta-fova[nnadir]
 grid=[grid1,grid2]
-print,grid
+;print,grid
 ;print,tb
 tbint=fltarr(nbgrid)
 for i=0,na-1 do begin
    ind=where(abs(grid-fova[i]) le pixatot[i]/2,nii)
-   print,i,nii
+  ; print,i,nii
    if nii ne 0 then begin
       tbint[ind]=tb[i]
    endif
+   if nii eq 0 then begin
+   endif
+
+      
 endfor
 latgrid=interpol(latlu,fova,grid)