Index

Rundeck: P1SoM40

Created: Wed May 8 03:30:18 EDT 2024

KPPE

File: OCNKPP.f

Summary: KPPE contains variables and routines for KPP mixing scheme
Author : NCAR (modifications by Gavin Schmidt)
Version:

---

Subroutines:

---

Global Variables:

lmo : used from ocean

Used by: | :bldepth |

fkph : used from oceanres

= vertical diffusion coefficient (cm**2/sec) (taken from OCEANRES)

fkpm : used from oceanres

= vertical viscosity coefficient (cm**2/sec) (taken from OCEANRES)

lsrpd : used from sw2ocean

Used by: | :bldepth |

fsr : used from sw2ocean

Used by: | :bldepth |

fsrz : used from sw2ocean

dfsrdz : used from sw2ocean

Used by: | :bldepth |

dfsrdzb : used from sw2ocean

Used by: | :bldepth |

km : integer, parameter

is max number of KPP layers
Initial Value = LMO

c0 : real*8, parameter

Initial Value = 0.0

c1 : real*8, parameter

Initial Value = 1.0

c2 : real*8, parameter

Initial Value = 2.0

c4 : real*8, parameter

Initial Value = 4.0

c5 : real*8, parameter

Initial Value = 5.0

c8 : real*8, parameter

Initial Value = 8.0

c16 : real*8, parameter

Initial Value = 16.0

c360 : real*8, parameter

Initial Value = 360.0

p25 : real*8, parameter

Initial Value = 0.25

p5 : real*8, parameter

Initial Value = 0.5

p75 : real*8, parameter

Initial Value = 0.75

epsln : real*8, parameter

Initial Value = 1.0e-20

c24 : real*8, parameter

Initial Value = 24.0

c60 : real*8, parameter

Initial Value = 60.0

c1440 : real*8, parameter

Initial Value = 1440.0

r24 : real*8, parameter

Initial Value = c1/c24

r60 : real*8, parameter

Initial Value = c1/c60

r1440 : real*8, parameter

Initial Value = c1/c1440

secday : real*8, parameter

Initial Value = c1/(c60*c1440)

c1p5 : real*8, parameter

Initial Value = 1.5

c3 : real*8, parameter

Initial Value = 3.

c35 : real*8, parameter

Initial Value = 35.

c10 : real*8, parameter

Initial Value = 10.

c100 : real*8, parameter

Initial Value = 100.

c1000 : real*8, parameter

Initial Value = 1000.

c10000 : real*8, parameter

Initial Value = 10000.

r3 : real*8, parameter

Initial Value = c1/c3

r10 : real*8, parameter

Initial Value = c1/c10

r100 : real*8, parameter

Initial Value = r10/c10

fkph : used from oceanres

= vertical diffusion coefficient (cm**2/sec) (taken from OCEANRES)

fkpm : used from oceanres

= vertical viscosity coefficient (cm**2/sec) (taken from OCEANRES)

bvdc :

= background vertical diffusion constant

bvvc :

= background vertical viscosity constant

vvclim : real*8, parameter

= vertical viscosity coefficient limit
Initial Value = 1000d0

vdclim : real*8, parameter

= vertical diffusion coefficient limit
Initial Value = 1000d0

vvcric : real*8, parameter

= maximum viscosity due to shear instability (cm**2/s)
Initial Value = 50d0

vdcric : real*8, parameter

= maximum diffusivity due to shear instability (cm**2/s)
Initial Value = 50d0

r1000 : real*8, parameter

Initial Value = c1/c1000

r10000 : real*8, parameter

Initial Value = r10*r1000

epsl : real*8, parameter

= epsln in "pconst.h" (set in "kmixinit") = 1.0e-20
Initial Value = epsln
Used by: | :bldepth |

epsilon : real*8, parameter

= nondimensional extent of the surface layer = 0.1
Initial Value = 0.1d0
Used by: | :bldepth |

vonk : real*8, parameter

= von Karman's constant = 0.4
Initial Value = 0.4d0
Used by: | :bldepth |

zetas : real*8, parameter

Initial Value = -1d0

conc3 : real*8, parameter

Initial Value = 16d0

concs : real*8, parameter

Initial Value = 98.96d0

conas : real*8, parameter

Initial Value = -28.86d0

zetam : real*8, parameter

Initial Value = -0.2d0

conc2 : real*8, parameter

Initial Value = 16d0

concm : real*8, parameter

Initial Value = 8.380d0

conam : real*8, parameter

Initial Value = 1.257d0

conc1 : real*8, parameter

Initial Value = 5d0
Used by: | :bldepth |

ricr : real*8, parameter

= critical bulk Richardson Number = 0.3 or 1.
Initial Value = 0.3d0
Used by: | :bldepth |

cekman : real*8, parameter

= coefficient for ekman depth = 0.7
Initial Value = 0.7d0

cmonob : real*8, parameter

= coefficient for Monin-Obukhov depth = 1.0
Initial Value = 1d0

concv : real*8, parameter

= ratio of interior buoyancy frequency to
Initial Value = 1.8d0

hbf : real*8, parameter

= fraction of bounadry layer depth to
Initial Value = 1d0

vtc : real*8

= non-dimensional coefficient for velocity
Used by: | :bldepth |

nni : integer, parameter

= number of values for zehat in the look up table
Initial Value = 890
Used by: | :bldepth |

nnj : integer, parameter

= number of values for ustar in the look up table
Initial Value = 480
Used by: | :bldepth |

wmt : real*8, dimension(0:nni+1, 0:nnj+1)

= lookup table for wm, the turbulent velocity scale
Used by: | :bldepth |

wst : real*8, dimension(0:nni+1, 0:nnj+1)

= lookup table for ws, the turbulent velocity scale
Used by: | :bldepth |

deltaz : real*8, parameter

= delta zehat in table
Initial Value = (zmax-zmin)

deltau :

= delta ustar in table

rdeltaz :

= recipricol of delta zehat in table
Used by: | :bldepth |

rdeltau :

= recipricol of delta ustar in table
Used by: | :bldepth |

zmin : real*8, parameter

= minimum limit for zehat in table (m3/s3)
Initial Value = -4d-7
Used by: | :bldepth |

zmax : real*8, parameter

= maximum limit for zehat in table
Initial Value = 0.
Used by: | :bldepth |

umin : real*8, parameter

= minimum limit for ustar in table (m/s)
Initial Value = 0.
Used by: | :bldepth |

umax : real*8, parameter

= maximum limit for ustar in table
Initial Value = 4d-2

riinfty : real*8, parameter

= local Richardson Number limit
Initial Value = 0.7d0

rriinfty : real*8, parameter

= 1/Riinfty
Initial Value = 1./Riinfty

difm0 : real*8, parameter

= viscosity max due to shear instability = vvcric
Initial Value = vvcric * r10000

difs0 : real*8, parameter

= diffusivity .. = vdcric
Initial Value = vdcric * r10000

difmiw : real*8, parameter

= viscosity background due to internal waves = fkpm
Initial Value = fkpm * r10000

difsiw : real*8, parameter

= diffusivity .. = fkph
Initial Value = fkph * r10000

difmcon : real*8, parameter

= viscosity due to convective instability = vvclim
Initial Value = vvclim * r10000

difscon : real*8, parameter

= diffusivity .. = vdclim
Initial Value = vdclim * r10000

bvsqcon : real*8, parameter

= value of N^2 where the convection coefs first become max
Initial Value = -1d-7

rbvsqcon : real*8, parameter

= 1/BVSQcon
Initial Value = 1./BVSQcon

num_v_smooth_ri : integer, parameter

= number of vertical smoothings of Ri
Initial Value = 1

rrho0 : real*8, parameter

= limit for double diffusive density ratio
Initial Value = 1.9d0

dsfmax : real*8, parameter

= maximum diffusivity in case of salt fingering (m2/s)
Initial Value = 0.001d0

cstar : real*8, parameter

= proportionality coefficient for nonlocal transport
Initial Value = 10d0

cg : real*8

= non-dimensional coefficient for counter-gradient term

diftop : real*8, parameter

diffusion scale for topographic mixing
Initial Value = 0d0

fz500 : real*8, dimension(km, km)

diffusion scale as function of height above topography