search

MPAS-Dev / MPAS

Repository for private MPAS development prior to the MPAS v6.0 release.
Other
4 stars 0 forks source link

Initial commit of iceberg advection solver. #1446

Closed darincomeau closed 7 years ago

darincomeau commented 7 years ago

This PR is a work in progress of the iceberg advection solver using incremental remapping. This PR rearranges and adds iceberg variables in the Registry to be consistent with IR. A new tracer linked list, bergTracersHead, is defined in mpas_cice_advection_incremental_remap_tracers.F, where bergAreaCategory is the parent tracer, and bergHeight, bergLength, are child tracers. A logical advectBergs is used and passed through IR subroutines as an optional argument to indicate if sea ice or icebergs are to be used. The code is currently failing inside the incremental_remap_block subroutine when called for icebergs (line 3026 of mpas_cice_advection_incremental_remap.F sets advectBergs = .false. for now to avoid this for debugging). The code aborts from an MPI_ABORT when run in parallel, but runs to completion when run in serial, though icebergs are not moving.

BFB when config_use_bergs = false.

darincomeau commented 7 years ago

@akturner Can you add the Work-In-Progress and Don't Merge tags?