Changes between Version 3 and Version 4 of ticket/1611_is_oce_cpl


Ignore:
Timestamp:
2015-10-05T12:44:05+02:00 (5 years ago)
Author:
mathiot
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • ticket/1611_is_oce_cpl

    v3 v4  
    2020 
    2121==== Method ==== 
    22 notation: bt = barotropic transport, if not precise, U/V/T/S mean properties of the top wet cell. 
     22 notation: bt = barotropic transport, if not precise, U/V/T/S mean properties of the top wet cell. 
    2323 * Thin/enlarge a top cell : T/S/SSH are unchanged, U/V in the top cell are corrected to keep bt_b = bt_n (usefull to improve stability) 
    2424 * 'Dry a cell' case: set mask = 0, T/S = 0, U/V = 0, bt is corrected to keep bt_b = bt_n and ssh_b = ssh_n 
     
    2727 * 'Wet a column' case: set mask to 1, T/S = MEAN(T(ji+1,jj), T(ji-1,jj), T(ji,jj+1), T(ji,jj-1)) (if no neighbour along i/j axis check along the k axis), SSH unchanged, U/V = 0 
    2828 
    29 As the before and now fields are not really compatible (modification of the geometry), the restart time step is prescribed to be an euler time step instead of a leap frog and fields_b = fields_n (only if coupling between ice sheet and ocean). 
     29 As the before and now fields are not really compatible (modification of the geometry), the restart time step is prescribed to be an euler time step instead of a leap frog and fields_b = fields_n (only if coupling between ice sheet and ocean). 
    3030 
    3131 An option is available to keep trend due to unconservation of Vol./heat/salt to 0 using what is done for runoff with specified T/S. I.e the location and the amount of extra/loss vol./heat/salt is diagnosed, saved and used to come back to a acceptable level of conservation before the next restart time step. The vol./heat/salt is removed/added as close as possible to the source/sink location using the techniques used in sbcrnf to prescribed a source of volume/heat/salt in the runoff. This option is only suitable for global ocean modeling. It has only been tested in the ISOMIP+ 3 and 4 experiments during 5 coupling steps. Only 5, because afterward the sea level was going ridiculously high.   
     
    4040 * divcur.F90: apply volume correction 
    4141 * trasbc.F90: apply heat/aslt correction 
     42 * lib_fortran.F90: at this stage, need all of these to diagnose all is OK. HAVE to be removed before or during review process.  
     43 * lbclnk.F90: add a routine to add correction term over the halo and afterward set the halo to 0. (Need to set a flag because not done for north fold but as the north fold do nut cut the Greenland, it is OK) 
     44 
    4245 
    4346----