# Cookbook 6: Pseudo-Free-Surface Formulation [CitcomS] steps = 100 [CitcomS.controller] monitoringFrequency = 10 [CitcomS.solver] datafile = cookbook6 rayleigh = 4.312616e+08 [CitcomS.solver.mesher] nprocx = 2 nprocy = 2 nprocz = 1 nodex = 61 nodey = 61 nodez = 25 coor = 1 coor_file = ./coord.dat [CitcomS.solver.tsolver] fixed_timestep = 7.77e-10 [CitcomS.solver.vsolver] precond = on accuracy = 1e-5 vlowstep = 100000 piterations = 100000 [CitcomS.solver.ic] # To try a different hot blob, one can specify three parameters: # 1. blob_center: Coordinates(theta,fi,r) of the center. # 2. blob_radius: Radius as non-dimensional distance. # 3. blob_dT: The amount of temperature increase/decrease # within the blob as non-dimensional temperature. # # Note: It is required to keep "tic_method = 2" to have the initial # hot blob. tic_method = 2 blob_center = 1.570800e+00,1.570800e+00,9.246600e-01 blob_radius = 6.278334e-02 blob_dT = 0.18 half_space_age = 2500.0 mantle_temp = 0.82 [CitcomS.solver.bc] topvbc = 2 pseudo_free_surf = on bottbcval = 0.82 # To convert pseudo-free-surface B.C. to free-slip B.C., do the # following: #topvbc = 0 #pseudo_free_surf = off [CitcomS.solver.param] start_age = 60 [CitcomS.solver.visc] TDEPV = on visc0 = 1,1,1,1 viscE = 9.50614,9.50614,9.50614,9.50614 viscT = 1.02126,1.02126,1.02126,1.02126 [CitcomS.solver.const] radius = 6.371e+06 density = 3270.0 thermdiff = 1.0e-06 gravacc = 10.0 thermexp = 3.0e-05 refvisc = 1.0e+21