source:
NEMO/branches/2019/dev_r11842_SI3-10_EAP/tests/ICE_RHEO/EXPREF/README
@
12263
Last change on this file since 12263 was 12263, checked in by stefryn, 4 years ago | |
---|---|
File size: 1.7 KB |
Rev | Line | |
---|---|---|
[12263] | 1 | |
2 | ------- | |
3 | purpose | |
4 | ------- | |
5 | This demonstration case allows for examining the sea ice rheology under idealised forcing. | |
6 | The set-up follows the paper "Stress and deformation characteristics of sea ice in a | |
7 | high-resolution, anisotropic sea ice model" by Heorton et al, Phil. Trans, 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 by different rheologies. The default set-up reproduces the leftmost panels of figure 3a (ln_rhg_EAP true in namelist_ice_cfg) and 3b (ln_rhg_EVP true in namelist_ice_cfg). Increase the last time step in namelist_cfg appropriately to generate the other panels. | |
15 | ||
16 | - The resolution of the box can be changed analogous to figure 2 of the paper (change rn_dx and rn_dy in namelist_cfg, default is 2000m) | |
17 | ||
18 | - The ratio of the windstress component can be changed to cause different modes of failure (Rwind in usrdef_sbc.F90, default is -0.8) | |
19 | ||
20 | It takes about 10 minutes to run a simulation on 16 processors | |
21 | ||
22 | ---------- | |
23 | How to run | |
24 | ---------- | |
25 | ||
26 | a) Compile and run the model once to get mesh_mask_xxxx.nc files | |
27 | ../../../makenemo -a ICE_RHEO -m X64_MOBILIS -j 4 | |
28 | mpirun ./nemo -np 1 | |
29 | ||
30 | b) Combine the mesh mask files into one | |
31 | ../../../tools/REBUILD_NEMO/rebuild_nemo mesh_mask 16 | |
32 | ||
33 | c) Create the initial condition file for sea-ice (initice.nc) by running this python script: | |
34 | python ./make_INITICE.py | |
35 | ||
36 | d) Run the model a second time | |
37 | mpirun ./nemo -np 16 | |
38 | ||
39 | --------------- | |
40 | What to look at | |
41 | --------------- | |
42 | EVP rheology has less defined shear strain maxima (linear kinematic features). If the model is run longer one can see that the LKFs appearing with EVP rheology intersect at a different angle than with EAP rheology | |
43 |
Note: See TracBrowser
for help on using the repository browser.