Changeset 455 for trunk

Timestamp:

12/11/11 21:53:53 (12 years ago)

Author:

lelod

Message:

MAJ interpolcarrec

File:

• trunk/src/interpol_correc.pro (modified) (1 diff)

trunk/src/interpol_correc.pro

@@ -110 +110 @@
 ;
 @cm_project
-
+; mettre au carre!!!!!!!!!
+nomcanal=canal(0)
+numcanal=str2num(canal,1)
 
 pixelsize,pixatot,pixbtot,alongatot,alongbtot
 
-
-fovb=fltarr(45)
-fova=fltarr(15)
-alongb=fltarr(45)
-alongab=alongbtot[45:89] ; verifier l'intervalle!!!!!
-pixb=pixbtot[45:89]
-alonga=alongatot[15:29]
-pixa=pixatot[15:29]
-
-;ifov nadir : diametre
+; parametres AMSU et calcul de la fauchee
+nfova=30
+nfovb=90
+ndfova=nfova/2
+ndfovb=nfovb/2
+fovb=fltarr(ndfovb)
+fova=fltarr(ndfova)
+alongb=fltarr(ndfovb)
+alongab=alongbtot[ndfovb:nfovb-1] ; verifier l'intervalle!!!!!
+pixb=pixbtot[ndfovb:nfovb]
+alonga=alongatot[ndfova:nfova]
+pixa=pixatot[ndfova:nfova]
+;ifov nadir : diametre (km)
 ifova=48.05 & ifovb=16
-
 ; position dans la fauchee: premier point au nadir a ifov/2
 ;
 coefb=0.8
 fovb(0)=ifovb/2.
-for i=1,44 do fovb(i)=fovb(i-1)+coefb*(pixb(i)/2 +pixb(i)/2)
-fov_b=fltarr(45)
+for i=1,ndfovb-1 do fovb(i)=fovb(i-1)+coefb*(pixb(i-1)/2 +pixb(i)/2)
+fov_b=fltarr(ndfovb)
 ;fov_b(44)=fovb(0)
 ;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,44 do fov_b(i)=-fovb(44-i)
-
+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,14 do fova(i)=fova(i-1)+coefa*(pixa(i)/2 +pixa(i)/2)
-fov_a=fltarr(15)
+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,14 do fov_a(i)=-fova(14-i)
-
-;lecture des fichiers de correction amsub
+for i=0,ndfova-1 do fov_a(i)=-fova(ndfova-1-i)
+
+
 path=project_id_env
 ; lecture des fichiers d'ajustement au nadir pour AMSUA et B
-cor_landa=project_id_env+'AMSU/CORR_LAND_AMSUA_JUIL2006.DAT'
-cor_seaa=project_id_env+'AMSU/CORR_SEA_AMSUA_JUIL2006.DAT'
-cor_la=fltarr(9,30)
-cor_sa=fltarr(9,30)
-openr, lun1,cor_landa, /get_lun
-for i=0, 8 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]
-
-    cor_la[i,*]=a
-endfor
-free_lun, lun1
-openr, lun1,cor_seaa, /get_lun
-for i=0, 8 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]
-    cor_sa[i,*]=a
-endfor
-free_lun, lun1
-
-cor_landb=path+'AMSU/CORR_LAND_AMSUB_JUIL2006.DAT'
-cor_seab=path+'AMSU/CORR_SEA_AMSUB_JUIL2006.DAT'
-    cor_lb=fltarr(5,90)
-    cor_sb=fltarr(5,90)
-    openr, lun1,cor_landb, /get_lun
-    for i=0, 4 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]
-     ;interpolation
-        cor_lb(i,45:89)=interpol(a(15:29),fova,fovb,/spline)
-        cor_lb(i,0:44)=interpol(a(0:14),fov_a,fov_b,/spline)
-;print,'i',i,a
-;print,cor_lb(i,*)
-    endfor
-    free_lun, lun1
-    openr, lun1,cor_seab, /get_lun
-    for i=0, 4 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]
-        ;interpolation
-        cor_sb(i,45:89)=interpol(a(15:29),fova,fovb,/spline)
-        cor_sb(i,0:44)=interpol(a(0:14),fov_a,fov_b,/spline)
-;print,'i',i,a
-;print,cor_sb(i,*)
-    endfor
-    free_lun, lun1
+; cas amsua
+if nomcanal eq 'a' then begin
+   cor_landa=project_id_env+'AMSU/CORR_LAND_AMSUA_JUIL2006.DAT'
+   cor_seaa=project_id_env+'AMSU/CORR_SEA_AMSUA_JUIL2006.DAT'
+   cor_l=fltarr(nfova)
+   cor_s=fltarr(nfova)
+   openr, lun1,cor_landa, /get_lun
+   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    ; normalement le dernier canal lu
+   free_lun, lun1
+   openr, lun1,cor_seaa, /get_lun
+   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
+   cor_s=a    ; normalement le dernier canal lu
+endif
+
+;cas amsub
+if nomcanal eq 'b' then begin
+   cor_landb=path+'AMSU/CORR_LAND_AMSUB_JUIL2006.DAT'
+   cor_seab=path+'AMSU/CORR_SEA_AMSUB_JUIL2006.DAT'
+   cor_l=fltarr(nfovb)
+   cor_s=fltarr(nfovb)
+   openr, lun1,cor_landb, /get_lun
+   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)
+
+   free_lun, lun1
+   openr, lun1,cor_seab, /get_lun
+   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)
+   free_lun, lun1
+endif
 
 end