New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 13794 for NEMO/branches/2019 – NEMO

Changeset 13794 for NEMO/branches/2019


Ignore:
Timestamp:
2020-11-13T18:39:50+01:00 (3 years ago)
Author:
acc
Message:

Branch dev_r11842_SI3-10_EAP. Tweaks to ensure ICE_RHEO test case can actually be built and run as suggested in its README file

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  
    1818   USE sbc_ice         ! Surface boundary condition: ice fields 
    1919   USE phycst          ! physical constants 
    20    USE ice, ONLY       : at_i_b, a_i_b 
     20   USE ice, ONLY       : jpl, at_i_b, a_i_b 
    2121   USE icethd_dh       ! for CALL ice_thd_snwblow 
    2222   ! 

  • NEMO/branches/2019/dev_r11842_SI3-10_EAP/tests/ICE_RHEO/EXPREF/README

    r12263 r13794  
    1  
    21------- 
    32purpose 
    43------- 
    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. 
     4This demonstration case allows for examining the sea ice rheology under idealised 
     5forcing.  The set-up follows the paper "Stress and deformation characteristics of sea 
     6ice in a high-resolution, anisotropic sea ice model" by Heorton et al, Phil. Trans, 
     7R. Soc. 2018. 
    88 
    99----------- 
    … …  
    1212A range of experiments can be configured: 
    1313 
    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. 
    1519 
    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) 
    1722 
    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)  
    1925 
    2026It takes about 10 minutes to run a simulation on 16 processors 
    … …  
    2430---------- 
    2531 
    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 
     32a) Compile and run the model once to get mesh_mask_xxxx.nc files. 
     33   E.g., starting in the root NEMO directory: 
    2934 
    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 
    3238 
    33 c) Create the initial condition file for sea-ice (initice.nc) by running this python script: 
    34 python ./make_INITICE.py 
     39c) 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 
    3543 
    3644d) Run the model a second time 
    37 mpirun ./nemo -np 16 
     45 
     46 mpirun -np 16 ./nemo 
    3847 
    3948--------------- 
    4049What to look at 
    4150--------------- 
    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 
     51EVP rheology has less defined shear strain maxima (linear kinematic 
     52features). If the model is run longer one can see that the LKFs appearing 
     53with EVP rheology intersect at a different angle than with EAP rheology 
    4354 

  • NEMO/branches/2019/dev_r11842_SI3-10_EAP/tests/demo_cfgs.txt

    r10516 r13794  
    99WAD OCE 
    1010BENCH OCE ICE TOP 
     11ICE_RHEO OCE SAS ICE 
Note: See TracChangeset for help on using the changeset viewer.