jmmp-group / severn-swot

Severn estuary 500m ocean model

NEMO regional configuration of the Severn Estuary

This model configuration has been developed in order to test the NEMO model in a tidal coastal area at 500 m resolution (in preparation for the UK wide NEMO 500 m model development).

Model Summary

Severn Estuary in the UK (-5E to -2E and 50.1N to 51.8N)

The model grid target resolution 500m; with 31 sigma vertical levels. Featuring:

• FES2014 tides
• ERA5 atmospheric forcing
• CMEMS open boundaries & Initial conditions
• Freshwater forcing
• Wetting and Drying

Model Setup

The following process is followed to build and get started with this configuration

git clone https://github.com/JMMP-Group/SEVERN-SWOT.git

Then follow descriptions in: https://github.com/JMMP-Group/SEVERN-SWOT/wiki

Experiment Summary

In the master branch:

• EXP_unforced_constTS Constant initial T and S. No tides. No met.

• EXP_barotropicTide Constant initial T and S. Only tidal forcing. No met.

• EXP_Tide_ERA5 Constant initial T and S. Tidal forcing. ERA5 forcing (sea-level pressure, 10m winds, fluxes etc.).

• EXP_Tide_ERA5_BDY Tidal forcing. ERA5 forcing (sea-level pressure, 10m winds, fluxes etc.), 3D boundary conditions and initial conditions.

In the WAD branch:

• EXP_barotropicTide_WAD Constant initial T and S. Only tidal forcing. No met. With Wetting And Drying.

• EXP_Tide_ERA5_WAD Constant initial T and S. Tidal forcing. ERA5 forcing (sea-level pressure, 10m winds, fluxes etc.). With Wetting And Drying.

• EXP_Tide_ERA5_River_IC_BDY_WAD Tidal forcing. ERA5 forcing (sea-level pressure, 10m winds, fluxes etc.), 3D boundary conditions and initial conditions. River forcing. With Wetting And Drying.

severn-swot's Issues

Add open boundaries (SSH + 2d velocities)

Tidal forcing is harmonic whereas SSH and 2D velocity forcing is time varying. All boundary files can be made with the PyNEMO tool.

This step does not require temperature and salinity structure in the simulation so is a natural progression in forcing from tide only, though will be more interesting when atm forcing is also added.

Generate tidal boundaries

Generate tidal boundary forcing as harmonic files using FES2014 data and the PyNEMO tool.

1. Compile PyNEMO (ideally on Archer2, otherwise livljobs)
2. Set up namelists
3. get necessary input files together
4. run and generate netcdf files of boundary harmonics

Acceptance criteria:

• Description in the wiki
• Any new scripts in SCRIPTS

Tools/DOMAINcfg/src files don't download

For some reason the src code for DOMAINcfg doesn't download but the src code for other tools does...
lines 10-16 in https://github.com/JMMP-Group/SEVERN-SWOT/blob/master/SCRIPTS/make_tools.sh

Documentation to build NEMO and XIOS execs

A walk through on how to build the NEMO and XIOS executables for Severn domain

Build domain_cfg file

Generate static files that define the domain. Work flow documentation is in wiki and SCRIPTS folder

Port file prep to ARCHER2

Some pre processing happens on livljobs8 because of nco command availability.
It would be easy to rapidly cycle through the build process if everything was on ARCHER2.

Action:

• Migrate wiki notes on preparing bathymetry and other things(?) to work on ARCHER2

Run tide only

Run the Severn domain with tides only forcing.

Instructions:

• Wiki recipe
• SCRIPTS

Run tide + met

Experiment set up with example met forcing (ERA5: pressure and winds) + tides

tide + met + WAD

WAD makes the waterline move back and forth. It also allows to relax the hardwired condition for minimum depth (currently 10m).

1st create new branch: feature/metWAD

Steps:

1. switch on WAD. This is the namelist namwad ! Wetting and Drying (WaD) in the namelist_cfg file
2. test run with WAD switched on
3. regenerate bathymetry to relax the hardware condition of 10m minimum depth
4. run

Add wetting and drying

Wetting and drying permits the water line to move. It also relaxes the requirement for a "minimum depth" being hardwired into the bathymetry (typically about 10m)

Run unforced

Does it timestep without forcing? Not a bad test to make sure ARCHER2 accounts are set up correctly.

Need coordinates.bdy.nc using PyNEMO. Then this runs https://github.com/JMMP-Group/SEVERN-SWOT/blob/master/SCRIPTS/run_unforced.sh

Run tide + met + wave forcing

Waves can be added as an external forcing (which is probably as sensible step prior to coupling)

Run tide + barotropic forcing

Create EXPeriment with barotropic forcing (SSH and depth averaged currents) + tidal forcing

NEMO-OASIS-WW coupling

Replicate NEMO wave coupling from EP model

Model hangs

Model seems to run but hangs without terminating properly

Check/update NEMO tools to r4.0-HEAD

I needed domzgr.f90.melange to get the vertical coord working. This might actually already be in r4.0-HEAD (James said). Probably need to review how the tools where checked out. (Though it seems to work for now with old namelist_cfg files)

tide boundary generation

From @micdom: @jpolton @mpayopayo I'm back to this:
what are the files you're listing in namelist*.bdy for sn_src_hgr, sn_src_zgr, sn_src_msk?
I don't have access to the /projectsa/accord folder (and probably I shouldn't to run a different configuration).

I think those files contain the domain from which the tides are picked up to be generated into our domain, maybe @jpolton can confim? Maybe there could be a copy of those files in projectsa/SWOT and pick them from there?

Run tide + met + waves + rivers