examples file to generate plots in Gibson et al 2016 Frontiers FORWARD paper


wx=3
wy=1
!p.multi=[0,wx,wy]


; Figure 1: examples of FORWARD output

for_drive,'cavmorph',inst='wl',line='pb',thcs=45,cavlength=150,rfilter='NRGF_FILTER',nwinx=wx*512,nwiny=wy*512,winnum=0
for_drive,'cavmorph',inst='wl',line='pb',thcs=45,cavlength=150,gridtype='Carrmap',cmer=0,charsize=.85,nwinx=wx*512,nwiny=wy*512,winnum=0,/noerase
for_drive,'cavmorph',inst='wl',line='pb',thcs=45,cavlength=150,gridtype='user',ruser=dblarr(201)+1.05,thuser=dindgen(201)*.15+30,phuser=dblarr(201)-30.,quantmap=quantmap,/noerase
plot,dindgen(201)*.15+30,alog10(quantmap.data),yrange=[1.3,1.7],xrange=[30.,60.],title='log(pB) vs. colatitude',/noerase


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

; Figure 4: PSI at different wavelengths

for_drive,'psimas',date='2012-01-04',xxmin=-2.5,xxmax=2.5,yymin=-2.5,yymax=2.5,units='PPM',nwinx=wx*512,nwiny=wy*512,winnum=1,colortable=0
for_drive,'psimas',date='2012-01-04',xxmin=-2.5,xxmax=2.5,yymin=-2.5,yymax=2.5,/xrt,line='AL-MESH',usecolor=3,nwinx=wx*512,nwiny=wy*512,winnum=1,/noerase,colortable=0
;
; note if this is run once, it will save a map and then can used commented version
; also note-- there are some pixels that are dark that were not in the paper figure
;  this is because we have changed our way of dealing with situations where there is no
;  good spectral line fit for the synthetic CoMP data. This may occur, for example, 
;  in the very low density/temperature regions in the far field (or, if velocity shifts the line 
;  out of the CoMP spectral band).  We now identify such points as NaN, and they show up as dark
;  or one can change the keyword "NULLDATACOLOR" (in the plot widget) to a color to more easily 
;  identify them.
for_drive,'psimas',date='2012-01-04',xxmin=-2.5,xxmax=2.5,yymin=-2.5,yymax=2.5,units='PPM',/comp,ngrid=256,ngy=256,nwinx=wx*512,nwiny=wy*512,winnum=1,/noerase,colortable=0,/savemap,mapname='psimas_10747_01042012_BIG_FOV'
;for_drive,'psimas',readmap='psimas_10747_01042012_BIG_FOV'
for_drive,'psimas',date='2012-01-04',xxmin=-2.5,xxmax=2.5,yymin=-2.5,yymax=2.5,/aia,nwinx=wx*512,nwiny=wy*512,winnum=1,/noerase,colortable=0


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

; Figure 5: Cavity with cool dense prominence (two populations)

for_drive,'cavmorph',/nougat,thcs=45,cdens=1e10,cff_noug=[.2,8.,0,0,0,0],nougwidth=.025,nougtop_r=.9,pop2t=2,xxmin=0.4,yymin=0.4,yymax=1.1,xxmax=1.1,nwinx=wx*512,nwiny=wy*512,winnum=2
for_drive,'cavmorph',/nougat,thcs=45,cdens=1e10,cff_noug=[.2,8.,0,0,0,0],nougwidth=.025,nougtop_r=.9,pop2t=2,xxmin=0.4,yymin=0.4,yymax=1.1,xxmax=1.1,nwinx=wx*512,nwiny=wy*512,winnum=2,/aia,/noerase


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

; Figure 6: MAS model polarization

;for_drive,'psimas',date='2012-01-04',ngrid=256,ngy=256,/comp,/fieldlines,/stklines,/savemap,mapname='psimas_10747_01042012',nwinx=wx*512,nwiny=wy*512,winnum=3,colortable=0
for_drive,readmap='psimas_10747_01042012',line='StokesI',nwinx=wx*512,nwiny=wy*512,winnum=3,/noerase,colortable=0
for_drive,readmap='psimas_10747_01042012',line='VoI',nwinx=wx*512,nwiny=wy*512,winnum=3,/noerase,colortable=0
for_drive,readmap='psimas_10747_01042012',line='StokesAz',nwinx=wx*512,nwiny=wy*512,winnum=3,/noerase,colortable=0
for_drive,readmap='psimas_10747_01042012',line='LoI',/fieldlines,bscale=-.6,/stklines,pscale=-.2,nwinx=wx*512,nwiny=wy*512,winnum=3,/noerase,colortable=0


; Figure 6: PFSS model polarization


for_drive,'pfssmod',date='2012-01-04',ngrid=256,ngy=256,/comp,/fieldlines,/stklines,/savemap,mapname='pfssmod_10747_01042012',nwinx=wx*512,nwiny=wy*512,winnum=4,colortable=0
for_drive,readmap='pfssmod_10747_01042012',line='VoI',nwinx=wx*512,nwiny=wy*512,winnum=4,/noerase,colortable=0
for_drive,readmap='pfssmod_10747_01042012',line='Az',nwinx=wx*512,nwiny=wy*512,winnum=4,/noerase,colortable=0
for_drive,readmap='pfssmod_10747_01042012',line='LoI',/fieldlines,bscale=-.6,/stklines,pscale=-.2,nwinx=wx*512,nwiny=wy*512,winnum=4,/noerase,colortable=0


; Figure 7: multiwavelength V/I comparison 

for_drive,readmap='psimas_10747_01042012',line='VoI',imin=-0.0001,imax=0.0001,nwinx=wx*512,nwiny=wy*512,winnum=5,colortable=0
;
; The slight difference between this and the figure in the paper is because the default now is to allow
; the PSI MAS velocities to affect the results.  If you add "velimpose=0.0001" to the above, you will get the figure in the paper.
;
; note after running the next command once, can avoid the recalculation by using the commented line below it
;
for_drive,'psimas',date='2012-01-04',ngrid=256,ngy=256,/greencomp,line='VoI',imin=-0.0001,imax=0.0001,nwinx=wx*512,nwiny=wy*512,winnum=5,/noerase,colortable=0,mapname='psimas_10747_01042012_green',/savemap
;for_drive,readmap='psimas_10747_01042012_green',line='VoI',imin=-0.0001,imax=0.0001,nwinx=wx*512,nwiny=wy*512,winnum=5,colortable=0
for_drive,'psimas',date='2012-01-04',ngrid=256,ngy=256,/radio,line='VoI',frequency_MHz=100,imin=-0.01,imax=0.01,nwinx=wx*512,nwiny=wy*512,winnum=5,/noerase,colortable=0
for_drive,'psimas',date='2012-01-04',ngrid=256,ngy=256,/radio,line='VoI',frequency_MHz=1000,imin=-0.01,imax=0.01,nwinx=wx*512,nwiny=wy*512,winnum=5,/noerase,colortable=0

; Figure 10: observations of cavity

for_plotfits,date='2012-01-04',/aia,xxmin=.2,xxmax=.8,yymin=.7,yymax=1.3,occult=-1.05,upoccult=1.29,nwinx=wx*512,nwiny=wy*512,winnum=6,colortable=0
for_plotfits,date='2012-01-04',/comp,xxmin=.2,xxmax=.8,yymin=.7,yymax=1.3,line='LoI',imin=-2.,imax=-1.,nwinx=wx*512,nwiny=wy*512,winnum=6,/noerase,colortable=0
;
; note this next line will only work if the file has been downloaded
; to the current directory
; so it is commented out and the previous line repeated
;
;for_plotfits,/comp,xxmin=.2,xxmax=.8,yymin=.7,yymax=1.3,filename='20120104.194037.comp.1074.dynamics.3.fts',line='DOPPLERVLOS',nwinx=wx*512,nwiny=wy*512,winnum=6,/noerase,colortable=0
for_plotfits,date='2012-01-04',/comp,xxmin=.2,xxmax=.8,yymin=.7,yymax=1.3,line='LoI',imin=-2.,imax=-1.,nwinx=wx*512,nwiny=wy*512,winnum=6,/noerase,colortable=0
for_plotfits,date='2012-01-04',/comp,line='Az',xxmin=.2,xxmax=.8,yymin=.7,yymax=1.3,occult=1.05,upoccult=1.29,nwinx=wx*512,nwiny=wy*512,winnum=6,/noerase,colortable=0,/stklines,stkcolor=1,nstarrow=30,pscale=.5

; Figure 11: model vs data full corona

for_plotfits,date='2012-01-04',/comp,line='LoI',imin=-2.,imax=-1.,nwinx=wx*512,nwiny=wy*512,winnum=7,colortable=0
for_plotfits,date='2012-01-04',/comp,line='Az',occult=-1.05,upoccult=1.29,nwinx=wx*512,nwiny=wy*512,winnum=7,/noerase,colortable=0
for_drive,readmap='psimas_10747_01042012',line='LoI',occult=1.05,upoccult=1.29,nwinx=wx*512,nwiny=wy*512,winnum=7,/noerase,colortable=0
for_drive,readmap='psimas_10747_01042012',line='Az',occult=1.05,upoccult=1.29,nwinx=wx*512,nwiny=wy*512,winnum=7,/noerase,colortable=0

; Figure 12:  model vs data southeast quadrant

for_plotfits,date='2012-01-04',/comp,line='LoI',imin=-2.,imax=-1.,xxmax=0.,yymax=0.,nwinx=wx*512,nwiny=wy*512,winnum=8,colortable=0
for_plotfits,date='2012-01-04',/comp,line='Az',occult=-1.05,upoccult=1.29,xxmax=0.,yymax=0.,nwinx=wx*512,nwiny=wy*512,winnum=8,/noerase,colortable=0
for_drive,readmap='psimas_10747_01042012',line='LoI',occult=1.05,upoccult=1.29,xxmax=0.,yymax=0.,nwinx=wx*512,nwiny=wy*512,winnum=8.,/noerase,colortable=0
for_drive,readmap='psimas_10747_01042012',line='Az',occult=1.05,upoccult=1.29,xxmax=0.,yymax=0.,nwinx=wx*512,nwiny=wy*512,winnum=8.,/noerase,colortable=0

;
; not in paper but interesting: northwest quadrant

for_plotfits,date='2012-01-04',/comp,line='LoI',imin=-2.,imax=-1.,xxmin=0.,yymin=0.,nwinx=wx*512,nwiny=wy*512,winnum=13,colortable=0
for_plotfits,date='2012-01-04',/comp,line='Az',occult=-1.05,upoccult=1.29,xxmin=0.,yymin=0.,nwinx=wx*512,nwiny=wy*512,winnum=13,/noerase,colortable=0
for_drive,readmap='psimas_10747_01042012',line='LoI',occult=1.05,upoccult=1.29,xxmin=0.,yymin=0.,nwinx=wx*512,nwiny=wy*512,winnum=13.,/noerase,colortable=0
for_drive,readmap='psimas_10747_01042012',line='Az',occult=1.05,upoccult=1.29,xxmin=0.,yymin=0.,nwinx=wx*512,nwiny=wy*512,winnum=13.,/noerase,colortable=0

; Figure 13; Low&Hundhausen effects of flux rope on linear polarization

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


for_drive,'lowhund',line='bx',thetao=45.,x_oinput=1.,xxmin=0.6,xxmax=0.9,yymin=0.6,yymax=0.9,/fieldlines,imin=-14,imax=3,nwinx=wx*512,nwiny=wy*512,winnum=9,colortable=0
for_drive,'lowhund',line='bx',thetao=45.,x_oinput=0.,xxmin=0.6,xxmax=0.9,yymin=0.6,yymax=0.9,/fieldlines,imin=-14,imax=3,nwinx=wx*512,nwiny=wy*512,winnum=9,colortable=0,/noerase
for_drive,'lowhund',line='bx',thetao=45.,x_oinput=-.5,xxmin=0.6,xxmax=0.9,yymin=0.6,yymax=0.9,/fieldlines,imin=-14,imax=3,nwinx=wx*512,nwiny=wy*512,winnum=9,colortable=0,/noerase
for_drive,'lowhund',/comp,line='LoI',thetao=45.,x_oinput=1.,xxmin=0.6,xxmax=0.9,yymin=0.6,yymax=0.9,/fieldlines,nwinx=wx*512,nwiny=wy*512,winnum=9,/noerase,colortable=0,imax=-.5,imin=-2.
for_drive,'lowhund',/comp,line='LoI',thetao=45.,x_oinput=0.,xxmin=0.6,xxmax=0.9,yymin=0.6,yymax=0.9,/fieldlines,nwinx=wx*512,nwiny=wy*512,winnum=9,/noerase,colortable=0,imax=-.5,imin=-2.
for_drive,'lowhund',/comp,line='LoI',thetao=45.,x_oinput=-0.5,xxmin=0.6,xxmax=0.9,yymin=0.6,yymax=0.9,/fieldlines,nwinx=wx*512,nwiny=wy*512,winnum=9,colortable=0,imax=-.5,imin=-2.,/noerase
for_drive,'lowhund',/comp,line='Az',thetao=45.,x_oinput=1.,xxmin=0.6,xxmax=0.9,yymin=0.6,yymax=0.9,/fieldlines,nwinx=wx*512,nwiny=wy*512,winnum=9,/noerase,colortable=0
for_drive,'lowhund',/comp,line='Az',thetao=45.,x_oinput=0.,xxmin=0.6,xxmax=0.9,yymin=0.6,yymax=0.9,/fieldlines,nwinx=wx*512,nwiny=wy*512,winnum=9,/noerase,colortable=0
for_drive,'lowhund',/comp,line='Az',thetao=45.,x_oinput=-.5,xxmin=0.6,xxmax=0.9,yymin=0.6,yymax=0.9,/fieldlines,nwinx=wx*512,nwiny=wy*512,winnum=9,/noerase,colortable=0


; Figure 15: Fan synthetic test bed
; 

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

for_drive,'numcube',cubename='$FORWARD_DB/TESTBEDS/fullthermodynamic_erupting_qp_mhd_Fan',xxmin=0.8,xxmax=1.4,yymin=-.5,yymax=.5,units='PPM',nwinx=wx*512,nwiny=wy*512,winnum=10,cuberot=-25.,colortable=0
for_drive,'numcube',cubename='$FORWARD_DB/TESTBEDS/fullthermodynamic_erupting_qp_mhd_Fan',xxmin=0.8,xxmax=1.4,yymin=-.5,yymax=.5,/aia,/occult,imin=-.41,imax=3.5,cuberot=-25,nwinx=wx*512,nwiny=wy*512,winnum=10,/noerase,colortable=0
for_drive,'numcube',cubename='$FORWARD_DB/TESTBEDS/fullthermodynamic_erupting_qp_mhd_Fan',xxmin=0.8,xxmax=1.4,yymin=-.5,yymax=.5,/comp,line='VOI',nwinx=wx*512,nwiny=wy*512,winnum=10,/noerase,cuberot=-25.,colortable=0,ngrid=256,ngy=256
for_drive,'numcube',cubename='$FORWARD_DB/TESTBEDS/fullthermodynamic_erupting_qp_mhd_Fan',xxmin=0.8,xxmax=1.4,yymin=-.5,yymax=.5,/radio,line='VOI',nwinx=wx*512,nwiny=wy*512,winnum=10,/noerase,cuberot=-25.,colortable=0

; Figure 16: noise 20cm vs 150cm

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

for_drive,readmap='psimas_10747_01042012',line='V',/donoise,nwinx=wx*512,nwiny=wy*512,winnum=11
for_drive,readmap='psimas_10747_01042012',line='V',/donoise,aperture=150,nwinx=wx*512,nwiny=wy*512,winnum=11,/noerase

end
;