Changeset 572 for trunk

Timestamp:

06/04/12 12:35:35 (12 years ago)

Author:

lelod

Message:

maj extract_amsua et cresamsu

Location:

trunk/src

Files:

• cresamsu.pro (modified) (4 diffs)
• extract_amsua.pro (modified) (7 diffs)

trunk/src/cresamsu.pro

@@ -1 @@
-; docformat = 'rst'
+ ; docformat = 'rst'
 ;+
 ;
@@ -314 +314 @@
 msg = report(['iii : ' + routine + ' : derniers points de la grille (temps) = ' + string(max(hrfin))]);-----------------------------------
 t=indgen(nbjour)*pasthor+hrdeb ; en heures
-print,'t',t
+;print,'t',t
 
 ; initialisation des tableaux de sortie (s valeurs, z poids, kont nb points)
@@ -332 +332 @@
 geomax = geolocation_to_string_idl(lonmax, latmax, look,1)
 
-fovmax=46
-openr,1,project_id_env+'AMSU/AMSU'+numch+'_correcbis_t_JJAS2009.dat'
-moytb=fltarr(fovmax)
-readf,1,moytb,format='(46(f6.2,1X))'
-close,1
-moytbliss=smooth(moytb,5)
-correction_t=moytbliss-mean(moytbliss[21:26])
-;print, correction
-openr,2,project_id_env+'AMSU/AMSU'+numch+'_correcbis_m_JJAS2009.dat'
-moytb=fltarr(fovmax)
-readf,2,moytb,format='(46(f6.2,1X))'
-close,2
-close,/all
-moytbliss=smooth(moytb,5)
-correction_m=moytbliss-mean(moytbliss[22:24])
+; correction supplementaire si necessaire
+correcbis='non'
+if correcbis eq 'oui' then begin
+    fovmax=46
+    openr,1,project_id_env+'AMSU/AMSU'+numch+'_correcbis_t_JJAS2009.dat'
+    moytb=fltarr(fovmax)
+    readf,1,moytb,format='(46(f6.2,1X))'
+    close,1
+    moytbliss=smooth(moytb,5)
+    correction_t=moytbliss-mean(moytbliss[21:26])
+    openr,2,project_id_env+'AMSU/AMSU'+numch+'_correcbis_m_JJAS2009.dat'
+    moytb=fltarr(fovmax)
+    readf,2,moytb,format='(46(f6.2,1X))'
+    close,2
+    close,/all
+    moytbliss=smooth(moytb,5)
+    correction_m=moytbliss-mean(moytbliss[22:24])
+endif
 
 
@@ -393 +396 @@
 
 ; boucle sur les points du fichier
+         tb1=tbinit
+         if correcbis eq "oui" then begin
+             print,"ATTENTION: correction supplementaire activee"
 ; application de la correction bis
-         tb1=fltarr(nn)
-         print,'fauchee',min(fov),max(fov)
-         for ifov=min(fov)+1,fovmax-1 do begin
-            ind=where(fov eq ifov and mask eq 1,nbfov)
-            if nbfov ne 0 then begin
-               tb1[ind]=tbinit[ind] -correction_t[ifov]
-            endif
-         endfor
-         ind=where(fov eq 1,smallfov)
-         if smallfov ne 0 then begin
-            tb1[ind]=tbinit[ind] -correction_t[0]
+             tb1=fltarr(nn)
+             print,'fauchee',min(fov),max(fov)
+             for ifov=min(fov)+1,fovmax-1 do begin
+                 ind=where(fov eq ifov and mask eq 1,nbfov)
+                 if nbfov ne 0 then begin
+                     tb1[ind]=tbinit[ind] -correction_t[ifov]
+                 endif
+             endfor
+             ind=where(fov eq 1,smallfov)
+             if smallfov ne 0 then begin
+                 tb1[ind]=tbinit[ind] -correction_t[0]
+             endif
+             ind=where(fov ge fovmax,bigfov)
+             if bigfov ne 0 then begin
+                 tb1[ind]=tbinit[ind] -correction_t[fovmax-1]
+             endif
+             for ifov=min(fov)+1,fovmax-1 do begin
+                 ind=where(fov eq ifov and mask eq 0,nbfov)
+                 if nbfov ne 0 then begin
+                     tb1[ind]=tbinit[ind]-correction_m[ifov]
+                 endif
+             endfor
+             ind=where(fov eq 1,smallfov)
+             if smallfov ne 0 then begin
+                 tb1[ind]=tbinit[ind] -correction_m[0]
+             endif
+             ind=where(fov ge fovmax,bigfov)
+             if bigfov ne 0 then begin
+                 tb1[ind]=tbinit[ind] -correction_m[fovmax-1]
+             endif
          endif
-         ind=where(fov ge fovmax,bigfov)
-         if bigfov ne 0 then begin
-            tb1[ind]=tbinit[ind] -correction_t[fovmax-1]
-        endif
-         for ifov=min(fov)+1,fovmax-1 do begin
-            ind=where(fov eq ifov and mask eq 0,nbfov)
-            if nbfov ne 0 then begin
-               tb1[ind]=tbinit[ind] -correction_m[ifov]
-            endif
-         endfor
-         ind=where(fov eq 1,smallfov)
-         if smallfov ne 0 then begin
-            tb1[ind]=tbinit[ind] -correction_m[0]
-         endif
-         ind=where(fov ge fovmax,bigfov)
-         if bigfov ne 0 then begin
-            tb1[ind]=tbinit[ind] -correction_m[fovmax-1]
-        endif
-
-
-;         print,'effet de la correction'
-;         print,'moyenne, min, max',mean(tb1-tbinit),min(tb1-tbinit),max(tb1-tbinit)
 
          for i=0,nn-1 do begin

trunk/src/extract_amsua.pro

@@ -332 +332 @@
 s = size(SCOPE_TRACEBACK(/STRUCTURE),/DIMENSION)
 routine = (SCOPE_TRACEBACK(/STRUCTURE))[s-1].Routine
-key_performance=0
+key_performance=1
 
 if key_performance EQ 1 THEN BEGIN
@@ -377 +377 @@
 
 IF key_performance EQ 1 THEN BEGIN
-    msg = report(['ppp : ' + routine + ' :  durï¿œe intermediaire demarrage boucle sur fichier ' $
+    msg = report(['ppp : ' + routine + ' :  duree intermediaire demarrage boucle sur fichier ' $
                    + string(SYSTIME(1)-time1,format='(F12.6)')])
 ENDIF
@@ -386 +386 @@
    COMMON amsua_data  ,ama_scan
    IF key_performance EQ 1 THEN BEGIN
-       msg = report(['ppp : ' + routine + ' : durï¿œe intermediaire ouverture et lecture du fichier' $
+       msg = report(['ppp : ' + routine + ' : duree intermediaire ouverture et lecture du fichier' $
                    + string(SYSTIME(1)-time1,format='(F12.6)')])
    ENDIF
@@ -395 +395 @@
    na_fov=na[1]                 ; nb points dans la fauchee
    nbpix= na_fov
-   nosat=ama_head.h_satid       ; numero du satellite NOAA (15 a 19)
+   nosat=ama_head.h_satid       ; numero du satellite NOAA (15 a 19) ou METOP 2
    print,'NUMERO DU SATELLITE ',nosat
    print,'INSTRUMENT',ama_head.h_instrument
@@ -419 +419 @@
    midpix=fix(nbpix/2)
    amalong=REFORM(ama_scan.latlon[1,*,*]/1.E4)
-   print, 'amalong',amalong
+  ; print, 'amalong',amalong
    amalati=REFORM(ama_scan.latlon[0,*,*]/1.E4)
    amcha=REFORM(ama_scan.btemps[nocanal-1,*,*]/100.)
@@ -452 +452 @@
 ;      endelse
 ;
-; correction nadir des donnees
+; correction nadir des donnees. En principe le 17 est elimine des listes
          isat=0
          if nosat eq 15 then begin
@@ -507 +507 @@
                    + string(SYSTIME(1)-time1,format='(F12.6)')])
         ENDIF
-        ; moych=fltarr(nbpix)
-;nfov=nbpix
-
-;moytb=fltarr(nfov)
-;sigtb=fltarr(nfov)
-;for ifov=0,nfov-1 do begin 
-;   ind=where(landseamask[ifov,*] eq 1 ,nifov) 
-;      moytb[ifov]=mean(ch_nadir[ifov,ind]) 
-;      sigtb[ifov]=stdev(ch_nadir[ifov,ind]) 
-;   endfor 
-;window,2
-;plot,indgen(nfov)+1,moytb,psym=7,xstyle=1,ystyle=1,yrange=[min(moytb)-max(sigtb),max(moytb)+max(sigtb)]
-
-;errplot,moytb-sigtb,moytb+sigtb