Changeset 580 for trunk

Timestamp:

06/20/12 15:30:15 (12 years ago)

Author:

lelod

Message:

extension a amsub

File:

• trunk/src/interpol_correc.pro (modified) (5 diffs)

trunk/src/interpol_correc.pro

@@ -60 +60 @@
 ; ====
 ;
+;
 ; on peut changer le codage en dur des fauchees ou demi-fauchees, mais
 ; programme adapte a AMSUA et AMSUB, pas MHS (90 points, mais fauchee
@@ -76 +77 @@
 ; EVOLUTIONS
 ; ==========
+;
+; - lelod 20120620
+;   amelioration du code et introduction des corrections par
+;   satellite pour amsub
+; impasse sur difference amsub - mhs a ce niveau
+; remis les donnees de correction dans data_raf_amsu pour update VARAMMA
 ;
 ; - lelod 20120525
@@ -129 +136 @@
 reads, strmid(numch,1,1),numcanal,format='(i1.1)'
 
-if nomcanal eq 'a' then nbpixel=30
-if nomcanal eq 'b' then nbpixel=90
 pixelsize,pixatot,pixbtot,alongatot,alongbtot
 if nomcanal eq 'a' then begin
@@ -140 +145 @@
 endelse
 
-; parametres AMSU et calcul de la fauchee
+; parametres AMSU et calcul de la geometrie de la fauchee
+; pourinterpolation de la correction (calculee sur geometrie AMSUA) a AMSUB
 nfova=30
 nfovb=90
-ndfova=nfova/2
-ndfovb=nfovb/2
-print,'nb points',numch,ndfova,ndfovb
-fovb=fltarr(ndfovb)
+ndfova=nfova/2 ; 15
+ndfovb=nfovb/2 ; 45
+; cas AMSUA
 fova=fltarr(ndfova)
-alongb=fltarr(ndfovb)
-alongab=alongbtot[ndfovb:nfovb-1] ; verifier l'intervalle!!!!!
-pixb=pixbtot[ndfovb:nfovb-1]
 alonga=alongatot[ndfova:nfova-1]
 pixa=pixatot[ndfova:nfova-1]
 ;ifov nadir : diametre (km)
-ifova=48.05 & ifovb=16
+ifova=48.05
 ; position dans la fauchee: premier point au nadir a ifov/2
-;
+; zone pixel : ellipse de grand axe pixb et petit axe alongb
+coefa=1 ; recouvrement des pixels dans la fauchee
+fova(0)=ifova/2.
+for i=1,ndfova-1 do fova(i)=fova(i-1)+coefa*(pixa(i-1)/2 +pixa(i)/2) ; verifier!!!
+fov_a=fltarr(ndfova)
+;fov_a(14)=fova(0)
+;for i=0,13 do fov_a(13-i)=fov_a(14-i)+coefa*(pixa(13-i)/2 +pixa(14-i)/2)
+for i=0,ndfova-1 do fov_a(i)=-fova(ndfova-1-i)
+
+fovb=fltarr(ndfovb)
+alongb=fltarr(ndfovb)
+alongb=alongbtot[ndfovb:nfovb-1] ; verifier l'intervalle!!!!!
+pixb=pixbtot[ndfovb:nfovb-1]
+;ifov nadir : diametre (km)
+ifovb=16.
 coefb=0.8
 fovb(0)=ifovb/2.
@@ -164 +180 @@
 ;for i=0,43 do  fov_b(43-i)=fov_b(44-i)+coefb*(pixb(43-i)/2 +pixb(44-i)/2)
 for i=0,ndfovb-1 do fov_b(i)=-fovb(ndfovb-1-i)
-; zone pixel : ellipse de grand axe pixb et petit axe alongb
-;
-coefa=1
-fova(0)=ifova/2.
-for i=1,ndfova-1 do fova(i)=fova(i-1)+coefa*(pixa(i-1)/2 +pixa(i)/2) ; verifier!!!
-fov_a=fltarr(ndfova)
-;fov_a(14)=fova(0)
-;for i=0,13 do fov_a(13-i)=fov_a(14-i)+coefa*(pixa(13-i)/2 +pixa(14-i)/2)
-for i=0,ndfova-1 do fov_a(i)=-fova(ndfova-1-i)
 
 
 path = project_env + '/src/dataref_amsu/'
 ; lecture des fichiers d'ajustement au nadir pour AMSUA et B
-; cas amsua: noaa17 saute - lecture pour les 5 satellites utilises N15, 16,
-;                           18, 19, M02
+; noaa17 elimine - lecture pour les 5 satellites utilises N15, 16, 18, 19, M02
 if nomcanal eq 'a' then begin
-   cor_l=fltarr(nfova,5)
-   cor_s=fltarr(nfova,5)
-      cor_seaa15=project_id_env+'AMSU/CORR_SEA_AMSUA'+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA15.DAT'
-      cor_landa15=project_id_env+'AMSU/CORR_LAND_AMSUA'+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA15.DAT'
-      cor_seaa16=project_id_env+'AMSU/CORR_SEA_AMSUA'+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA16.DAT'
-      cor_landa16=project_id_env+'AMSU/CORR_LAND_AMSUA'+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA16.DAT'
-      cor_seaa18=project_id_env+'AMSU/CORR_SEA_AMSUA'+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA18.DAT'
-      cor_landa18=project_id_env+'AMSU/CORR_LAND_AMSUA'+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA18.DAT'
-      cor_seaa19=project_id_env+'AMSU/CORR_SEA_AMSUA'+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA19.DAT'
-      cor_landa19=project_id_env+'AMSU/CORR_LAND_AMSUA'+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA19.DAT'
-      cor_seaM2=project_id_env+'AMSU/CORR_SEA_AMSUA'+strmid(numch,1,1)+'_DataFromJune2010_Afrique_METOP.DAT'
-      cor_landM2=project_id_env+'AMSU/CORR_LAND_AMSUA'+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
-          corrland=listcorr_l[nosat]
-         openr, lun1,corrland, /get_lun, ERROR = error
-         IF (error NE 0) then begin
-            ras = report(['eee : can not open for reading '$
-                          + '!C' $
-                          + 'code : ' + !ERROR_STATE.MSG $
-                          + corrland])
-            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
-          corrsea=listcorr_s[nosat]
-          openr, lun1,corrsea, /get_lun, ERROR = error
-          IF (error NE 0) then begin
-            ras = report(['eee : can not open for reading '$
-                          + '!C' $
-                          + corrsea])
-            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]
-      free_lun, lun1
-         cor_s[*,nosat]=interpol(a,nfova,/lsquadratic) 
-         free_lun, lun1
-      endfor
-
-endif
-
-;cas amsub
-if nomcanal eq 'b' then begin
-   cor_landb = path + 'CORR_LAND_AMSUB_JUIL2006.DAT'
-   cor_seab = path + 'CORR_SEA_AMSUB_JUIL2006.DAT'
-   cor_l=fltarr(nfovb)
-   cor_s=fltarr(nfovb)
-   openr, lun1,cor_landb, /get_lun, ERROR = error
+   instr='A'
+   nfov=nfova
+endif else begin
+   instr='B'
+   nfov=nfovb
+endelse
+cor_l=fltarr(nfov,5)
+cor_s=fltarr(nfov,5)
+cor_sea15=path+'AMSU/CORR_SEA_AMSU'+instr+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA15.DAT'
+cor_land15=path+'AMSU/CORR_LAND_AMSU'+instr+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA15.DAT'
+cor_sea16=path+'AMSU/CORR_SEA_AMSU'+instr+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA16.DAT'
+cor_land16=path+'AMSU/CORR_LAND_AMSU'+instr+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA16.DAT'
+cor_sea18=path+'AMSU/CORR_SEA_AMSU'+instr+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA18.DAT'
+cor_land18=path+'AMSU/CORR_LAND_AMSU'+instr+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA18.DAT'
+cor_sea19=path+'AMSU/CORR_SEA_AMSU'+instr+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA19.DAT'
+cor_land19=path+'AMSU/CORR_LAND_AMSU'+instr+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA19.DAT'
+cor_seaM2=path+'AMSU/CORR_SEA_AMSU'+instr+strmid(numch,1,1)+'_DataFromJune2010_Afrique_METOP.DAT'
+cor_landM2=path+'AMSU/CORR_LAND_AMSU'+instr+strmid(numch,1,1)+'_DataFromJune2010_Afrique_METOP.DAT'
+listcorr_l=[cor_land15,cor_land16,cor_land18,cor_land19,cor_landM2]
+for nosat=0,4 do begin
+   corrland=listcorr_l[nosat]
+   openr, lun1,corrland, /get_lun, ERROR = error
    IF (error NE 0) then begin
       ras = report(['eee : can not open for reading '$
-            + '!C' $
-            + cor_landb])
+                    + '!C' $
+                    + 'code : ' + !ERROR_STATE.MSG $
+                    + corrland])
       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
-                       ;interpolation
-   cor_l[ndfovb:nfovb-1]=interpol(a[ndfova:nfova-1],fova,fovb,/spline)
-   cor_l[0:ndfovb-1]=interpol(a[0:ndfova-1],fov_a,fov_b,/spline)
-
+   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]
+   if nomcanal='a' then begin
+      cor_l[*,nosat]=interpol(a,nfova,/lsquadratic) 
+   endif else begin
+      cor_l[ndfovb:nfovb-1,nosat]=interpol(a[ndfova:nfova-1],fova,fovb,/spline)
+      cor_l[0:ndfovb-1,nosat]=interpol(a[0:ndfova-1],fov_a,fov_b,/spline)
+   endelse
    free_lun, lun1
-   openr, lun1,cor_seab, /get_lun, ERROR = error
+endfor
+
+listcorr_s=[cor_sea15,cor_sea16,cor_sea18,cor_sea19,cor_seaM2]
+for nosat=0,4 do begin
+   corrsea=listcorr_s[nosat]
+   openr, lun1,corrsea, /get_lun, ERROR = error
    IF (error NE 0) then begin
       ras = report(['eee : can not open for reading '$
-            + '!C' $
-            + cor_seab])
+                    + '!C' $
+                    + corrsea])
       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
-                       ;interpolation
-   cor_s[ndfovb:nfovb]=interpol(a[ndfova:nfova-1],fova,fovb,/spline)
-   cor_s[0:ndfovb-1]=interpol(a[0:ndfova-1],fov_a,fov_b,/spline)
+   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]
    free_lun, lun1
-endif
+   
+   if nomcanal='a' then begin
+      cor_s[*,nosat]=interpol(a,nfova,/lsquadratic) 
+   endif else begin
+      cor_s[ndfovb:nfovb-1,nosat]=interpol(a[ndfova:nfova-1],fova,fovb,/spline)
+      cor_s[0:ndfovb-1,nosat]=interpol(a[0:ndfova-1],fov_a,fov_b,/spline)
+   endelse
+   free_lun, lun1
+endfor
 
 end