summary of lnd=>glc exchange for CLM=>CISM and CLM=>DGLC%NOEVOLVE

[mvertens]

for CLM and (CISM or DGLC%NOEVOLVE) the following hold:

is_local%wrap%lnd2glc_coupling = .true.
is_local%wrap%accum_lnd2glc = .true.

• Accumulate land input
  • input for glc is accumulated every land coupling period on the land grid
    by the call to med_phases_prep_glc_accum_lnd in med_phases_post_lnd
  • is done on the land grid for the following fields:
    fldnames_fr_lnd(3) = (/'Flgl_qice_elev','Sl_tsrf_elev ','Sl_topo_elev '/)
  • is averaged in med_phases_prep_glc_avg which is called from med_phases_prep_glc in the run sequence
• Map accumulated land input
  • field bundle from land grid (with elevation classes) to glc grid (without elevation classes)
    based on the setting of the xml variable GLC_AVG_PERIOD which can either be glc_coupling_period or yearly(default)
  • done in med_phases_prep_glc_map_lnd2glc (which is called from med_phases_prep_glc_avg)
  • map the accumulate land field from the land grid (in multiple elevation classes)
    to the glc grid (in multiple elevation classes) using bilinear interpolation
    • initialize accumulated field bundle on the glc grid to zero before doing the mapping
    • map accumlated land fields to each ice sheet (normalize by the land fraction in the mapping)
    • determine elevation class of each glc grid gridcell (elevclass_g)
  • loop over fields in export field bundle to glc for ice sheet and
    perform vertical interpolation of data onto ice sheet topography -
    this maps all of the input elevation classes into an export to glc without elevation classes
  • renormalize surface mass balance (i.e. dataexp_g) so that the global
    integral on the glc grid is equal to the global integral on the land grid.
    this is done in the call to med_phases_prep_glc_renormalize_smb.
• Set FBExp(compglc(ns)) data
• Zero land accumulator and accumulated field bundles on land grid