Changeset 584

Timestamp:

07/10/12 11:58:21 (12 years ago)

Author:

lelod

Message:

amelioration interpol_correc pour amsub

File:

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

trunk/src/interpol_correc.pro

@@ -194 +194 @@
 cor_l=fltarr(nfov,5)
 cor_s=fltarr(nfov,5)
+corl=fltarr(nfova,5)
+cors=fltarr(nfova,5)
 cor_sea15=path+'CORR_SEA_AMSU'+instr+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA15.DAT'
 cor_land15=path+'CORR_LAND_AMSU'+instr+strmid(numch,1,1)+'_DataFromJune2010_Afrique_NOAA15.DAT'
@@ -224 +226 @@
       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)
+      corl[*,nosat]=interpol(a,nfova,/lsquadratic) 
+      cor_l[ndfovb:nfovb-1,nosat]=interpol(corl[ndfova:nfova-1,nosat],fova,fovb)
+      cor_l[0:ndfovb-1,nosat]=interpol(corl[0:ndfova-1,nosat],fov_a,fov_b)
    endelse
    free_lun, lun1
@@ -258 +261 @@
       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)
+      cors[*,nosat]=interpol(a,nfova,/lsquadratic) 
+      cor_s[ndfovb:nfovb-1,nosat]=interpol(cors[ndfova:nfova-1,nosat],fova,fovb)
+      cor_s[0:ndfovb-1,nosat]=interpol(cors[0:ndfova-1,nosat],fov_a,fov_b)
    endelse
    free_lun, lun1