Changeset 560

Timestamp:

05/25/12 18:10:58 (12 years ago)

Author:

lelod

Message:

update correction nadir a5

Location:

trunk/src

Files:

• extract_amsua.pro (modified) (3 diffs)
• interpol_correc.pro (modified) (4 diffs)

trunk/src/extract_amsua.pro

@@ -155 +155 @@
 ; ==========
 ;
+; - lelod 20120525
+;
+;  revisite de la correction nadir, calculee par satellite
+;  utilisation du no de satellite lu dans le fichier pour choix de la
+;  correction
+;  pb pour Metop: pas lu par read_amsua1c, et pas de donnee AQUA dans
+;  la bdd
+
 ; $Id$
 ;
@@ -411 +419 @@
 ;
 ; correction nadir des donnees
+         if nosat eq 15 then isat=0
+         if nosat eq 16 then isat=1
+         if nosat eq 18 then isat=2
+         if nosat eq 19 then isat=3
          ch_nadir=fltarr(nbpix,nzon)
          landseamask=intarr(nbpix,nzon)+2 ; valeur hors zone selectionnee
@@ -440 +452 @@
                         bate2=bate1[*,yind]
                         if (mean(bate2) le 0.5) then begin
-                           ch_nadir[ifo,isc]=amch[ifo,isc]-cor_s[ifo]
+                           ch_nadir[ifo,isc]=amch[ifo,isc]-cor_s[ifo,isat]
                            landseamask[ifo,isc]=0 ; mer
                         endif else begin
-                           ch_nadir[ifo,isc]=amch[ifo,isc]-cor_l[ifo]
+                           ch_nadir[ifo,isc]=amch[ifo,isc]-cor_l[ifo,isat]
                            landseamask[ifo,isc]=1 ; terre
                         endelse

trunk/src/interpol_correc.pro

@@ -77 +77 @@
 ; ==========
 ;
+; - lelod 20120525
+;
+; analyse de la correction nadir : refonte et calcul par satellite
+; avec adaptation à la region choisie (Afrique de l'ouest)
+; changement a faire dans extract_amsua.pro
 ; $Id$
 ;
@@ -172 +177 @@
 path = project_env + '/src/dataref_amsu/'
 ; lecture des fichiers d'ajustement au nadir pour AMSUA et B
-; cas amsua
+; cas amsua: noaa17 saute - lecture pour les 5 satellites utilises N15, 16,
+;                           18, 19, M02
 if nomcanal eq 'a' then begin
-   cor_landa = path + 'CORR_LAND_AMSUA_JUIL2006.DAT'
-   cor_seaa = path +'CORR_SEA_AMSUA_JUIL2006.DAT'
-   cor_l=fltarr(nfova)
-   cor_s=fltarr(nfova)
-   ;test du nouveau fichier de correction pour a5 (FK avril 2012)
-   if numcanal eq 4 or numcanal eq 5 then begin
-      cor_landa=project_id_env+'AMSU/CORR_LAND_AMSUA'+strmid(numch,1,1)+'_DataFromJune2010_40deg_OrogInf500.DAT'
-      cor_seaa=project_id_env+'AMSU/CORR_SEA_AMSUA'+strmid(numch,1,1)+'_DataFromJune2010_40deg.DAT'
-      openr, lun1,cor_landa, /get_lun, ERROR = error
-      IF (error NE 0) then begin
-         ras = report(['eee : can not open for reading '$
-               + '!C' $
-               + 'code : ' + !ERROR_STATE.MSG $
-               + cor_landa])
-         STOP
-      ENDIF
-      readf, lun1, t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,$
-             t11,t12,t13,t14,t15,t16,t17,t18,t19,t20,$
-             t21,t22,t23,t24,t25,t26,t27,t28,t29,t30
-      a=[t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,$
-         t11,t12,t13,t14,t15,t16,t17,t18,t19,t20,$
-         t21,t22,t23,t24,t25,t26,t27,t28,t29,t30]
-      cor_l=interpol(a,nfova,/lsquadratic) 
-      free_lun, lun1
-      openr, lun1,cor_seaa, /get_lun, ERROR = error
-      IF (error NE 0) then begin
-         ras = report(['eee : can not open for reading '$
-               + '!C' $
-               + cor_seaa])
-         STOP
-      ENDIF
+   cor_l=fltarr(nfova,5)
+   cor_s=fltarr(nfova,5)
+      cor_seaa15=project_id_env+'AMSU/CORR_SEA_AMSU'+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA15.DAT'
+      cor_landa15=project_id_env+'AMSU/CORR_LAND_AMSU'+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA15.DAT'
+      cor_seaa16=project_id_env+'AMSU/CORR_SEA_AMSU'+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA16.DAT'
+      cor_landa16=project_id_env+'AMSU/CORR_LAND_AMSU'+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA16.DAT'
+      cor_seaa18=project_id_env+'AMSU/CORR_SEA_AMSU'+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA18.DAT'
+      cor_landa18=project_id_env+'AMSU/CORR_LAND_AMSU'+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA18.DAT'
+      cor_seaa19=project_id_env+'AMSU/CORR_SEA_AMSU'+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA19.DAT'
+      cor_landa19=project_id_env+'AMSU/CORR_LAND_AMSU'+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA19.DAT'
+      cor_seaM2=project_id_env+'AMSU/CORR_SEA_AMSU'+strmid(numch,1,1)+'_DataFromJune2010_Afrique_METOP.DAT'
+      cor_landM2=project_id_env+'AMSU/CORR_LAND_AMSU'+strmid(numch,1,1)+'_DataFromJune2010_Afrique_METOP.DAT'
+      listcorr_l=[cor_landa15,cor_landa16,cor_landa18,cor_landa19,cor_landM2]
+      for nosat=0,4 do begin
+         openr, lun1,listcorr_l[nosat], /get_lun, ERROR = error
+         IF (error NE 0) then begin
+            ras = report(['eee : can not open for reading '$
+                          + '!C' $
+                          + 'code : ' + !ERROR_STATE.MSG $
+                          + cor_landa])
+            STOP
+         ENDIF
+         readf, lun1, t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,$
+                t11,t12,t13,t14,t15,t16,t17,t18,t19,t20,$
+                t21,t22,t23,t24,t25,t26,t27,t28,t29,t30
+         a=[t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,$
+            t11,t12,t13,t14,t15,t16,t17,t18,t19,t20,$
+            t21,t22,t23,t24,t25,t26,t27,t28,t29,t30]
+         cor_l[nosat,*]=interpol(a,nfova,/lsquadratic) 
+         free_lun, lun1
+      endfor
+     listcorr_s=[cor_seaa15,cor_seaa16,cor_seaa18,cor_seaa19,cor_seaM2]
+      for nosat=0,4 do begin
+         openr, lun1,listcorr_s[nosat], /get_lun, ERROR = error
+         IF (error NE 0) then begin
+            ras = report(['eee : can not open for reading '$
+                          + '!C' $
+                          + cor_seaa])
+            STOP
+         ENDIF
       readf, lun1, t1,t2,t3,t4,t5,$
              t6,t7,t8,t9,t10,$
@@ -219 +234 @@
          t26,t27,t28,t29,t30]
       free_lun, lun1
-      ;_s=a
-      cor_s=interpol(a,nfova,/lsquadratic) ; normalement le dernier canal lu
-    endif else begin
-      openr, lun1,cor_landa, /get_lun, ERROR = error
-      IF (error NE 0) then begin
-         ras = report(['eee : can not open for reading '$
-               + '!C' $
-               + cor_landa])
-         STOP
-      ENDIF
-      for i=0, numcanal-1 do begin
-          readf, lun1, t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,$
-                 t11,t12,t13,t14,t15,t16,t17,t18,t19,t20,$
-                 t21,t22,t23,t24,t25,t26,t27,t28,t29,t30
-          a=[t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,$
-             t11,t12,t13,t14,t15,t16,t17,t18,t19,t20,$
-             t21,t22,t23,t24,t25,t26,t27,t28,t29,t30]
-       endfor
-       ;cor_l=a
-       cor_l=interpol(a,nfova,/lsquadratic)    ; normalement le dernier canal lu
-       free_lun, lun1
-       openr, lun1,cor_seaa, /get_lun, ERROR = error
-      IF (error NE 0) then begin
-         ras = report(['eee : can not open for reading '$
-               + '!C' $
-               + cor_seaa])
-         STOP
-      ENDIF
-       for i=0, numcanal-1 do begin
-          readf, lun1, t1,t2,t3,t4,t5,$
-                 t6,t7,t8,t9,t10,$
-                 t11,t12,t13,t14,t15,$
-                 t16,t17,t18,t19,t20,$
-                 t21,t22,t23,t24,t25,$
-                 t26,t27,t28,t29,t30
-    
-          a=[t1,t2,t3,t4,t5,$
-             t6,t7,t8,t9,t10,$
-             t11,t12,t13,t14,t15,$
-             t16,t17,t18,t19,t20,$
-             t21,t22,t23,t24,t25,$
-             t26,t27,t28,t29,t30]
-       endfor
-       free_lun, lun1
-       ;_s=a
-       cor_s=interpol(a,nfova,/lsquadratic) ; normalement le dernier canal lu
-    endelse
+         cor_s[nosat,*]=interpol(a,nfova,/lsquadratic) 
+         free_lun, lun1
+      endfor
+
 endif
 
@@ -281 +253 @@
       STOP
    ENDIF
-for i=0, numcanal-1 do begin
+   for i=0, numcanal-1 do begin
       readf, lun1, t1,t2,t3,t4,t5,$
              t6,t7,t8,t9,t10,$