The 3D forward models


The interest in a new 3D forward model is to compare the solutions of various codes among each other but also to investigate how well the coastal and oceanic effects are dealt with in the current codes.


Before the workshop the participants had to calculate responses from two forward models DTM3.0 and DTM3.1 (shown below) based on work by Ichihara & Mogi (2009). Periods to compute: from 100Hz to 10.000s with four periods per decade (use the rounded periods: 0.01 s, 0.018 s, 0.032 s, 0.056 s, 0.1 s, 0.18 s, 0.32 s, ...., 5600 s, 10000 s). Sites positions can be downloaded here.


Based on the workshop discussion we prepare updated DTM3.2 model. The number of periods to compute is smaller: from 10s to 10.000s with four periods per decade (use the rounded periods: 10 s, 18 s, 32 s, 56 s, 100 s, 180 s, 320 s, 560 s, 1000 s, 1800 s, 3200 s, 5600 s, 10000 s). We chose more denser sites distribution, where positions of totally 225 sites can be downloaded here (first column is east-west oriented axis). The dimension of the local and regional conductor are same as in DTM3.1, so it should be easy to modified DTM3.1 to DTM3.2. Don't forget to add 10Ohm.m conductor from 100km. Also here you can download matlab script for new bathymetry.


Please email the results to vozar@cp.dias.ie as a plain ASCII format file with the following columns: x (in km), y (in km), period (in s), ρxx, Φxx, ρxy, Φxy, ρyx, Φyx, ρyy and Φyy. It is your responsibility to send the data using the correct coordinate system for the impedance (right-handed Cartesian with x pointing North, y pointing East and z positive downwards) and using the e+iωt sign convention for the time dependency. The kind of code (FD, FE, integral equation,...) you used, mesh dimensions and the computing time together with some properties (CPU speed, memory,...) of the used computer/cluster might also be interesting for comparison. It would be great if you could include a few words about that in your email.