Regarding Inundation Modelling with Mirone (tsun2)

[Patanjali]

Hi Joaquim luis,

I am new to this group and let me thank you for sharing a fantastic
stand alone program built around using matlab and gmt.

I am new to Tsunami modelling and I have explored the tsunami
inundation tutorial provided by you., Thanks a lot for the tutorial
and I have some doubts as follows:

Can you please let me know the grid preparation (especially in UTM
coordinate system) or grid conversion directly into UTM from GEBCO for
my area of interest to use in tsunami inundation modeling ? As the
bathymetric data from GEBCO is in Cartesian coordinate system (please
correct me, if I am wrong).

Is there any in-built routine in mirone package to convert the GEBCO
data to UTM coordinate system (similar to data format provided in
inundation tutorial)?

Is it necessary to convert the bathymetric grids into UTM coordinate
system to do the inundation modelling?

Can we do the inundation modelling in Cartesian coordinate system
also?


Please clarify my doubts, thanks in advance

with regards

Patanjali

[jluis]

Hi Patanjali

> I am new to this group and let me thank you for sharing a fantastic
> stand alone program built around using matlab and gmt.

Thanks

Mirone has one tool to easily convert grids in geographic into UTM.
Access it via "Projections -> GDAL project"
First box should have
+proj=latlong
and second
+proj=utm +zone=29 +datum=WGS84
(in above replace the zone for your data)

However, this does no provide a fine control on the grid spacing and
region. When I need to do that I follow the GMT recipe, that I
illustrate below. For exact commands you must consult GMT manual, and
they depend obviously on your data
grd2xyz grid_geog.grd > xyz_geo.dat
mapproject xyz_geo.dat -I -Ju??/1:1 -F -C > xyz_utm.dat
surface xyz_utm.dat -R.... -I.. -Ggrid_out.grd

This said, Mirone's solution for inundation is quite old (it dates
pre-2004 event) and it's not even very correct as it doesn't transmit
the momentum.
On the other hand, I'm sorry but GEBCO data is no good to do
inundation studies. For one side it's not very good at all and for
other it's resolution is way too coarse to be used for inundations
where you want to have at least 50 m grid spacings.
What I recommend you to learn/use is ANUGA. Mirone is still very
usefull here as you can compute the oceanic progation with the SWAN
model and output the data in .sww or most formats that you will later
use to feed the boundary conditions of the ANUGA run.
I'm sorry but I don't do that for some time now and so I don't have
the details fresh in mind to explain more. Perhaps Rjaraman has.

Good luck

Joaquim Luis