Skip to content

Commit 1a3f760

Browse files
vanroekelvanroekel
authored
Merge pull request #23 from BrodiePearson/ocean/addADCMixing
Add Stokes terms to ADC scheme
2 parents c546898 + f4c494d commit 1a3f760

File tree

4 files changed

+271
-74
lines changed

4 files changed

+271
-74
lines changed

‎src/core_ocean/adc_mixing/Registry_adc_mixing_fields_opts.xml‎

Lines changed: 68 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -19,6 +19,22 @@
1919
description="number of decimal places to keep, inverse is position of decimal rounding"
2020
possible_values="large positive reals"
2121
/>
22+
<nml_option name="config_adc_langmuir" type="logical" default_value=".false." units="unitless"
23+
description="Switch to turn on wave effects (Langmuir turbulence) in ADC scheme. If false the Stokes drift of waves is set to zero. "
24+
possible_values=".true. or .false."
25+
/>
26+
<nml_option name="config_adc_La_t" type="real" default_value="0.3" units="unitless"
27+
description="Turbulent Langmuir number La_t = SQRT(friction velocity / surface Stokes drift magnitude). Default value is that of an 'equilibrium wind-sea'."
28+
possible_values="positive real numbers, increasing values indicate weaker waves"
29+
/>
30+
<nml_option name="config_adc_stokes_depth" type="real" default_value="4.8" units="m"
31+
description="Exponential decay depth of Stokes drift profile. Default value is common in LES studies."
32+
possible_values="positive real numbers"
33+
/>
34+
<nml_option name="config_adc_stokes_angle" type="real" default_value="0.0" units="radians"
35+
description="Angle between wind and waves (surface stress and surface Stokes drift vector respectively). Positive values denote waves to the left of the wind."
36+
possible_values="positive or negative real numbers"
37+
/>
2238
<nml_option name="config_adc_tau_o" type="real" default_value="1800" units="s^{-1}"
2339
description="characterstic eddy turnover timescale"
2440
possible_values="positive real numbers"
@@ -63,13 +79,33 @@
6379
description="alpha_2 third (vorticity) rapid pressure coefficient"
6480
possible_values="values near 0.6"
6581
/>
82+
<nml_option name="config_adc_alpha_st_0" type="real" default_value="2.2" units="unitless"
83+
description="alpha_st_0 first (TKE) Stokes pressure coefficient for second-moment budgets"
84+
possible_values="values near 2.2 (Pearson 2019)"
85+
/>
86+
<nml_option name="config_adc_alpha_st_1" type="real" default_value="1.4" units="unitless"
87+
description="alpha_st_1 second (strain) Stokes pressure coefficient"
88+
possible_values="values near 1.4 (Pearson 2019)"
89+
/>
90+
<nml_option name="config_adc_alpha_st_2" type="real" default_value="-0.5" units="unitless"
91+
description="alpha_st_2 third (vorticity) Stokes pressure coefficient"
92+
possible_values="values near -0.5 (Pearson 2019)"
93+
/>
6694
<nml_option name="config_adc_alpha_tracer1" type="real" default_value="0.6" units="unitless"
6795
description="alpha_scalar1 first rapid pressure coefficient in tracer flux budgets (Mironov 2001)"
68-
possible_values="suggested values: 0.6 (Canuto 2001) or equal to alpha_scalar2, >0, <1 (Harcourt 2013)"
96+
possible_values="suggested values: 0.6 (Canuto 2001) or 0.7 (Harcourt 2013) or positive and equal to alpha_tracer2 (production cancellation)"
6997
/>
7098
<nml_option name="config_adc_alpha_tracer2" type="real" default_value="1.0" units="unitless"
7199
description="alpha_scalar2 rapid pressure coefficient in tracer flux budgets (Mironov 2001)"
72-
possible_values="suggested values: 1.0 (Canuto 2001) or equal to alpha_scalar2, >0, <1 (Harcourt 2013)"
100+
possible_values="suggested values: 1.0 (Canuto 2001) or 0.7 (Harcourt 2013) or positive and equal to alpha_tracer1 (production cancellation)"
101+
/>
102+
<nml_option name="config_adc_alpha_st_tracer1" type="real" default_value="0.6" units="unitless"
103+
description="alpha_st_scalar1 first Stokes pressure coefficient in tracer flux budgets (analogous to Mironov 2001 shear terms)"
104+
possible_values="suggested values: 0.6 (Canuto analogy modified for Stokes production of wscalar) or 0.7 (Harcourt 2013) or POSITIVE and equal to MINUS alpha_tracer2 (production cancellation)"
105+
/>
106+
<nml_option name="config_adc_alpha_st_tracer2" type="real" default_value="-1.0" units="unitless"
107+
description="alpha_st_scalar2 second Stokes pressure coefficient in tracer flux budgets (analogous to Mironov 2001 shear terms)"
108+
possible_values="suggested values: -1.0 (Canuto analogy modified for Stokes production of wscalar) or 0.7 (Harcourt 2013) or NEGATIVE and equal to MINUS alpha_tracer1 (production cancellation)"
73109
/>
74110
<nml_option name="config_adc_c11" type="real" default_value="0.1" units="unitless"
75111
description="c_11 is the buoyancy contribution in the pressure correlation in w3 eqn"
@@ -292,24 +328,30 @@
292328
<var name="w2tend6" type="real" dimensions="nVertLevelsP1 nCells Time" units="m^2/s^3"
293329
description="conversion of w2 to u2/v2 due to splat effects"
294330
/>
295-
<var name="w3tend1" type="real" dimensions="nVertLevels nCells Time" units="m^2/s^3"
331+
<var name="w2tend7" type="real" dimensions="nVertLevelsP1 nCells Time" units="m^2/s^3"
332+
description="Stokes production of w2 [-2*[uw*d(ustokes)/dz + vw*d(vstokes)/dz]"
333+
/>
334+
<var name="w3tend1" type="real" dimensions="nVertLevels nCells Time" units="m^3/s^4"
296335
description="entrainment/detrainment term -- similar to dissipation"
297336
/>
298-
<var name="w3tend2" type="real" dimensions="nVertLevels nCells Time" units="m^2/s^3"
337+
<var name="w3tend2" type="real" dimensions="nVertLevels nCells Time" units="m^3/s^4"
299338
description="turbulent transport of w3 [ -d(w^4)/dz ]"
300339
/>
301-
<var name="w3tend3" type="real" dimensions="nVertLevels nCells Time" units="m^2/s^3"
340+
<var name="w3tend3" type="real" dimensions="nVertLevels nCells Time" units="m^3/s^4"
302341
description="gradient production of w3 [ (3/2)*d( [w2]^2 )/dz ]"
303342
/>
304-
<var name="w3tend4" type="real" dimensions="nVertLevels nCells Time" units="m^2/s^3"
343+
<var name="w3tend4" type="real" dimensions="nVertLevels nCells Time" units="m^3/s^4"
305344
description="parameterized pressure terms (Rotta, buoyancy)"
306345
/>
307-
<var name="w3tend5" type="real" dimensions="nVertLevels nCells Time" units="m^2/s^3"
346+
<var name="w3tend5" type="real" dimensions="nVertLevels nCells Time" units="m^3/s^4"
308347
description="buoyancy production and sub-plume scale effects"
309348
/>
310-
<var name="w3tend6" type="real" dimensions="nVertLevels nCells Time" units="m^2/s^3"
349+
<var name="w3tend6" type="real" dimensions="nVertLevels nCells Time" units="m^3/s^4"
311350
description="destruction due to splat effects"
312351
/>
352+
<var name="w3tend7" type="real" dimensions="nVertLevels nCells Time" units="m^3/s^4"
353+
description="Stokes production of w3 -3*[uw2*d(ustokes)/dz + vw2*d(ustokes)/dz]"
354+
/>
313355
<var name="wttend1" type="real" dimensions="nVertLevelsP1 nCells Time" units="m C/s^2"
314356
description="entrainment/detrainment term -- similar to dissipation"
315357
/>
@@ -328,6 +370,9 @@
328370
<var name="wttend6" type="real" dimensions="nVertLevelsP1 nCells Time" units="m C/s^2"
329371
description="subplume scale production term"
330372
/>
373+
<var name="wttend7" type="real" dimensions="nVertLevelsP1 nCells Time" units="m C/s^2"
374+
description="Stokes production of wt -[ut*d(ustokes)/dz + vt*d(vstokes)/dz]"
375+
/>
331376
<var name="wstend1" type="real" dimensions="nVertLevelsP1 nCells Time" units="m PSU/s^2"
332377
description="entrainment/detrainment term -- similar to dissipation"
333378
/>
@@ -346,6 +391,9 @@
346391
<var name="wstend6" type="real" dimensions="nVertLevelsP1 nCells Time" units="m PSU/s^2"
347392
description="subplume scale production term"
348393
/>
394+
<var name="wstend7" type="real" dimensions="nVertLevelsP1 nCells Time" units="m PSU/s^2"
395+
description="Stokes production of ws -[us*d(ustokes)/dz + vs*d(vstokes)/dz]"
396+
/>
349397
<var name="uwtend1" type="real" dimensions="nVertLevelsP1 nCells Time" units="m^2/s^3"
350398
description="turbulent transport of uw [ -d(uww)/dz ]"
351399
/>
@@ -361,6 +409,9 @@
361409
<var name="uwtend5" type="real" dimensions="nVertLevelsP1 nCells Time" units="m^2/s^3"
362410
description="dissipation of uw [-kappa_FL*[d(uw)/dz]^2] "
363411
/>
412+
<var name="uwtend6" type="real" dimensions="nVertLevelsP1 nCells Time" units="m^2/s^3"
413+
description="Stokes production of uw -[u2*d(ustokes)/dz + uv*d(vstokes)/dz]"
414+
/>
364415
<var name="vwtend1" type="real" dimensions="nVertLevelsP1 nCells Time" units="m^2/s^3"
365416
description="turbulent transport of vw [ -d(vww)/dz ]"
366417
/>
@@ -376,6 +427,9 @@
376427
<var name="vwtend5" type="real" dimensions="nVertLevelsP1 nCells Time" units="m^2/s^3"
377428
description="dissipation of uw [-kappa_FL*[d(uw)/dz]^2]"
378429
/>
430+
<var name="vwtend6" type="real" dimensions="nVertLevelsP1 nCells Time" units="m^2/s^3"
431+
description="Stokes production of vw -[uv*d(ustokes)/dz + v2*d(vstokes)/dz]"
432+
/>
379433
<var name="u2tend1" type="real" dimensions="nVertLevelsP1 nCells Time" units="m^2/s^3"
380434
description="turbulent transport of u2 [ -d(uuw)/dz ]"
381435
/>
@@ -586,6 +640,7 @@
586640
<var name="w2tend4"/>
587641
<var name="w2tend5"/>
588642
<var name="w2tend6"/>
643+
<var name="w2tend7"/>
589644
<var name="u2tend1"/>
590645
<var name="u2tend2"/>
591646
<var name="u2tend3"/>
@@ -601,26 +656,31 @@
601656
<var name="uwtend3"/>
602657
<var name="uwtend4"/>
603658
<var name="uwtend5"/>
659+
<var name="uwtend6"/>
604660
<var name="vwtend1"/>
605661
<var name="vwtend2"/>
606662
<var name="vwtend3"/>
607663
<var name="vwtend4"/>
608664
<var name="vwtend5"/>
665+
<var name="vwtend6"/>
609666
<var name="w3tend1"/>
610667
<var name="w3tend2"/>
611668
<var name="w3tend3"/>
612669
<var name="w3tend4"/>
613670
<var name="w3tend5"/>
671+
<var name="w3tend6"/>
614672
<var name="wttend1"/>
615673
<var name="wttend2"/>
616674
<var name="wttend3"/>
617675
<var name="wttend4"/>
618676
<var name="wttend5"/>
677+
<var name="wttend6"/>
619678
<var name="wstend1"/>
620679
<var name="wstend2"/>
621680
<var name="wstend3"/>
622681
<var name="wstend4"/>
623682
<var name="wstend5"/>
683+
<var name="wstend6"/>
624684
<var name="u2"/>
625685
<var name="v2"/>
626686
<var name="w2"/>

0 commit comments

Comments
 (0)