;<pre><i> example file for the CAVMORPH temperature and density model. 

	;<a href="http://www.hao.ucar.edu/FORWARD/FOR_SSW/docs/MODELPAPERS/streammorph.pdf">Paper </a>describing 3D streamer model underlying CAVMORPH
	;<a href="http://www.hao.ucar.edu/FORWARD/FOR_SSW/docs/MODELPAPERS/cavmorph.pdf">Paper </a>describing addition of cavity 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>

	;<i>*** Streamer without cavity</i>

for_drive,'cavmorph',line='Dens',/plotlog,ngrid=512,photolength=1000.

	;<i>*** Streamer with cavity --  August 2007 fit 
	; note we choose cmer=60 to position the cavity at the West limb
	; the cavity position is defined with the save file aug07fitnew.genx
	; as about -30</i>

cmer=60.

for_drive,'cavmorph',line='Dens',/plotlog,ngrid=512,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',maxheight=1.3,cmer=cmer

	;
	; <i>try putting it in different quadrants</i>
	;

for_drive,'cavmorph',line='Dens',/plotlog,ngrid=512,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',maxheight=1.3,thcs=135,winnum=3,cavtop_th=135,cmer=cmer

for_drive,'cavmorph',line='Dens',/plotlog,ngrid=512,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',maxheight=1.3,thcs=45,winnum=4,cavtop_th=45,cmer=cmer-180.

for_drive,'cavmorph',line='Dens',/plotlog,ngrid=512,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',maxheight=1.3,thcs=135,winnum=5,cavtop_th=135,cmer=cmer-180.

for_drive,'cavmorph',line='Dens',/plotlog,ngrid=512,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmin=-1.2,xxmax=0.,yymin=0,yymax=1.2,winnum=3,cmer=cmer

	;<i>*** Vary nonradiality (keeping cavity top height and size fixed)</i>

for ia = 0,15 do for_drive,'cavmorph',line='Dens',/plotlog,ngrid=512,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmax=0,xxmin=-1.2,yymin=0,yymax=1.2,cavtop_th=45.+ia,extratitle='cavtop_th='+strtrim(string(45+ia),2),cmer=cmer

	;<i>*** Vary cavity bottom location (keeping cavity top position fixed)</i>

for ch = .1,.5,.05 do for_drive,'cavmorph',line='Dens',/plotlog,ngrid=512,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmax=0,xxmin=-1.2,yymin=0,yymax=1.2,cavheight=ch,extratitle='cavity height='+strtrim(string(ch),2),cmer=cmer

	<i>;*** Now for LOS-integrated observables</i>

!p.multi=0

for_drive,'cavmorph',/xrt,/plotlog,ngrid=128,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmax=0,xxmin=-1.2,yymin=0,yymax=1.2,winnum=1,nwinx=512,nwiny=512,cmer=cmer

for_drive,'cavmorph',/trace,/plotlog,ngrid=128,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmax=0,xxmin=-1.2,yymin=0,yymax=1.2,winnum=2,nwinx=512,nwiny=512,cmer=cmer

for_drive,'cavmorph',/eis,line='fe12_195_119',/plotlog,ngrid=128,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmax=0,xxmin=-1.2,yymin=0,yymax=1.2,winnum=3,nwinx=512,nwiny=512,cmer=cmer

	;
	;<i>*** Add nougat </i>
	;	

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

for_drive,'cavmorph',line='Dens',/plotlog,ngrid=ngr,nwinx=wx*ngr,nwiny=wy*ngr,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmin=-1.,xxmax=-.4,yymin=0.4,yymax=1.,winnum=2,cmer=cmer,/nougat,nougmult=3.,title='enhanced density',charsize=charsize

for_drive,'cavmorph',line='Temp',/plotlog,ngrid=ngr,nwinx=wx*ngr,nwiny=wy*ngr,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmin=-1.,xxmax=-.4,yymin=0.4,yymax=1.,winnum=2,cmer=cmer,/nougat,nougmult=3.,Trim=1.6e6,Tnoug=2.e6,/noerase,title='enhanced temperature',charsize=charsize

	;<i>*** Nougat with XRT lines</i>

wx=3
wy=1
ngr=300
!p.multi=[0,wx,wy]
charsize=1.5

	;***no slope or bang

for_drive,'cavmorph',/xrt,/plotlog,ngrid=ngr,nwinx=wx*ngr,nwiny=wy*ngr,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmin=-1.,xxmax=-.4,yymin=0.4,yymax=1.,winnum=3,cmer=cmer,/nougat,nougmult=3.,title='XRT Be-thin, enhanced density, no slope',mang=0,charsize=charsize

for_drive,'cavmorph',/xrt,/plotlog,ngrid=ngr,nwinx=wx*ngr,nwiny=wy*ngr,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmin=-1.,xxmax=-.4,yymin=0.4,yymax=1.,winnum=3,cmer=cmer,/nougat,Trim=1.6e6,Tnoug=2.e6,/noerase,title='XRT Be-thin, enhanced temperature, no slope',mang=0,charsize=charsize

for_drive,'cavmorph',/xrt,/plotlog,ngrid=ngr,nwinx=wx*ngr,nwiny=wy*ngr,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmin=-1.,xxmax=-.4,yymin=0.4,yymax=1.,winnum=3,cmer=cmer,/nougat,Trim=1.6e6,Tnoug=2.e6,nougmult=3.,title='XRT Be-thin, enhanced dens/temp, no slope',/noerase,mang=0,charsize=charsize

	;***model slope 

for_drive,'cavmorph',/xrt,/plotlog,ngrid=ngr,nwinx=wx*ngr,nwiny=wy*ngr,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmin=-1.,xxmax=-.4,yymin=0.4,yymax=1.,winnum=5,cmer=cmer,/nougat,nougmult=3.,title='XRT Be-thin, enhanced density, with slope',charsize=charsize

for_drive,'cavmorph',/xrt,/plotlog,ngrid=ngr,nwinx=wx*ngr,nwiny=wy*ngr,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmin=-1.,xxmax=-.4,yymin=0.4,yymax=1.,winnum=5,cmer=cmer,/nougat,Trim=1.6e6,Tnoug=2.e6,/noerase,title='XRT Be-thin, enhanced temperature, with slope',charsize=charsize

for_drive,'cavmorph',/xrt,/plotlog,ngrid=ngr,nwinx=wx*ngr,nwiny=wy*ngr,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmin=-1.,xxmax=-.4,yymin=0.4,yymax=1.,winnum=5,cmer=cmer,/nougat,Trim=1.6e6,Tnoug=2.e6,nougmult=3.,title='XRT Be-thin, enhanced dens/temp, with slope',/noerase,charsize=charsize

	;***model slope and Earth's bang for August 2007

for_drive,'cavmorph',/xrt,/plotlog,ngrid=ngr,nwinx=wx*ngr,nwiny=wy*ngr,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmin=-1.,xxmax=-.4,yymin=0.4,yymax=1.,winnum=4,cmer=cmer,/nougat,nougmult=3.,title='XRT Be-thin, enhanced density, with slope, bang',bang=-6.341,charsize=charsize

for_drive,'cavmorph',/xrt,/plotlog,ngrid=ngr,nwinx=wx*ngr,nwiny=wy*ngr,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmin=-1.,xxmax=-.4,yymin=0.4,yymax=1.,winnum=4,cmer=cmer,/nougat,Trim=1.6e6,Tnoug=2.e6,/noerase,title='XRT Be-thin, enhanced temperature, with slope, bang',bang=-6.341,charsize=charsize

for_drive,'cavmorph',/xrt,/plotlog,ngrid=ngr,nwinx=wx*ngr,nwiny=wy*ngr,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmin=-1.,xxmax=-.4,yymin=0.4,yymax=1.,winnum=4,cmer=cmer,/nougat,Trim=1.6e6,Tnoug=2.e6,nougmult=3.,title='XRT Be-thin, enhanced dens/temp, with slope, bang',/noerase,bang=-6.341,charsize=charsize

	;***model slope and Earth's bang for August 2007; vary length

for_drive,'cavmorph',/xrt,/plotlog,ngrid=ngr,nwinx=wx*ngr,nwiny=wy*ngr,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmin=-1.,xxmax=-.4,yymin=0.4,yymax=1.,winnum=6,cmer=cmer,/nougat,Trim=1.6e6,Tnoug=2.e6,nougmult=3.,nouglength=25,title='XRT Be-thin, 25 deg length',bang=-6.341,charsize=charsize

for_drive,'cavmorph',/xrt,/plotlog,ngrid=ngr,nwinx=wx*ngr,nwiny=wy*ngr,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmin=-1.,xxmax=-.4,yymin=0.4,yymax=1.,winnum=6,cmer=cmer,/nougat,Trim=1.6e6,Tnoug=2.e6,nougmult=3.,nouglength=15,title='XRT Be-thin, 15 deg length',/noerase,bang=-6.341,charsize=charsize

for_drive,'cavmorph',/xrt,/plotlog,ngrid=ngr,nwinx=wx*ngr,nwiny=wy*ngr,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmin=-1.,xxmax=-.4,yymin=0.4,yymax=1.,winnum=6,cmer=cmer,/nougat,Trim=1.6e6,Tnoug=2.e6,nougmult=3.,nouglength=5,title='XRT Be-thin, 5 deg length',/noerase,bang=-6.341,charsize=charsize


	;***Nougat with AIA lines

for_drive,'cavmorph',/aia,/plotlog,ngrid=ngr,nwinx=wx*ngr,nwiny=wy*ngr,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmin=-1.,xxmax=-.4,yymin=0.4,yymax=1.,winnum=7,cmer=cmer,/nougat,Trim=1.6e6,Tnoug=2.e6,Tcav=1.2e6,nougmult=3.,title='AIA 171: Tr=1.6,Tn=2,Tc=1.2; Nn=3*Nc=.3*Nr',bang=-6.341,line=171,charsize=charsize

for_drive,'cavmorph',/aia,/plotlog,ngrid=ngr,nwinx=wx*ngr,nwiny=wy*ngr,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmin=-1.,xxmax=-.4,yymin=0.4,yymax=1.,winnum=7,cmer=cmer,/nougat,Trim=1.6e6,Tnoug=2.e6,Tcav=1.2e6,nougmult=3.,title='AIA 193',bang=-6.341,/noerase,line=193,charsize=charsize

for_drive,'cavmorph',/aia,/plotlog,ngrid=ngr,nwinx=wx*ngr,nwiny=wy*ngr,readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',xxmin=-1.,xxmax=-.4,yymin=0.4,yymax=1.,winnum=7,cmer=cmer,/nougat,Trim=1.6e6,Tnoug=2.e6,Tcav=1.2e6,nougmult=3.,title='AIA 335',bang=-6.341,/noerase,line=335,charsize=charsize

        ;<i>can also set a second population of e.g. cool material$
         (so far only set up for analytic model CavMorph, but could be applied to$
           a numerical model -- instructions will be put within $
            http://www.hao.ucar.edu/FORWARD/FOR_SSW/idl/MODELS/NUMCUBE/make_my_cube.pro)  </i>

!p.multi=0

for_drive,'cavmorph',/nougat,thcs=45,pop2T=2,cdens=1e10,cff_noug=[.2,8.,0,0,0,0],nougtop_r=.9,nougwidth=.025,xxmin=.4,yymin=.4,xxmax=1.1,yymax=1.1,winnum=0,/sunedge
for_drive,'cavmorph',/aia,/nougat,thcs=45,pop2T=2,cdens=1e10,cff_noug=[.2,8.,0,0,0,0],nougtop_r=.9,nougwidth=.025,xxmin=.4,yymin=.4,xxmax=1.1,yymax=1.1,winnum=1,/sunedge

for_drive,'cavmorph',/losem,/nougat,thcs=45,nougtop_r=.9,nougwidth=.025,xxmin=.4,yymin=.4,xxmax=1.1,yymax=1.1,winnum=2,/sunedge
for_drive,'cavmorph',/colden,/nougat,thcs=45,nougtop_r=.9,nougwidth=.025,xxmin=.4,yymin=.4,xxmax=1.1,yymax=1.1,winnum=3,/sunedge
for_drive,'cavmorph',/pop2losem,/nougat,thcs=45,pop2T=2,cdens=1e10,cff_noug=[.2,8.,0,0,0,0],nougtop_r=.9,nougwidth=.025,xxmin=.4,yymin=.4,xxmax=1.1,yymax=1.1,winnum=4,/sunedge
for_drive,'cavmorph',/pop2colden,/nougat,thcs=45,pop2T=2,cdens=1e10,cff_noug=[.2,8.,0,0,0,0],nougtop_r=.9,nougwidth=.025,xxmin=.4,yymin=.4,xxmax=1.1,yymax=1.1,winnum=5,/sunedge

	;*** Carrington map examples

	;<i>Shows how streamer (without cavity) varies when positioned at different colatitudes. Note, to avoid discontinuities at the pole, a DC offset is subtracted (the value of the streamer delta-gaussian at the pole -- fmax in cavmorph.pro).  For this reason, when it is actually centered on the pole, the streamer delta completely disappears.</i>

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

ith=0.

for_drive,'cavmorph',line='Dens',ngrid=ngr,gridt='CarrMap',limb='East',charsize=1.5,thcs=ith,rheight=1.1,extratitle='colat='+strtrim(string(ith),2),nwinx=wx*ngr,nwiny=wy*ngr/2.,/plotlog,imin=8.,imax=9.,phcs=70,winnum=3

for ith=30.,90.,30. do for_drive,'cavmorph',line='Dens',ngrid=ngr,gridt='CarrMap',limb='East',charsize=1.5,thcs=ith,rheight=1.1,extratitle='colat='+strtrim(string(ith),2),nwinx=wx*ngr,nwiny=wy*ngr/2.,/plotlog,imin=8.,imax=9.,phcs=70,/noerase

	;<i>Now show how vary with m angle</i>

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

im=0.

for_drive,'cavmorph',line='Dens',ngrid=ngr,gridt='CarrMap',nwinx=wx*ngr,nwiny=wy*ngr/2.,/plotlog,imin=8.,imax=9.,length=50.,extratitle='mang ='+strtrim(String(im),2),limb='East',phcs=110,mang=0.,rheight=1.1

for im=30.,90.,30. do for_drive,'cavmorph',line='Dens',ngrid=ngr,gridt='CarrMap',nwinx=wx*ngr,nwiny=wy*ngr/2.,/plotlog,imin=8.,imax=9.,length=50.,mang=im,/noerase,extratitle='mang ='+strtrim(String(im),2),limb='East',phcs=110,rheight=1.1
 
	;<i>Now show how vary with B angle</i>

!p.multi=0

for ib=-10,10 do for_drive,'cavmorph',line='dens',nwinx=400,nwiny=800,/plotlog,cavtop_r=1.2,yymin=.5,yymax=2.,xxmin=0,xxmax=1.3,thcs=30,bang=ib

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

ib=0.

for_drive,'cavmorph',line='Dens',ngrid=ngr,gridt='CarrMap',nwinx=wx*ngr,nwiny=wy*ngr/2.,/plotlog,imin=8.,imax=9.,length=50.,extratitle='bang ='+strtrim(String(ib),2),limb='West',phcs=110,rheight=1.1

for ib=30.,90.,30. do for_drive,'cavmorph',line='Dens',ngrid=ngr,gridt='CarrMap',nwinx=wx*ngr,nwiny=wy*ngr/2.,/plotlog,imin=8.,imax=9.,length=50.,bang=ib,/noerase,extratitle='bang ='+strtrim(String(ib),2),limb='West',phcs=110,rheight=1.1
 
ib=0.

for_drive,'cavmorph',line='Dens',ngrid=ngr,gridt='CarrMap',nwinx=wx*ngr,nwiny=wy*ngr/2.,/plotlog,imin=8.,imax=9.,length=50.,extratitle='bang ='+strtrim(String(ib),2),limb='East',winnum=3,phcs=110,rheight=1.1

for ib=30.,90.,30. do for_drive,'cavmorph',line='Dens',ngrid=ngr,gridt='CarrMap',nwinx=wx*ngr,nwiny=wy*ngr/2.,/plotlog,imin=8.,imax=9.,length=50.,bang=ib,/noerase,extratitle='bang ='+strtrim(String(ib),2),limb='East',winnum=3,phcs=110,rheight=1.1

	;<i>SPARTAN fit shown in Gibson et al 2003  JGR 108, A12, 1444   (density and pB)
	;Reproduces Figure 4 a-b (top streamer only)</i>

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

for_drive,'cavmorph',line='Dens',/plotlog,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',maxheight=6.,ngrid=ngr,nwinx=wx*ngr,nwiny=wy*ngr,charsize=.5

for_drive,'cavmorph',/plotlog,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',maxheight=6.,ngrid=ngr,nwinx=wx*ngr,nwiny=wy*ngr,/noerase,charsize=.5

	;<i>Reproduces Figure 5  (top streamer only -- note it does not match exactly because it doesnt include the large equatorial streamer)</i>

wx=2
wy=3
ngr=90
!p.multi=[0,wx,wy]

for_drive,'cavmorph',line='Dens',ngrid=ngr,gridt='CarrMap',limb='West',charsize=1.5,rheight=1.25,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,phcs=110,imin=6.5,imax=8.25

for_drive,'cavmorph',ngrid=ngr,gridt='CarrMap',limb='West',charsize=1.5,rheight=1.25,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,/noerase,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,phcs=110,imin=.3,imax=1.4

for_drive,'cavmorph',line='Dens',ngrid=ngr,gridt='CarrMap',limb='West',charsize=1.5,rheight=1.75,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,/noerase,phcs=110,imin=5.,imax=7.2

for_drive,'cavmorph',ngrid=ngr,gridt='CarrMap',limb='West',charsize=1.5,rheight=1.75,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,/noerase,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,phcs=110,imin=-1.,imax=0.

for_drive,'cavmorph',line='Dens',ngrid=ngr,gridt='CarrMap',limb='West',charsize=1.5,rheight=2.5,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,/noerase,phcs=110,imin=4.5,imax=6.3

for_drive,'cavmorph',ngrid=ngr,gridt='CarrMap',limb='West',charsize=1.5,rheight=2.5,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,/noerase,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,phcs=110,imin=-2.5,imax=-.4665

	;<i><Reproduces Figure 6  (top streamer only)</i>

for_drive,'cavmorph',line='Dens',ngrid=ngr,gridt='CarrMap',limb='West',charsize=1.5,rheight=1.25,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,mang=10.,phcs=110,winnum=3,imin=6.5,imax=8.25

for_drive,'cavmorph',ngrid=ngr,gridt='CarrMap',limb='West',charsize=1.5,rheight=1.25,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,/noerase,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,mang=10.,phcs=110,winnum=3,imin=.3,imax=1.4

for_drive,'cavmorph',line='Dens',ngrid=ngr,gridt='CarrMap',limb='West',charsize=1.5,rheight=1.75,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,/noerase,mang=10.,phcs=110,winnum=3,imin=5.,imax=7.2

for_drive,'cavmorph',ngrid=ngr,gridt='CarrMap',limb='West',charsize=1.5,rheight=1.75,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,/noerase,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,mang=10.,phcs=110,winnum=3,imin=-1.,imax=0.

for_drive,'cavmorph',line='Dens',ngrid=ngr,gridt='CarrMap',limb='West',charsize=1.5,rheight=2.5,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50,/noerase,mang=10.,phcs=110,winnum=3,imin=4.5,imax=6.3

for_drive,'cavmorph',ngrid=ngr,gridt='CarrMap',limb='West',charsize=1.5,rheight=2.5,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,/noerase,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,mang=10.,phcs=110,winnum=3,imin=-2.5,imax=-.4665

	;<i>Reproduces Figure 7  (top streamer only)</i>

et=' Bang=7.25 '

for_drive,'cavmorph',ngrid=ngr,gridt='CarrMap',limb='East',charsize=1.5,rheight=1.25,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,mang=10,Bang=7.25,extratitle=et,winnum=4,imin=.3,imax=1.4,phcs=110

for_drive,'cavmorph',ngrid=ngr,gridt='CarrMap',limb='West',charsize=1.5,rheight=1.25,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,/noerase,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,mang=10,Bang=7.25,extratitle=et,winnum=4,imin=.3,imax=1.4,phcs=110

for_drive,'cavmorph',ngrid=ngr,gridt='CarrMap',limb='East',charsize=1.5,rheight=1.75,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,/noerase,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,mang=10,Bang=7.25,extratitle=et,winnum=4,imin=-1.,imax=0.,phcs=110

for_drive,'cavmorph',ngrid=ngr,gridt='CarrMap',limb='West',charsize=1.5,rheight=1.75,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,/noerase,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,mang=10,Bang=7.25,extratitle=et,winnum=4,imin=-1.,imax=0.,phcs=110

for_drive,'cavmorph',ngrid=ngr,gridt='CarrMap',limb='East',charsize=1.5,rheight=2.5,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,/noerase,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,mang=10,Bang=7.25,extratitle=et,winnum=4,imin=-2.,imax=-.4665,phcs=110

for_drive,'cavmorph',ngrid=ngr,gridt='CarrMap',limb='West',charsize=1.5,rheight=2.5,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,/noerase,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,mang=10,Bang=7.25,extratitle=et,winnum=4,imin=-2.,imax=-.4665,phcs=110

	;<i>and for reference</i>

for_drive,'cavmorph',line='Dens',ngrid=ngr,gridt='CarrMap',limb='East',charsize=1.5,rheight=1.25,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,mang=10,Bang=7.25,extratitle=et,winnum=5,imin=6.5,imax=8.25,phcs=110

for_drive,'cavmorph',line='Dens',ngrid=ngr,gridt='CarrMap',limb='West',charsize=1.5,rheight=1.25,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,/noerase,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,mang=10,Bang=7.25,extratitle=et,winnum=5,imin=6.5,imax=8.25,phcs=110

for_drive,'cavmorph',line='Dens',ngrid=ngr,gridt='CarrMap',limb='East',charsize=1.5,rheight=1.75,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,/noerase,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,mang=10,Bang=7.25,extratitle=et,winnum=5,imin=5.,imax=7.2,phcs=110

for_drive,'cavmorph',line='Dens',ngrid=ngr,gridt='CarrMap',limb='West',charsize=1.5,rheight=1.75,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,/noerase,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,mang=10,Bang=7.25,extratitle=et,winnum=5,imin=5.,imax=7.2,phcs=110

for_drive,'cavmorph',line='Dens',ngrid=ngr,gridt='CarrMap',limb='East',charsize=1.5,rheight=2.5,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,/noerase,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,mang=10,Bang=7.25,extratitle=et,winnum=5,imin=4.5,imax=6.3,phcs=110

for_drive,'cavmorph',line='Dens',ngrid=ngr,gridt='CarrMap',limb='West',charsize=1.5,rheight=2.5,nwinx=wx*360,nwiny=wy*360/2.,/plotlog,/noerase,readprams='$FORWARD/MODELS/CAVMORPH/spartanprams.genx',photolength=50.,mang=10,Bang=7.25,extratitle=et,winnum=5,imin=4.5,imax=6.3,phcs=110

end
;</pre>