paramsf2p.dat

#############################
### Simulation parameters ###
#############################

directory	        out1		# directory's name
physical_units	    Yes			# convert field in physical units? (Yes or No)
override_units	    No			# do we change code units? If so, the two lines below must be specified
new_unit_length	    4.5e12	    # new unit of length in meters
new_unit_mass	    4e30		# new unit of mass in kg

movie		        No			# make a movie? requires on to take the form onmin,onmax (eg: 1,10)
on	   	            20		    # fargo's output number(s), separated by a comma (,) if several
                                # can also be 'all' to display all available frames

verbose             Yes         # displays more information on screen


###############
### Roadmap ###
###############

### plot time evolution of disc torque or power on planets?
# type No, torque, indtorque or power
plot_tqwk	           No
running_time_average   No
normalize_torque       Yes		# normalize by reference Gamma_0 torque?	

### plot planets orbital elements?
# type x,y to plot y vs. x: x can be t (time) or a (semi-major axis);
# y can be a, e (eccentricity), r (orbital radius), m (planet mass),
# p (orbital period ratio, if several planets) or mmr (mean-motion resonant angles)
plot_planet  	No			 

# if plot_planet is set to Yes and y to mmr,
# type fo,fi to compute the fo:fi mean-motion resonant angles (eg, 3,2 or 2,1)
mmr_integers	5,3

### plot dust quantities for Lagrangian particles (Dusty FARGO-ADSG only)?
# type x,y to plot y vs. x: x can be r (radius) or phi (azimuth);
# y can be phi, size (dust size), Stokes (Stokes number), vr (radial velocity) or vphi (azimuthal velocity)
plot_dust	    No


### if all above plot booleans are set to No, gas quantities are displayed:

fluid		    gas			# what fluid? gas, dust or all (FARGO3D only)
whatfield	    dens		# full list of available fields to be displayed at the bottom of this file
onedprofile	    No			# default is to plot (1D) azimuthally-averaged radial profile. Other options: Cut, Median
nodiff	        Yes			# relative difference wrt initial field? (Yes, No, nao or normnao)
                            # nao = X - <X>_phi, normnao = X / <X>_phi - 1

fieldmin        		    # if specified, minimum value of field to display (for color bar)
fieldmax				    # if specified, maximum value of field to display (for color bar)

fieldofview	    polar		# field of view: cartesian (cart), polar (polar, default: phi vs. r), vertical or latitudinal?
rvsphi		    No			# if field of view is polar, do we plot radius vs phi instead?
slice           midplane    # can be: 'midplane', 'upper' or 'lower' if field of view is polar or cart
                            # or can be 'average' or simply left undefined if field of view is vertical or latitudinal
z_average       No          # 3D cylindrical runs w/o vertical stratification: vertical average of disc fields before plot?

rotate_angle			    # if set, rotate 2D field by angle specified in degrees
flip_xaxis	    No	        # if Yes, flip x axis to -x axis

showdust	    No			# show Lagrangian dust particles? (Dusty FARGO-ADSG only)
sizemin		    1e-5		# if specified, minimum particle size (in meters)
sizemax		    1e-2		# if specified, maximum particle size (in meters)

showplanet	    No			# spot planet position with a filled circle

streamlines	    No			# show gas streamlines?
nstreamlines	20		    # number of streamlines to show


######################
### Common options ###
######################

log_colorscale    Yes	    # logarithmic color scale for 2D plots?
log_xyplots_y	  No		# logarithmic y-scale for y=f(x) plots?
log_xyplots_x	  No		# logarithmic x-scale for y=f(x) plots?

use_legend				    # manually enter legend to be used for 1D plots if several directories 
                            # ex: 'First_run','Second_run','Third_run' (the '_' will be replaced by white spaces)

myrmin					    # if specified, minimum radius to display field
myrmax	        			# if specified, maximum radius to display field
rbox                        # if specified, radius range about planet's orbital radius to display
                            # specify 'planet_to_track' below if there are several planets
planet_to_track             # index of planet to display squared box around
box_around_planet           # size of squared box around tracked planet

myphimin				    # if specified, minimum azimuth to display field (between 0 and 2pi)
myphimax				    # if specified, maximum azimuth to display field (between 0 and 2pi)
myzmin                      # if specified, minimum z to display field
myzmax                      # if specified, maximum z to display field

mytmin					    # if specified, minimum time for y=f(time) plots
mytmax					    # if specified, maximum time for y=f(time) plots
myymin					    # if specified, minimum y-value for y=f(time or x) plots
myymax	        			# if specified, maximum y-value for y=f(time or x) plots

time_in_inner_planet_orbits No   # by default time is expressed in code's units except if  
                                 # it is specified to be in orbital periods at innermost planet

take_one_point_every        # if specified, take one point every N in y=f(x) plots

saveaspdf	     No			# save image as pdf file? (Yes or No)
saveaspng	     Yes	    # save image as png file? (Yes or No)
mycolormap	     nipy_spectral 	    # color scale for 2D contour plots

display_label               # add user-defined label in top-right corner of 2D plots

dynamical_colorscale        # if an animation is requested, colorscale varies from one frame to another
                            # instead of being set throughout all frames

################################
### Other options in roadmap ###
################################

### plot disc's centre-of-mass?
# type 'xy' to plot x vs. y (time in colour) or 'tr' to show time evolution of radial distance 
# of centre-of-mass from frame's origin (ie the central star), or 'tpa' to show time evolution 
# of the position angle of the disc's centre of mass
plot_disccom	No
central_binary  No          

### plot time evolution of disc mass?
plot_discmass	No

### plot time evolution of disc's averaged eccentricity?
plot_discecc	No

### plot time evolution of disc's averaged pericentre argument?
plot_discperarg	No

### plot time evolution of librating, orbit-crossing flow vortensities?
plot_libcross   No

### plot gas quantities for turbulent disc simulations (eg MRI-driven MHD turbulence 
# with FARGO3D) -- plot_turb can be set to: No, power_spectrum, auto_correlation
plot_turb       No


################################
### List of available fields ###
################################

# dens: surface density (2D runs) or volume density (3D runs)
# surfacedens: vertically-integrated volume density (3D runs)
# rtadens: running-time averaged particle density (2D runs with Lagrangian particles)
# vrad: radial velocity
# vtheta: azimuthal velocity
# vcol: latitudinal velocity
# brad: radial component of magnetic field
# btheta: azimuthal component of magnetic field
# bcol: latitudinal component of magnetic field
# pressure: thermal pressure
# temp: temperature
# toomre: Toomre-Q parameter (2D runs)
# entropy: specific entropy
# vorticity: vertical component of velocity curl (2D runs only so far!)
# vortensity: vorticity divided by gas density (2D runs only so far!)
# invvortensity: inverse vortensity
# normvorticity: ratio of horizontal epicyclic frequency squared over angular frequency squared
# drl: radial derivative of specific angular momentum
# ecc: disc eccentricity assuming each grid cell interacts with the star as a two-body problem (2D)
# stokes: dust Stokes number (2D like 3D)
# label: passive scalar (original 2D)
# mdot: mass accretion rate (2D)
# dgratio: dust-to-gas density ratio
# sgaccr: radial self-gravitating acceleration
# sgacctheta: azimuthal self-gravitating acceleration
# bettacooling: cooling-timescale-to-orbital-period ratio
# direct_torque: direct disc torque on planet
# cumulative_direct_torque: cumulative direct disc torque on planet (int_rmin^r Gamma(x)dx)
# indirect_torque: indirect disc torque on planet
# cumulative_indirect_torque: cumulative indirect disc torque on planet (int_rmin^r Gamma(x)dx)
# alpha_reynolds: vertically and azimuthally averaged radial profile of Reynolds alpha parameter
# alpha_maxwell: vertically and azimuthally averaged radial profile of Maxwell alpha parameter
# bc: beta cooling if computed by FARGO3D
# betacooling: computed via density and temperature outputs, assuming BL94 opacities
# kappa: Rossland opacities if computed by FARGO3D
# rossby: Rossby number (~normalized relative vorticity)