#------------------------------------------------------------------ #tbgrav060tau140.sjob 1.04.2003 #------------------------------------------------------------------ #self-gravitating runs for the TB #using code transferred from halla 080403 #individual eroma #gravitational viscosity #tau=0.5 #4*4lcr #epsn=0.5 #epst=1.0 #------------------------------------------------------------------ #------------------------------------------------------------------ #1) create directory for results rm -r ~hsalo/SUNBOX/RESULTS/tbgrav060tau140 mkdir ~hsalo/SUNBOX/RESULTS/tbgrav060tau140 #------------------------------------------------------------------ #------------------------------------------------------------------ #2) go to this directory + copy job-file there cd ~hsalo/SUNBOX/RESULTS/tbgrav060tau140 cp ~hsalo/SUNBOX/JOBS/tbgrav060tau140.sjob ./ #------------------------------------------------------------------ #------------------------------------------------------------------ #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 #individual eromax #dorb_eropart,safe,coll update-interval, safe-factor, # minimum ncoll for using eromax #------------------------------------------------ # 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-> fixed 0-> free -1-> free+gravity # 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 ../emn_t000/emn_t000.outcross ./incross #on way of giving particle indices to be followed: via file (norbit neg) echo 1 > orbit.index echo 2 >> orbit.index #spin-versio #control-file #eromax # #------------------------------------------------------------------ ~hsalo/HALLA_CODE/HALLA150403/halla150403.exe << EOF > tbgrav060tau140.log 'tbgrav060tau140',99 RUNID, DUMMY 900.,900. ROO1,R002 900., 0.,-1.,1. rho_moon,rad(0_no),cen,rot 1000.,100000 orb_cross,n_cross 0,0 i_boundary(0-old,5-fresh),n_boundary 0.,0.,0.,0. xyvel,zvel,z,spin 1 isearch_order 1-x 2-y 5,0. nmodes,afidl .02 dorb_fidl 0.1 qzinit/(r*omega) 0 startup (0) or cont (1) -> RINIT 5, .1d0 norbit orbtest 0 nomix(0) mix(1) 4542,60.D6,-1 LKM ETA(M) SEED 1.000d0 ,1.d0 kappa/omega wz/omega 1 radial groups =lkmr 4542, -1977.2, 1977.2, 6.,9. NC(lkm) CX1 CX2 CI CE (lkmr lines) 0.,1.4 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,4542, 1.,1. SIG,LKM,if,ef (q smin smax) <- RINIT 9990000,1000 collmax colCOL 30.1, 0.01, .1 ORBMAX orbd, orbe (CIDL) .002d0 DT 3., 2. , 25. raja1,raja1min,raja1max 1000., 1. DGIDL, PIDL 0.,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,128,128 ENCMIN (lcr/2),ngrav,NXG,NYG 8,3.00 nzg,fzrange 1 ncoll 0,0 mx,my 1 0.50,0.000077, 0.0 ,0.0 IMODEL LAMDA V0 beta vcrita -.0,0 lfric (neg) fricmodel(0,3,3,3) 20, 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 .001 ctul_frac 1000 ntauphot 0,20.,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,128,128 vglim1(-1-grav),vglim2,nx,ny 2. cic (+ngp, 2-cic) 1.,1 myyspinfac1,myyspinfac2 0.00000 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 1., 1. gsheet, gcorr 0.1,2.,10 dorb_eropart,safe,coll EOF #------------------------------------------------------------------ #------------------------------------------------------------------ #4) convert fortran output-files to idl-savefiles ~hsalo/SUNBOX/boxread.sh tbgrav060tau140 & #------------------------------------------------------------------ #------------------------------------------------------------------ #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