Changeset 498 for trunk

Timestamp:

03/19/12 15:24:05 (12 years ago)

Author:

lelod

Message:

implementation nouvelle correction et trace

Location:

trunk/src

Files:

• interpol_correc.pro (modified) (1 diff)
• plot_amsu.pro (modified) (1 diff)

trunk/src/interpol_correc.pro

@@ -172 +172 @@
    cor_landa = path + 'CORR_LAND_AMSUA_JUIL2006.DAT'
    cor_seaa = path +'CORR_SEA_AMSUA_JUIL2006.DAT'
-  ;test de l'autre fichier de correction: quel est le bon???
-  ; cor_landa=project_id_env+'AMSU/CORR_LAND_AMSUA.DAT'
-  ; cor_seaa=project_id_env+'AMSU/CORR_SEA_AMSUA.DAT'
+  ;test du nouveau fichier de correction (FK mars 2012)
+  cor_landa=project_id_env+'AMSU/CORR_LAND_AMSUA5_DataFromJune2010.DAT'
+  cor_seaa=project_id_env+'AMSU/CORR_SEA_AMSUA5_DataFromJune2010.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
+  ;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
+   ;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

trunk/src/plot_amsu.pro

@@ -71 +71 @@
          nn=n_elements(tb1)
 
-tbmin=100. & tbmax=300.
-jnd=where(tb1 gt tbmin and tb1 lt tbmax)
+; plot sur continent
+tbmin=241. & tbmax=280.
+jnd=where(tb1 gt tbmin and tb1 lt tbmax and mask eq 1)
 tb=tb1[jnd]
 lolo=lon[jnd]
 lala=lat[jnd]
-
+print,'continent ',min(tb),max(tb)
 window,0
 loadct,39
 plot,lon,lat,xrange=[lon_min,lon_max],yrange=[lat_min,lat_max],xstyle=1,ystyle=1,/nodata
-;ind=where(mask eq 1)
-;oplot,lon[ind],lat[ind],psym=1
-
-
 colscale=intarr(n_elements(tb))
 for j=0,n_elements(tb)-1 do colscale[j]=(tb[j]-min(tb))*255/(max(tb)-min(tb))
- ;endfor
-print,min(colscale),max(colscale)
+plots,lolo,lala,psym=5,symsize=1,color=colscale
+; plot sur ocean
+tbmin=200. & tbmax=280.
+jnd=where(tb1 gt tbmin and tb1 lt tbmax and mask eq 0)
+tb=tb1[jnd]
+lolo=lon[jnd]
+lala=lat[jnd]
+print,'ocean ',min(tb),max(tb)
+window,1
+loadct,39
+plot,lon,lat,xrange=[lon_min,lon_max],yrange=[lat_min,lat_max],xstyle=1,ystyle=1,/nodata
+colscale=intarr(n_elements(tb))
+for j=0,n_elements(tb)-1 do colscale[j]=(tb[j]-min(tb))*255/(max(tb)-min(tb))
 plots,lolo,lala,psym=5,symsize=1,color=colscale
 
+;bas=where (tb le min(tb)+0.5)
+;plots,lolo[bas],lala[bas],color=0,symsize=2,psym=5
+;print,lolo[bas],lala[bas]
 end