Sediment setup
Peter Holterman
I'm currently helping to setup a GETM model of a part of the West coast
of Sri Lanka. One of the objectives is to study the sediment
distribution of a river discharging into the Indian Ocean. Does have
GETM sediment implemented? Is there an example test case I can have a
look at?
Any help is appreciated!
Best regards,
Peter Holtermann
Manolo Ruiz
Dear Peter,
GETM has indeed a sediment set-up that has been applied to different
set-ups:
- an spm module was integrated into the GETM developers version and was
used to study the Elbe Estuarine Turbidity Maximum. Two configurations
were set-up and tested: an idealised 2d estuary and a realistic 3d Elbe
application (Burchard et al. 2004)
- to the Sylt-Romo Bight, Burchard et al. JPO 2008
- to some Wadden Sea set-ups by Emil Stanev and other people in Oldenburg.
I think the Sylt-Romo test case must be up and running in the last GETM
version.
These test-cases are mainly estuarine applications to domains with
extensive intertidal flats and if your domain is a shelf configuration,
there might be other issues to concentrate in: SPM fluxes through the
open boundary, sediment classes distribution on the shelf...
hope it helps,
Manuel
References
Burchard, H., K. Bolding, and M. R. Villarreal, Three-dimensional
modelling of estuarine turbidity maxima in a tidal estuary, Ocean
Dynamics, 54, 250-265, 2004.
Burchard, H., G. Flöser, J. V. Staneva, T. H. Badewien, and R.
Riethmüller, Impact of density gradients on net sediment transport into
the Wadden Sea, J. Phys. Oceanogr., 38, 566-587, 2008.
Stanev EV, Brink-Spalink G, Wolff J-O (2007) Sediment dynamics in
tidally dominated environments controlled by transport and turbulence. A
case study for the East Frisian Wadden Sea. J Geophys Res 112, C04018,
doi:10.1029/2005JC003045
Peter Holterman escribió:
==================================
Dr. Manuel Ruiz Villarreal
Instituto Español de Oceanografía
Centro Oceanográfico de A Coruña
Muelle de Animas s/n
15001 A Coruña, Galicia, SPAIN
ALTERNATIVE E-MAIL: mrui...@gmail.com
(use if delivery problems and cc in urgent or important mails)
Tel.: (+34) 981 20 53 62 Ext. 23252
Direct: (+34) 981 21 82 52
Fax: (+34) 981 22 90 77
http://www.ieo.es
==================================
La información contenida en este e-mail y sus ficheros adjuntos es totalmente confidencial y no debería ser usado si no fuera usted alguno de los destinatarios. Si ha recibido este e-mail por error, por favor avise al remitente y bórrelo de su buzón o de cualquier otro medio de almacenamiento. This email is confidential and should not be used by anyone who is not the original intended recipient. If you have received this e-mail in error please inform the sender and delete it from your mailbox or any other storage mechanism.
Richard Hofmeister
It is time for a "Howto SPM"...
The implemented spm module has a seperate namelist and can be activated
with the compiler option GETM_SPM=true
Try to take the namelist "spm.inp" from the interleaving test case.
The code is located in $GETMDIR/src/3d/spm.F90 - you might take a look
into that, because there not very much documentation beside that. The
spm method has two state variables: the spm_pool in a virtual bed layer
(2D variable) and the suspended matter itself (3d variable in the model
grid). You can set sinking velocities and critical stresses for the
exchange between the potential spm bed pool and the spm in the water.
At the moment, there is a bug in the developers code in the write and
read routines for netCDF restart files. The spm_pool variable name is
incorrect and the spm variables are not imported, i can send a bug fix
to Karsten.
Richard
Peter Holterman schrieb:
Richard Hofmeister
Leibniz Institute for Baltic Sea Research
Seestr. 15, 18119 Warnemünde
Tel: 0381-5197103 Fax: 0381-5197440
Richard Hofmeister
The compiler option is "-DSPM" for the code itself, but you have to set
an environment variable called "GETM_SPM=true" for the Makefile.
Richard
Richard Hofmeister schrieb:
Karsten Bolding
On Mon, Jul 13, 2009 at 11:03:48 +0200, Richard Hofmeister wrote:
>
> Hello Peter,
>
> It is time for a "Howto SPM"...
Sounds like a good idea.
> The implemented spm module has a seperate namelist and can be activated
> with the compiler option GETM_SPM=true
> Try to take the namelist "spm.inp" from the interleaving test case.
The code has not been changed for a long time so I assume it will work
OK - also with the spm.inp you point to.
>
> The code is located in $GETMDIR/src/3d/spm.F90 - you might take a look
> into that, because there not very much documentation beside that. The
> spm method has two state variables: the spm_pool in a virtual bed layer
> (2D variable) and the suspended matter itself (3d variable in the model
> grid). You can set sinking velocities and critical stresses for the
> exchange between the potential spm bed pool and the spm in the water.
Johan van der Molen has implemented a much more advanced sediment module
in GOTM (not in CVS yet). When available in GOTM it should be included
in GETM as well. It allows for up to 10 (correct me if I'm wrong Johan)
sediment classes with different characterstics.
>
> At the moment, there is a bug in the developers code in the write and
> read routines for netCDF restart files. The spm_pool variable name is
> incorrect and the spm variables are not imported, i can send a bug fix
> to Karsten.
I'll update CVS according to Richards patch.
>
> Richard
>
Karsten
Johan Van Der Molen (Cefas)
<snip>
> >
> > The code is located in $GETMDIR/src/3d/spm.F90 - you might take a
> > look into that, because there not very much documentation
> beside that.
> > The spm method has two state variables: the spm_pool in a
> virtual bed
> > layer (2D variable) and the suspended matter itself (3d variable in
> > the model grid). You can set sinking velocities and
> critical stresses
> > for the exchange between the potential spm bed pool and the
> spm in the water.
>
> Johan van der Molen has implemented a much more advanced
> sediment module in GOTM (not in CVS yet). When available in
> GOTM it should be included in GETM as well. It allows for up
> to 10 (correct me if I'm wrong Johan) sediment classes with
> different characterstics.
>
Not just 10: as many as you want (within the memory/computing time
limits of your system :-) )!
<snip>
Johan
***********************************************************************************
This email and any attachments are intended for the named recipient only. Its unauthorised use, distribution, disclosure, storage or copying is not permitted. If you have received it in error, please destroy all copies and notify the sender. In messages of a non-business nature, the views and opinions expressed are the author's own and do not necessarily reflect those of the organisation from which it is sent. All emails may be subject to monitoring.
***********************************************************************************
Peter Holterman
this would be the perfect feature for our setup, we want to look at the
sedimentation patterns of different sediment types from a river source.
-Peter
Johan Van Der Molen (Cefas)
> -----Original Message-----
> From: getm-...@googlegroups.com
> [mailto:getm-...@googlegroups.com] On Behalf Of Peter Holterman
> Sent: 20 July 2009 10:54
> To: getm-...@googlegroups.com
> Subject: [getm-users: 313] Re: Sediment setup
>
>
> Hi,
>
> this would be the perfect feature for our setup, we want to
> look at the sedimentation patterns of different sediment
> types from a river source.
>
-clean up and release 1D GOTM version (there are some variables that
need to be re-defined/re-directed)
-implement some kind of wave model/parameterisation to run alongside
GETM; wave data are required input for the SPM model (for the 1D version
I used observations)
-implement, test and release 3D GETM version of the SPM model
Realistically, I don't think I'm going to have much time for this in the
next 6 months or so. But if someone from the getm developers group would
like to have a go at it I would be happy to help him/her along.
For anyone interested, here's a ref to the pre-release 1D version:
van der Molen, J., K. Bolding, N. Greenwood, and D. K. Mills (2009), A
1-D vertical multiple grain size model of suspended particulate matter
in combined currents and waves in shelf seas, J. Geophys. Res., 114,
F01030, doi:10.1029/2008JF001150.
Johan