1 |
Seaice-only verification experiment in idealized periodic channel |
2 |
----------------------------------------------------------------- |
3 |
|
4 |
1) main forward experiment (code, input) |
5 |
|
6 |
Re-entrant zonally periodic channel (80x42 grid points) with just level (Nr=1) |
7 |
uniform resolution (5.km, 10m), solid Southern boundary with triangular shape |
8 |
coastline ("bathy_3c.bin") |
9 |
|
10 |
Use seaice (dynamics & thermodynamics from pkg/thsice) with EXF (see data.pkg) |
11 |
with initial ice thickness of 0.2 m (but no snow) |
12 |
(thSIceThick_InitFile='const+20.bin', in "input/data.ice") |
13 |
Initial seaice concentration is 100 % everywhere |
14 |
(thSIceFract_InitFile='const100.bin', in "input/data.ice") |
15 |
and seaice is initially at rest. |
16 |
|
17 |
At runtime turn off time-stepping in 'data', PARM01, using: |
18 |
momStepping = .FALSE., |
19 |
saltStepping = .FALSE., |
20 |
tempAdvection=.FALSE., |
21 |
And just keep surface temp relaxation (tauRelax = 1 month) toward fixed SST: |
22 |
in data.exf : |
23 |
> climsstperiod = 0.0, |
24 |
> climsstTauRelax = 2592000., |
25 |
> climsstfile = 'tocn.bin', |
26 |
|
27 |
Forcing: |
28 |
None of the forcing vary with time; Most of the input files have been |
29 |
generated using matlab script "input/gendata.m". |
30 |
SST relaxation field is uniform in X, parabolic function of Y with |
31 |
maximum close to Southern boundary. |
32 |
|
33 |
Atmospheric air temp is uniform in Y, and only vary with X (~sin(2.pi.x/Lx)) |
34 |
with an amplitude of 4.K ('tair_4x.bin'); |
35 |
Uses constant Relative Humidity (70%, file 'qa70_4x.bin') |
36 |
constant and uniform downward shortwave (100.W/m2, 'dsw_100.bin'), |
37 |
downward longwave (250.W/m^2, 'dlw_250.bin'), |
38 |
zonal wind (10.m/s, 'windx.bin'), |
39 |
no meridional wind, no precip. |
40 |
|
41 |
Ocean surface currents comes from a 3 levels ocean-only run (without seaice) |
42 |
using the same wind forcing (uVel_3c0.bin, vVel_3c0.bin) (matlab script: |
43 |
"input/getdata.m") |
44 |
|
45 |
Two matlab scripts ("input/grph_res.m" and "input/grph_diag.m") |
46 |
are provided to make some basic plots (snap-shot and time-averaged |
47 |
diagnostics). |
48 |
|
49 |
2) other (secondary) experiments (using the same executable) |
50 |
|
51 |
a) seaice-dynamics only with LSR solver (input.dyn_lsr/) |
52 |
use same forcing as main forward experiment, |
53 |
( link input files from "input.dyn_lsr/" dir 1rst, then from "input/" dir ) |
54 |
but does not use at all pkg/thsice (advect pkg/seaice properties using |
55 |
pkg/seaice advection S/R.); turn off seaice thermodynamics |
56 |
(usePW79thermodynamics=.FALSE., in input.dyn_lsr/data.seaice). |
57 |
|
58 |
b) seaice-dynamics only with JNFK solver (input.dyn_jfnk/) |
59 |
use same forcing as main forward experiment, |
60 |
( link input files from "input.dyn_jfnk/" dir 1rst, then from "input/" dir ) |
61 |
turn off seaice thermodynamics (thSIce_skipThermo=.TRUE., in file |
62 |
"input.dyn_jfnk/data.ice") and advect pkg/thsice properties using |
63 |
pkg/thsice advection S/R. |
64 |
|
65 |
c) seaice-thermodynamics only from pkg/seaice (input.thermo/) |
66 |
use same forcing as main forward experiment, |
67 |
( link input files from "input.thermo/" dir 1rst, then from "input/" dir ) |
68 |
Except a different bathymetry (no triangular shape coast line in the South, |
69 |
just a simple channel of constant width, with constant and uniform |
70 |
zonal current of 0.2 m/s): |
71 |
in "input.thermo/data": |
72 |
bathyFile = 'channel.bin', |
73 |
uVelInitFile = 'const+20.bin', |
74 |
vVelInitFile = 'const_00.bin', |
75 |
And use different initial seaice concentration from the main forward experiment: |
76 |
(in "input.thermo/data.seaice", AreaFile = 'ice0_area.bin', |
77 |
and HeffFile = 'ice0_heff.bin',) |
78 |
100 % everywhere except in the 1rst and last 2 rows: 0., 10.% in the South |
79 |
and 1.%, 0. in the North: |
80 |
Turn off seaice dynamics (SEAICEuseDYNAMICS =.FALSE., in |
81 |
"input.thermo/data.seaice") which implies no advection (zero ice velocity). |
82 |
|
83 |
d) seaice-thermodynamics only from pkg/thsice (input.thsice/) |
84 |
use same forcing as main forward experiment, |
85 |
( link input files from "input.thsice/" dir 1rst, then from "input/" dir ) |
86 |
Except for bathymetry and initial seaice concentration which are identical |
87 |
to the one used in input.thermo/ : |
88 |
(thSIceFract_InitFile='ice0_area.bin', in "input.thsice/data.ice") |
89 |
No seaice dynamics and no seaice advection, so that none of pkg/seaice S/R |
90 |
are used (useSEAICE = .TRUE., commented out in "input.thsice/data.pkg"), |
91 |
|
92 |
3) main adjoint experiment (code_ad, input_ad) |
93 |
Use similar forward set-up as "input.thermo/". |
94 |
(to be completed ...) |