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
mandli
commented
1 year ago
The Coriolis parameter
omegais set ingeoclaw_module.f90based on the rotational velocity of the earth. We should allow this parameter to be specified insetrun.pylikegravityandearth_radiusare, 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 == 1but doesn't need to be. The factor111d3appearing in the code is2*pi*earth_radius/360in units of meters / degree latitude, so it seems like we should replaceby