LIGGGHTS WWW Site - LAMMPS WWW Site - LIGGGHTS Documentation - LIGGGHTS Commands

fix seaice/coriolis command

Syntax:

fix ID group-ID seaice/coriolis phi keyword value ... 

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).