Changeset 563

Timestamp:

05/29/12 10:15:03 (12 years ago)

Author:

lelod

Message:

modifs chaine

Location:

trunk/src

Files:

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

trunk/src/correc_amsu.pro

@@ -6 +6 @@
 ; ============
 ;
-;  correc_amsu,numch, yyyyb, mmb, ddb, yyyye, mme, dde, lonmin, lonmax, latmin, latmax
+;  correc_amsu,numch, yyyyb, mmb, ddb, yyyye, mme, dde, lonmin, lonmax, latmin, latmax,surf
 ; DESCRIPTION
 ; ===========
 ;
-; programme d'interpolation avec algorithme de Cressmann'ajustement
-;des donnees nadir (en complément de la correction de FK)
+; programme calculant une correction additionnelle 
+; des donnees nadir (en complément de la correction de FK) sur terre ou
+; mer (surf = 0)
 ;
-;; grille de sortie en long / lat
-; recoit en entree les jours de l'annee jdeb et jfin
-;correspondant au debut et la fin du mois
+; recoit en entree les dates de debut et fin, la zone et surf
 ;
 
@@ -34 +33 @@
 ;    latmin=-5.
 ;    latmax=20.
+;    surf="m"
 
 ;    cresamsu,numch, yyyyb, mmb, ddb $
 ;            , yyyye, mme, dde $
-;            , lonmin, lonmax, latmin, latmax
+;            , lonmin, lonmax, latmin, latmax,surf
 ; 
 ;
@@ -57 +57 @@
 ;
 
-pro correc_amsu, numch, yyyyb, mmb, ddb, yyyye, mme, dde, lonmin, lonmax, latmin, latmax
+pro correc_amsu, numch, yyyyb, mmb, ddb, yyyye, mme, dde, lonmin, lonmax, latmin, latmax,surf
 
 compile_opt idl2, strictarrsubs
@@ -68 +68 @@
 ;
 nparam = N_PARAMS()
-IF (nparam NE 11) THEN BEGIN
+IF (nparam NE 12) THEN BEGIN
     ras = report(['Incorrect number of arguments.' $
           + '!C' $
@@ -77 +77 @@
 tbmin=100
 tbmax=350
-
+isurf=1
+if surf eq 'm' then isurf=0
 prefix=numch
 look = 'filename'
@@ -145 +146 @@
                if nfovorb eq nfov then begin
                   for ifov=0,nfov-1 do begin 
-                     ind=where(fovo eq ifov+1 and masko eq 1 ,nifov) 
+                     ind=where(fovo eq ifov+1 and masko eq isurf ,nifov) 
                      if (nifov ne 0) then begin 
                         sumtb[ifov]=sumtb[ifov]+total(tb[ind])
@@ -169 +170 @@
 sigtb=sqrt((vartb-kont*moytb*moytb)/kont)
 print,'moytb',moytb
+fic_correc=project_id_env+'/AMSU/AMSU'+numch+'_correcbis_'+surf+'.dat'
+openw,1,fic_correc
+printf,1,moytb[*],format='(46(f6.2,1X))'
+
 window,0
 plot,fovv,moytb,psym=7,xstyle=1,ystyle=1,yrange=[min(moytb)-max(sigtb),max(moytb)+max(sigtb)]

trunk/src/cresamsu.pro

@@ -102 +102 @@
 ; ==========
 ;
+;  lelod, fin mai 2012
+;  ajout de correction supplementaire pour ameliorer les champs
+; lecture et application de la correction avant interpolation
 ;  $Id: 
 ;
@@ -330 +333 @@
 
 fovmax=46
-openr,1,project_id_env+'AMSU/AMSU'+numch+'_correcbis.dat'
+openr,1,project_id_env+'AMSU/AMSU'+numch+'_correcbis_t_MJJAS2006.dat'
 moytb=fltarr(fovmax)
 readf,1,moytb,format='(46(f6.2,1X))'
+close,1
 moytbliss=smooth(moytb,5)
-correction=moytbliss-mean(moytbliss[21:30])
+correction_t=moytbliss-mean(moytbliss[21:30])
 ;print, correction
+openr,2,project_id_env+'AMSU/AMSU'+numch+'_correcbis_m_MJJAS2006.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:23])
 
 
@@ -386 +397 @@
          print,'fauchee',min(fov),max(fov)
          for ifov=min(fov)+1,fovmax-1 do begin
-            ind=where(fov eq ifov,nbfov)
+            ind=where(fov eq ifov and mask eq 1,nbfov)
             if nbfov ne 0 then begin
-               tb1[ind]=tbinit[ind] ;-correction[ifov]
+               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[0]
+            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[fovmax-1]
+            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)

trunk/src/extract_amsua.pro

@@ -304 +304 @@
 s = size(SCOPE_TRACEBACK(/STRUCTURE),/DIMENSION)
 routine = (SCOPE_TRACEBACK(/STRUCTURE))[s-1].Routine
-;
+key_performance=0
+
 if key_performance EQ 1 THEN BEGIN
     time1 = SYSTIME(1)
@@ -366 +367 @@
    na_fov=na[1]                 ; nb points dans la fauchee
    nbpix= na_fov
-   nosat=ama_head.h_satid
+   nosat=ama_head.h_satid       ; numero du satellite NOAA (15 a 19)
+   print,'NUMERO DU SATLLITE ',nosat
                                 ; nb pixels dans la fauchee AMSUA
    fov=indgen(nbpix)
@@ -419 +421 @@
 ;
 ; correction nadir des donnees
-         if nosat eq 15 then isat=0
+         isat=0
+         if nosat eq 15 then begin
+             isat=0
+             print,'no satellite',nosat, isat
+         endif
          if nosat eq 16 then isat=1
          if nosat eq 18 then isat=2
@@ -468 +474 @@
                    + 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
+        ; 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