#------------------------------------------------------------------ #------------------------------------------------------------------ #1) create directory for results rm -r ~hsalo/SUNBOX/RESULTS/ide2007_e50_r300n mkdir ~hsalo/SUNBOX/RESULTS/ide2007_e50_r300n #------------------------------------------------------------------ #------------------------------------------------------------------ #2) go to this directory + copy job-file there cd ~hsalo/SUNBOX/RESULTS/ide2007_e50_r300n cp ~hsalo/SUNBOX/JOBS/ide2007_e50_r300n.job ./ #------------------------------------------------------------------ #------------------------------------------------------------------ #3) run parameters in the order they are read in: #RUNID = run-identification #NLEN = number of characters in 'RUNID' #ROO1,ROO2 = internal density of particles if ROO1=ROO2 # For study of hysteresis: roo increases linearly # from roo1(at orb=0) to roo2(at ORB=ORBMAX/2) and then back #moonlet related (see below): #------------------------------------ #RHO_MOON = moonlet density #RAD_MOON = moomlet radius #CEN_MOON = treatment of moonlet center #ORB_CROSS = collect boundary crossings for ORB>ORB_CROSS #N_CROSS = total number of crossing collected #I_BOUNDARY how to treat boundaries #N_BOUNDARY number of fresh particles used #------------------------------------ #ISEARCH_ORDER close pair searched either in x(1) or y(2) # direction: use 1 for xl>yl, 2 for yl> xl # related to fourier-mode tabulation # NMODES, AIDL # DORB_FIDL # #LKM = number of particles #ETA = Saturnocentric distance (meters) #NC = LKM (must be set to same value as above) #CX1, CX2 = initial x-limits of the distribution #CI, CE = initial rms inclination and eccentricity #XL, YL = half-size of the calculation box #SIG,LKM = radius of particles; LKM must again be same value as above) #ORBMAX = maximum number of orbits simulated #DT = timestep (in orbital periods) #PIDL = interval of runid.pidl output (position data) #MODEL = 1 -> eps= constant = lamda # 2 -> eps= eps-bridges (v_n/v0) #LAMDA,V0 #ENCLIM = outer limiting distance for inclusion of self-gravity via FFT # <0 -> in units of XL #ENCMIN = inner limiting distance for FFT3 # inside this use pairwise calculation #NGX,NYG = number of subdivision in FFT #------------------------------------------------ # MODIFICATIONS IN THE CASE OF SIZE DISTRIBUTION #IDE: TAU=0.5, R=1M, RHO=900 #SIGMA=4*RHO*R/3*TAU=600 #SIZE: TAU=0.256, RMAX/RMIN=10, RMAX=10, RHO=900 #-> SIGMA=600 #INSTEAD OF BOXCOMP_F USE BOXCOMPS_F #IDENTICAL #1 dist: 1-> discrete sizes 2-> power-law #1 number of size groups #1.,2000, 1.,1. SIG,LKM,ifac,efac (q smin smax) #SIZE #2 dist: 1-> discrete sizes 2-> power-law #6 number of size groups #3.,.5, 5. q smin smax #------------------------------------------------ #MOONLET RUNS #------------------------------------------------ #parameters: # rho_moon moonlet density # rad_moon moonlet radius 0-> no moonlet # cen_moon 1-> moonlet gravity subtracted 0-> free # -1 --> moonlet does not feel the gravity of particles # (either use cen=0, rho zero, or cen=-1, rho>0 ?) # orb_cross, n_cross start tabulation at, number of tabulated crossing # i_boundary 0-> old 1-> fresh 2-> fresh/no overlaps # i_boundary 5 - correct treatment # n_boundary number of crossings used from incross #same program is used both for sampling of tabulated crossing #and for using the tabulated crossings: #------------------------------------------------ #RUNS FOR STORING PARTICLE CROSSINGS #set parameters: # orb_cross 10. start sampling at (example) # n_cross 100000 maximum number sampled (example) # i_boundary=0 use old periodic boundaries (FIXED) # rad_moon =0. do not use satellite (FIXED) # note: optical depth, elastic properties must be the same # as in the actual run. #------------------------------------------------ #RUNS FOR USING TABULATED PARTICLE CROSSINGS: #copy correct outcross-file to ./incross #set parameters: # orb_cross > ORBMAX sampling probably not wanted # n_boundary 10000 make sure there are enough orbits in incross-file! # i_boundary=5 correct moonlet boundary treatment # rad_moon = correct values #------------------------------------------------ # storing orbits: #------------------------------------------------ # two different ways of specifying particles # EXAMPLE 1) # 2, .0025d0 norbit orbtest2_1000g # takes norbit-particles from each size group (now only one group) # moonlet index=1 # orbtest2_1000g=stroing interval # EXAMPLE 2) # norbit negative -> # reads -norbit indexes form file 'orbit.index' # in the current example this is created by echo # stored orbits are read with: emorbit_read.pro #------------------------------------------------ # tabulated density and velocity fields #------------------------------------------------ #1.0,100,250,0.02,1. orb2d,n2dx,n2dy,dorb2d,forb2d #orb2d = output interval #n2dx,n2dy = number of cells #dorb2d = collection interval #forb2d = 1 - all data between outputs, (i-1)*orb2d - i*orb2d, is stored #forb2d = 0 - only last output, at i*orb2d, is stored # can have any value between these, (i-forb2d)*orb2d - i*orb2d # stored tables read with emtab_read.pro #------------------------------------------------ #cp ~hsalo/SUNBOX/RESULTS/nattest003.inipos ./inipos #on way of giving particle indices to be followed: via file (norbit neg) echo 1 > orbit.index echo 2 >> orbit.index rm alftab.dat echo 0.1 0.1 > alftab.dat echo 0.1 0.1 >> alftab.dat #copy correct outcross-file to ./incross #cp ~hsalo/SUNBOX/RESULTS/lnattest_5m_00/*.outcross ./incross #rho_moon=410 a=100 000 --> rhill/rmoon=1 for test particles #------------------------------------------------------------------ ~/SSBOX/CODE/SSBOX040405/ssbox_030206.exe << EOF > ide2007_e50_r300n.log 'ide2007_e50_r300n',99 RUNID, DUMMY 1.,1.,1. NEW:fh1,fh2 corio,epicyclic norm=1,1 300.,300.,300. ROO1,R002,ROO3 10.,10.,10.,10. troo1,troo2,troo3,troo3 1800., 0.,3.,1. rho_moon,rad(0_no),cen,rot 1000.,100000 orb_cross,n_cross 0,10000 i_boundary(0-old,5-fresh),n_boundary 1.,1,1,0. xyvel,zvel,z,spin 0., 0., 0., 0., 0., 0. em: cx,cy,eks,ink,delta,gamma 0., 4., 8. brange(0-old,1-y,2-pos,3-neg),bmin,bmax 2 isearch_order 1-x 2-y 40,0.,0.,1,1 NEW:nmodes,afidl,afidly,m1-m2 .1 dorb_fidl 0.001 qzinit/(r*omega) 0.0, 2., 0., 40, 1 NEW:tamp, vamp, eamp,namp,nwave 5000.,50. NEW:start_amp,dstart_amp 9999.002,999.04 NEW: dorbk1,dorbk2 (kcollect) 0 startup (0) or cont (1) -> RINIT 5, .1d0 norbit orbtest 0 nomix(0) mix(1) 8000,130.D6,-1 LKM ETA(M) SEED 1.000d0 ,1.d0 kappa/omega wz/omega 1 radial groups =lkmr 8000, -850., 850., 8., 12. NC(lkm) CX1 CX2 CI CE (lkmr lines) 152.508, 152.508 XL YL (XL=0 -> TAU=YL) 40, 5, -1 prowid sig4 def (>0 ->CYMAX) 1 dist: 1-> discrete sizes 2-> power-law 1 number of size groups 1.6667,8000.,1.,1. q smin smax) 9990000,100000 collmax colCOL 50.1, 0.01, .1 ORBMAX orbd, orbe (CIDL) .0025d0 DT 1., .5 , 25. raja1,raja1min,raja1max 1000., 0.5 DGIDL, PIDL 10000.,0.50, 20 orbpos,dorpos,norbpos (MIDL) 6 gmet 8-gravi 2-nongravi 0.5 glim 0. encz (0-all 1-only zcomp) 0., -1. encsoft ENCLIM(neg-rel to xl) 20.,1,64,64 ENCMIN (lcr/2),ngrav,NXG,NYG 8,3.00 nzg,fzrange 1 ncoll 0,0 mx,my 1, 0.50,0.000077,0.10, 0.0 ,0.0 IMODEL LAMDA V0 beta vcrita -.0,0 lfric (neg) fricmodel(0,3,3,3) 40, 400, 0 alfstep alfome (in OMEGAs) 9 dmet 0, 1,1,0,0 ihertz (0-no),rho,a,d,b2 (JM-params) -1.,-1. erovak1,erovak2 .0001 ctul_frac 10000 ntauphot 0,200.,0.05 animin,animax,dani (ANIDIL) 100.0,100,100,0.02,1. orb2d,n2dx,n2dy,dorb2d,forb2d (2DTAB) 000,.01 n_tracer,dorb_tracer -1.,100. orb-visco (-1-cdata),dorb_visco -1.,30,32,32 vglim1(-1-grav),vglim2,nx,ny 2. cic (+ngp, 2-cic) 1.,1 myyspinfac1,myyspinfac2 0.0 gamma_adh (Choksi 0.370) 0 shear_nzone 1.5,1.5,1.,1. sh1,sh2,wt1,wt2 0.,0.,1,2 debug:t1,t2,ii,jj 0., 0. gsheet, gcorr 0.05,2.,-10 dorb_eropart,safe,coll_eropart(neg-no) 1.,1000.,1. ovisco111,ovisco2,dvisc 0, 0 calc_egrav, calc_vg (1,1 -> calc) 0 nob_imp_rho600_size (0,1-no_xy,2-no_x,3-no_y,ig-i) 10000.1,0.2,40,40,40,80 pp: start,dorb,nrad,nx,ndx,ndy 9999.002,999.04 NEW: dorbn1,dorbn2 (ncollect) 20,9. NEW:neigh,raja 0.,30.,60.,1. RHO_AXI,m_AXI,m_axi2,met 0 ssbox_fft_versio 1-xyz no 0 is_addgrav (0-no, 1 -all) 0. tk_friction 0. slim_nograv (0-no effect) 40, 0.01, 1. EOF #------------------------------------------------------------------ #------------------------------------------------------------------ #4) convert fortran output-files to idl-savefiles ~hsalo/SUNBOX/boxread.sh ide2007_e50_r300n & #------------------------------------------------------------------ #------------------------------------------------------------------ #MISC INFO: #erovak1,erovak2 #dero=EROVAK1*sqrt(dvx**2+dvy**2+dvz**2)*dt+EROVAK2*raja1 #------------------------------------------------------------------ # IMODEL=1 ALFA=LAMDA # IMODEL=2 VC=V0 # IMODEL=3 ALFA=LAMDA/(1+C.V/V0) # BETA=0 NO FRICTION # BETA=1 VTAN SUPPRESSES #---------------------------------------------------------------------- # new models in ssbox_d3.exe: # v in cm/sec # IMODEL=-11 BRIDGES eps1=0.32*v^(-0.24) FROSTED T=210 K r=2.5 # IMODEL=-12 HATZES eps2=0.48*v^(-0.20) FROSTED T=123 K r=2.5 # IMODEL=-13 HATZES compacted frost r=2.5 cm # eps3=0.90*exp(-0.395*v) +0.002*v^(-0.9) # IMODEL=-14 HATZES compacted frost r=20 cm # eps4=0.90*exp(-0.22*v) +0.01*v^(-0.6) #---------------------------------------------------------------------- #------------------------------------------------------------------------ #CALCULATION OF IMPACTS (USE 4 with moonlets) # DMET=1 original dimpacts # DMET=2 dimpacts2: POTB faster' # DMET=3 dimpacts3: POTB faster + 2D' # DMET=4 dimpactsrk: RK4' # DMET=5 dimpactsrk: RK4 single' # DMET=6 dimpactsrk: RK4 single 2D' # DMET=7 dimpactsrk: RK4 linear or hertzian or tabulated # DMET=8 dimpactsrk: RK4 +friction #------------------------------------------------------------------------ #SELF-GRAVITY (USE 5): # GMET=0 GRAVIOLD # GMET=1 GRAVI # GMET=2 NONGRAVI # GMET=3 GRAVIGRID # GMET=4 GRAVIGRID (FFT2) # GMET=5 GRAVIGRID (FFT3) #--------------------------------------------------------------------- # DIST < ENCMIN -> calculate force directly # DIST > ENCMIN and DIST < ENCLIM -> calculate force via GRID # NGRAV: update forces every NGRAV steps # NXG,NYG,NZG: number of grid divisions # FZRANGE: z-extent of GRID = 2*fzrange*STDEV(ZZ) #------------------------------------------------------------------ # SIZE_DISTRIBUTION: # 2 1->DISCRETE 2-> CONT # 6 SIZE GROUPS # 3.,0.64,6.4 qmin smin smax) # IDENTICAL: # 1 1->DISCRETE 2-> CONT # 1 SIZE GROUPS # 1.0 100, 1., 1. sig,lkm,ifac,efac #----------------------------------------------------------------- unalias rm alias rmf rm -fr mv fort.99 debug.dat rmf *midl* rmf *.radi rmf *.enet rmf *.zidl rmf *.proj #rmf fort* rmf *junk* rmf *zprof #rmf *ctim rmf *tidl rmf *dgidl rmf *dggidl rmf *logi #rmf *.log #rmf *.2dtab #rmf *.alfa #rmf *.aniidl #rmf *.iter #rmf *.oidl #rmf *.tauphot #rmf *.ufidl