## 'nfast' and 'ndtfast' on Barotropic time step

Report or discuss software problems and other woes

Moderators: arango, robertson

Message
Author
yj7054
Posts: 9
Joined: Mon Sep 24, 2018 7:40 pm
Location: CSIRO - Hobart Site

### 'nfast' and 'ndtfast' on Barotropic time step

Dear all,

nfast: Number of barotropic timesteps needed to compute time-averaged barotropic variables centered at time level n+1
ndtfast: Number of barotropic timesteps between each baroclinic timestep.

I am confused with the time-stepping in ROMS.

For example, I set dt = 150s, ndtfast = 30 in ocean.in, so I got dtfast = 5s, and I also got nfast = 42 computed by set_weights.F. In main3d.F, there are 42 barotropic time-steps for each dt, rather than 30 steps. SO my question is that are there 12 time steps(M*-M) double counted for each dt? For the next baroclinic step n+2, the barotropic step starts from m=M or from m=M* on the figure attached?

Could anyone give me a further explain? Thank you very much!
Attachments
The split time stepping used in the model.
The split time stepping used in the model.png (15.72 KiB) Viewed 1037 times

wilkin
Posts: 511
Joined: Mon Apr 28, 2003 5:44 pm
Location: Rutgers University
Contact:

### Re: 'nfast' and 'ndtfast' on Barotropic time step

https://www.myroms.org/wiki/Numerical_S ... g_Overview
and the Shchepetkin papers.
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu

yj7054
Posts: 9
Joined: Mon Sep 24, 2018 7:40 pm
Location: CSIRO - Hobart Site

### Re: 'nfast' and 'ndtfast' on Barotropic time step

wilkin wrote:Read this explanation on WikiROMS:
https://www.myroms.org/wiki/Numerical_S ... g_Overview
and the Shchepetkin papers.
Hi John,

I read the explanation before but it was complex for me. Now I think I got the main idea roughly:

There are M* barotropic steps for each baroclinic step, and the fast-time averages among these M* steps are transmitted to the baroclinic mode, as well as used as the initial values for barotropic mode during the next baroclinic step.

arango
Posts: 1091
Joined: Wed Feb 26, 2003 4:41 pm
Location: IMCS, Rutgers University
Contact:

### Re: 'nfast' and 'ndtfast' on Barotropic time step

Nope, there are actually M* barotropic timesteps The reason for it is that the cosine-square shape filter needs to provide time-averaged values of ubar, vbar, and zeta centered at the n+1 baroclinic timestep, which gives us second-order temporal accuracy for time-averaged barotropic motions. You may check routine set_weights.F to see how the time weights are computed. The above diagram is not clear enough.

yj7054
Posts: 9
Joined: Mon Sep 24, 2018 7:40 pm
Location: CSIRO - Hobart Site

### Re: 'nfast' and 'ndtfast' on Barotropic time step

arango wrote:Nope, there are actually M* barotropic timesteps The reason for it is that the cosine-square shape filter needs to provide time-averaged values of ubar, vbar, and zeta centered at the n+1 baroclinic timestep, which gives us second-order temporal accuracy for time-averaged barotropic motions.
Hi Hernan,

Thank you very much!

So the time-averaged values of ubar, vbar, and zeta you mentioned above will be used as initial value for the barotropic modes from m=M to m=M+M* during next baroclinic step n+2, right?

If I am wrong, where does the initial value come from during the next baroclinic step n+2?

arango