[MITgcm-devel] pkg/icefront

[Dimitris Menemenlis]

Martin and Jean-Michel, for application to Greenland glaciers, I am thinking of
extending pkg/shelfice so that it can deal with vertical ice fronts.  This would involve:

1) a routine (similar to dwnslp_init_fixed.F) which locates 3 things:
   a) horizontal locations of wet grid points adjacent to glacier calving faces
   b) vertical location of points in a)
   c) ratio of vertical face area to horizontal grid area for points in a)

2) a routine similar to shelfice_thermodynamics which applies the ice shelf melt/freeze
    equations to locations detected in 1) above and which scales them by 1c

3) a routine similar to shelfice_forcing which applies the forcing computed
     in 2) to locations computed in 1)

The question I have is whether above computations should be incorporated in
pkg/shelfice or whether they should be carried out in a separate package, e.g.,
pkg/icefront, that reproduces most of the functionality of pkg/shelfice ?

Let me know preferred way to proceed.  D.

Dimitris Menemenlis <menemenlis at jpl.nasa.gov>
Jet Propulsion Lab, California Institute of Technology
MS 300-323, 4800 Oak Grove Dr, Pasadena CA 91109-8099, USA
tel: 818-354-1656;  cell: 818-625-6498;  fax: 818-393-6720