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

Timestamp:

2015-10-05T12:44:05+02:00 (9 years ago)

Author:

mathiot

Comment:

--

ticket/1611_is_oce_cpl

@@ -20 +20 @@
 
 ==== Method ====
-notation: bt = barotropic transport, if not precise, U/V/T/S mean properties of the top wet cell.
+ notation: bt = barotropic transport, if not precise, U/V/T/S mean properties of the top wet cell.
  * 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)
  * '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
@@ -27 +27 @@
  * '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
 
-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).
+ 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).
 
  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.  
@@ -40 +40 @@
  * divcur.F90: apply volume correction
  * trasbc.F90: apply heat/aslt correction
+ * lib_fortran.F90: at this stage, need all of these to diagnose all is OK. HAVE to be removed before or during review process. 
+ * 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)
+
 
 ----