;<pre><i> Point of view

	;In these examples we will use the CAVMORPH model to give use something to look at.
	;We will also only be plotting density in the POS becasue it is fast.</i>

	;<i>plot using a grid of 512 </i>

ngrid=512

	;<i>By default the models are centered on Carrington Longitude 0, 
	;so if you to see a feature there on the west limb set</i>

cmer=-90 ;<i>(this is actually the default value for CMER)</i>

	;<i>plotting density ('dens') in the plane of the sky </i>

for_drive,'cavmorph',line='Dens',/plotlog,ngrid=ngrid,cmer=cmer

	;<i>We can also plot on the solar disk using the keyword DODISK.</i>

for_drive,'cavmorph',line='Dens',/plotlog,ngrid=ngrid,cmer=cmer,/DODISK

	;<i>This is just showing the density at r=1 R_sun. this makes more 
	;sense for other observables, which we will get to latter.</i>

	;<i>Plot a Cavity
	;Kind of dull, so let's define the cavity a bit more. We could 
	;do this with keywords, ;but for the moment let's read in a preset 
	;parameter file:</i>

for_drive,'cavmorph',line='Dens',/plotlog,ngrid=ngrid,cmer=cmer,$
		readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx'
 		
	;<i>Lets look at this from a different point of view, 
	;+60 degrees from the central meridian, rather than -90:</i>

cmer=60.

for_drive,'cavmorph',line='Dens',/plotlog,ngrid=ngrid,cmer=cmer,$
		readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx'
 		
	;<i>let's leave the plane of the ecliptic and look at it from an angle 30	
	;above the ecliptic: </i>

Bang=-40

for_drive,'cavmorph',line='Dens',/plotlog,ngrid=ngrid,cmer=cmer,Bang=Bang,$
		readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx'
 
	;<i>Looking at the Sun from a particular date from the Earth, say</i>

date='9-Aug-2007'		

for_drive,'cavmorph',line='Dens',/plotlog,ngrid=ngrid,date=date,$
		readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx'
 
	;<i>Plot only a section of the sun</i>

Bang=5
xxmin=-1.3 & xxmax=-.2
yymin=.2 & yymax=1.3

for_drive,'cavmorph',line='Dens',/plotlog,ngrid=ngrid,cmer=cmer,BAng=BAng,$
		readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',$
		xxmin=xxmin,xxmax=xxmax,yymin=yymin,yymax=yymax
 
	;<i> Carrington map examples
	;The default grid is a simple plane of the sky grid, 
	;but another option is a carrington projection based on 
	;data from one of the limbs.</i>

GridType='CarrMap'

	;<i>use data 0.1 Rsun from the east limb:</i>

rheight=1.1
limb='east'

for_drive,'cavmorph',line='Dens',/plotlog,ngrid=ngrid,cmer=cmer,BAng=BAng,$
		readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',$
		GridType=GridType,rheight=rheight,limb=limb
		
	;<i>try a different point of view</i>

cmer=-30

for_drive,'cavmorph',line='Dens',/plotlog,ngrid=ngrid,cmer=cmer,BAng=BAng,$
		readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',$
		GridType=GridType,rheight=rheight,limb=limb

	;<i>out of the ecliptic:</i>

Bang=25
for_drive,'cavmorph',line='Dens',/plotlog,ngrid=ngrid,cmer=cmer,BAng=BAng,$
		readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',$
		GridType=GridType,rheight=rheight,limb=limb

	;<i>If you want you can obtain the model data at particular
	;locations by specifying the 3D coordinates.
	;altitude in solar radii</i>

ruser=findgen(10)/100.+1.01   

	;<i>colatitude/heliographic latitude (not polar angle) in degrees</i>

thuser=replicate(1.,10)		

	;<i>longitude in spherical/heliographic coordinates (not central meridian)</i>

phuser=replicate(45.,10)
gridt='USERINPUT'

for_drive,'cavmorph',line='Dens',$
		readprams='$FORWARD/MODELS/CAVMORPH/aug07fitnew.genx',$
		gridt=gridt,ruser=ruser,thuser=thuser,phuser=phuser,quant=density

plot,ruser,density.data,title=density.id,xtitle='R (R_Sun)',color=0

end
;</pre>