Syntax:
fix ID group-ID seaice/coriolis phi keyword value ...
phi can be a variable (see below)
region value = region-ID region-ID = ID of region atoms must be in to have added force
Examples:
fix f1 all seaice/coriolis 60.0 fix f1 all seaice/coriolis -70.0 region regID
Description:
This fix is part of the USER-SEAICE package; it is only enabled if LIGGGHTS was built with that package.
It adds Coriolis force to the x and y components of force for each disk-shaped atom in the group.
It is assumed that the Coriolis force is horizontal, i.e., it does not contribute to the force along the z-axis. It is also assumed that the Coriolis force does not contribute to the torque acting on the disks.
This fix requires a global property, rhoIce, that can be set with the fix property/global command:
fix pg1 all property/global rhoIce scalar 910.0
No default value for this constant is predefined and the lack of the above definition in the input script will produce an error.
The values of phi (in degrees) are positive in the northern hemisphere and negative in the southern hemisphere. Phi can be specified as an equal-style or atom-style variable. If the value is a variable, it should be specified as v_name, where name is the variable name. In this case, the variable will be evaluated at each timestep. For more details concerning specifying variable forcing, see the documentation of the fix addforce command.
If the region keyword is used, the atom must also be in the specified geometric region in order to have the Coriolis force added to it.
Restart, fix_modify, output, run start/stop, minimize info:
See the documentation of the fix addforce command.
Restrictions: none
Related commands:
fix addforce, fix seaice/current, fix seaice/waves, fix seaice/wind
Default: none
(Herman) Agnieszka Herman, Geosci. Model Devel., submitted (2015).