lauren

ADCIRC outputs for PTM

13 posts in this topic

I'm a new SMS user, and am using ADCIRC to generate the input files for a PTM model.

I've got the coastlines, mesh, and bathymetry data working, but my question is: how do I generate the XMDF data files containing water surface elevations and currents? From what I understand ADCIRC is supposed to generate fort.63 (global elevation) and fort.64 (global velocity) files when run, based on tidal forcing. I'm modeling a canal off of a lake, however, so I don't want velocities generated from tides. What I'd like to do is input SWE and velocity at several points (from gage data), extend that to the rest of the mesh, and then save that in XMDF format for use in PTM.

Is this possible?

Share this post


Link to post
Share on other sites

First create a scatter set of the WSE and velocity data you want to use. You can then:

  1. Interpolate the scatter set values to the mesh
  2. Use the data set toolbox (Temporal | Sample time steps) to create a transient data set if your data is steady state
  3. In the PTM model control, on the files tab use "Create input file(s) from data" to create a mesh file and XMDF flow (hydrodynamic) file

Feel free to post follow up questions if this isn't clear.

- John

Share this post


Link to post
Share on other sites

Another option would be to create an ASCII dataset file and then read that into SMS and export it as an XMDF file as described above.

Share this post


Link to post
Share on other sites

Thank you John. I am working on putting the velocity and SWE scatter data together. When I load those files in, do I do it in the same project where I have the bathymetry data, or do I start a new one?

Also: when I open the .sms file where the bathymetry mesh is, I get two error messages.

"NETA does not match the number of open boundary nodes." and

"There is a forcing boundary but no constituents specified and the fort.19 cannot be read."

What is NETA?

How can I turn off the forcing boundary? I don't think tidal constituents make sense in my case.

Thanks to you again!

Share this post


Link to post
Share on other sites

I prefer to keep everything in one project. This will be useful when you do post processing after you run PTM so you can see the particles move through the domain and show velocity vectors, contours, and maybe even a functional surface. PTM animations are usually pretty interesting.

The best place to learn about NETA and the other ADCIRC file parameters is at http://www.adcirc.org/

PTM won't use all the parameters ADCIRC requires - it basically just needs the mesh file so it knows the shape of the domain and where particles can leave the domain (ocean boundary). Then of course the hydrodynamics also.

Share this post


Link to post
Share on other sites

I have created a scatter set of the velocity data I want to use, but have been unsuccessful at interpolating the vectors to the mesh. I can interpolate the scalar values to the mesh, but how do I apply direction to them after-the-fact? Thanks.

Share this post


Link to post
Share on other sites

To interpolate vectors:

  1. In the Dataset toolbox, split the vector into two scalar quantities (Vx and Vy or Magnitude and Direction)
  2. Interpolate the two scalar datasets to the mesh
  3. In the Dataset toolbox, recombine the two scalar quantities into a vector

Share this post


Link to post
Share on other sites

Hi John. Thanks again for all your help. I have another question.

You said:

"Use the data set toolbox (Temporal | Sample time steps) to create a transient data set if your data is steady state."

How do I do this? I'd like to add time steps to my velocity and water surface elevation data. I don't see a way to create a "Temporal" data set. Under the "data" menu the options are "Geometry" and "Coastal," neither of which seem to apply.

Thanks again.

Share this post


Link to post
Share on other sites

A couple of questions:

  • Are you using SMS 10.1?
  • Which module is the data in?

Share this post


Link to post
Share on other sites

No, I am using SMS 10.0.

I imported the velocity and water surface data as scatter sets, then interpolated them to the mesh.

Share this post


Link to post
Share on other sites

Lauren,

The ability to create a transient dataset from a steady state dataset is a new feature in SMS 10.1. Unfortunately, SMS 10.0 does not have this functionality.

Share this post


Link to post
Share on other sites

John,

Thanks again for your assistance. Right now my PTM model will not run because it says that my hydrodynamic data is temporally shorter than the simulation run. I understand that I cannot create transient data sets from within SMS 10.0. But is it possible to "add time" to steady-state velocity data, so it applies for a longer timeframe?

ERROR: *** FLOW FILE DOES NOT RUN BEYOND STOP TIME ***

Thanks John.

Share this post


Link to post
Share on other sites

Lauren,

Unfortunately the model PTM does not have an option to loop or extrapolate hydrodynamic data. Depending on what model your hydrodynamic output is from, it may be possible to manually alter the file to add additional timesteps.

- John

1 person likes this

Share this post


Link to post
Share on other sites

Create an account or sign in to comment

You need to be a member in order to leave a comment

Create an account

Sign up for a new account in our community. It's easy!


Register a new account

Sign in

Already have an account? Sign in here.


Sign In Now