source: NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/ICEDYN/EXPREF @ 10419

Name Size Rev Age Author Last Change
../
1_context_nemo.xml 21 bytes 7821   4 years clem continue debugging test cases
1_namelist_cfg 15.8 KB 10073   2 years nicolasmartin Switch to test cases folder and apply previous modifications
1_namelist_ice_cfg 5.6 KB 10419   23 months smasson dev_r10164_HPC09_ESIWACE_PREP_MERGE: merge with trunk@10418, see #2133
1_namelist_ice_ref 42 bytes 9740   2 years clem repair broken links
1_namelist_ref 38 bytes 9740   2 years clem repair broken links
AGRIF_FixedGrids.in 22 bytes 9159   3 years clem change size of the zoom (enlarge)
AGRIF_FixedGrids.in_1ghost 22 bytes 9159   3 years clem change size of the zoom (enlarge)
AGRIF_FixedGrids.in_3ghosts 22 bytes 9159   3 years clem change size of the zoom (enlarge)
context_nemo.xml 2.4 KB 9930   2 years smasson trunk: reactivate predefined zoom with xios2, see #2115
domain_def_nemo.xml 45 bytes 9684   2 years smasson trunk: fix symbolic links
field_def_nemo-ice.xml 48 bytes 9726   2 years clem add missing xml file in SAS_BIPER configuration
field_def_nemo-oce.xml 48 bytes 9684   2 years smasson trunk: fix symbolic links
file_def_nemo-ice.xml 5.8 KB 10419   23 months smasson dev_r10164_HPC09_ESIWACE_PREP_MERGE: merge with trunk@10418, see #2133
file_def_nemo-oce.xml 2.2 KB 9572   2 years nicolasmartin Global renaming for reference configurations and test cases (./CONFIG): 2 …
grid_def_nemo.xml 43 bytes 9930   2 years smasson trunk: reactivate predefined zoom with xios2, see #2115
iodef.xml 1.4 KB 9509   3 years clem SASBIPER test case formally introduced in sette tests (for now …
make_INITICE.py 4.6 KB 9019   3 years timgraham Merge of dev_CNRS_2017 into branch
namelist_cfg 15.8 KB 10073   2 years nicolasmartin Switch to test cases folder and apply previous modifications
namelist_ice_cfg 6.3 KB 10419   23 months smasson dev_r10164_HPC09_ESIWACE_PREP_MERGE: merge with trunk@10418, see #2133
namelist_ice_ref 42 bytes 9684   2 years smasson trunk: fix symbolic links
namelist_ref 38 bytes 9684   2 years smasson trunk: fix symbolic links
README 2.1 KB 9789   2 years clem change the name of SAS_BIPER configuration to a more explicit one: ICEDYN

README

-------
Purpose
-------
This demonstration case can serve to different purposes:

1) The main one: study of the advection of a patch of ice in a bi-periodic channel with a slab ocean (i.e. 1 ocean layer)
              ==> set ln_icethd=false in namelist_ice_cfg

2) Secondary: study of ice thermodynamics in the same basin
              ==> set ln_icedyn=false in namelist_ice_cfg
              ==> eventually change usrdef_sbc.F90 in MY_SRC to fit your needs

-----------
Experiments
-----------
Two experiments can be configured:

a) A simple channel at 3km horizontal resolution (slightly variable => +- 10%)
              ==> remove key_agrif in cpp keys

b) Same as a) but with an Agrif zoom in the center 1:3, and 3 ghost cells
              ==> add key_agrif in cpp keys

----------
How to run
----------

a) Compile and run the model once to get a mesh_mask.nc file with the following command:
../../../makenemo -a TEST_CASES -n ICEDYN -m X64_ADA -j 4
poe ./opa -procs 1

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

c) Run the model a second time
poe ./opa -procs 4

---------------
What to look at
---------------
In case of purpose 1, 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

In case of purpose 2, one can test conservation of properties:
---
c) ice should not change at all if surface fluxes = 0 and SST = freezing temperature

--------------------------------------
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)
Note: See TracBrowser for help on using the repository browser.