Changes between Version 3 and Version 4 of ticket/1299

Timestamp:

2014-10-29T11:46:13+01:00 (9 years ago)

Author:

jamesharle

Comment:

--

ticket/1299

@@ -15 +15 @@
 Implement an on-line vertical interpolation of BDY inputs needed to allow easy switching between vertical coordinate systems without the need to alter input boundary files.
 
-Files affected
-:
+Files affected:
+
 * fldread.F90
 * bdydta.F90
@@ -24 +24 @@
 * iom.F90
 
+Testing carried out on the AMM12 configuration in the trunk code. Six experiments were carried out to check the consistency of the code. These were:
 
+# Constant T/S boundary conditions on original vertical grid, with nb_jpk_bdy = 51 (first sanity test)
+# Constant T/S boundary conditions on staggered vertical grid, with nb_jpk_bdy = 51 (gdept(jpk = [1 51] ) same as native grid all depths in-between offset from native grid, second sanity test)
+# Realistic boundary conditions on original vertical grid, with nb_jpk_bdy = -1 (i.e. should be the identical to trunk code)
+# Realistic boundary conditions on original vertical grid, with nb_jpk_bdy = 51 (i.e. should be the identical to trunk code within rounding errors in interpolation routine)
+# Realistic boundary conditions on staggered vertical grid, with nb_jpk_bdy = 51 (gdept(jpk = [1 51] ) same as native grid all depths in-between offset from native grid)
+# Realistic boundary conditions on a reduced z-level staggered vertical grid, with nb_jpk_bdy = 31
+
+The datasets for the offset and reduced level input data where generated using the interp routine in matlab.
+
+Results from the tests:
+
+Outstanding issues:
+
+How best to handle the baroclinic velocities in the case of a bathymetry mismatch between input data and native grid:
+# Set velocities to zero in depths on native grid deeper than bathymetry in input file
+# Down fill velocities for missing data
+
+in all cases a barotropic correction has to applied. At present a conservative remapping option does not exist. This should be a straight forward addition of a subroutine in fldread.F90
 
 ----