Changeset 456 for trunk

Timestamp:

12/11/11 22:01:00 (13 years ago)

Author:

lelod

Message:

preparation vf de extract

File:

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

trunk/src/extract_amsua.pro

@@ -77 +77 @@
 ;   ${PROJECT_ID}/AMSU/2006/08/a5_20060801_060w30s_050e45n.dat must have been created.
 ;
+; TODO
+;
+; lelod 20111211
+; attention, il faut transformer les donnees 2D en 1D
+; reprendre methode "extract_amsuab", qui utilisait un "where"?
+; ou passer directement dans l'ecriture???
+;
 ; :History:
 ;
 ; $Id$
 ;
+; lelod 20111211
+; mis en commentaires la partie interpolation
+; prevu decodage de numch pour connaitre le nom instrument et le
+;numero du canal (pas fini, au debut)
 ; $URL: svn+ssh://lelod@forge.ipsl.jussieu.fr/ipsl/forge/projets/varamma/svn/trunk/src/extract_amsuab.pro $
 ;
-; evolution (lelod, dec 2011) pour extraire amsua seul avec resolution
-; amelioree
 ;
 ; - fplod 20111209T125312Z aedon.locean-ipsl.upmc.fr (Darwin)
@@ -184 +193 @@
 ;
 @cm_project
-
-;===================================
-; fichier de sortie comportant les observations amsu
-;
-;====================================
+nomcanal=strmid(numch,0,1)
+;nocanal=strmid(numch,1,1)
+;??????? passage en numerique
+nocanal=5
 
 
@@ -263 +271 @@
 ; latitude, et les Tbs des differents canaux (15 pour AMSUA)
 
-    ;amch1=REFORM(ama_scan.btemps[0,*,*]/100.)
-   ; amch2=REFORM(ama_scan.btemps[1,*,*]/100.)
-    ;amch3=REFORM(ama_scan.btemps[2,*,*]/100.)
-   ; amch4=REFORM(ama_scan.btemps[3,*,*]/100.)
-   ; amch5=REFORM(ama_scan.btemps[4,*,*]/100.)
-   ; amch6=REFORM(ama_scan.btemps[5,*,*]/100.)
-   ; amch7=REFORM(ama_scan.btemps[6,*,*]/100.)
-   ; amch8=REFORM(ama_scan.btemps[7,*,*]/100.)
-   ; amch9=REFORM(ama_scan.btemps[8,*,*]/100.)
-   ; amch10=REFORM(ama_scan.btemps[9,*,*]/100.)
-   ; amch11=REFORM(ama_scan.btemps[10,*,*]/100.)
-   ; amch12=REFORM(ama_scan.btemps[11,*,*]/100.)
-   ; amch13=REFORM(ama_scan.btemps[12,*,*]/100.)
-   ; amch14=REFORM(ama_scan.btemps[13,*,*]/100.)
-   ; amch15=REFORM(ama_scan.btemps[14,*,*]/100.)
 ; extraction du canal traite
-    nocanal=5
+
     amch=REFORM(ama_scan.btemps[nocanal,*,*]/100.)
     amzen=REFORM(ama_scan.angles[0,*,*]/100.)
@@ -296 +289 @@
         desc=0
     endelse
+; attention tableaux a deux dimensions
+
+
 
 ; appel a interpolswatw pour ajuster les pixels amsua sur une grille
@@ -301 +297 @@
 ; et selection de la zone conservee
 ;tek_color
-resol=1
+;resol=1
 ;nbgrid=0
-chint=fltarr(100)
-latint=fltarr(100)
-lonint=fltarr(100)
-fovint=fltarr(100)
-timeint=fltarr(100)
-time=fltarr(100)
+;chint=fltarr(100)
+;latint=fltarr(100)
+;lonint=fltarr(100)
+;fovint=fltarr(100)
+;timeint=fltarr(100)
+;time=fltarr(100)
 ;n_scan-1L
-for i=0L,10L do begin
-   tb=amch[*,i]
-   lon=amalon[*,i]
-   lat=amalat[*,i]
+;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
+;   interpolswath,tb,lat,lon,resol,nbgrid,tbgrid,latgrid,longrid
+;     if i eq 0L then begin
        ; plot,lon,tb,psym=1
-       ; oplot,longrid,tbgrid,psym=4
-        fovgrid=indgen(nbgrid)+1
-     endif
+       ; 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)
-print, 'npt', npt
-  if npt ne 0 then begin
-     time[zone]=tt[i]
-     chint=[chint,tbgrid[zone]]
-     latint=[latint,latgrid[zone]]
-     lonint=[lonint,longrid[zone]]
-     fovint=[fovint,fovgrid[zone]]
-     timeint=[timeint,time[zone]]
-  endif else begin
-     print, 'www : aucun point dans la zone'
-  endelse
-
- endfor
-
- if (npt ne 0) then begin
- ; partie a remplacer par module ecriture standardise FP
-;
+;  zone=where((longrid ge lon_min) and (longrid le lon_max) $
+;             and (latgrid ge lat_min) and (latgrid le lat_max),npt)
+;print, 'npt', npt
+;  if npt ne 0 then begin
+;     time[zone]=tt[i]
+;     chint=[chint,tbgrid[zone]]
+;     latint=[latint,latgrid[zone]]
+;     lonint=[lonint,longrid[zone]]
+;     fovint=[fovint,fovgrid[zone]]
+;     timeint=[timeint,time[zone]]
+;  endif else begin
+;     print, 'www : aucun point dans la zone'
+;  endelse
+
+; endfor
+
+
+
       ; replace missing values -9999.99 by NaN
   index_missing = where(chint EQ -9999.99, nb_missing)
   if (nb_missing NE 0) THEN BEGIN
-     chint[index_missing]=!values.f_nan
+     tb[index_missing]=!values.f_nan
   endif
 ; param landseamask mis fictivement pour unification du format
         landseamask=!values.f_nan
-nn=max(size(chint))
+
 print,'nb points dans la zone',nn
 
@@ -376 +371 @@
                           , landseamask: 9$
                           , tb : chint[i] $
-                       }]
-         endelse
-     endfor
-  endif else begin
-     print, 'www : 2 aucun point dans la zone'
-  endelse
+
 
 ; fin boucle sur les fichiers lus