Use the Mader's SWAN code to model tsunami propagation. To use this program you need to have two grids that are EXACTLY compatible. This means they must cover the same region and have the same number of rows/columns. One of the grids represents the bathymetry (z positive up) in meters and other contains the initial deformation as produced by the Okada deformation module. The next thing you need is a parameter file. In the "data" directory of the root Mirone installation you will find a file named swan.par. Use this as an example for your case. I don't know the purpose of most of the parameters, only of those who seems to be important. The swan.par has short comments in front of each parameter.
Now that you are ready to do the simulation you have to choose between computing the simulation on specific points - called maregraphs here - or on the whole grid region (but see the NOTE below). If you want the first option, than you either plot those stations graphically (using the Plot Stations option) or import a x,y file with the maregraphs location. This later option is also available in the swan_options window. If you want to compute the simulation in the whole region, then proceed with the Compute option. In the swan_options window you will finally select what to do. The Bat and Source edit boxes must contain either the name of your grids or a gray line
with the In memory array string, which appears when the Swan module was called from a Mirone window in which a valid grid was previously loaded. In the Params edit box you load the swan.par parameter file mentioned above. If you already plotted stations in the bathymetry window grid the Maregraphs section will be automatically set, otherwise check the check box and give the necessary file names to compute the maregraphs (in case you want to do it - not mandatory). Maregraphs are nothing but a record of the water height vs time at a certain location(s). The Output grids check box tells the program that you want to write grids on disk containing the water height at particular time intervals of the simulation cycle. For example if in the swan.par dt = 1 and cumint = 30 that means the simulation is computed at increments of one second and grids (as well as maregraphs if that option was selected) are written at every 30 seconds of the simulation. Grids are automatically named after the name template provided in the edit box appended with the cycle number. As an alternative to writing grids on disk, you can compute an AVI movie. In Nș of cycles you select for how many cycles will the simulation is run. Again, using the same example of dt = 1 the default value of 1010 means the simulation will extend for 1010 seconds.
NOTE: if a rectangle exists on the window that contains the bathymetry file and you have selected to save grids, than that rectangle works as a sub-region selection. That is, the grids which are written cover only the region delimited by the rectangle. This option is original meant for creating input files to run the swan2sww command line utility which creates a .sww netCDF file to use with the ANUGA tsunami modeling code. But this is still a very fresh story.
I left behind on purpose the Plot stations on grid borders Swan option. This option's purpose is meant only for using in the Tsun2 module.
Copyright © <2010>, <Joaquim Luis>
Created with the Freeware Edition of HelpNDoc: Single source CHM, PDF, DOC and HTML Help creation