forked from ACCESS-NRI/access-om3-configs
-
Notifications
You must be signed in to change notification settings - Fork 0
/
MOM_input
1036 lines (868 loc) · 57.6 KB
/
MOM_input
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
/* WARNING: DO NOT EDIT this file. Any changes you make will be
overriden. To make changes in MOM6 parameters within CESM
framework, use SourceMods or user_nl_mom mechanisms.
This input file provides the adjustable run-time parameters
for version 6 of the Modular Ocean Model (MOM6). By default,
this file contains the out-of-the-box CESM configuration. A
full list of parameters for this case can be found in the
corresponding MOM_parameter_doc.all file which is generated
by the model at runtime. */
INPUTDIR = "./input/"
! "The directory in which input files are found."
TRIPOLAR_N = False
! "[Boolean] default = False
! Use tripolar connectivity at the northern edge of the
! domain. With TRIPOLAR_N, NIGLOBAL must be even."
NIHALO = 4
! "default = 2
! The number of halo points on each side in the
! x-direction. With STATIC_MEMORY_ this is set as NIHALO_
! in MOM_memory.h at compile time; without STATIC_MEMORY_
! the default is NIHALO_ in MOM_memory.h (if defined) or 2."
NJHALO = 4
! "default = 2
! The number of halo points on each side in the
! y-direction. With STATIC_MEMORY_ this is set as NJHALO_
! in MOM_memory.h at compile time; without STATIC_MEMORY_
! the default is NJHALO_ in MOM_memory.h (if defined) or 2."
NIGLOBAL = 360
! "The total number of thickness grid points in the
! x-direction in the physical domain. With STATIC_MEMORY_
! this is set in MOM_memory.h at compile time."
NJGLOBAL = 300
! "The total number of thickness grid points in the
! y-direction in the physical domain. With STATIC_MEMORY_
! this is set in MOM_memory.h at compile time.
! The number of processors in the x-direction. With
! STATIC_MEMORY_ this is set in MOM_memory.h at compile time.
! The number of processors in the x-direction. With
! STATIC_MEMORY_ this is set in MOM_memory.h at compile time.
! The processor layout that was acutally used."
IO_LAYOUT = 1, 1
! "default = 0
! The processor layout to be used, or 0,0 to automatically
! set the io_layout to be the same as the layout."
NK = 101
! "[nondim]
! The number of model layers."
USE_LEGACY_DIABATIC_DRIVER = False
! "[Boolean] default = True
! If true, use the a legacy version of the diabatic subroutine.
! This is temporary and is needed avoid change in answers."
DIABATIC_FIRST = True
! "[Boolean] default = False
! If true, apply diabatic and thermodynamic processes,
! including buoyancy forcing and mass gain or loss,
! before stepping the dynamics forward."
USE_REGRIDDING = True
! "[Boolean] default = False
! If True, use the ALE algorithm (regridding/remapping).
! If False, use the layered isopycnal algorithm."
THICKNESSDIFFUSE = True
! "[Boolean] default = False
! If true, interfaces or isopycnal surfaces are diffused,
! depending on the value of FULL_THICKNESSDIFFUSE."
THICKNESSDIFFUSE_FIRST = True
! "[Boolean] default = False
! If true, do thickness diffusion before dynamics.
! This is only used if THICKNESSDIFFUSE is true."
DT = 1800.0
! "[s]
! The (baroclinic) dynamics time step. The time-step that
! is actually used will be an integer fraction of the
! forcing time-step (DT_FORCING in ocean-only mode or the
! coupling timestep in coupled mode.)"
DT_THERM = 3600.0
! "[s] default = 3600.0
! The thermodynamic and tracer advection time step.
! Ideally DT_THERM should be an integer multiple of DT
! and less than the forcing or coupling time-step.
! By default DT_THERM is set to DT.
! The minimum amount of time in seconds between
! calculations of depth-space diagnostics. Making this
! larger than DT_THERM reduces the performance penalty
! of regridding to depth online."
HFREEZE = 10.0
! "[m] default = -1.0
! If HFREEZE > 0, melt potential will be computed. The actual depth
! over which melt potential is computed will be min(HFREEZE, OBLD),
! where OBLD is the boundary layer depth. If HFREEZE <= 0 (default),
! melt potential will not be computed."
DTBT_RESET_PERIOD = 0.0
! "[s] default = 7200.0
! The period between recalculations of DTBT (if DTBT <= 0).
! If DTBT_RESET_PERIOD is negative, DTBT is set based
! only on information available at initialization. If 0,
! DTBT will be set every dynamics time step. The default
! is set by DT_THERM. This is only used if SPLIT is true."
FRAZIL = True
! "[Boolean] default = False
! If true, water freezes if it gets too cold, and the
! the accumulated heat deficit is returned in the
! surface state. FRAZIL is only used if
! ENABLE_THERMODYNAMICS is true.
! If true, apply geothermal heating."
MIN_SALINITY = 0.0
! "[PPT] default = 0.01
! The minimum value of salinity when BOUND_SALINITY=True. The default is 0.01
! for backward compatibility but ideally should be 0."
BOUND_SALINITY = True
! "[Boolean] default = False
! If true, limit salinity to being positive. (The sea-ice
! model may ask for more salt than is available and
! drive the salinity negative otherwise.)"
C_P = 3992.0
! "[J kg-1 K-1] default = 3991.86795711963
! The heat capacity of sea water, approximated as a
! constant. This is only used if ENABLE_THERMODYNAMICS is
! true. The default value is from the TEOS-10 definition
! of conservative temperature."
USE_PSURF_IN_EOS = True
! "[Boolean] default = True
! If true, always include the surface pressure contributions in equation of
! state calculations."
CHECK_BAD_SURFACE_VALS = True
! "[Boolean] default = False
! If true, check the surface state for ridiculous values."
SAVE_INITIAL_CONDS = True
! "[Boolean] default = False
! If true, write the initial conditions to a file given
! by IC_OUTPUT_FILE."
GRID_ROTATION_ANGLE_BUGS = False
! "[Boolean] default = True
! If true, use an older algorithm to calculate the sine and
! cosines needed rotate between grid-oriented directions and
! true north and east. Differences arise at the tripolar fold."
DTFREEZE_DP = -7.75e-08
! "[deg C Pa-1] default = 0.0
! When TFREEZE_FORM=LINEAR,
! this is the derivative of the freezing potential
! temperature with pressure."
USE_IDEAL_AGE_TRACER = True
! "[Boolean] default = False
! If true, use the ideal_age_example tracer package."
USE_CFC_CAP = False
! "[Boolean] default = False
! If true, use the MOM_CFC_cap tracer package."
COORD_CONFIG = none
! "This specifies how layers are to be defined:
! file - read coordinate information from the file
! specified by (COORD_FILE).
! linear - linear based on interfaces not layesrs.
! ts_ref - use reference temperature and salinity
! ts_range - use range of temperature and salinity
! (T_REF and S_REF) to determine surface density
! and GINT calculate internal densities.
! gprime - use reference density (RHO_0) for surface
! density and GINT calculate internal densities.
! ts_profile - use temperature and salinity profiles
! (read from COORD_FILE) to set layer densities.
! USER - call a user modified routine.
! The file from which the coordinate densities are read."
REMAP_UV_USING_OLD_ALG = False
! "[Boolean] default = True
! If true, uses the old remapping-via-a-delta-z method for remapping u and v. If
! false, uses the new method that remaps between grids described by an old and
! new thickness."
GRID_CONFIG = mosaic
! "A character string that determines the method for
! defining the horizontal grid. Current options are:
! mosaic - read the grid from a mosaic (supergrid)
! file set by GRID_FILE.
! cartesian - use a (flat) Cartesian grid.
! spherical - use a simple spherical grid.
! mercator - use a Mercator spherical grid."
GRID_FILE = ocean_hgrid.nc
! "Name of the file from which to read horizontal grid data."
TOPO_CONFIG = file
! "This specifies how bathymetry is specified:
! file - read bathymetric information from the file
! specified by (TOPO_FILE).
! flat - flat bottom set to MAXIMUM_DEPTH.
! bowl - an analytically specified bowl-shaped basin
! ranging between MAXIMUM_DEPTH and MINIMUM_DEPTH.
! spoon - a similar shape to 'bowl', but with an vertical
! wall at the southern face.
! halfpipe - a zonally uniform channel with a half-sine
! profile in the meridional direction.
! benchmark - use the benchmark test case topography.
! DOME - use a slope and channel configuration for the
! DOME sill-overflow test case.
! DOME2D - use a shelf and slope configuration for the
! DOME2D gravity current/overflow test case.
! seamount - Gaussian bump for spontaneous motion test case.
! USER - call a user modified routine."
TOPO_FILE = topog.nc
! "default = 'topog.nc'
! The file from which the bathymetry is read."
MAXIMUM_DEPTH = 6000.0
! "[m]
! The maximum depth of the ocean."
MINIMUM_DEPTH = 0.5
! "[m] default = 0.0
! If MASKING_DEPTH is unspecified, then anything shallower than
! MINIMUM_DEPTH is assumed to be land and all fluxes are masked out.
! If MASKING_DEPTH is specified, then all depths shallower than
! MINIMUM_DEPTH but depper than MASKING_DEPTH are rounded to
! MINIMUM_DEPTH."
REMAPPING_SCHEME = PPM_H4
! "default = 'PLM'
! This sets the reconstruction scheme used
! for vertical remapping for all variables.
! It can be one of the following schemes:
! PCM (1st-order accurate)
! PLM (2nd-order accurate)
! PPM_H4 (3rd-order accurate)
! PPM_IH4 (3rd-order accurate)
! PQM_IH4IH3 (4th-order accurate)
! PQM_IH6IH5 (5th-order accurate)"
INIT_LAYERS_FROM_Z_FILE = True
! "[Boolean] default = False
! If true, intialize the layer thicknesses, temperatures,
! and salnities from a Z-space file on a latitude-
! longitude grid."
TEMP_SALT_Z_INIT_FILE = ocean_temp_salt.res.nc
! "default = 'temp_salt_z.nc'
! The name of the z-space input file used to initialize
! the layer thicknesses, temperatures and salinities."
Z_INIT_FILE_PTEMP_VAR = temp
! "default = 'ptemp'
! The name of the potential temperature variable in
! TEMP_SALT_Z_INIT_FILE. Conservative temperature in ocean_temp_salt.res.nc"
Z_INIT_FILE_SALT_VAR = salt
! "default = 'salt'
! The name of the salinity variable in
! TEMP_SALT_Z_INIT_FILE."
Z_INIT_REMAP_OLD_ALG = False
! "[Boolean] default = True
! If false, uses the preferred remapping algorithm for initialization. If true,
! use an older, less robust algorithm for remapping."
USE_VARIABLE_MIXING = True
! "[Boolean] default = False
! If true, the variable mixing code will be called. This
! allows diagnostics to be created even if the scheme is
! not used. If KHTR_SLOPE_CFF>0 or KhTh_Slope_Cff>0,
! this is set to true regardless of what is in the
! parameter file."
RESOLN_SCALED_KH = True
! "[Boolean] default = False
! If true, the Laplacian lateral viscosity is scaled away
! when the first baroclinic deformation radius is well
! resolved."
KHTH_SLOPE_CFF = 0.01
! "[nondim] default = 0.0
! The nondimensional coefficient in the Visbeck formula for the interface depth
! diffusivity"
DEPTH_SCALED_KHTH = False
! "[Boolean] default = False
! If true, KHTH is scaled away when the depth is shallower
! than a reference depth: KHTH = MIN(1,H/H0)**N * KHTH,
! where H0 is a reference depth, controlled via
! DEPTH_SCALED_KHTH_H0, and theexponent (N) is
! controlled via DEPTH_SCALED_KHTH_EXP."
KHTH_USE_EBT_STRUCT = True
! "[Boolean] default = False
! If true, uses the equivalent barotropic structure
! as the vertical structure of thickness diffusivity."
RESOLN_SCALED_KHTH = True
! "[Boolean] default = False
! If true, the interface depth diffusivity is scaled away
! when the first baroclinic deformation radius is well
! resolved."
USE_GM_WORK_BUG = False
! "[Boolean] default = True
! If true, compute the top-layer work tendency on the u-grid with the incorrect
! sign, for legacy reproducibility."
USE_STORED_SLOPES = True
! "[Boolean] default = False
! If true, the isopycnal slopes are calculated once and
! stored for re-use. This uses more memory but avoids calling
! the equation of state more times than should be necessary."
ETA_TOLERANCE = 1e-06
! "[m] default = 3.15E-09
! The tolerance for the differences between the
! barotropic and baroclinic estimates of the sea surface
! height due to the fluxes through each face. The total
! tolerance for SSH is 4 times this value. The default
! is 0.5*NK*ANGSTROM, and this should not be set less x
! than about 10^-15*MAXIMUM_DEPTH."
USE_NEUTRAL_DIFFUSION = True
! "[Boolean] default = False
! If true, enables the neutral diffusion module."
NDIFF_INTERIOR_ONLY = True
! "[Boolean] default = False
! If true, only applies neutral diffusion in the ocean interior. That is, the
! algorithm will exclude the surface and bottomboundary layers."
USE_LATERAL_BOUNDARY_DIFFUSION = True
! "[Boolean] default = False
! If true, enables the lateral boundary tracer's diffusion module."
LBD_LINEAR_TRANSITION = True
! "[Boolean] default = False
! If True, apply a linear transition at the base/top of the boundary.
! The flux will be fully applied at k=k_min and zero at k=k_max."
SIMPLE_TKE_TO_KD = True
! "[Boolean] default = False
! If true, uses a simple estimate of Kd/TKE that will
! work for arbitrary vertical coordinates. If false,
! calculates Kd/TKE and bounds based on exact
! energetics/nfor an isopycnal layer-formulation."
REGRIDDING_COORDINATE_MODE = Z*
! "default = 'LAYER'
! Coordinate mode for vertical regridding.
! Choose among the following possibilities:
! LAYER - Isopycnal or stacked shallow water layers
! ZSTAR, Z* - stetched geopotential z*
! SIGMA_SHELF_ZSTAR - stetched geopotential z* ignoring shelf
! SIGMA - terrain following coordinates
! RHO - continuous isopycnal
! HYCOM1 - HyCOM-like hybrid coordinate
! SLIGHT - stretched coordinates above continuous isopycnal
! ADAPTIVE - optimize for smooth neutral density surfaces"
ALE_COORDINATE_CONFIG = "FILE:ocean_vgrid.nc,nzv"
! "default = 'UNIFORM'
! Determines how to specify the coordinate
! resolution. Valid options are:
! PARAM - use the vector-parameter ALE_RESOLUTION
! UNIFORM[:N] - uniformly distributed
! FILE:string - read from a file. The string specifies
! the filename and variable name, separated
! by a comma or space, e.g. FILE:lev.nc,dz
! or FILE:lev.nc,interfaces=zw
! WOA09[:N] - the WOA09 vertical grid (approximately)
! FNC1:string - FNC1:dz_min,H_total,power,precision
! HYBRID:string - read from a file. The string specifies
! the filename and two variable names, separated
! by a comma or space, for sigma-2 and dz. e.g.
! HYBRID:vgrid.nc,sigma2,dz"
BOUND_CORIOLIS = True
! "[Boolean] default = False
! If true, the Coriolis terms at u-points are bounded by
! the four estimates of (f+rv)v from the four neighboring
! v-points, and similarly at v-points. This option would
! have no effect on the SADOURNY Coriolis scheme if it
! were possible to use centered difference thickness fluxes."
TIDE_M2 = True
! "[Boolean] default = False
! If true, apply tidal momentum forcing at the M2
! frequency. This is only used if TIDES is true."
TIDE_SAL_SCALAR_VALUE = 0.094
! "[m m-1]
! The constant of proportionality between sea surface
! height (really it should be bottom pressure) anomalies
! and bottom geopotential anomalies. This is only used if
! TIDES and TIDE_USE_SAL_SCALAR are true."
MASS_WEIGHT_IN_PRESSURE_GRADIENT = True
! "[Boolean] default = False
! If true, use mass weighting when interpolation T/S for
! top/bottom integrals in AFV pressure gradient calculation."
LAPLACIAN = True
! "[Boolean] default = False
! If true, use a Laplacian horizontal viscosity."
KH_VEL_SCALE = 0.0
! "[m s-1] default = 0.0
! The velocity scale which is multiplied by the grid
! spacing to calculate the Laplacian viscosity.
! The final viscosity is the largest of this scaled
! viscosity, the Smagorinsky viscosity and KH."
AH = 1000000000000.0
! "[m4 s-1] default = 0.0
! The background biharmonic horizontal viscosity."
AH_VEL_SCALE = 0.0
! "[m s-1] default = 0.0
! The velocity scale which is multiplied by the cube of
! the grid spacing to calculate the biharmonic viscosity.
! The final viscosity is the largest of this scaled
! viscosity, the Smagorinsky and Leith viscosities, and AH."
LEITH_AH = True
! "[Boolean] default = False
! If true, use a biharmonic Leith nonlinear eddy viscosity."
LEITH_BI_CONST = 128.0
! "[nondim] default = 0.0
! The nondimensional biharmonic Leith constant, typical values are thus far
! undetermined."
USE_LAND_MASK_FOR_HVISC = True
! "[Boolean] default = False
! If true, use Use the land mask for the computation of thicknesses
! at velocity locations. This eliminates the dependence on arbitrary
! values over land or outside of the domain. Default is False in order to
! maintain answers with legacy experiments but should be changed to True
! for new experiments."
HMIX_FIXED = 0.5
! "[m]
! The prescribed depth over which the near-surface
! viscosity and diffusivity are elevated when the bulk
! mixed layer is not used."
CHANNEL_DRAG = True
! "[Boolean] default = False
! If true, the bottom drag is exerted directly on each
! layer proportional to the fraction of the bottom it
! overlies."
HBBL = 10.0
! "[m]
! The thickness of a bottom boundary layer with a viscosity of KVBBL if
! BOTTOMDRAGLAW is not defined, or the thickness over which near-bottom
! velocities are averaged for the drag law if BOTTOMDRAGLAW is defined but
! LINEAR_DRAG is not."
PRANDTL_TURB = 1.0
! "[nondim] default = 0.0
! The turbulent Prandtl number applied to shear
! instability."
U_TRUNC_FILE = U_velocity_truncations
! "default = ''
! The absolute path to a file into which the accelerations
! leading to zonal velocity truncations are written.
! Undefine this for efficiency if this diagnostic is not
! needed."
V_TRUNC_FILE = V_velocity_truncations
! "default = ''
! The absolute path to a file into which the accelerations
! leading to meridional velocity truncations are written.
! Undefine this for efficiency if this diagnostic is not
! needed."
KV = 0.0001
! "[m2 s-1]
! The background kinematic viscosity in the interior.
! The molecular value, ~1e-6 m2 s-1, may be used."
MAXVEL = 6.0
! "[m s-1] default = 3.0E+08
! The maximum velocity allowed before the velocity
! components are truncated."
CFL_TRUNCATE_RAMP_TIME = 7200.0
! "[s] default = 0.0
! The time over which the CFL trunction value is ramped
! up at the beginning of the run."
Z_INIT_ALE_REMAPPING = True
! "[Boolean] default = False
! If True, then remap straight to model coordinate from file."
INTERPOLATE_RES_FN = False
! "[Boolean] default = True
! If true, interpolate the resolution function to the
! velocity points from the thickness points; otherwise
! interpolate the wave speed and calculate the resolution
! function independently at each point."
GILL_EQUATORIAL_LD = True
! "[Boolean] default = False
! If true, uses Gill's definition of the baroclinic
! equatorial deformation radius, otherwise, if false, use
! Pedlosky's definition. These definitions differ by a factor
! of 2 infront of the beta term in the denominator. Gill'sis the more appropriate definition."
DRAG_BG_VEL = 0.1
! "[m s-1] default = 0.0
! DRAG_BG_VEL is either the assumed bottom velocity (with
! LINEAR_DRAG) or an unresolved velocity that is
! combined with the resolved velocity to estimate the
! velocity magnitude. DRAG_BG_VEL is only used when
! BOTTOMDRAGLAW is defined."
BBL_USE_EOS = True
! "[Boolean] default = False
! If true, use the equation of state in determining the
! properties of the bottom boundary layer. Otherwise use
! the layer target potential densities."
BBL_THICK_MIN = 0.1
! "[m] default = 0.0
! The minimum bottom boundary layer thickness that can be
! used with BOTTOMDRAGLAW. This might be
! Kv / (cdrag * drag_bg_vel) to give Kv as the minimum
! near-bottom viscosity."
BOUND_BT_CORRECTION = True
! "[Boolean] default = False
! If true, the corrective pseudo mass-fluxes into the
! barotropic solver are limited to values that require
! less than 0.1*MAXVEL to be accommodated.
! The barotropic x-halo size that is actually used.
! The barotropic y-halo size that is actually used."
BT_PROJECT_VELOCITY = True
! "[Boolean] default = False
! If true, step the barotropic velocity first and project
! out the velocity tendancy by 1+BEBT when calculating the
! transport. The default (false) is to use a predictor
! continuity step to find the pressure field, and then
! to do a corrector continuity step using a weighted
! average of the old and new velocities, with weights
! of (1-BEBT) and BEBT."
BT_THICK_SCHEME = FROM_BT_CONT
! "default = 'HYBRID'
! A string describing the scheme that is used to set the
! open face areas used for barotropic transport and the
! relative weights of the accelerations. Valid values are:
! ARITHMETIC - arithmetic mean layer thicknesses
! HARMONIC - harmonic mean layer thicknesses
! HYBRID (the default) - use arithmetic means for
! layers above the shallowest bottom, the harmonic
! mean for layers below, and a weighted average for
! layers that straddle that depth
! FROM_BT_CONT - use the average thicknesses kept
! in the h_u and h_v fields of the BT_cont_type"
BEBT = 0.2
! "[nondim] default = 0.1
! BEBT determines whether the barotropic time stepping
! uses the forward-backward time-stepping scheme or a
! backward Euler scheme. BEBT is valid in the range from
! 0 (for a forward-backward treatment of nonrotating
! gravity waves) to 1 (for a backward Euler treatment).
! In practice, BEBT must be greater than about 0.05."
DTBT = -0.95
! "[s or nondim] default = -0.98
! The barotropic time step, in s. DTBT is only used with
! the split explicit time stepping. To set the time step
! automatically based the maximum stable value use 0, or
! a negative value gives the fraction of the stable value.
! Setting DTBT to 0 is the same as setting it to -0.98.
! The value of DTBT that will actually be used is an
! integer fraction of DT, rounding down."
KHTH = 0.0
! "[m2 s-1] default = 0.0
! The background horizontal thickness diffusivity."
KHTH_MAX = 0.0
! "[m2 s-1] default = 0.0
! The maximum horizontal thickness diffusivity."
USE_KH_IN_MEKE = True
! "[Boolean] default = False
! If true, uses the thickness diffusivity calculated here to diffuse MEKE."
KHTH_USE_FGNV_STREAMFUNCTION = True
! "[Boolean] default = False
! If true, use the streamfunction formulation of
! Ferrari et al., 2010, which effectively emphasizes
! graver vertical modes by smoothing in the vertical."
FGNV_C_MIN = 0.01
! "[m s-1] default = 0.0
! A minium wave speed used in the Ferrari et al., 2010,
! streamfunction formulation."
MIXEDLAYER_RESTRAT = True
! "[Boolean] default = False
! If true, a density-gradient dependent re-stratifying
! flow is imposed in the mixed layer.
! This is only used if BULKMIXEDLAYER is true."
FOX_KEMPER_ML_RESTRAT_COEF = 1.0
! "[nondim] default = 0.0
! A nondimensional coefficient that is proportional to
! the ratio of the deformation radius to the dominant
! lengthscale of the submesoscale mixed layer
! instabilities, times the minimum of the ratio of the
! mesoscale eddy kinetic energy to the large-scale
! geostrophic kinetic energy or 1 plus the square of the
! grid spacing over the deformation radius, as detailed
! by Fox-Kemper et al. (2010)
! The file that specifies the vertical grid for
! depth-space diagnostics, or blank to disable
! depth-space output.
! The number of depth-space levels. This is determined
! from the size of the variable zw in the output grid file."
MLE_FRONT_LENGTH = 1000.0
! "[m] default = 0.0
! If non-zero, is the frontal-length scale used to calculate the
! upscaling of buoyancy gradients that is otherwise represented
! by the parameter FOX_KEMPER_ML_RESTRAT_COEF. If MLE_FRONT_LENGTH is
! non-zero, it is recommended to set FOX_KEMPER_ML_RESTRAT_COEF=1."
MLE_MLD_DECAY_TIME = 345600.0
! "[s] default = 0.0
! The time-scale for a running-mean filter applied to the mixed-layer
! depth used in the MLE restratification parameterization. When
! the MLD deepens below the current running-mean the running-mean
! is instantaneously set to the current MLD."
USE_CVMix_CONVECTION = True
! "[Boolean] default = False
! If true, turns on the enhanced mixing due to convection
! via CVMix. This scheme increases diapycnal diffs./viscs.
! at statically unstable interfaces. Relevant parameters are
! contained in the CVMix_CONVECTION% parameter block."
BBL_MIXING_AS_MAX = False
! "[Boolean] default = True
! If true, take the maximum of the diffusivity from the
! BBL mixing and the other diffusivities. Otherwise,
! diffusiviy from the BBL_mixing is simply added."
USE_LOTW_BBL_DIFFUSIVITY = True
! "[Boolean] default = False
! If true, uses a simple, imprecise but non-coordinate dependent, model
! of BBL mixing diffusivity based on Law of the Wall. Otherwise, uses
! the original BBL scheme."
KD_MAX = 0.1
! "[m2 s-1] default = -1.0
! The maximum permitted increment for the diapycnal
! diffusivity from TKE-based parameterizations, or a
! negative value for no limit."
HORIZ_VARYING_BACKGROUND = True
! "[Boolean] default = False
! If true, apply vertically uniform, latitude-dependent background diffusivity,
! as described in Danabasoglu et al., 2012."
KD = 2e-05
! "[m2 s-1]
! The background diapycnal diffusivity of density in the
! interior. Zero or the molecular value, ~1e-7 m2 s-1,
! may be used."
PRANDTL_BKGND = 5.0
! "[nondim] default = 1.0
! Turbulent Prandtl number used to convert vertical background diffusivities
! into viscosities."
KD_MIN = 2e-06
! "[m2 s-1] default = 2.0E-07
! The minimum diapycnal diffusivity.
! If true, use an internal tidal dissipation scheme to
! drive diapycnal mixing, along the lines of St. Laurent
! et al. (2002) and Simmons et al. (2004)."
INT_TIDE_DECAY_SCALE = 300.3003003003003
! "[m] default = 500.0
! The decay scale away from the bottom for tidal TKE with
! the new coding when INT_TIDE_DISSIPATION is used."
USE_KPP = True
! "[Boolean] default = False
! If true, turns on the [CVMix] KPP scheme of Large et al., 1994,
! to calculate diffusivities and non-local transport in the OBL."
USE_LMD94 = True
! "[Boolean] default = False
! If true, use the Large-McWilliams-Doney (JGR 1994)
! shear mixing parameterization."
SMOOTH_RI = True
! "[Boolean] default = False
! If true, vertically smooth the Richardsonnumber by applying a 1-2-1 filter once."
USE_CVMIX_DDIFF = True
! "[Boolean] default = False
! If true, turns on double diffusive processes via CVMix.
! Note that double diffusive processes on viscosity are ignored
! in CVMix, see http://cvmix.github.io/ for justification."
MAX_ENT_IT = 20
! "default = 5
! The maximum number of iterations that may be used to
! calculate the interior diapycnal entrainment."
TOLERANCE_ENT = 1e-05
! "[m] default = 2.683281572999748E-05
! The tolerance with which to solve for entrainment values."
HMIX_MIN = 2.0
! "[m] default = 0.0
! The minimum mixed layer depth if the mixed layer depth
! is determined dynamically."
PEN_SW_SCALE = 15.0
! "[m] default = 0.0
! The vertical absorption e-folding depth of the
! penetrating shortwave radiation."
PEN_SW_FRAC = 0.42
! "[nondim] default = 0.0
! The fraction of the shortwave radiation that penetrates
! below the surface."
PRESSURE_DEPENDENT_FRAZIL = True
! "[Boolean] default = False
! If true, use a pressure dependent freezing temperature
! when making frazil. The default is false, which will be
! faster but is inappropriate with ice-shelf cavities."
TRACER_ADVECTION_SCHEME = "PPM:H3"
! "default = 'PLM'
! The horizontal transport scheme for tracers:
! PLM - Piecewise Linear Method
! PPM:H3 - Piecewise Parabolic Method (Huyhn 3rd order)"
KHTR_MIN = 50.0
! "[m2 s-1] default = 0.0
! The minimum along-isopycnal tracer diffusivity."
KHTR_MAX = 0.0
! "[m2 s-1] default = 0.0
! The maximum along-isopycnal tracer diffusivity."
DEBUG = False
! "If true, write out verbose debugging data."
CHECK_DIFFUSIVE_CFL = True
! "[Boolean] default = False
! If true, use enough iterations the diffusion to ensure
! that the diffusive equivalent of the CFL limit is not
! violated. If false, always use the greater of 1 or
! MAX_TR_DIFFUSION_CFL iteration."
MAX_TR_DIFFUSION_CFL = 2.0
! "[nondim] default = -1.0
! If positive, locally limit the along-isopycnal tracer
! diffusivity to keep the diffusive CFL locally at or
! below this value. The number of diffusive iterations
! is often this value or the next greater integer."
MAXTRUNC = 0
! "[truncations save_interval-1] default = 0
! The run will be stopped, and the day set to a very
! large value if the velocity is truncated more than
! MAXTRUNC times between energy saves. Set MAXTRUNC to 0
! to stop if there is any truncation of velocities."
OCEAN_SURFACE_STAGGER = A
! "default = 'C'
! A case-insensitive character string to indicate the
! staggering of the surface velocity field that is
! returned to the coupler. Valid values include
! 'A', 'B', or 'C'."
RESTORE_SALINITY = False
! "[Boolean] default = False
! If true, the coupled driver will add a globally-balanced
! fresh-water flux that drives sea-surface salinity
! toward specified values."
FLUXCONST = 0.5
! "[m day-1]
! The constant that relates the restoring surface fluxes
! to the relative surface anomalies (akin to a piston
! velocity). Note the non-MKS units."
ADJUST_NET_FRESH_WATER_TO_ZERO = True
! "[Boolean] default = False
! If true, adjusts the net fresh-water forcing seen
! by the ocean (including restoring) to zero."
GUST_CONST = 0.02
! "[Pa] default = 0.02
! The background gustiness in the winds."
FIX_USTAR_GUSTLESS_BUG = True
! "[Boolean] default = False
! If true correct a bug in the time-averaging of the gustless wind
! friction velocity."
RESTART_CONTROL = 3
! "default = 1
! An integer whose bits encode which restart files are
! written. Add 2 (bit 1) for a time-stamped file, and odd
! (bit 0) for a non-time-stamped file. A non-time-stamped
! restart file is saved at the end of the run segment
! for any non-negative value."
ENERGYSAVEDAYS = 1.0
! "[days] default = 1.0
! The interval in units of TIMEUNIT between saves of the
! energies of the run and other globally summed diagnostics."
TIDES = True
! "[Boolean] default = False
! If true, apply tidal momentum forcing."
USE_MEKE = True
! "[Boolean] default = False
! If true, turns on the MEKE scheme which calculates
! a sub-grid mesoscale eddy kinetic energy budget."
MEKE_GMCOEFF = 0.0
! "[nondim] default = -1.0
! The efficiency of the conversion of potential energy
! into MEKE by the thickness mixing parameterization.
! If MEKE_GMCOEFF is negative, this conversion is not
! used or calculated."
MEKE_KHMEKE_FAC = 0.5
! "[nondim] default = 0.0
! A factor that maps MEKE%Kh to Kh for MEKE itself."
MEKE_VISCOSITY_COEFF_KU = 0.2
! "[nondim] default = 0.0
! If non-zero, is the scaling coefficient in the expression forviscosity used to
! parameterize harmonic lateral momentum mixing byunresolved eddies represented
! by MEKE. Can be negative torepresent backscatter from the unresolved eddies."
MEKE_MIN_LSCALE = True
! "[Boolean] default = False
! If true, use a strict minimum of provided length scales rather than harmonic
! mean."
MEKE_ALPHA_RHINES = 1.0
! "[nondim] default = 0.05
! If positive, is a coefficient weighting the Rhines scale
! in the expression for mixing length used in MEKE-derived diffusiviity."
MEKE_ALPHA_EADY = 1.0
! "[nondim] default = 0.05
! If positive, is a coefficient weighting the Eady length scale
! in the expression for mixing length used in MEKE-derived diffusiviity."
MEKE_ALPHA_DEFORM = 1.0
! "[nondim] default = 0.0
! If positive, is a coefficient weighting the deformation scale in the
! expression for mixing length used in MEKE-derived diffusivity."
MEKE_ALPHA_FRICT = 1.0
! "[nondim] default = 0.0
! If positive, is a coefficient weighting the frictional arrest scale in the
! expression for mixing length used in MEKE-derived diffusivity."
MEKE_ALPHA_GRID = 1.0
! "[nondim] default = 0.0
! If positive, is a coefficient weighting the grid-spacing as a scale in the
! expression for mixing length used in MEKE-derived diffusivity."
KH_RES_SCALE_COEF = 0.4
! "[nondim] default = 1.0
! A coefficient that determines how KhTh is scaled away if RESOLN_SCALED_... is
! true, as F = 1 / (1 + (KH_RES_SCALE_COEF*Rd/dx)^KH_RES_FN_POWER)."
VISC_RES_SCALE_COEF = 0.4
! "[nondim] default = 1.0
! A coefficient that determines how Kh is scaled away if RESOLN_SCALED_... is
! true, as F = 1 / (1 + (KH_RES_SCALE_COEF*Rd/dx)^KH_RES_FN_POWER). This
! function affects lateral viscosity, Kh, and not KhTh."
MEKE_GEOMETRIC = True
! "[Boolean] default = False
! If MEKE_GEOMETRIC is true, uses the GM coefficient formulation from the
! GEOMETRIC framework (Marshall et al., 2012)."
MEKE_KHTH_FAC = 1.0
! "[nondim] default = 0.0
! A factor that maps MEKE%Kh to KhTh."
MEKE_KHTR_FAC = 1.0
! "[nondim] default = 0.0
! A factor that maps MEKE%Kh to KhTr."
MEKE_EQUILIBRIUM_ALT = True
! "[Boolean] default = False
! If true, use an alternative formula for computing the (equilibrium)initial
! value of MEKE."
MEKE_VISC_DRAG = False
! "[Boolean] default = True
! If true, use the vertvisc_type to calculate the bottom drag acting on MEKE."
MEKE_EQUILIBRIUM_RESTORING = True
! "[Boolean] default = False
! If true, restore MEKE back to its equilibrium value, which is calculated
! at each time step."
MEKE_RESTORING_TIMESCALE = 10000000.0
! "[s] default = 1.0E+06
! The timescale used to nudge MEKE toward its equilibrium value."
MEKE_ADVECTION_FACTOR = 1.0
! "[nondim] default = 0.0
! A scale factor in front of advection of eddy energy. Zero turns advection off.
! Using unity would be normal but other values could accommodate a mismatch
! between the advecting barotropic flow and the vertical structure of MEKE."
DIAG_COORD_DEF_Z = WOA09
! "default = WOA09
! Determines how to specify the coordinate
! resolution. Valid options are:
! PARAM - use the vector-parameter DIAG_COORD_RES_Z
! UNIFORM[:N] - uniformly distributed
! FILE:string - read from a file. The string specifies
! the filename and variable name, separated
! by a comma or space, e.g. FILE:lev.nc,dz
! or FILE:lev.nc,interfaces=zw
! WOA09[:N] - the WOA09 vertical grid (approximately)
! FNC1:string - FNC1:dz_min,H_total,power,precision
! HYBRID:string - read from a file. The string specifies
! the filename and two variable names, separated
! by a comma or space, for sigma-2 and dz. e.g.
! HYBRID:vgrid.nc,sigma2,dz"
MAX_RINO_IT = 25
! "[nondim] default = 50
! The maximum number of iterations that may be used to
! estimate the Richardson number driven mixing."
CHANNEL_CONFIG = global_1deg
! "default = 'none'
! A parameter that determines which set of channels are
! restricted to specific widths. Options are:
! none - All channels have the grid width.
! global_1deg - Sets 16 specific channels appropriate
! for a 1-degree model, as used in CM2G.
! list - Read the channel locations and widths from a
! text file, like MOM_channel_list in the MOM_SIS
! test case.
! file - Read open face widths everywhere from a
! NetCDF file on the model grid."
SMAG_BI_CONST = 0.06
! "[nondim] default = 0.0
! The nondimensional biharmonic Smagorinsky constant,
! typically 0.015 - 0.06."
KAPPA_SHEAR_ITER_BUG = False
! [Boolean] default = True If true, use an older, dimensionally inconsistent estimate of the derivative of diffusivity with energy in the Newton's method iteration. The bug causes undercorrections when dz > 1 m.
KAPPA_SHEAR_ALL_LAYER_TKE_BUG = False
! [Boolean] default = True If true, report back the latest estimate of TKE instead of the time average TKE when there is mass in all layers. Otherwise always report the time averaged TKE, as is currently done when there are some massless layers.