;
examples file for the LITESLOW model.
;Paper describing Lites&Low model
;If this is your first time using for_drive, then read
;Setup instructions
;*** slices in planes parallel to plane of sky
ymn=-.6
ymx=.6
xmn=.8
xmx=2.
for ixoffset=0,13 do for_drive,'liteslow',ngrid=128,sig1=90,nwinx=512,nwiny=512,xoffset=-.3+float(ixoffset)*0.05,extratitle='X-offset='+strtrim(float(-.3+float(ixoffset)*0.05),2)+' out=10',mapname='LitesLow_BX_X-offset='+strtrim(float(ixoffset),2),/gif,charsize=0.8,line='Bx',imin=-3.,imax=3.,/fieldlines,bscale=0.3,out=10,yymin=ymn,yymax=ymx,xxmin=xmn,xxmax=xmx,/pos,narrow=20,nstarrow=20
;***COMP****
;***FIRST SET THINGS UP***
ymn=-.6
ymx=.6
xmn=.8
xmx=2.
useout=0.01
usedens=0.001
;***sig90***
;**planar slices
for ixoffset=0,6 do for_drive,'liteslow',ngrid=128,sig1=90,nwinx=512,nwiny=512,xoffset=-.3+float(ixoffset)*0.05,extratitle='X-offset='+strtrim(float(-.3+float(ixoffset)*0.05),2)+'sig=90',mapname='LitesLow_sig90_COMP_LOI_X-offset='+strtrim(float(ixoffset),2),/gif,charsize=0.8,line='LOI',/comp,/pos,/fieldlines,bscale=0.,/stklines,/savemap,/plotlog,pscale=0.,out=useout,densscale=usedens,yymin=ymn,yymax=ymx,xxmin=xmn,xxmax=xmx,narrow=20,nstarrow=20
;**full integration
; warning it takes a while
for_drive,'liteslow',ngrid=32,line='LOI',nwinx=512,nwiny=512,/fieldlines,/comp,/verbose,/stklines,/savemap,mapname='liteslowcomp_sig90_LOI_int',sig1=90,losuse='xlos',extratitle='sigma=90',/gif,out=useout,densscale=usedens,yymin=ymn,yymax=ymx,xxmin=xmn,xxmax=xmx,narrow=20,nstarrow=20
;***sig0***
;**planar slices
for ixoffset=0,6 do for_drive,'liteslow',ngrid=128,sig1=0.,nwinx=512,nwiny=512,xoffset=-.6+float(ixoffset)*0.1,extratitle='X-offset='+strtrim(float(-.6+float(ixoffset)*0.1),2)+'sig=0',mapname='LitesLow_sig0_COMP_LOI_X-offset='+strtrim(float(ixoffset),2),/gif,charsize=0.8,line='LOI',/comp,/pos,/fieldlines,bscale=0,/stklines,/savemap,/plotlog,pscale=0,out=useout,densscale=usedens,yymin=ymn,yymax=ymx,xxmin=xmn,xxmax=xmx,narrow=20,nstarrow=20
;**full integration
for_drive,'liteslow',ngrid=32,line='LOI',nwinx=512,nwiny=512,/fieldlines,/comp,/verbose,/stklines,/savemap,mapname='liteslowcomp_sig0_LOI_int',sig1=0,losuse='xlos',extratitle='sigma=0',/gif,out=useout,densscale=usedens,yymin=ymn,yymax=ymx,xxmin=xmn,xxmax=xmx,narrow=20,nstarrow=20
;*** PLAY WITH SLICES****
;***sig90***
;***POS***
for_drive,'liteslow',nwinx=512,nwiny=512,readmap='LitesLow_sig90_COMP_LOI_X-offset=6.00000.sav',/gif,charsize=0.8,/comp,line='LOI',/stklines,pscale=0,/plotlog,/fieldlines,bscale=0
for_drive,'liteslow',nwinx=512,nwiny=512,readmap='LitesLow_sig0_COMP_LOI_X-offset=6.00000.sav',/gif,charsize=0.8,/comp,line='L',/stklines,pscale=0,/fieldlines,bscale=0.
;***ALL SLICES***
for ixoffset=0,6 do for_drive,'liteslow',nwinx=512,nwiny=512,readmap='LitesLow_sig90_COMP_LOI_X-offset='+strtrim(float(ixoffset),2)+'.sav',mapname='LitesLow_sig90_COMP_LOI_X-offset='+strtrim(float(ixoffset),2),/gif,charsize=0.8,/comp,line='LOI',/stklines,pscale=0,/fieldlines,bscale=0,/plotlog,narrow=20,nstarrow=20
;***sig0***
;***POS***
for_drive,'liteslow',nwinx=512,nwiny=512,readmap='LitesLow_sig0_COMP_LOI_X-offset=6.00000.sav',mapname='LitesLow_POS_sig0_COMP_LOI_int',/gif,charsize=0.8,/comp,line='LOI',/stklines,/plotlog,/fieldlines,bscale=0,pscale=0,narrow=20,nstarrow=20
for_drive,'liteslow',nwinx=512,nwiny=512,readmap='LitesLow_sig0_COMP_LOI_X-offset=6.00000.sav',mapname='LitesLow_POS_sig0_COMP_L',/gif,charsize=0.8,/comp,line='L',/stklines,pscale=0,/fieldlines,bscale=0.,narrow=20,nstarrow=20
for_drive,'liteslow',nwinx=512,nwiny=512,readmap='LitesLow_sig0_COMP_LOI_X-offset=6.00000.sav',mapname='LitesLow_POS_sig0_COMP_LoI',/gif,charsize=0.8,/comp,line='LoI',/stklines,pscale=0.,/plotlog,/fieldlines,bscale=0,narrow=20,nstarrow=20
;***ALL SLICES***
;****PLAY WITH full integration
;***sig 90
for_drive,'liteslow',line='L',nwinx=512,nwiny=512,/fieldlines,/comp,/stklines,readmap='liteslowcomp_sig90_LOI_int.sav',/gif,winnum=6,mapname='liteslowcomp_sig90_L_xlos_fieldlines',pscale=0,bscale=0,narrow=20,nstarrow=20
for_drive,'liteslow',line='LoI',nwinx=512,nwiny=512,/fieldlines,/comp,/stklines,readmap='liteslowcomp_sig90_LOI_int.sav',/gif,winnum=6,/plotlog,mapname='liteslowcomp_sig90_LoI_xlos_fieldlines',pscale=0,bscale=0,narrow=20,nstarrow=20
;***sig 0
for_drive,'liteslow',line='L',nwinx=512,nwiny=512,/fieldlines,/comp,/stklines,readmap='liteslowcomp_sig0_LOI_int.sav',/gif,winnum=6,mapname='liteslowcomp_sig90_L_xlos_fieldlines',pscale=0,bscale=0,narrow=20,nstarrow=20
for_drive,'liteslow',line='LoI',nwinx=512,nwiny=512,/fieldlines,/comp,/stklines,readmap='liteslowcomp_sig0_LOI_int.sav',/gif,winnum=6,/plotlog,mapname='liteslowcomp_sig90_LoI_xlos_fieldlines',pscale=0,bscale=0,narrow=20,nstarrow=20
end
;