[MITgcm-devel] Add sea ice surface forcing to pkg/seaice

[Menemenlis, Dimitris (3248)]

Jean-Michel, a couple of comments and questions about verification experiments:

1. How do I get permission to testreport adjoint verification experiments on baudelaire?

2. I was not able to test the adjoint seaice verification experiments.
    (I think that I am using an outdated version of staf and have contacted Thomas about it.)
    So there is a small possiblity that some adjoint seaice experiments will fail tonight.

3. I commented out snowprecipfile definition in
   global_ocean.cs32x15/input.seaice/data.exf
   global_ocean.cs32x15/input_ad.seaice/data.exf
   global_ocean.cs32x15/input_ad.seaice_dynmix/data.exf
because snowprecip is now active in pkg/seaice and this will change results.

4. Would it be OK to modify global_ocean.cs32x15/input.seaice
    as a test for snowprecip in pkg/seaice ?

Dimitris Menemenlis

On Apr 21, 2012, at 6:31 PM, Dimitris Menemenlis wrote:

> I have added snowPrecip forcing for pkg/seaice:
> 
>            HSNOW(I,J,bi,bj) = HSNOW(I,J,bi,bj) +
>      &          convertPRECIP2HI*ICE2SNOW*snowPrecip(i,j,bi,bj)
> 
> This is a quick-and-dirty way to simulate collapse of tidewater glaciers.
> Contrary to Precip, snowPrecip is not multiplied by sea ice
> concentration (AREA).  It falls as snow even on top of open water,
> helps cool SST, and eventually through "SEAICEuseFlooding"
> it generates some sea ice.  
> 
> #define SEAICE_VARIABLE_SALINITY
> is needed so that the useFlooding sea ice is added as fresh-water sea ice.