Changeset 472 for trunk

Timestamp:

12/15/11 23:16:23 (12 years ago)

Author:

lelod

Message:

suite correc nadir

Location:

trunk/src

Files:

• interpol_correc.pro (modified) (1 diff)
• test_correc.pro (modified) (5 diffs)

trunk/src/interpol_correc.pro

@@ -211 +211 @@
 ;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_landb=path+'AMSU/CORR_LAND_AMSUB_JUIL2006.DAT'
+   ;cor_seab=path+'AMSU/CORR_SEA_AMSUB_JUIL2006.DAT'
+   cor_landb=path+'AMSU/CORR_LAND_AMSUB.DAT'
+   cor_seab=path+'AMSU/CORR_SEA_AMSUB.DAT'
    cor_l=fltarr(nfovb)
    cor_s=fltarr(nfovb)

trunk/src/test_correc.pro

@@ -7 +7 @@
 reads, strmid(numch,1,1),numcanal,format='(i1.1)'
 
-if nomcanal eq 'a' then nbpix=30
-if nomcanal eq 'b' then nbpix=90
+nbpix=30
 
 path=project_id_env
@@ -55 +54 @@
    cor_s2006=a    ; normalement le dernier canal lu
 
-   cor_landa=project_id_env+'AMSU/CORR_LAND_AMSUA_AMSUB.DAT'
-   cor_seaa=project_id_env+'AMSU/CORR_SEA_AMSUA_AMSUB.DAT'
+   cor_landa=project_id_env+'AMSU/CORR_LAND_AMSUA.DAT'
+   cor_seaa=project_id_env+'AMSU/CORR_SEA_AMSUA.DAT'
    cor_l=fltarr(nbpix)
    cor_s=fltarr(nbpix)
@@ -95 +94 @@
    free_lun, lun1
    cor_sold=a    ; normalement le dernier canal lu
-window,10
-tek_color
-plot, indgen(nbpix),cor_l2006,yrange=[-15,5],psym=2,color=0
-oplot,indgen(nbpix),cor_lold,psym=2,color=3
-oplot,indgen(nbpix),cor_s2006,psym=4,color=5
-oplot,indgen(nbpix),cor_sold,psym=4,color=8
 
 endif
 
-
-
-
-
 ;cas amsub
 if nomcanal eq 'b' then begin
+   print,'canal',numch
    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
@@ -130 +118 @@
          t26,t27,t28,t29,t30]
    endfor
+   cor_l2006=a 
    free_lun, lun1
    openr, lun1,cor_seab, /get_lun
@@ -146 +135 @@
          t26,t27,t28,t29,t30]
    endfor
+   cor_s2006=a 
+
+   cor_landb=path+'AMSU/CORR_LAND_AMSUB.DAT'
+   cor_seab=path+'AMSU/CORR_SEA_AMSUB.DAT'
+    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
+   cor_lold=a 
+
+   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
+   cor_sold=a 
 endif
+window,10
+tek_color
+plot, indgen(nbpix),cor_l2006,yrange=[-15,15],psym=2,color=0
+oplot,indgen(nbpix),cor_lold,psym=2,color=3
+oplot,indgen(nbpix),cor_s2006,psym=4,color=5
+oplot,indgen(nbpix),cor_sold,psym=4,color=8
 
 end