Changeset 592

Timestamp:

07/20/12 10:18:27 (12 years ago)

Author:

lelod

Message:

modif plot amsu

File:

• trunk/src/plot_amsu.pro (modified) (1 diff)

trunk/src/plot_amsu.pro

@@ -80 +80 @@
 print,'nb orbites',max(desc)
 
-; plot sur continent orbites du matin
+; plot 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)
+rep=''
+descmin=0 & descmax=max(desc)
+
+print,'no d orbite min et max?'
+read,descmin,descmax
+
+lab_orbite:
+window,0
+loadct,39
+jnd=where(tb1 gt tbmin and tb1 lt tbmax and desc ge descmin and desc le descmax)
 tb=tb1[jnd]
 lolo=lon[jnd]
 lala=lat[jnd]
-;print,'continent ',min(tb),max(tb)
-window,0
-loadct,39
+for ide=descmin,descmax do print,'heure des orbites tracees',mean(hour[ide])
 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
+print,'autres orbites? (oui ou ...)'
+read,rep
+if rep eq 'oui' then goto,lab_orbite
+
+
 ;mer=0
 ;if mer ne 0 then begin