source: branches/2017/dev_CNRS_2017/NEMOGCM/CONFIG/TEST_CASES/SAS_BIPER/EXP00/README @ 8900

-------
Context
-------
The aim of this demonstration case is to study the advection of a patch of ice in a bi-periodic channel with a slab ocean (i.e. 1 ocean layer which does not interact with sea-ice).

Two experiments can be configured:
a) A simple channel at 3km horizontal resolution (slightly variable => +- 10%)
b) Same as a) but with an Agrif zoom in the center 1:3, and 3 ghost cells (defined in par_oce.F90)

----------
How to run
----------
Be careful, the model can only run on 1 processor in the N-S dimension

Compile and run the model once to get a mesh_mask.nc file with the following command (note; remove key_agrif in cpp keys for experiment a):
../../../makenemo -a TEST_CASES -n SAS_BIPER -m X64_ADA -j 4
poe ./opa -procs 1

Create the initial condition file for sea-ice (initice.nc) by running this python script: 
python ./make_INITICE.py

Run the model a second time
poe ./opa -procs 1

---------------
What to look at
---------------
One can test

a) the advection scheme: Ultimate-Macho (ln_adv_UMx=T) versus Prather (ln_adv_Pra=T)
   for a square (ice concentration) or a gaussian (ice volume)
      with either a constant velocity (ln_dynADV=T)
               or a constant ice-atm. stress, thus velocity is calculated by rheology (ln_dynRHGADV=T)
      with 1 or 5 ice categories (jpl=1 or 5)

b) the advection through an agrif zoom 1:3

--------------------------------------
Interpretation of the results (remarks)
--------------------------------------
- Prather conserves the max values but also creates side lobes
- UM does not conserve the max but does not create side lobes
- The "unmoving" ice (if any) is due to ice rheology which states that ice mass below a certain thresold (1kg/m2) is considered to move at the ocean velocity (thus 0 m/s)