<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hi all,<div><br></div><div>I've been running some 2D and 3D spherical models and just noticed what's probably a small bug somewhere in the post-processing section of the code.</div><div><br></div><div>In the example input file simple_2d_shell.prm, the radial gravity value is specified as follows:</div><div><font class="Apple-style-span" color="#da180f"> subsection Radial constant</font></div><div><font class="Apple-style-span" color="#da180f"> # Magnitude of the gravity vector in $m/s^2$. The direction is always</font></div><div><font class="Apple-style-span" color="#da180f"> # radially outward from the center of the earth.</font></div><div><font class="Apple-style-span" color="#da180f"> set Magnitude = 9.81</font></div><div><font class="Apple-style-span" color="#da180f"> end </font></div><div><br></div><div>In the .o file output during the run, this is not the listed value for gravity:</div><div><div><font class="Apple-style-span" color="#da180f"> Postprocessing:</font></div><div><font class="Apple-style-span" color="#da180f"><br></font></div><div><font class="Apple-style-span" color="#da180f"> Reference density (kg/m^3): 3300 </font></div><div><font class="Apple-style-span" color="#da180f"> Reference gravity (m/s^2): 7.714e+13</font></div><div><font class="Apple-style-span" color="#da180f"> Reference thermal expansion (1/K): 4e-05 </font></div><div><font class="Apple-style-span" color="#da180f"> Temperature contrast accross model domain (K): 3300</font></div><div><font class="Apple-style-span" color="#da180f"> Model domain depth (m): 2.855e+06</font></div><div><font class="Apple-style-span" color="#da180f"> Reference thermal diffusivity (m^2/s): 1.13939e-06</font></div><div><font class="Apple-style-span" color="#da180f"> Reference viscosity (Pas): 1e+22 </font></div><div><font class="Apple-style-span" color="#da180f"> Ra number: 6.86297e+19</font></div><div><font class="Apple-style-span" color="#da180f"> k_value: 4.7 </font></div><div><font class="Apple-style-span" color="#da180f"> reference_cp: 1250 </font></div><div><font class="Apple-style-span" color="#da180f"> reference_thermal_diffusivity: 1.13939e-06</font></div></div><div><br></div><div>So it seems the correct gravity value is not being accessed in the post-processing section of the code. The model would not run with Ra = 1e19, so it almost certainly seems to be restricted to post-processing. </div><div><br></div><div>The calls to access the gravity value in velocity_statistics.cc and table_velocity_statistics.cc </div><div>seem to work fine for the box models, so I'm not sure where the source of the error lies.</div><div><br></div><div>Let me know if I can pass on any more info.</div><div> <br><div apple-content-edited="true">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>Cheers,</div><div>John</div><div><br></div><div>John Naliboff</div><div>Postdoctoral Scholar</div><div>Geology Department</div><div>University of California, Davis</div><div>Earth and Physical Sciences Bldg., Rm. 2119</div><div>Davis, CA 95616</div><br class="Apple-interchange-newline"><br class="Apple-interchange-newline"></div>
</div>
<br></div></body></html>