Changeset 566

Timestamp:

05/30/12 09:45:28 (12 years ago)

Author:

lelod

Message:

plot_amsu avec correcbis

File:

• trunk/src/plot_amsu.pro (modified) (2 diffs)

trunk/src/plot_amsu.pro

@@ -63 +63 @@
 prefix=numch
 look = 'filename'
+openr,1,project_id_env+'AMSU/AMSU'+numch+'_correcbis_t_JJAS2009.dat'
+moytb=fltarr(fovmax)
+readf,1,moytb,format='(46(f6.2,1X))'
+close,1
+moytbliss=smooth(moytb,5)
+correction_t=moytbliss-mean(moytbliss[21:26])
+;print, correction
+openr,2,project_id_env+'AMSU/AMSU'+numch+'_correcbis_m_JJAS2009.dat'
+moytb=fltarr(fovmax)
+readf,2,moytb,format='(46(f6.2,1X))'
+close,2
+close,/all
+moytbliss=smooth(moytb,5)
+correction_m=moytbliss-mean(moytbliss[22:24])
 
 geomin = geolocation_to_string_idl(lonmin, latmin, look,1)
@@ -74 +88 @@
          lat=result.data.lat
          mask=result.data.landseamask
-         tb1=result.data.tb
+         tbinit=result.data.tb
          ; decodage du nb de points (nn)
          nn=n_elements(tb1)
+tb1=fltarr(nn)
+print,'fauchee',min(fov),max(fov)
+for ifov=min(fov)+1,fovmax-1 do begin
+   ind=where(fov eq ifov and mask eq 1,nbfov)
+   if nbfov ne 0 then begin
+      tb1[ind]=tbinit[ind] -correction_t[ifov]
+   endif
+endfor
+         
+ind=where(fov eq 1,smallfov)
+if smallfov ne 0 then begin
+   tb1[ind]=tbinit[ind] -correction_t[0]
+endif
+ind=where(fov ge fovmax,bigfov)
+if bigfov ne 0 then begin
+   tb1[ind]=tbinit[ind] -correction_t[fovmax-1]
+endif
+for ifov=min(fov)+1,fovmax-1 do begin
+   ind=where(fov eq ifov and mask eq 0,nbfov)
+   if nbfov ne 0 then begin
+      tb1[ind]=tbinit[ind] -correction_m[ifov]
+   endif
+endfor
+ind=where(fov eq 1,smallfov)
+if smallfov ne 0 then begin
+   tb1[ind]=tbinit[ind] -correction_m[0]
+endif
+ind=where(fov ge fovmax,bigfov)
+if bigfov ne 0 then begin
+   tb1[ind]=tbinit[ind] -correction_m[fovmax-1]
+endif
 
 ; plot sur continent