;<pre><i>  examples file for the GIBLOW model.

    	;<a href="http://www.hao.ucar.edu/FORWARD/FOR_SSW/docs/MODELPAPERS/giblow98.pdf">Paper </a>describing Gibson&Low model

	;If this is your first time using for_drive, then read
        ;<a href = "http://www.hao.ucar.edu/FORWARD/FOR_SSW/docs/INSTRUCTIONS/examples0_setup.html">Setup instructions</a> </i>

	;<i>**Plane of the sky slice examples</i>

!p.multi=0
for_drive,'giblow',line='Dens',/plotlog,ngrid=512,out=1.
for_drive,'giblow',line='Dens',/plotlog,ngrid=512,out=1.,/bonly,winnum=3

	;<i>* note bonly doesnt solve full mhd eqns, but puts in isothermal temperature
	;* and density depletion in bubble</i>

for_drive,'giblow',line='Temp',/plotlog,limb='East',out=1.,ngrid=512,cmer=90.,extratitle='cmer=90',chars=.8

for_drive,'giblow',line='Pres',/plotlog,limb='West',out=1.,ngrid=512,thetao=50.,chars=.8

	;<i>**shows what the bubble looks like rotating into and the out of view (plane of sky slice of magnetic field strength)</i>

for ic=-70.,-110.,-5. do for_drive,'giblow',line='Bmag',out=1.,ngrid=512,cmer=ic,extratitle='Central Meridian='+strtrim(string(fix(ic)),2),imin=.0,imax=1.5,chars=.8,/dodisk,plotlog=0,colortable=0

	;<i>**Example of varying a model parameter  (in this case, time)</i>

for itime=0.,5.,1. do for_drive,'giblow',line='Bmag',out=1.,ngrid=512,extratitle='time='+strtrim(string(itime),2),time=itime,velmult=.1

	;<i>**Magnetic field overplotted with plane of sky direction, oriented along line of sight (x-axis)</i>

wx=4
wy=1
ngr=512
!p.multi=[0,wx,wy]

for_drive,'giblow',line='Bmag',imin=0.,imax=1.5,yymin=-1,yymax=1.,xxmin=.8,xxmax=1.8,ngrid=ngr,nwinx=wx*ngr/2,nwiny=wy*ngr,/fieldlines,narrow=20,bscale=0,out=0,sig=90,cc=8,plotlog=0,colortable=0

for_drive,'giblow',line='Bx',imin=-0.5,imax=0.5,yymin=-1,yymax=1.,xxmin=.8,xxmax=1.8,ngrid=ngr,nwinx=wx*ngr/2,nwiny=wy*ngr,/fieldlines,/noerase,bscale=0,narrow=20,sig=90,cc=8,plotlog=0,colortable=0

for_drive,'giblow',line='By',imin=-0.5,imax=0.5,yymin=-1,yymax=1.,xxmin=.8,xxmax=1.8,ngrid=ngr,nwinx=wx*ngr/2,nwiny=wy*ngr,/fieldlines,/noerase,bscale=0,narrow=20,sig=90,cc=8,plotlog=0,colortable=0

for_drive,'giblow',line='Bz',imin=-0.5,imax=0.5,yymin=-1,yymax=1.,xxmin=.8,xxmax=1.8,ngrid=ngr,nwinx=wx*ngr/2,nwiny=wy*ngr,/fieldlines,/noerase,bscale=0,narrow=20,sig=90,cc=8,plotlog=0,colortable=0

	;<i>**shows what the bubble looks like as bubble neutral line rotates from parallel to limb to perpendicular and back to parallel (rotation about y axis)</i>

wx=4
wy=1
ngr=512
!p.multi=[0,wx,wy]

ib = 0.01

for_drive,'giblow',line='Bx',sig=ib,extratitle='Sigma='+strtrim(string(ib),2),imin=-0.5,imax=.5,yymin=-1,yymax=1.,xxmin=.8,xxmax=1.8,ngrid=ngr,nwinx=wx*ngr/2.,nwiny=wy*ngr,/fieldlines,out=0,winnum=5,cc=8,plotlog=0,colortable=0,narrow=20

for ib=30.01,90.01,30. do for_drive,'giblow',line='Bx',sig=ib,extratitle='Sigma='+strtrim(string(ib),2),imin=-0.5,imax=.5,charsize=1.5,yymin=-1,yymax=1.,xxmin=.8,xxmax=1.8,ngrid=ngr,nwinx=wx*ngr/2.,nwiny=wy*ngr,/noerase,/fieldlines,out=0,cc=8,plotlog=0,colortable=0,narrow=20

ib = 0.01

for_drive,'giblow',line='Bz',sig=ib,extratitle='Sigma='+strtrim(string(ib),2),imin=-0.5,imax=.5,yymin=-1,yymax=1.,xxmin=.8,xxmax=1.8,ngrid=ngr,nwinx=wx*ngr/2.,nwiny=wy*ngr,/fieldlines,out=0,winnum=6,cc=8,plotlog=0,colortable=0,narrow=20

for ib=30.01,90.01,30. do for_drive,'giblow',line='Bz',sig=ib,extratitle='Sigma='+strtrim(string(ib),2),imin=-0.5,imax=.5,charsize=1.5,yymin=-1,yymax=1.,xxmin=.8,xxmax=1.8,ngrid=ngr,nwinx=wx*ngr/2.,nwiny=wy*ngr,/noerase,/fieldlines,out=0,cc=8,plotlog=0,colortable=0,narrow=20

	;<i>** Dens shows core structure</i>

wx=4
wy=1
ngr=512
!p.multi=[0,wx,wy]

ib = 0.01

for_drive,'giblow',line='Dens',out=1.,sig=ib,extratitle='Sigma='+strtrim(string(ib),2),imin=5d7,imax=5d8,yymin=-1,yymax=1.,xxmin=.8,xxmax=1.8,ngrid=ngr,nwinx=wx*ngr/2.,nwiny=wy*ngr,/fieldlines,narrow=20,winnum=6,plotlog=0,colortable=0,/bonly

for ib=30.01,90.01,30. do for_drive,'giblow',line='Dens',out=1.,sig=ib,extratitle='Sigma='+strtrim(string(ib),2),imin=5d7,imax=5d8,charsize=1.5,yymin=-1,yymax=1.,xxmin=.8,xxmax=1.8,ngrid=ngr,nwinx=wx*ngr/2.,nwiny=wy*ngr,/noerase,narrow=20,/fieldlines,plotlog=0,colortable=0,/bonly


	;<i>**shows what the bubble looks like as viewing angle changes from ecliptic to pole - note extrema (B=0,90) correspond to sigma (0,90)</i>

wx=4
wy=1
ngr=512
!p.multi=[0,wx,wy]

ib = 0.01

for_drive,'giblow',line='Bx',bang=ib,extratitle='B-angle='+strtrim(string(ib),2),imin=-0.5,imax=.5,yymin=-1,yymax=1.,xxmin=.8,xxmax=1.8,ngrid=ngr,nwinx=wx*ngr/2.,nwiny=wy*ngr,narrow=20,/fieldlines,bscale=0,winnum=3,cc=9,colortable=0,plotlog=0,docont=0

for ib=30.01,90.01,30. do for_drive,'giblow',line='Bx',bang=ib,extratitle='B-angle='+strtrim(string(ib),2),imin=-0.5,imax=.5,charsize=1.5,yymin=-1,yymax=1.,xxmin=.8,xxmax=1.8,ngrid=ngr,nwinx=wx*ngr/2.,nwiny=wy*ngr,/noerase,narrow=20,/fieldlines,bscale=0,cc=9,colortable=0,plotlog=0,docont=0

ib = 0.01

for_drive,'giblow',line='Bz',bang=ib,extratitle='B-angle='+strtrim(string(ib),2),imin=-0.5,imax=.5,yymin=-1,yymax=1.,xxmin=.8,xxmax=1.8,ngrid=ngr,nwinx=wx*ngr/2.,nwiny=wy*ngr,narrow=20,/fieldlines,bscale=0,winnum=1,cc=9,colortable=0,plotlog=0,docont=0

for ib=30.01,90.01,30. do for_drive,'giblow',line='Bz',bang=ib,extratitle='B-angle='+strtrim(string(ib),2),imin=-0.5,imax=.5,charsize=1.5,yymin=-1,yymax=1.,xxmin=.8,xxmax=1.8,ngrid=ngr,nwinx=wx*ngr/2.,nwiny=wy*ngr,/noerase,/fieldlines,bscale=0,cc=9,colortable=0,plotlog=0,docont=0,narrow=20

	;<i>***view from North pole- see how changing carrington longitude affects location of structure (Carrington longitude fixed coord)</i>

!p.multi=0

for ic=360.,180.,-30. do for_drive,'giblow',line='Br',ngrid=512,out=.2,cmer=ic,bang=90.001,extratitle='cmer='+strtrim(string(fix(ic)),2)+' Bang=90',/gif

	;<i>(*note gif files will be saved locally and can be made into movie)</i>
	;<i>**see how thetao, Bang, and cmer all work together 
	;** can be compared to tecplot 3D field line projections</i>

sig=90
bang=30
thetao=90
cmer=-90

for_drive,'giblow',line='bmag',/fieldlines,cc=8,ngrid=512,yymin=-.6,yymax=.6,xxmin=0.8,xxmax=2.,narrow=30,bscale=.5,nwinx=650,nwiny=650,charsize=.8,extratitle='sig='+strtrim(string(sig),2)+' bang='+strtrim(string(bang),2)+' thetao='+strtrim(string(thetao),2)+' cmer='+strtrim(string(cmer),2),sig=sig,bang=bang,thetao=thetao,cmer=cmer

sig=90
bang=10
thetao=80
cmer=-80

for_drive,'giblow',line='bmag',/fieldlines,cc=8,ngrid=512,yymin=-.4,yymax=.8,xxmin=0.8,xxmax=2.,narrow=30,bscale=.5,nwinx=650,nwiny=650,charsize=.8,extratitle='sig='+strtrim(string(sig),2)+' bang='+strtrim(string(bang),2)+' thetao='+strtrim(string(thetao),2)+' cmer='+strtrim(string(cmer),2),sig=sig,bang=bang,thetao=thetao,cmer=cmer,winnum=3

	;<i>***phio is a somewhat redundant parameter with cmer, its uses is more clear in next example, but this shows from north pole how it affects structure location</i>

ip=73

for_drive,'giblow',line='Br',ngrid=512,out=.2,cmer=0,phio=ip,bang=90.001,extratitle='cmer=0, phio='+strtrim(string(fix(ip)),2)+' Bang=90'

	;<i>***Carrington map with bubble positioned in long/colat</i>

ip=73
it=133

for_drive,'giblow',line='Dens',ngrid=512,phio=ip,thetao=it,extratitle='colat='+strtrim(string(fix(it)),2)+' phio='+strtrim(string(fix(ip)),2),gridt='CarrMap',out=1,/bonly

	;<i>***Carrington map with changing colatitude</i>

wx=1
wy=4
ngr=360
!p.multi=[0,wx,wy]

ith=0.

for_drive,'giblow',line='Dens',ngrid=ngr,gridt='CarrMap',limb='East',charsize=1.5,thetao=ith,rheight=1.1,extratitle='colat='+strtrim(string(ith),2),nwinx=wx*ngr,nwiny=wy*ngr/2.,imin=5d7,imax=5d8,phio=55,out=1,plotlog=0,/Bonly

for ith=30.,90.,30. do for_drive,'giblow',line='Dens',ngrid=ngr,gridt='CarrMap',limb='East',charsize=1.5,thetao=ith,rheight=1.1,extratitle='colat='+strtrim(string(ith),2),/noerase,nwinx=wx*ngr,nwiny=wy*ngr/2.,imin=5d7,imax=5d8,phio=55,out=1,plotlog=0,/Bonly

	;<i>***Carrington map with changing B angle</i>

wx=1
wy=4
ngr=360
!p.multi=[0,wx,wy]

ith=0.

for_drive,'giblow',line='Dens',ngrid=ngr,gridt='CarrMap',limb='East',charsize=1.5,thetao=80.,rheight=1.1,extratitle='colat=80.',nwinx=wx*ngr,nwiny=wy*ngr/2.,imin=5d7,imax=5d8,winnum=2,phio=155,out=1,/Bonly,plotlog=0

for ith=30.,90.,30. do for_drive,'giblow',line='Dens',ngrid=ngr,gridt='CarrMap',limb='East',charsize=1.5,thetao=80.,rheight=1.1,extratitle='colat=80.'+'bang='+strtrim(string(fix(ith)),2),/noerase,nwinx=wx*ngr,nwiny=wy*ngr/2.,bang=ith,imin=5d7,imax=5d8,phio=155,out=1,/Bonly,plotlog=0

ith=0.

for_drive,'giblow',line='Dens',ngrid=ngr,gridt='CarrMap',limb='West',charsize=1.5,thetao=80.,rheight=1.1,extratitle='colat=80.',nwinx=wx*ngr,nwiny=wy*ngr/2.,imin=5d7,imax=5d8,winnum=3,phio=155,out=1,/Bonly,plotlog=0

for ith=30.,90.,30. do for_drive,'giblow',line='Dens',ngrid=ngr,gridt='CarrMap',limb='West',charsize=1.5,thetao=80.,rheight=1.1,extratitle='colat=80.'+'bang='+strtrim(string(fix(ith)),2),/noerase,nwinx=wx*ngr,nwiny=wy*ngr/2.,bang=ith,imin=5d7,imax=5d8,phio=155,out=1,/Bonly,plotlog=0

	;<i>***finally, integration along los (pB is default)</i>

	;<i>***disk</i>

wx=2
wy=1
ngr=128
!p.multi=[0,wx,wy]

ib = 0.01

for_drive,'giblow',sig=ib,extratitle='Sigma='+strtrim(string(fix(ib)),2),yymin=-1,yymax=1.,xxmin=.8,xxmax=1.8,ngrid=ngr/2,nwinx=wx*512,nwiny=wy*512,winnum=5,/plotlog,charsize=.7,cc=8,out=1,/bonly

ib = 90.01

for_drive,'giblow',sig=ib,extratitle='Sigma='+strtrim(string(fix(ib)),2),yymin=-1,yymax=1.,xxmin=.8,xxmax=1.8,ngrid=ngr/2,nwinx=wx*512,nwiny=wy*512,winnum=5,/plotlog,charsize=.7,/overplot,/noerase,cc=8,mapname='savepB',/savemap,out=1,/bonly

	;<i>***test the reading back in of a saved file</i>

for_drive,'giblow',readmap='savepB.sav',winnum=4,nwinx=512,nwiny=512,/plotlog

	
	;<i>***integration of Comp observable </i>

	;<i>** to make from scratch:</i>

for_drive,'giblow',/comp,nwinx=512,nwiny=512,xxmin=0.7,xxmax=2.5,yymin=-1.,yymax=1.,sig1=90.,line='I',charsize=0.7,winnum=1,/fieldlines,/savemap,mapname='giblowcomp',/Bonly,/verbose,ngr=64

	;<i>*** to read from file plot stoke lines:</i>
!p.multi=0

for_drive,'giblow',charsize=1.2,line='Q',/comp,/Bonly,bscale=0,/fieldlines,pscale=0.,readmap='giblowcomp.sav',/stklines,winnum=2,nwinx=512,nwiny=512,colortable=0

for_drive,'giblow',charsize=1.2,line='U',/comp,/Bonly,bscale=0,/fieldlines,pscale=0.,readmap='giblowcomp.sav',/stklines,winnum=3,nwinx=512,nwiny=512,colortable=0

for_drive,'giblow',charsize=1.2,line='V',/comp,/Bonly,bscale=0,/fieldlines,pscale=0.,readmap='giblowcomp.sav',/stklines,winnum=4,nwinx=512,nwiny=512,colortable=0

for_drive,'giblow',charsize=1.7,line='LoI',/comp,/Bonly,bscale=0,/fieldlines,pscale=0.,readmap='giblowcomp.sav',/stklines,winnum=6,nwinx=512,nwiny=512

;</pre>
end