To couple or not to couple
Download
1 / 33

To Couple or Not To Couple - PowerPoint PPT Presentation


  • 246 Views
  • Uploaded on

To Couple or Not To Couple. John Warner, USGS. Outline. Overview of some recent advancements to ROMS sediment transport components wave/current interactions model coupling Focus on Model Coupling MCT, ROMS, SWAN How to develop a coupled ROMS-SWAN application (inlet test example)

loader
I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
capcha
Download Presentation

PowerPoint Slideshow about 'To Couple or Not To Couple' - akando


An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
To couple or not to couple

To CoupleorNot To Couple

John Warner, USGS


Outline
Outline

  • Overview of some recent advancements to ROMS

    • sediment transport components

    • wave/current interactions

    • model coupling

  • Focus on Model Coupling

    • MCT, ROMS, SWAN

  • How to develop a coupled ROMS-SWAN application (inlet test example)

  • Summary /where are we going?


What is model coupling
What is model coupling?

Model A running on M nodes.

Model B running on N nodes.

Model C ………

………

……

with communications

between all models.

(it also works here)

78 processor cluster -

WHOI shared


Why couple models
Why Couple models?

  • Your model doesn't have all the physics !

    • For example, waves are needed from a different source, such as another model or data

  • Provide dynamic feedback

    • nonlinear interactions between forcing

      (i.e. wave-currents-morphology).

http://www.cofc.edu


Where when is model coupling important
Where / when is model coupling important ?

Groundwater

Surface winds

Waves

Sea ice

http://www.cofc.edu/CGOInquiry/tidalinlets.htm


To couple or not to couple

How are we coupling models?

  • Model Coupling Toolkit

  • Mathematics and Computer Science Division Argonne National Laboratory

  • http://www-unix.mcs.anl.gov/mct/

  • R. Jacob, J. Larson, E. Ong, “M×N Communication and Parallel Interpolation in CCSM Using the Model Coupling Toolkit”, (Preprint) ANL/MCSP1225-0205, Mathematics and Computer Science Division, Argonne National Laboratory, Feb 2005. Submitted to International Journal for High Performance Computing Applications.

  • J. Larson, R. Jacob, E. Ong, “The Model Coupling Toolkit: A New Fortran90 Toolkit for Building Multiphysics Parallel Coupled Models”, (Preprint) ANL/MCS-P1208-1204, Mathematics and Computer Science Division, Argonne National Laboratory, Dec 2004. Submitted to International Journal for High Performance Computing Applications.

Model B

"N" processors

Model A

"M" processors

MCT is an open-source package that provides MPI based communications between all nodes of a distributed memory modeling component system. Download and compile as libraries that are linked to. (like netcdf)


Plan to change coupler to esmf in the near future
Plan to change coupler to ESMF in the near future

Earth System Modeling Frameworkhttp://www.esmf.ucar.edu/

Partners:NOAA Geophysical Fluid Dynamics Laboratory NOAA National Centers for Environmental Prediction

NSF National Center for Atmospheric Research NASA Goddard Global Modeling and Assimilation Office

NASA Goddard Institute for Space Studies NASA Jet Propulsion Laboratory NASA Goddard Land Information Systems project DOD Naval Research Laboratory

DOD Air Force Weather Agency DOD Army Engineer Research and Development Center DOE Los Alamos National Laboratory DOE Argonne National Laboratory

University of Michigan Princeton University Massachusetts Institute of Technology UCLA

Center for Ocean-Land-Atmosphere Studies Programme for Integrated Earth System Modeling (PRISM)

Common Component Architecture (CCA)


What do we need from swan and what does swan want from roms
What do we need from SWAN, and what does SWAN want from ROMS?

u, v, h, h

Dwave, Hwave, Lwave,

Pwave_top, Pwave_bot,

Ub_swan, Wave_dissip



3d momentum equations
3D momentum equations

(wave - current interaction adds new terms)

x:

y:

s:

continuity:

scalar

transport:


Radiation stress terms
Radiation stress terms

horizontal

radiation stresses:

vertical

radiation stresses:

where

roller function:

vertical

shape

functions:

wave energy:

stokes velocities:

Mellor, G. L. (2003). The three-dimensional current and surface wave equations. Journal of Physical Oceanography 33, 1978-1989.

Mellor, G. L. (2005). Some consequences of the three-dimensional current and surface wave equations. Journal of Physical Oceanography 35, 2291-2298.

Roller model based on:

Svendsen, I. A. (1984). Wave heights and set-up in a surf zone, Coastal Engineering, 8, 303-329.

Svendsen, I, Haas, K and Zhao, Q (2002). Quasi-3D Nearshore Circulation Model SHORECIRC, CACR Report 2002-01.


Depth integrated equations
Depth-integrated equations

x:

y:

continuity:

where

horizontal

depth-integrated

radiation stresses:

group speed:

stokes velocities:

Velocities are Lagrangian !! So we are recoding to make them all Eulerian.


Swan the wave model
SWAN (the wave model)

N = wave action density

(energy density / relative frequency)

  • cx, cy = propagation velocities (x- and y- directions)

  • s = relative frequency

  • = wave direction

    S = source/sink term for:

    - wind-wave generation

    - wave breaking

    - bottom dissipation

    - nonlinear wave-wave interactions

    SWAN accounts for shoaling, diffraction, partial transmission, and reflection.

Booij, N., R.C. Ris and L.H. Holthuijsen, 1999, A third-generation wave model for coastal regions, Part I, Model description and validation,

J.Geoph.Research, 104, C4, 7649-7666.

Booij, N., R.C. Ris and L.H. Holthuijsen, 1999, A third-generation wave model for coastal regions, Part II, Model description and validation,

J.Geoph.Research, 104, C4, 7649-7666.

Booij, N., Haagsma, IJ.G., Holthuijsen, L.H., Kieftenburg, A.T.M.M., Ris, R.C., van der Westhuysen, A.J., and Zijlema, M. (2004).

SWAN Cycle III version 40.41 User Manual, Delft University of Technology.

pretty good user’s manual !!!


Simple example of an uncoupled application to demonstrate wave driven flows
Simple example of an uncoupled application to demonstrate wave driven flows ….


Shoreface test case
Shoreface test case

test case setup

1) use SWAN to predict

wave heights, wave dissip, wave direction,

wave length, wave period

2) results from SWAN were used as forcing

to ROMS and Shorecirc.

Thanks to Kevin Haas at Ga Tech Savannah



Model organization

master.F

Model organization

mpi_init

init_file (# procs/model)

{

init

run

finalize

SWAN

{

init

run

finalize

ROMS


Init run and finalize
init, run, and finalize

ROMS

SWAN

init_param

init_parallel

init_scaclars

init_coupling

MPI_INIT

init

(grid decomp)

roms_init

SWINIT

SWREAD (grid)

init_coupling

SWINITMPI

run

(sync. point)

main3d

.....

waves_coupling

...

swanmain

.....

ocean_coupling

...

roms_run

SWMAIN

mpi_finalize

close_io

roms_finalize

finalize

SWEXITMPI

mpi_finalize

close_io


Grid decomposition during initialization
Grid decomposition (during initialization)

SWAN

ROMS

  • Each tile is on a

  • separate processor.

  • Each tile registers

  • with MCT.


Init coupling
init_coupling

ROMS- init_coupling

SWAN- init_coupling

1

1

2

2

3

3

processed by each ROMS tile

processed by each SWAN tile


Synchronization run phase
Synchronization (run phase)

ROMS- cocean_oupling

SWAN- waves_coupling

MCT

MCT

processed by each ROMS tile

processed by each SWAN tile


How to set up a coupled roms swan application
How to set up a Coupled ROMS / SWAN application

Use INLET_TEST as an example.

This application is provided in ROMS 2.3

You need to compile and install MCT (only once).

You need all the input files for a stand alone ROMS simulation.

You need all the input files for a stand alone SWAN simulation.

We have some m files to make ROMS grids, SWAN grids from ROMS grids,

SWAN wave forcing files, SWAN wind forcing files, etc.


Step 1 cppdefs h
Step 1: cppdefs.h

INLET_TEST

is already set up to

have model coupling

WAVES_OCEAN activates coupling to the wave model.

NEARSHORE_MELLOR activates radiation stress terms

in momentum equations.


Step 2 coupling input file
Step 2: coupling input file

coupling_inlet_test.in

Set number of processors

for each model.

Set synchronization

interval.

File names.

(SWAN = INPUT)


Step 3 roms input file
Step 3: ROMS input file

ocean_inlet_test.in

Make sure

number of tiles

= number of processors

specified in

coupling input file


Step 4 swan input file
Step 4: SWAN input file

INPUT

Keep these 2 sets of

commands for

CURRENT and

WLEV. These activate

internal SWAN flags.

MANY MORE commands further in this file. RTFSM.


Step 5 makefile
Step 5: makefile

Needs MPI

SWAN_COUPLE

activates compiler

directives for

MCT libraries.


Step 6 compile and run
Step 6: compile and run

make

mpirun –np M*N ./oceanM ROMS/External/coupling_inlet_test.in


Inlet test
Inlet Test

1200 m

2

ubar = 0.5 m/s

depth (m)

Hs = 2.0 m

T = 10 s

16

1200 m

Use the inlet test case as an example. This test case will be distributed with ROMS 2.3


Inlet test preliminary look
Inlet test (preliminary look)

ROMS

SWAN

ROMS

SWAN

1) Two-way coupled simulation shows wave-current interaction:

reduced magnitude of currents and increased wave height.

2) Coupled simulation has lower velocity and less erosion offshore of inlet.

3) Lateral BC effects of wave model.



New directory structure for roms 2 3
New directory structure for ROMS 2.3

(Natalie Perlin OSU, but not formalized into Rutgers version yet)

COAMPS

WRF

Others …


Summary
Summary

  • Model enhancements include :

    • sediment transport components

    • wave/current interaction

    • model coupling

  • Developed a coupled ocean-wave modeling system.

  • New directory structure to allow coupling to other models.

  • ROMS 2.3 release – in a few weeks after finalize Lagrangian to Eulerian issue.