Coriolis parameters should be more flexible

[rjleveque]

The Coriolis parameter omega is set in geoclaw_module.f90 based on the rotational velocity of the earth. We should allow this parameter to be specified in setrun.py like gravity and earth_radius are, to make it easier to change to different conditions (e.g. for modeling tsunamis on Mars, as requested by a user).

Also note that the radius of the earth is effectively hardwired into the beta-plane approximation when icoord_system == 1 but doesn't need to be. The factor 111d3 appearing in the code is 2*pi*earth_radius/360 in units of meters / degree latitude, so it seems like we should replace

theta = y / 111d3 * DEG2RAD + theta_0

by

theta = y/earth_radius + theta_0

[mandli]

👍