<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p style="margin-top:0;margin-bottom:0">Hi everyone,</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">I am experimenting with timesteps and have discovered something strange. I have deltaT=600, and for the first few months of simulation I need to have deltaTmom=60 so that some areas with an unstable water column in the
initial conditions don't cause the model to blow up. After the first few months, I can remove deltaTmom=60 (so it defaults back to deltaT=600) and the model is stable.</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">But when I do this, the model slows down! <span style="font-size: 12pt;">To confirm this was the case, </span><span style="font-size: 12pt;">I ran two simulations starting from the same timestep (1 year in). The simulation
with deltaTmom=60 completed 22464 timesteps in 20 minutes. The simulation with deltaTmom=deltaT=600 only completed 8496 timesteps. </span></p>
<p style="margin-top:0;margin-bottom:0"><span style="font-size: 12pt;"><br>
</span></p>
<p style="margin-top:0;margin-bottom:0"><span style="font-size: 12pt;">I have two questions about this:</span></p>
<p style="margin-top:0;margin-bottom:0"><span style="font-size: 12pt;"><br>
</span></p>
<p style="margin-top:0;margin-bottom:0">1) Is this behaviour as counterintuitive as it seems? I'm wondering if larger momentum timesteps are somehow causing the elliptic solver to converge more slowly, which more than compensates for the reduced number of timesteps?</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">2) Since a smaller deltaTmom is so much faster, is there anything wrong with having a different value for deltaTmom and deltaT? Some of my colleagues thought that having different timesteps for ocean tracers and momentum
might be dangerously unphysical.</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Many thanks,</p>
<p style="margin-top:0;margin-bottom:0">Kaitlin Naughten</p>
</div>
<hr>
<small>This message (and any attachments) is for the recipient only. NERC is subject to the Freedom of Information Act 2000 and the contents of this email and any reply you make may be disclosed by NERC unless it is exempt from release under the Act. Any material
supplied to NERC may be stored in an electronic records management system</small>.
<hr>
</body>
</html>