restore,'waccmx_20170821eclipse.sav',/verbose ; date,ut,z,tn_b,tn_e,wi_b,wi_e,diffu_b,diffu_e,opdt_b,opdt_e, ;dfield_b,dfield_e,dwind_b,dwind_e,prod_b,prod_e,loss_b,loss_e, ; ; _b variables: base case; _e variabiles: eclipase case ;plot z=120,200,275,350 nut=n_elements(ut) ;millstone partial, maximum, and end time tut1=17.+27./60. tut2=18.+46./60. tut3=19.+59./60. ; index=where(abs(ut-tut1) eq min(abs(ut-tut1))) iut=index[0] ; tz=[120,200,275,300] ntz=n_elements(tz) color8 s=1.5 tmp_loss=fltarr(nut) tmp_prod=fltarr(nut) tmp=fltarr(nut) !p.multi=[0,2,2] for i=0,ntz-1 do begin index=where(abs(z-tz[i]) eq min(abs(z-tz[i]))) iz=index[0] for m=0,nut-1 do begin if m eq 0 then tmp_loss[m]=loss_e[iz,m]-loss_e[iz,m] $ else tmp_loss[m]=(loss_e[iz,m]-loss_e[iz,m-1])-(loss_b[iz,m]-loss_b[iz,m-1]) endfor for m=0,nut-1 do begin if m eq 0 then tmp_prod[m]=prod_e[iz,m]-prod_e[iz,m] $ else tmp_prod[m]=(prod_e[iz,m]-prod_e[iz,m-1])-(prod_b[iz,m]-prod_b[iz,m-1]) endfor tmp=tmp_prod-tmp_loss if i le 1 then begin plot,ut,tmp,xtitle='UT(hours)',ytitle='cm!u-3!ns!u-1!n',$ title=string(tz[i],format='(i3)')+' km', $ /xstyle,xrange=[14,22], $ /ystyle,yrange=[-20,20], $ charsize=s,/nodata endif else begin plot,ut,tmp,xtitle='UT(hours)',ytitle='cm!u-3!ns!u-1!n',$ title=string(tz[i],format='(i3)')+' km', $ /xstyle,xrange=[14,22], $ /ystyle,yrange=[-10,10], $ charsize=s,/nodata endelse oplot,ut,tmp,color=5 for m=0,nut-1 do begin if m eq 0 then tmp[m]=dwind_e[iz,m]-dwind_e[iz,m] $ ;differential difference else tmp[m]=(dwind_e[iz,m]-dwind_e[iz,m-1])-(dwind_b[iz,m]-dwind_b[iz,m-1]) endfor oplot,ut,tmp for m=0,nut-1 do begin if m eq 0 then tmp[m]=dfield_e[iz,m]-dfield_e[iz,m] $ else tmp[m]=(dfield_e[iz,m]-dfield_e[iz,m-1])-(dfield_b[iz,m]-dfield_b[iz,m-1]) endfor oplot,ut,tmp,color=1 for m=0,nut-1 do begin if m eq 0 then tmp[m]=diffu_e[iz,m]-diffu_e[iz,m] $ else tmp[m]=(diffu_e[iz,m]-diffu_e[iz,m-1])-(diffu_b[iz,m]-diffu_b[iz,m-1]) endfor oplot,ut,tmp,color=3 oplot,[tut1,tut1],[-400,100],line=1 oplot,[tut2,tut2],[-400,100],line=1 oplot,[tut3,tut3],[-400,100],line=1 endfor xyouts,14.5,7,'prod-loss',color=5 xyouts,14.5,5,'trans_wind' xyouts,14.5,3,'trans_efield',color=1 xyouts,14.5,1,'amb_diff',color=3 end