Changeset 13794
- Timestamp:
- 2020-11-13T18:39:50+01:00 (4 years ago)
- Location:
- NEMO/branches/2019/dev_r11842_SI3-10_EAP/tests
- Files:
-
- 2 deleted
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11842_SI3-10_EAP/tests/ICE_ADV2D/MY_SRC/usrdef_sbc.F90
r13662 r13794 18 18 USE sbc_ice ! Surface boundary condition: ice fields 19 19 USE phycst ! physical constants 20 USE ice, ONLY : at_i_b, a_i_b20 USE ice, ONLY : jpl, at_i_b, a_i_b 21 21 USE icethd_dh ! for CALL ice_thd_snwblow 22 22 ! -
NEMO/branches/2019/dev_r11842_SI3-10_EAP/tests/ICE_RHEO/EXPREF/README
r12263 r13794 1 2 1 ------- 3 2 purpose 4 3 ------- 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. 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 8 9 9 ----------- … … 12 12 A range of experiments can be configured: 13 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. 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. 15 19 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) 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) 17 22 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) 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) 19 25 20 26 It takes about 10 minutes to run a simulation on 16 processors … … 24 30 ---------- 25 31 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 32 a) Compile and run the model once to get mesh_mask_xxxx.nc files. 33 E.g., starting in the root NEMO directory: 29 34 30 b) Combine the mesh mask files into one 31 ../../../tools/REBUILD_NEMO/rebuild_nemo mesh_mask 16 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 32 38 33 c) Create the initial condition file for sea-ice (initice.nc) by running this python script: 34 python ./make_INITICE.py 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 35 43 36 44 d) Run the model a second time 37 mpirun ./nemo -np 16 45 46 mpirun -np 16 ./nemo 38 47 39 48 --------------- 40 49 What to look at 41 50 --------------- 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 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 43 54 -
NEMO/branches/2019/dev_r11842_SI3-10_EAP/tests/demo_cfgs.txt
r10516 r13794 9 9 WAD OCE 10 10 BENCH OCE ICE TOP 11 ICE_RHEO OCE SAS ICE
Note: See TracChangeset
for help on using the changeset viewer.