------------------------------------------------ EXAMPLES/README - Example job-files for different MC-code drivers 4.1.2002/18.1.2002 HS ------------------------------------------------ 1) Asymmetry vs longitude 2) Opposition-curve 3) Tilt-curve cd to ~/CODE_TO_DONES/EXAMPLE-directory and give the following source-commands: -------------------------------------------- 1) Azimuthal asymmetry - mcnewd_asyobs.exe -------------------------------------------- Particle positions read from files in subdirectory ../nr450a130n4_100, contains ascii position files nr450a130n4_100.posdatNN, with NN=10 to 49 These examples use posdat-files 40-49 Nphot = 1000 37 azimuth-values For the geometry of Fig. 5 in Franklin et al. transmitted Solar light (270 cpusecs) source asytest1.job & For the geometry of Fig. 5 in Franklin et al. Saturn-shine, assuming blue planet (200 cpusecs) source asytest2.job & Both examples write two output files (to EXAMPLE-directory), e.g. asytest1.res - logfile asytest1.all - data to be read by idl To read asytest1.all file run the idl-program (need to give tek_color command first to define color-palettes) tek_color mc_asyobs_read,'asytest1.all',/plot,print=1 -> plots and stores results to asytest1.all.save file three plots/page (in ccw order) normalized intensity vs. longitude relative to Eastern ansa I/F vs azimuth, total + different scattering orders I/F vs. azimuth, total = average of separate posfiles (all plots for albedo=1) call mc_asyobs_read without parameters -> outputs brief info of keywords variables in save-file explained in MC_ASYOBS_READ.SAVEINFO contains also some examples of plotting I/F vs azimuth for an assumed albedo All necessary idl-subroutines should also be found in EXAMPLES directory - if not please let me know! -------------------------------------------- 2) Opposition effect - mcnewd_oppo.exe -------------------------------------------- Particle positions read from files in subdirectory ../SIZEDIST_DATA This example uses tmodel020.posdat19 size-distribution from 0.20 to 5 meters q=3, N=15548 TAU_DYN=1.26 XL=YL=50m (semi-widths) Bridges epsilon same data is stored into idl-savefile tmodel020.pos19 restore,'../SIZEDIST_DATA/tmodel020.pos19' slice,xx,yy,zz,ss ; shows slice through z=0 2.1) Photometric run BSUN=BOBS=26 degrees ALPHA = 0 - 6 degrees (NA=17 values) Nphot = 5000 HG-phasefunction with g=-0.35 160 cpusecs source oppotest1.job & writes to files oppotest1.res - logfile oppotest1.www - data to be read by idl 2.2) Other photometric run size-distribution from 0.20 to 5 meters q=3, N=15548 BSUN=BOBS=6 degrees instead of 26 290 cpusecs source oppotest2.job & 2.3 Third photometric run size-distribution from 0.50 to 5 meters q=3, N=15548 BSUN=BOBS=6 degrees cpusecs source oppotest3.job & To read oppotest1.www file run the idl-program mc_wread_f,'oppotest1' which stores into idl-file oppotest1.outwsave mc_wread_f,'oppotest2' mc_wread_f,'oppotest3' ;To plot I/F vs alpha: restore,'oppotest1.outwsave' ; testm na*10 table ; testm(*,0) = 1st order-scattering for na-values of alpha ; .... ; testm(*,8) = 9th order-scattering -"- ; iperf = conversion factor alb=0.5 testa=testm(*,0)*0. for i=0,8 do begin & testa=testa+alb^(i+1.)*testm(*,i) & endfor testa=testa*iperf nwin plot,alpha*!radeg,testa,xtitle='phase-angle',ytitle='I/F' restore,'oppotest2.outwsave' alb=0.5 testa=testm(*,0)*0. for i=0,8 do begin & testa=testa+alb^(i+1.)*testm(*,i) & endfor testa=testa*iperf oplot,alpha*!radeg,testa -------------------------------------------- 3) Tilt effect - mcnewd_tilt.exe -------------------------------------------- Particle positions read from files in subdirectory ../SIZEDIST_DATA This eaxmple uses tmodel020.posdat19 size-distribution from 0.20 to 5 meters q=3, N=15548 TAU_DYN=1.26 XL=YL=50m (semi-widths) Bridges epsilon 3.1 Photometric run ALPHA = 6 degrees (one value only) BSUN=BOBS= 1. 2. 4. 6. 9. 12. 15. 18. 21. 24. 26.74 (11 values) Nphot = 1000 Callisto phase-function 85 cpusecs source tilttest1.job & writes to files tilttest1.res - logfile tilttest1.ttt - data to be read by idl To read tilttest1.ttt file run the idl-program mc_tiltread_f,'tilttest1',/plot which stores into idl-file oppotest1.tiltsave and plots the I/F vs B curve for A=1 ;To plot I/F vs B=B': restore,'tilttest1.tiltsave' ; sunlat_tab = B ; testw_tab = contains 11 I/F values now (different B=B', A=1) ; testm_tab = (11,1,10) = (B, ialpha, order) ; these already multiplied by I/F-conversion alb=0.5 testm=reform(testm_tab(*,0,*)) testa=testm*0. for i=0,8 do begin & testa=testa+alb^(i+1.)*testm(*,i) & endfor nwin plot,sunlat_tab,testa,xtitle="B=B'",ytitle='I/F for A=0.5' 3.2 Another example for tilt-curve: This example uses tmodel050.posdat19 size-distribution from 0.50 to 5 meters q=3 N=3444 TAU_DYN=1.26 XL=YL=50m (semi-widths) Bridges epsilon same data is stored into idl-savefile tmodel020.pos19 restore,'../SIZEDIST_DATA/tmodel020.pos19' slice,xx,yy,zz,ss ; shows slice through z=0 Photometric run ALPHA = 6 degrees (one value only) BSUN=BOBS= 1. 2. 4. 6. 9. 12. 15. 18. 21. 24. 26.74 (11 values) Nphot = 1000 Callisto phase-function 85 cpusecs source tilttest1.job & mc_tiltread_f,'tilttest2',/plot restore,'tilttest2.tiltsave' alb=0.5 testm=reform(testm_tab(*,0,*)) testa=testm*0. for i=0,8 do begin & testa=testa+alb^(i+1.)*testm(*,i) & endfor oplot,sunlat_tab,testa,col=2 *************************************************** The commands for oppo and tilt-plots are also collected to oppo_tilt.pro ***************************************************