Changes between Version 20 and Version 21 of ticket/1611_is_oce_cpl


Ignore:
Timestamp:
2015-10-21T18:33:06+02:00 (5 years ago)
Author:
mathiot
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • ticket/1611_is_oce_cpl

    v20 v21  
    1313=== Description === 
    1414==== Purpose ==== 
    15 Allow NEMO to restart with a new geometry beneath the ice shelf (grounding line migration, i.e advance or retreat as well as calving front migration). This branch allow to test ice sheet/ocean coupling at each restart time step through netcdf file exchange.  
     15Allow NEMO to restart with a new geometry beneath the ice shelf (grounding line migration, i.e advance or retreat as well as calving front migration). This branch allows to test ice sheet/ocean coupling at each restart time step through netcdf file exchange.  
    1616 
    1717==== Limitation ==== 
    1818 * It is not suitable if high frequency coupling is needed. Because the model need to be stop each time you want to couple ocean/ice sheet. 
    19  * It not a conservative method (NEMO is creating/losing volume/mass/salt/heat at the coupling time step if the geometry change). I.e., under a pure dynamics ice sheet advance, your mean ssh will stay constant. 
     19 * It not a conservative method (NEMO is creating/losing volume/mass/salt/heat at the coupling time step if the geometry change). I.e., under a pure dynamics ice sheet advance, your mean ssh will stay constant. An option has been added to reduce the unbalance to 0 after a certain amount of time step (ie we apply the total correction over X time step).  
    2020 
    2121==== Method ==== 
     
    4444 * domngb.F90: I add an optional argument to be able to do this beneath the ice shelf (i.e. at level jk instead of 1). 
    4545 * domain.F90: Call to domwri modified, because if coupling, mask could be change in iscpl routine (closed pool or other weird configuration due to the extrapolation process. (3 solutions: let as it is now, simplify it and move it to the end of istate.F90 or build a proper flood filling algorithm to apply in domzgr and keep the call to dom_wri as it is). 
    46  * domzgr.F90: modification of zgr_isf, because in case of coupling, we do not dig the bathymetry to ensure at least 2 wet cells in the water column. (if no coupling, isf draft and bathy could be dig to ensure this condition). '''NOT YET CLEAN''' 
     46 * domzgr.F90: modification of zgr_isf, because in case of coupling, we do not want to dig the bathymetry to ensure at least 2 wet cells in the water column. (if no coupling, isf draft and bathy could be dig to ensure this condition). 
    4747 * + small modifications in order to have dev_r5589_is_oce_cpl branch compatible to dev_r5151_UKMO_ISF. Could leads to some issue during the merge of these 2 branches (modification of umask_i by ssumask ...). 
    4848 
     
    7272 
    7373=== Bit Comparability === 
    74 ||Does this change preserve answers in your tested standard configurations (to the last bit) ?||'''YES'''|| 
    75 ||Does this change bit compare across various processor configurations. (1xM, Nx1 and MxN are recommended)||!'''YES/NO!'''|| 
     74||Does this change preserve answers in your tested standard configurations (to the last bit) ?||'''YES (ISOMIP+ EXP4'''|| 
     75||Does this change bit compare across various processor configurations. (1xM, Nx1 and MxN are recommended)||'''YES (8x4 and 4x8)'''|| 
    7676||Is this change expected to preserve answers in all possible model configurations?||'''YES'''|| 
    7777||Is this change expected to preserve all diagnostics? [[BR]]!,,!''Preserving answers in model runs does not necessarily imply preserved diagnostics. !''||'''YES'''||