Hi Forum,
Till now, I have been running models with moderate resolution, and these run reasonably fast on the medium-range cluster. A 200 Myr run takes around 1.5 weeks to complete using 160 cores using the gridding below.
subsection Mesh refinement
set Initial adaptive refinement = 1
set Initial global refinement = 4
set Minimum refinement level = 3
set Refinement fraction = 0.4
set Coarsening fraction = 0.05
set Strategy = nonadiabatic temperature,minimum refinement function, composition, vis
cosity
set Time steps between mesh refinement = 3
subsection Minimum refinement function
set Function constants =
set Function expression = if (depth > 2500000,6,if(depth > 660000,4,if(depth > 80000,5,6)))
set Variable names = depth,y,z
end
end
#################################################
Now, I am trying to move toward high-resolution models, and this is the one I got from a recently published paper: https://agupubs.onlinelibrary.wiley.com/doi/full/10.1029/2025JB032510. I was running this model with 960 cores on Stampede3, but it is very slow—I can only run about 7 Myr in 18 hours. Does anyone have any idea how many cores I should use to complete these runs within one week on Stampede3?
subsection Mesh refinement
set Initial adaptive refinement = 3
set Initial global refinement = 5
set Minimum refinement level = 3
set Refinement fraction = 0.35
set Coarsening fraction = 0.15
set Strategy = strain rate, maximum refinement function, minimum refinement function, temperature, viscosity, velocity
set Time steps between mesh refinement = 8
subsection Minimum refinement function
set Function constants =
set Function expression = if (depth < 400e6, 4, if (depth < 660e3, 4,3))
set Variable names = depth,y,z
end
subsection Maximum refinement function
set Function constants =
set Function expression = if (depth < 300e3, 8, if (depth < 450e3, 7, if (depth < 800e3, 7, 6)))set Variable names = depth,y,z
end
end