← Back to branch summary

~fluidity-core/fluidity/ice-melt-mesh-movement

« back to all changes in this revision

Viewing changes to main/Fluids.F90

  • Committer: Stephan Kramer
  • Date: 2013-05-18 13:41:19 UTC
  • Revision ID: s.kramer@imperial.ac.uk-20130518134119-zn5e6pqem6p02yys
Mesh movement according to ice melt.

To use choose ice_melt as mesh movement scheme and add a prescribed field 
"ThermalConductivity" with the conductivity values for water and ice in the
corresponding regions. As the ice melt depends on the normal temperature
gradient at the interface you'll need enough resolution to resolve the 
diffusive boundary layer.

Show diffs side-by-side

added added

removed removed

Lines of Context:
main/Fluids.F90
517
517
       ! to ensure those properties are satisfied on the new mesh not the old one.
518
518
       call move_mesh_imposed_velocity(state)
519
519
       call move_mesh_pseudo_lagrangian(state)
 
520
       call move_mesh_ice_melt(state)
520
521
 
521
522
       call enforce_discrete_properties(state, only_prescribed=.true., &
522
523
            exclude_interpolated=.true., &