| 36 | |
| 37 | The restart is now working for the basic 3D kernel. The problem was that we were writing output after calling rhs3d. This overwrites the values of ru and rv. We moved output above of rhs3d and created a new routine set_zeta to assign the free-surface to its time averaged value. This affects the adjoint, so we are supporting both options for now. Many thanks to John Warner for his persistance. |
| 38 | |
| 39 | We now get zero error in the 3D upwelling test case for either ANA_VMIX or GLS_MIXING: |
| 40 | |
| 41 | Variable RMS Correlation |
| 42 | |
| 43 | zeta 0.00000000e+00 1.00000000e+00 |
| 44 | ubar 0.00000000e+00 1.00000000e-00 |
| 45 | vbar 0.00000000e+00 1.00000000e+00 |
| 46 | u 0.00000000e+00 1.00000000e+00 |
| 47 | v 0.00000000e+00 1.00000000e+00 |
| 48 | temp 0.00000000e+00 1.00000000e-00 |
| 49 | salt 0.00000000e+00 1.00000000e-00 |
| 50 | |
| 51 | We are still working on the perfect restart for sediment, biology, and other algorithms |