Changes between Version 2 and Version 3 of 2019WP/SI3-10_EAP_rheology
- Timestamp:
- 2020-05-28T16:23:00+02:00 (5 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
2019WP/SI3-10_EAP_rheology
v2 v3 23 23 24 24 Current EVP rheology assumes ice is isotropic, but observations show linear kinematic features (LKF) which are sign of anisotropic behaviour. 25 Elas itc 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]]25 Elastic 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]] 26 26 27 27 The IMMERSE project has a milestone and deliverable to implement EAP rheology in SI3. … … 171 171 === Test case === 172 172 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.173 We 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. 174 174 175 175 === Documentation updates … … 184 184 185 185 update: 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 reas sed later.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 reasessed later. 187 187 188 188 == Tests … … 192 192 }}} 193 193 194 ''...'' 194 The 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 196 Current 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 198 Convergence of the velocities can be checked with code from changeset 12741 which is not yet incorporated in this branch. 195 199 196 200 == Review