| 5 |
============================== |
============================== |
| 6 |
|
|
| 7 |
o pkg/seaice: |
o pkg/seaice: |
| 8 |
|
- replace ZMIN=4.e8 with ZMIN=0. in dynsolver.F as in C-grid code |
| 9 |
|
- fix bug in metric terms of u-equation in LSR solver (implicit part), |
| 10 |
|
changes verification/lab_sea |
| 11 |
|
- rearrange implicit part of u-equation in LSR solver to make it easier to |
| 12 |
|
read, understand, and debug. This causes changes at the truncation level. |
| 13 |
|
- fix (serious) bug in the boundary conditions (sign error in the implicit |
| 14 |
|
part) of the LSR solver. Changes all LSR results. |
| 15 |
|
- fixed problem with HSNOW initialization to avoid thick snow when AREA |
| 16 |
|
is initilized with small values - does not change lab_sea results. |
| 17 |
|
- added HsnowFile initialization capability |
| 18 |
|
o pkg/smooth: addition to model/src and pkg/ctrl so allow |
| 19 |
|
initialization of pkg/smooth and application to control vector |
| 20 |
|
o pkg/seaice: fix another bug in the LSR no-slip boundary conditons |
| 21 |
|
o dependencies: pkg/profiles requires pkg/cal |
| 22 |
|
o pkg/profiles: |
| 23 |
|
- added option ALLOW_PROFILES_GENERICGRID |
| 24 |
|
-> pkg/profiles can thus be used with any grid (e.g. cs32x15) assuming |
| 25 |
|
that interpolation information (grid points and coefficients) is provided |
| 26 |
|
within the netcdf input files. |
| 27 |
|
-> sample matlab scripts can be found in MITgcm_contrib/gael |
| 28 |
|
that allow you to prepare netcdf input files. |
| 29 |
|
- also modified: cleaner include statements, cleaner error print statements, |
| 30 |
|
cleaner stops when error, and more tests of the netcdf input files consistency. |
| 31 |
|
|
| 32 |
|
o additional test (in global_ocean.cs32x15) using seaice+thsice pkgs. |
| 33 |
|
o pkg/seaice: |
| 34 |
|
- move seaice-diagnostics from seaice_do_diags to do_statevars_diags; this |
| 35 |
|
requires a new routine: seaice_diagnostics_fill.F |
| 36 |
|
- make sure that seaice does not move if not forced (pressure replacement |
| 37 |
|
method, Hibler and Ip, 1995), changes results |
| 38 |
- and new damping criterion for EVP solver according to Hunke, JCP, 2001 |
- and new damping criterion for EVP solver according to Hunke, JCP, 2001 |
| 39 |
- SEAICE_evpDampC = > 0. turns it on (615. is the recommended value). |
- SEAICE_evpDampC = > 0. turns it on (615. is the recommended value). |
| 40 |
It is turned off by default (=-1.), because it changes the results. |
It is turned off by default (=-1.), because it changes the results. |
Parent Directory
| 
Revision Log
| 
Revision Graph
| 
Patch