New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
README in NEMO/trunk/tests/SAS_BIPER/EXPREF – NEMO

source: NEMO/trunk/tests/SAS_BIPER/EXPREF/README @ 9740

Last change on this file since 9740 was 9019, checked in by timgraham, 6 years ago

Merge of dev_CNRS_2017 into branch

File size: 2.2 KB
Line 
1-------
2Purpose
3-------
4This demonstration case can serve to different purposes:
5
61) 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)
7              ==> set ln_icethd=false in namelist_ice_cfg
8
92) Secondary: study of ice thermodynamics in the same basin
10              ==> set ln_icedyn=false in namelist_ice_cfg
11              ==> eventually change usrdef_sbc.F90 in MY_SRC to fit your needs
12
13-----------
14Experiments
15-----------
16Two experiments can be configured:
17
18a) A simple channel at 3km horizontal resolution (slightly variable => +- 10%)
19              ==> remove key_agrif in cpp keys
20
21b) Same as a) but with an Agrif zoom in the center 1:3, and 3 ghost cells (defined in par_oce.F90)
22              ==> add key_agrif in cpp keys
23
24----------
25How to run
26----------
27Be careful, the model can only run on 1 processor in the N-S dimension
28
29a) Compile and run the model once to get a mesh_mask.nc file with the following command:
30../../../makenemo -a TEST_CASES -n SAS_BIPER -m X64_ADA -j 4
31poe ./opa -procs 1
32
33b) Create the initial condition file for sea-ice (initice.nc) by running this python script:
34python ./make_INITICE.py
35
36c) Run the model a second time
37poe ./opa -procs 1
38
39---------------
40What to look at
41---------------
42In case of purpose 1, One can test
43---
44a) the advection scheme: Ultimate-Macho (ln_adv_UMx=T) versus Prather (ln_adv_Pra=T)
45   for a square (ice concentration) or a gaussian (ice volume)
46      with either a constant velocity (ln_dynADV=T)
47               or a constant ice-atm. stress, thus velocity is calculated by rheology (ln_dynRHGADV=T)
48      with 1 or 5 ice categories (jpl=1 or 5)
49
50b) the advection through an agrif zoom 1:3
51
52In case of purpose 2, one can test conservation of properties:
53---
54c) ice should not change at all if surface fluxes = 0 and SST = freezing temperature
55
56--------------------------------------
57Interpretation of the results (remarks)
58--------------------------------------
59- Prather conserves the max values but also creates side lobes
60- UM does not conserve the max but does not create side lobes
61- 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.