Warning, /verification/offline_exf_seaice/README is written in an unsupported language. File is not indexed.
view on githubraw file Latest commit 09461bf9 on 2022-04-12 23:30:40 UTC
6ee3c45c30 Jean*0001 Seaice-only verification experiment in idealized periodic channel
0002 -----------------------------------------------------------------
0003
0004 1) main forward experiment (code, input)
0005
0006 Re-entrant zonally periodic channel (80x42 grid points) with just level (Nr=1)
0007 uniform resolution (5.km, 10m), solid Southern boundary with triangular shape
0008 coastline ("bathy_3c.bin")
0009
0010 Use seaice (dynamics & thermodynamics from pkg/thsice) with EXF (see data.pkg)
0011 with initial ice thickness of 0.2 m (but no snow)
0012 (thSIceThick_InitFile='const+20.bin', in "input/data.ice")
0013 Initial seaice concentration is 100 % everywhere
0014 (thSIceFract_InitFile='const100.bin', in "input/data.ice")
0015 and seaice is initially at rest.
0016
0017 At runtime turn off time-stepping in 'data', PARM01, using:
0018 momStepping = .FALSE.,
0019 saltStepping = .FALSE.,
0020 tempAdvection=.FALSE.,
0021 And just keep surface temp relaxation (tauRelax = 1 month) toward fixed SST:
0022 in data.exf :
0023 > climsstperiod = 0.0,
0024 > climsstTauRelax = 2592000.,
0025 > climsstfile = 'tocn.bin',
0026
0027 Forcing:
0028 None of the forcing vary with time; Most of the input files have been
0029 generated using matlab script "input/gendata.m".
0030 SST relaxation field is uniform in X, parabolic function of Y with
0031 maximum close to Southern boundary.
0032
0033 Atmospheric air temp is uniform in Y, and only vary with X (~sin(2.pi.x/Lx))
0034 with an amplitude of 4.K ('tair_4x.bin');
0035 Uses constant Relative Humidity (70%, file 'qa70_4x.bin')
0036 constant and uniform downward shortwave (100.W/m2, 'dsw_100.bin'),
0037 downward longwave (250.W/m^2, 'dlw_250.bin'),
0038 zonal wind (10.m/s, 'windx.bin'),
0039 no meridional wind, no precip.
0040
0041 Ocean surface currents comes from a 3 levels ocean-only run (without seaice)
0042 using the same wind forcing (uVel_3c0.bin, vVel_3c0.bin) (matlab script:
0043 "input/getdata.m")
0044
0045 Two matlab scripts ("input/grph_res.m" and "input/grph_diag.m")
0046 are provided to make some basic plots (snap-shot and time-averaged
0047 diagnostics).
0048
0049 2) other (secondary) experiments (using the same executable)
0050
0051 a) seaice-dynamics only with LSR solver (input.dyn_lsr/)
0052 use same forcing as main forward experiment,
0053 ( link input files from "input.dyn_lsr/" dir 1rst, then from "input/" dir )
0054 but does not use at all pkg/thsice (advect pkg/seaice properties using
0055 pkg/seaice advection S/R.); turn off seaice thermodynamics
0056 (usePW79thermodynamics=.FALSE., in input.dyn_lsr/data.seaice).
0057
0058 b) seaice-dynamics only with JNFK solver (input.dyn_jfnk/)
0059 use same forcing as main forward experiment,
0060 ( link input files from "input.dyn_jfnk/" dir 1rst, then from "input/" dir )
0061 turn off seaice thermodynamics (thSIce_skipThermo=.TRUE., in file
0062 "input.dyn_jfnk/data.ice") and advect pkg/thsice properties using
0063 pkg/thsice advection S/R.
0064
09461bf9e9 Jean*0065 c) seaice-dynamics only with LSR solver and elliptical yieldcurve with
0066 non-normal flow rule (input.dyn_ellnnfr)
0067 use same forcing as main forward experiment,
0068 ( link input files from "input.dyn_ellnnfr/" dir 1rst, then from
0069 "input.dyn_lsr/" and "input/" dir )
0070 advect pkg/seaice properties using pkg/seaice advection S/R.;
0071 turn off seaice thermodynamics (usePW79thermodynamics=.FALSE., in
0072 input.dyn_ellnnfr/data.seaice).
0073
0074 d) seaice-dynamics only with LSR solver and Mohr-Coulomb
0075 yieldcurve with elliptical plastic potential (input.dyn_mce)
0076 use same forcing as main forward experiment,
0077 ( link input files from "input.dyn_mce/" dir 1rst, then from
0078 "input.dyn_lsr/" and "input/" dir )
0079 advect pkg/seaice properties using pkg/seaice advection S/R.;
0080 turn off seaice thermodynamics (usePW79thermodynamics=.FALSE., in
0081 input.dyn_mce/data.seaice).
0082
0083 e) seaice-dynamics only with KRYLOV solver and parabolic lens yieldcurve
0084 (input.dyn_paralens)
0085 use same forcing as main forward experiment,
0086 ( link input files from "input.dyn_paralens/" dir 1rst, then from
0087 "input.dyn_jfnk/" and "input/" dir )
0088 turn off seaice thermodynamics (thSIce_skipThermo=.TRUE., in file
0089 "input.dyn_mce/data.ice" linked from "input.dyn_jfnk/") and
0090 advect pkg/thsice properties using pkg/thsice advection S/R.
0091
0092 f) seaice-dynamics only with JFNK solver and teardrop yieldcurve
0093 (input.dyn_teardrop)
0094 use same forcing as main forward experiment,
0095 ( link input files from "input.dyn_teardrop/" dir 1rst, then from
0096 "input.dyn_jfnk/" and "input/" dir )
0097 turn off seaice thermodynamics (thSIce_skipThermo=.TRUE., in file
0098 "input.dyn_mce/data.ice" linked from "input.dyn_jfnk/") and
0099 advect pkg/thsice properties using pkg/thsice advection S/R.
0100
0101 g) seaice-thermodynamics only from pkg/seaice (input.thermo/)
6ee3c45c30 Jean*0102 use same forcing as main forward experiment,
0103 ( link input files from "input.thermo/" dir 1rst, then from "input/" dir )
0104 Except a different bathymetry (no triangular shape coast line in the South,
0105 just a simple channel of constant width, with constant and uniform
0106 zonal current of 0.2 m/s):
0107 in "input.thermo/data":
0108 bathyFile = 'channel.bin',
0109 uVelInitFile = 'const+20.bin',
0110 vVelInitFile = 'const_00.bin',
0111 And use different initial seaice concentration from the main forward experiment:
0112 (in "input.thermo/data.seaice", AreaFile = 'ice0_area.bin',
0113 and HeffFile = 'ice0_heff.bin',)
0114 100 % everywhere except in the 1rst and last 2 rows: 0., 10.% in the South
0115 and 1.%, 0. in the North:
0116 Turn off seaice dynamics (SEAICEuseDYNAMICS =.FALSE., in
0117 "input.thermo/data.seaice") which implies no advection (zero ice velocity).
0118
09461bf9e9 Jean*0119 h) seaice-thermodynamics only from pkg/thsice (input.thsice/)
6ee3c45c30 Jean*0120 use same forcing as main forward experiment,
0121 ( link input files from "input.thsice/" dir 1rst, then from "input/" dir )
0122 Except for bathymetry and initial seaice concentration which are identical
0123 to the one used in input.thermo/ :
0124 (thSIceFract_InitFile='ice0_area.bin', in "input.thsice/data.ice")
0125 No seaice dynamics and no seaice advection, so that none of pkg/seaice S/R
0126 are used (useSEAICE = .TRUE., commented out in "input.thsice/data.pkg"),
0127
0128 3) main adjoint experiment (code_ad, input_ad)
09461bf9e9 Jean*0129 Use similar forward set-up as "input.thermo/",
0130 no seaice dynamics and no seaice advection.
0131
0132 4) other (secondary) adjoint experiment (using the same executable)
0133
0134 a) seaice-thermodynamics form pkg/thsice (input_ad.thsice)
0135 use similar forward set-up as "input.thsice",
0136 no seaice dynamics and no seaice advection.
