| 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 | {{{ |
| 42 | Variable RMS Correlation |
| 43 | |
| 44 | zeta 0.00000000e+00 1.00000000e+00 |
| 45 | ubar 0.00000000e+00 1.00000000e-00 |
| 46 | vbar 0.00000000e+00 1.00000000e+00 |
| 47 | u 0.00000000e+00 1.00000000e+00 |
| 48 | v 0.00000000e+00 1.00000000e+00 |
| 49 | temp 0.00000000e+00 1.00000000e-00 |
| 50 | salt 0.00000000e+00 1.00000000e-00 |
| 51 | }}} |
| 52 | |
| 53 | We are still working on the perfect restart for sediment, biology, and other algorithms. |