Warning, /verification/seaice_itd/README is written in an unsupported language. File is not indexed.
view on githubraw file Latest commit 0ba3967d on 2014-06-24 08:00:16 UTC
0ba3967dec Mart*0001 Seaice-only verification experiment in idealized periodic channel with
0002 ice thickness distribution (otherwise very similar to
0003 offline_exf_seaice): CPP-flag SEAICE_ITD is defined
0004 -----------------------------------------------------------------
0005
0006 1) main forward experiment (code, input)
0007
0008 Re-entrant zonally periodic channel (80x42 grid points) with just level (Nr=1)
0009 uniform resolution (5.km, 10m), solid Southern boundary with triangular shape
0010 coastline ("bathy_3c.bin")
0011
0012 Use seaice (dynamics & thermodynamics from pkg/seaice) with EXF (see data.pkg)
0013 with initial ice thickness ranging from nearly 0 m in the "south"
0014 to over 7 m in the "north"(but no snow)
0015 (HeffFile = 'heff_quartic.bin', in "input/data.seaice")
0016 Initial seaice concentration is 100 % everywhere
0017 (AreaFile='const100.bin', in "input/data.seaice")
0018 and seaice is initially at rest.
0019
0020 Ridging is computed according to Thorndyke et al (1975) and Hibler
0021 (1980). Ice strength P is computed following Rothrock (1975)
0022
0023 At runtime turn off time-stepping in 'data', PARM01, using:
0024 momStepping = .FALSE.,
0025 saltStepping = .FALSE.,
0026 tempAdvection=.FALSE.,
0027
0028 Forcing:
0029 None of the forcing vary with time; the input files have been
0030 generated using the python script "input/gendata.py".
0031 SST relaxation field is uniform in X, parabolic function of Y with
0032 maximum close to Southern boundary.
0033
0034 Atmospheric air temp is uniform in Y, and only vary with X (~sin(2.pi.x/Lx))
0035 with an amplitude of 4.K ('tair_4x.bin');
0036 Uses constant Relative Humidity (70%, file 'qa70_4x.bin')
0037 constant and uniform downward shortwave (100.W/m2, 'dsw_100.bin'),
0038 downward longwave (250.W/m^2, 'dlw_250.bin'),
0039 zonal wind (10.m/s, 'windx.bin'),
0040 no meridional wind, no precip.
0041
0042 Ocean surface currents comes from a 3 levels ocean-only run (without seaice)
0043 using the same wind forcing (uVel_3c0.bin, vVel_3c0.bin) (matlab script:
0044 "input/getdata.m")
0045
0046 2) other (secondary) experiments (with the same executable)
0047
0048 a) input.lipscomb07: seaice-dynamics only with ridging scheme
0049 following Lipscomb et al (2007); uses same forcing as main forward
0050 experiment, without thermodynamic forcing
0051 (usePW79thermodynamics=.FALSE., SEAICEpartFunc = 1, SEAICEredistFunc = 1,
0052 in input.dyn_lsr/data.seaice).
0053
0054 b) input.thermo: seaice-dynamics and thermodynamics with ridging, but
0055 ice strength is computed following Hibler (1979):
0056 P = 27.5e4*(h*c)*exp(-20*(1-c))
0057
