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.
2019WP/SI3-10_EAP_rheology (diff) – NEMO

Changes between Version 2 and Version 3 of 2019WP/SI3-10_EAP_rheology


Ignore:
Timestamp:
2020-05-28T16:23:00+02:00 (4 years ago)
Author:
stefryn
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • 2019WP/SI3-10_EAP_rheology

    v2 v3  
    2323 
    2424Current EVP rheology assumes ice is isotropic, but observations show linear kinematic features (LKF) which are sign of anisotropic behaviour. 
    25 Elasitc anisotropic plastic (EAP) rheology has been developed by UCL to address this shortcoming (refs 1-3 theorical development, refs 4-5 for implemation in CICE). It is routinely used in CICE simulations and shows more realistic LKF patterns.[[BR]] 
     25Elastic anisotropic plastic (EAP) rheology has been developed by UCL to address this shortcoming (refs 1-3 theoretical development, refs 4-5 for implementation in CICE). It is routinely used in CICE simulations and shows more realistic LKF patterns.[[BR]] 
    2626 
    2727The IMMERSE project has a milestone and deliverable to implement EAP rheology in SI3. 
     
    171171=== Test case === 
    172172 
    173 We wil use the ice_adv2d testcase modified to examine the effect of the rheology by imposing diffent controlled boudary conditions. The approach is inspired by the idealised simulations in Heorton et al. 2018. 
     173We will use the ice_adv2d testcase modified to examine the effect of the rheology by imposing different controlled boundary conditions. The approach is inspired by the idealised simulations in Heorton et al. 2018. 
    174174 
    175175=== Documentation updates 
     
    184184 
    185185update: EAP rheology has an impact on the closing rate for ridging/rafting. Currently this is hardcoded appropriate for EVP in icedyn_rdgrft.F90. An option will be added in icedyn_rdgrft.F90 to use EAP rate instead. 
    186 update 2: Current implementation in CICE ignores the impact of advection on the anisotropy. Initial implementation in SI3 will do the same, so no extra tracer is needed for now. This can be reassed later. 
     186update 2: Current implementation in CICE ignores the impact of advection on the anisotropy. Initial implementation in SI3 will do the same, so no extra tracer is needed for now. This can be reasessed later. 
    187187 
    188188== Tests 
     
    192192}}} 
    193193 
    194 ''...'' 
     194The test case is called ICE_RHEO, included in the development branch in the tests folder. It is based on the ice_adv2d test case, but with some differences in the forcing and boundary conditions, see the usrdef files in the MY_SRC directory. There is no single verification value, rather the intention is to examine the deformation patterns caused by different rheologies, which can be compared to those in Heorton et al. 2018 (CICE-based). The rheologies are EVP vs EAP, but it could be used for other future rheology options too. 
     195 
     196Current status: there is still an issue preventing to reproduce Harry's results, probably in the boundary conditions, which I need to recheck with him. 
     197 
     198Convergence of the velocities can be checked with code from changeset 12741 which is not yet incorporated in this branch. 
    195199 
    196200== Review