Changeset 522

Timestamp:

04/16/12 18:04:38 (12 years ago)

Author:

pinsard

Message:

correction for correction reading if a5 vs a

File:

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

trunk/src/interpol_correc.pro

@@ -80 +80 @@
 ;
 ; $URL$
+;
+; - fplod 20120416
+;
+;   * correction for correction reading if a5 vs a
 ;
 ; - fplod 20120106
@@ -172 +176 @@
    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)
+   cor_l=fltarr(nfova)
+   cor_s=fltarr(nfova)
+   ;test du nouveau fichier de correction pour a5 (FK avril 2012)
    if numch eq 'a5' then begin
       cor_landa=project_id_env+'AMSU/CORR_LAND_AMSUA5_DataFromJune2010_40deg_OrogInf500.DAT'
@@ -186 +190 @@
          t21,t22,t23,t24,t25,t26,t27,t28,t29,t30]
       cor_l=interpol(a,nfova,/lsquadratic) 
-   endif
-      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
-   cor_l=interpol(a,nfova,/lsquadratic)    ; normalement le dernier canal lu
-   free_lun, lun1
-   openr, lun1,cor_seaa, /get_lun
-   ;for i=0, numcanal-1 do begin
+      free_lun, lun1
+      openr, lun1,cor_seaa, /get_lun
       readf, lun1, t1,t2,t3,t4,t5,$
              t6,t7,t8,t9,t10,$
@@ -214 +205 @@
          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
+      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
+      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
+       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
 endif