Ignore:
Timestamp:
2019-01-15T17:31:25+01:00 (19 months ago)
Author:
clem
Message:

change demo config name from ICEDYN to ICE_AGRIF

Location:
NEMO/trunk/tests/ICE_AGRIF
Files:
1 edited
1 moved

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/tests/ICE_AGRIF/EXPREF/README

    r9789 r10516  
    22Purpose 
    33------- 
    4 This demonstration case can serve to different purposes: 
    5  
    6 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) 
    7               ==> set ln_icethd=false in namelist_ice_cfg 
    8  
    9 2) 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 
     4This demonstration case can be used to study the advection of a patch of ice in a bi-periodic channel 
     5with a slab ocean (i.e. 1 ocean layer) and an agrif zoom in the center 
    126 
    137----------- 
     
    1711 
    1812a) A simple channel at 3km horizontal resolution (slightly variable => +- 10%) 
    19               ==> remove key_agrif in cpp keys 
     13              ==> in MY_SRC/usrdef_hgr.F90, uncomment the part "variable scale factors" and comment part "constant scale factors" 
    2014 
    21 b) Same as a) but with an Agrif zoom in the center 1:3, and 3 ghost cells 
    22               ==> add key_agrif in cpp keys 
     15b) A simple channel at 3km horizontal resolution (constant scale factors) 
     16              ==> in MY_SRC/usrdef_hgr.F90, comment the part "variable scale factors" and uncomment part "constant scale factors" 
    2317 
    2418---------- 
     
    2721 
    2822a) Compile and run the model once to get a mesh_mask.nc file with the following command: 
    29 ../../../makenemo -a TEST_CASES -n ICEDYN -m X64_ADA -j 4 
    30 poe ./opa -procs 1 
     23../../../makenemo -r ICE_AGRIF -n ICE_AGRIF -m X64_ADA -j 4 
     24mpirun ./nemo -np 1 
    3125 
    3226b) Create the initial condition file for sea-ice (initice.nc) by running this python script:  
     
    3428 
    3529c) Run the model a second time 
    36 poe ./opa -procs 4 
     30mpirun ./nemo -np 4 
    3731 
    3832--------------- 
    3933What to look at 
    4034--------------- 
    41 In case of purpose 1, One can test 
    42 --- 
    43 a) the advection scheme: Ultimate-Macho (ln_adv_UMx=T) versus Prather (ln_adv_Pra=T) 
     35the advection through an agrif zoom 1:3 
    4436   for a square (ice concentration) or a gaussian (ice volume) 
    45       with either a constant velocity (ln_dynADV=T) 
     37      with either a constant velocity (ln_dynADV2D=T) 
    4638               or a constant ice-atm. stress, thus velocity is calculated by rheology (ln_dynRHGADV=T) 
    4739      with 1 or 5 ice categories (jpl=1 or 5) 
    48  
    49 b) the advection through an agrif zoom 1:3 
    50  
    51 In case of purpose 2, one can test conservation of properties: 
    52 --- 
    53 c) ice should not change at all if surface fluxes = 0 and SST = freezing temperature 
    54  
    55 -------------------------------------- 
    56 Interpretation of the results (remarks) 
    57 -------------------------------------- 
    58 - Prather conserves the max values but also creates side lobes 
    59 - UM does not conserve the max but does not create side lobes 
    60 - 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) 
     40         (note that ln_dynADV2D=T only works with jpl=1) 
Note: See TracChangeset for help on using the changeset viewer.