4 |
Notes on tags used in MITgcmUV |
Notes on tags used in MITgcmUV |
5 |
============================== |
============================== |
6 |
|
|
7 |
|
o pkg/exf/exf_radiation.F: add the factor ocean_emissivity to lwdown, for |
8 |
|
within #ifdef EXF_LWDOWN_WITH_EMISSIVITY/#endif. Unless ocean_emissivity=1. |
9 |
|
this factor is absolutely required, because the lw exitance (= out-going |
10 |
|
long wave radiation) is emissivity*stefanBoltzmann*T^4 + rho*lwdown, |
11 |
|
where the reflectivity rho = 1-emissivity for conservation reasons. Once |
12 |
|
this is established, we can remove the CPP-Flag (if we want to). |
13 |
|
o pkg/seaice: |
14 |
|
- qualify a stop statement in S/R seaice_check: EXF_READ_EVAP and |
15 |
|
ALLOW_BULKFORMULAE are only required if SEAICE_EXTERNAL_FLUXES is defined |
16 |
|
- sort out emissivity and long wave radiation, step 1: make SEAICE_emissivity |
17 |
|
what the name implies (and not emissivity*BoltzmannConstant), this means |
18 |
|
the old default value of 5.5e-8 now means |
19 |
|
5.5e-8/5.670e-8 = 0.97001763668430343479 as ocean_emissivity in |
20 |
|
S/R exf_readparms |
21 |
|
- FV-discretization of the B-grid LSOR solver, turn on with |
22 |
|
SEAICE_LSRBNEW. This new solver also involves a slight re-organization |
23 |
|
of the calling routine dynsolver, so that B and C-grid code look a |
24 |
|
little more alike. |
25 |
|
- fix index bugs in ostres.F (B-grid code) |
26 |
|
o pkg/seaice/seaice_ocean_budget.F: change (unused) code to make it look better |
27 |
o pkg/seaice/seaice_readparms.F: make ocean_emissivity*stefanBoltzmann the |
o pkg/seaice/seaice_readparms.F: make ocean_emissivity*stefanBoltzmann the |
28 |
default for SEAICE_emissivity, if useEXF; requires setting a wrong |
default for SEAICE_emissivity, if useEXF; requires setting a wrong |
29 |
parameter value in global_ocean.cs32x15/input.seaice/data.seaice for |
parameter value in global_ocean.cs32x15/input.seaice/data.seaice for |