Changeset 10516 for NEMO/trunk/tests/ICE_AGRIF/EXPREF/README

Timestamp:

2019-01-15T17:31:25+01:00 (5 years ago)

Author:

clem

Message:

change demo config name from ICEDYN to ICE_AGRIF

Location:

NEMO/trunk/tests/ICE_AGRIF

Files:

• . (moved) (moved from NEMO/trunk/tests/ICEDYN)
• EXPREF/README (modified) (4 diffs)

NEMO/trunk/tests/ICE_AGRIF/EXPREF/README

@@ -2 +2 @@
 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
+This demonstration case can be used to study the advection of a patch of ice in a bi-periodic channel
+with a slab ocean (i.e. 1 ocean layer) and an agrif zoom in the center
 
 -----------
@@ -17 +11 @@
 
 a) A simple channel at 3km horizontal resolution (slightly variable => +- 10%)
-              ==> remove key_agrif in cpp keys
+              ==> in MY_SRC/usrdef_hgr.F90, uncomment the part "variable scale factors" and comment part "constant scale factors"
 
-b) Same as a) but with an Agrif zoom in the center 1:3, and 3 ghost cells
-              ==> add key_agrif in cpp keys
+b) A simple channel at 3km horizontal resolution (constant scale factors)
+              ==> in MY_SRC/usrdef_hgr.F90, comment the part "variable scale factors" and uncomment part "constant scale factors"
 
 ----------
@@ -27 +21 @@
 
 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
+../../../makenemo -r ICE_AGRIF -n ICE_AGRIF -m X64_ADA -j 4
+mpirun ./nemo -np 1
 
 b) Create the initial condition file for sea-ice (initice.nc) by running this python script: 
@@ -34 +28 @@
 
 c) Run the model a second time
-poe ./opa -procs 4
+mpirun ./nemo -np 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)
+the advection through an agrif zoom 1:3
    for a square (ice concentration) or a gaussian (ice volume)
-      with either a constant velocity (ln_dynADV=T)
+      with either a constant velocity (ln_dynADV2D=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 that ln_dynADV2D=T only works with jpl=1)