Changeset 591

Timestamp:

07/19/12 18:47:22 (12 years ago)

Author:

lelod

Message:

simplification plot_amsu pour amsub

File:

• trunk/src/plot_amsu.pro (modified) (3 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)
@@ -88 +74 @@
          lat=result.data.lat
          mask=result.data.landseamask
-         tbinit=result.data.tb
+         tb1=result.data.tb
          ; decodage du nb de points (nn)
-         nn=n_elements(tb1)
-tb1=fltarr(nn)
+         nn=n_elements(tbinit)
 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
+print,'nb orbites',max(desc)
 
-; plot sur continent
-tbmin=241.
-tbmax=280.
-jnd=where(tb1 gt tbmin and tb1 lt tbmax and mask eq 1)
+; plot sur continent orbites du matin
+tbmin=100.
+tbmax=300.
+;jnd=where(tb1 gt tbmin and tb1 lt tbmax and mask eq 1 and hour le 12)
+;plot partout
+jnd=where(tb1 gt tbmin and tb1 lt tbmax and hour le 12)
 tb=tb1[jnd]
 lolo=lon[jnd]
 lala=lat[jnd]
-print,'continent ',min(tb),max(tb)
+;print,'continent ',min(tb),max(tb)
 window,0
 loadct,39
@@ -137 +96 @@
 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
+;mer=0
+;if mer ne 0 then begin
 ; 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=[lonmin,lonmax],yrange=[latmin,latmax],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
+;    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=[lonmin,lonmax],yrange=[latmin,latmax],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
+;endif
 
-;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