Hi,
I computed open boundary conditions for my child domain (using VVL) from a parent model (also using VVL). I ran a first simulation (simulation 1) using ln_full_vel = .true. and providing bn_u3d and bn_v3d.
Then, for a second simulation (simulation 2), I decomposed the full BDY currents into their barotropic and baroclinic components, using water level to separate these two parts (because of VVL), to test whether I could reproduce the results of simulation 1.
- Simulation 1: ln_full_vel = .true., with bn_u3d and bn_v3d provided.
- Simulation 2: ln_full_vel = .false., with bn_u3d, bn_v3d, bn_u2d, and bn_v2d provided, as NEMO would normally separate them.
So far, the results differ significantly between the two simulations.
Does anyone have suggestions or insights on this? My goal is to eventually understand the impact of using currents coming from different baroclinic and barotropic sources at the boundaries.