[8900] | 1 | ------- |
---|
[8962] | 2 | Purpose |
---|
[8900] | 3 | ------- |
---|
[8962] | 4 | This demonstration case can serve to different purposes: |
---|
[8900] | 5 | |
---|
[8962] | 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 |
---|
| 12 | |
---|
| 13 | ----------- |
---|
| 14 | Experiments |
---|
| 15 | ----------- |
---|
[8900] | 16 | Two experiments can be configured: |
---|
[8962] | 17 | |
---|
[8900] | 18 | a) A simple channel at 3km horizontal resolution (slightly variable => +- 10%) |
---|
[8962] | 19 | ==> remove key_agrif in cpp keys |
---|
| 20 | |
---|
[8900] | 21 | b) Same as a) but with an Agrif zoom in the center 1:3, and 3 ghost cells (defined in par_oce.F90) |
---|
[8962] | 22 | ==> add key_agrif in cpp keys |
---|
[8900] | 23 | |
---|
| 24 | ---------- |
---|
| 25 | How to run |
---|
| 26 | ---------- |
---|
| 27 | Be careful, the model can only run on 1 processor in the N-S dimension |
---|
| 28 | |
---|
[8962] | 29 | a) Compile and run the model once to get a mesh_mask.nc file with the following command: |
---|
[8900] | 30 | ../../../makenemo -a TEST_CASES -n SAS_BIPER -m X64_ADA -j 4 |
---|
| 31 | poe ./opa -procs 1 |
---|
| 32 | |
---|
[8962] | 33 | b) Create the initial condition file for sea-ice (initice.nc) by running this python script: |
---|
[8900] | 34 | python ./make_INITICE.py |
---|
| 35 | |
---|
[8962] | 36 | c) Run the model a second time |
---|
[8900] | 37 | poe ./opa -procs 1 |
---|
| 38 | |
---|
| 39 | --------------- |
---|
| 40 | What to look at |
---|
| 41 | --------------- |
---|
[8962] | 42 | In case of purpose 1, One can test |
---|
| 43 | --- |
---|
[8900] | 44 | a) 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 | |
---|
| 50 | b) the advection through an agrif zoom 1:3 |
---|
| 51 | |
---|
[8962] | 52 | In case of purpose 2, one can test conservation of properties: |
---|
| 53 | --- |
---|
| 54 | c) ice should not change at all if surface fluxes = 0 and SST = freezing temperature |
---|
| 55 | |
---|
[8900] | 56 | -------------------------------------- |
---|
| 57 | Interpretation 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) |
---|