read_eg_ecmwf,lat,lon,t,q pos=[0.15,0.15,0.9,0.8] set_plot,'x' window,0,retain=2 map_set,/mercator , limit=[0,-30,40,40],position=pos map_continents, limit=[0,-30,40,40],position=pos wvmr_levs=indgen(11)*2. loadct,39 contour,q(*,*,20,0)*1000.,lon,lat,levels=wvmr_levs,/fill,/overplot map_continents, limit=[0,-30,40,40],position=pos,thick=3 t_levs=indgen(25)*5.+270. contour,t(*,*,20,0),lon,lat,/overplot,levels=t_levs,c_labels=t_levs*0.+1 contour,wvmr_levs#(intarr(11)+1),wvmr_levs,indgen(11),position=[0.15,0.88,0.9,0.93],levels=wvmr_levs,xrange=[0,20],/fill,/noerase,ystyle=4,title='WVMR (gkg!E-1!N)',xstyle=1 contour,wvmr_levs#(intarr(11)+1),wvmr_levs,indgen(11),levels=wvmr_levs,ystyle=4,/overplot write_png,'~lecjm/test_png_out.png',tvrd(true=1)