[MITgcm-devel] ICE workshop - recap and homework!

[Jean-Michel Campin]

Hi,

Some details below of the new fwd lab_sea tests.
If you have suggestions or if you see a bad combination of
params, let me know.
Also, I am not sure if this (lab_sea/input.salt_plume/data.seaice)
 is the right way  of doing SEAICE_VARIABLE_SALINITY with SItracers ? 
Gael, can you confirm ?
Right now, it seems that siHSALT from the pickup-file is not
used to initialise the salinity SItracer.

Cheers,
Jean-Michel

--------------------------------------------------------------------
Details:

lab_sea/code/SEAICE_OPTIONS.h
 I took the one that Martin (version.2) sent, 
  - remove SEAICE_VARIABLE_FREEZING_POINT
    (has been removed completely from the seaice code)
  - define SEAICE_ALLOW_FREEDRIFT (so that it can be tested in input.fd)
  - update the description of SEAICE_MULTICATEGORY
  - some minor editing

input/data.seaice:
 I took the one that Martin (version.2) sent, change the setting order
  (dynamics, seaice state & transport, thermodyn, IO params)
  set SEAICEheatConsFix=.TRUE., (as we discuss)
 and added 
  SEAICE_multDim  = 7, (but this is the current default) 
 and also kept the 3rd namelist with the same 2 SItracers as before.

input.fd/data.seaice: similar to input/data.seaice, except:
 - SEAICEuseFREEDRIFT = .TRUE.,
   SEAICEheatConsFix commented out (= default = F)
 - remove SItracers ; 
 - SEAICEadvSchHeff   = 33,
  # add small diffKh to test diffusion with multi-dim advect.
   SEAICEdiffKhHeff   = 20.,
 - SEAICE_frazilFrac  = 0., 
  so that we have 1 lab_sea test that does not freeze sea-water if below freezing.

input.hb87/data.seaice: similar to input/data.seaice, except:
 - SEAICE_deltaTevp   = 30.,
 - useHB87stressCoupling = .TRUE.,
 - no SItracers (as it used to be), no SEAICEheatConsFix
 - all 4 initial condition files (same as before)
 - use default Adv.Scheme with DIFF1 = 4.E-3,
   (so that it's tested in one of the lab-sea test)
 and also added:
 - SEAICE_multDim     = 1, (will run single category)
  # strong ice-ocean turb-flux coeff:
   SEAICE_mcPheePiston= 0.001,
   SEAICE_frazilFrac  = 0.5,

input.salt_plume/data.seaice: 
 has not been changed so much (kept old solve4temp)
 except: Adv-scheme=7 (since the ocean is also using 7), 
         const. freezing point, + salinity as SItracer :
    SEAICE_salinityTracer = .TRUE.,
  # SEAICE_saltFrac    = 0.3,
  #- constant seawater freezing point:
    SEAICE_tempFrz0    = -1.96,
    SEAICE_dTempFrz_dS = 0.,
with 2 SItracers:
    SItrNumInUse       = 2,
    SItrName(1)        = 'ridge',
    SItrMate(1)        = 'AREA',
    SItrName(2)        = 'salinity',
    SItrFromOceanFrac(2) = 0.3,

--------------------------------------------------------------------
And here is the fwd-lab_sea relevant part of seaice-verif-notes.txt
that Patrick sent:

* lab_sea
---------
  - give up legacy
  - try different setups with different options
    forw1: input: exact baseline parameters + SEAICEheatConsFix=.TRUE. + LSR
    forw2: input.fd: SEAICEheatConsFix=.FALSE., SEAICEuseFREEDRIFT=.TRUE., --- start June 1st
    forw3: input.salt_plume: const. freezing point, old solve4temp + salinity as SItracer
    forw4: input.hb87: forw1 + EVP + useHB87stressCoupling