1 | ------- |
---|
2 | purpose |
---|
3 | ------- |
---|
4 | This demonstration case allows for examining the sea ice rheology under idealised |
---|
5 | forcing. The set-up follows the paper "Stress and deformation characteristics of sea |
---|
6 | ice in a high-resolution, anisotropic sea ice model" by Heorton et al, Phil. Trans, |
---|
7 | R. Soc. 2018. |
---|
8 | |
---|
9 | ----------- |
---|
10 | Experiments |
---|
11 | ----------- |
---|
12 | A range of experiments can be configured: |
---|
13 | |
---|
14 | - The suggested experiment is to examine the different shear patterns formed |
---|
15 | by different rheologies. The default set-up reproduces the leftmost panels |
---|
16 | of figure 3a (ln_rhg_EAP true in namelist_ice_cfg) and 3b (ln_rhg_EVP true in |
---|
17 | namelist_ice_cfg). Increase the last time step in namelist_cfg appropriately |
---|
18 | to generate the other panels. |
---|
19 | |
---|
20 | - The resolution of the box can be changed analogous to figure 2 of the paper |
---|
21 | (change rn_dx and rn_dy in namelist_cfg, default is 2000m) |
---|
22 | |
---|
23 | - The ratio of the windstress component can be changed to cause different modes |
---|
24 | of failure (Rwind in usrdef_sbc.F90, default is -0.8) |
---|
25 | |
---|
26 | It takes about 10 minutes to run a simulation on 16 processors |
---|
27 | |
---|
28 | ---------- |
---|
29 | How to run |
---|
30 | ---------- |
---|
31 | |
---|
32 | a) Compile and run the model once to get mesh_mask_xxxx.nc files. |
---|
33 | E.g., starting in the root NEMO directory: |
---|
34 | |
---|
35 | ./makenemo -a ICE_RHEO -r ICE_RHEO_TEST -m X64_MOBILIS -j 4 |
---|
36 | cd tests/ICE_RHEO_TEST/EXP00 |
---|
37 | mpirun -np 1 ./nemo |
---|
38 | |
---|
39 | c) Create the initial condition file for sea-ice (initice.nc) by running |
---|
40 | this python script: |
---|
41 | |
---|
42 | python ../../ICE_RHEO/EXPREF/make_INITICE.py |
---|
43 | |
---|
44 | d) Run the model a second time |
---|
45 | |
---|
46 | mpirun -np 16 ./nemo |
---|
47 | |
---|
48 | --------------- |
---|
49 | What to look at |
---|
50 | --------------- |
---|
51 | EVP rheology has less defined shear strain maxima (linear kinematic |
---|
52 | features). If the model is run longer one can see that the LKFs appearing |
---|
53 | with EVP rheology intersect at a different angle than with EAP rheology |
---|
54 | |
---|