Index | Rundeck: P1SoM40 | Created: Fri May 10 03:30:18 EDT 2024 |
gissmix | Module: | File: OCNGISSVM.f |
Author : AHoward and YCheng
Version:
Declaration:
subroutine gissmix( n, ze, zg, db, dv2, adt, bds, rho, ustarb2, exy, fc, hbl, strait, ri, rr, bv2, km, kh, ks, kc, e)
- grav :
- acceleration of gravity m/s^2
- lmo :
- max. number of vertical layers
- taubx :
- (IM,J_0H:J_1H) x component of tau_b = kinematic bottom drag (m/s)^2
- tauby :
- (IM,J_0H:J_1H) y component of tau_b = kinematic bottom drag (m/s)^2
- exya :
- internal tidal energy (w/m^2)
- ut2a :
- unresolved bottom shear squared (m/s)^2
- n :
integer
- number of vert. layers on this column
- ze :
real*8, dimension(0:lmo)
- vertical grid-edge depth (m), > 0
- zg :
real*8, dimension(0:lmo+1)
- vertical grid depth (m), < 0
- db :
real*8, dimension(lmo)
- -grav/rho*d(rho) (m/s^2)
- dv2 :
real*8, dimension(lmo)
- vel. diff. squared btw layers (m/s)^2
- adt :
real*8, dimension(lmo)
- rho*alpha*DT (kg/m^3)
- bds :
real*8, dimension(lmo)
- rho*beta*DS (kg/m^3)
- rho :
real*8, dimension(lmo)
- density
- ustarb2 :
real*8
- velocity squared at zg(n)
- exy :
real*8
- tidal power input
- fc :
real*8
- Coriolis parameter=2*omega*sin(lat) (1/s)
- hbl :
real*8
- pbl depth (m)
- strait :
integer
- 0: not in strait; 1: in strait
- ri :
real*8, dimension(0:lmo+1)
- local richardson number
- rr :
real*8, dimension(0:lmo+1)
- local alpha*dTdz/(beta*dSdz)
- bv2 :
real*8, dimension(0:lmo+1)
- Brunt Vaisala frequency squared (1/s**2)
- km :
real*8, dimension(0:lmo+1)
- vertical momentun diffusivity (m**2/s)
- kh :
real*8, dimension(0:lmo+1)
- vertical heat diffusivity (m**2/s)
- ks :
real*8, dimension(0:lmo+1)
- vertical salinity diffusivity (m**2/s)
- kc :
real*8, dimension(0:lmo+1)
- vertical passive scalar diffusivity (m**2/s)
- e :
real*8, dimension(lmo)
- ocean turbulent kinetic energy (m/s)**2
- flag :
integer
- =0 if abs(rr)<=1; =1 if abs(rr)>1
- l0 :
real*8
- constant length scale within obl
- l1 :
real*8
- length scale before reduced by stable buoyancy
- l1min :
real*8, parameter
- minimum of l1 below obl
Initial Value= 3.d0
- lr :
real*8
- length scale reduction factor by stable buoyancy
- l2 :
real*8
- length scale after reduced by stable buoyancy
- l2min :
real*8, parameter
- minimum of l2
Initial Value= .05d0
- len :
real*8
- final length scale
- num_smooth :
integer, parameter
- Initial Value
= 0
- vs2 :
real*8, dimension(0:lmo+1)
- velocity shear squared (1/s**2)
- iter :
integer
- mr :
integer
- len :
real*8
- final length scale
- l :
integer
- jlo :
integer
- jhi :
integer
- klo :
integer
- khi :
integer
- a1 :
real*8
- a2 :
real*8
- b1 :
real*8
- b2 :
real*8
- c1 :
real*8
- c2 :
real*8
- c3 :
real*8
- c4 :
real*8
- ril :
real*8
- rrl :
real*8
- gm :
real*8
- sm :
real*8
- sh :
real*8
- ss :
real*8
- sc :
real*8
- kml :
real*8
- khl :
real*8
- ksl :
real*8
- kcl :
real*8
- etau :
real*8
- kz :
real*8
- zbyh :
real*8
- bydz :
real*8
- zl :
real*8
- tmp :
real*8
- tmp1 :
real*8
- fbyden :
real*8
- afc :
real*8
- ltn :
real*8
- bvbyf :
real*8
- fac :
real*8
- kmbg :
real*8
- khbg :
real*8
- ksbg :
real*8
- den :
real*8
- fz :
real*8
- epstd_byn2 :
real*8
- kmtd :
real*8
- khtd :
real*8
- kstd :
real*8
- phim2 :
real*8
- zb :
real*8
- unr20 :
real*8
- lr :
real*8
- length scale reduction factor by stable buoyancy
- l0 :
real*8
- constant length scale within obl
- l1 :
real*8
- length scale before reduced by stable buoyancy
- l1min :
real*8, parameter
- minimum of l1 below obl
Initial Value= 3.d0
- l2 :
real*8
- length scale after reduced by stable buoyancy
- l2min :
real*8, parameter
- minimum of l2
Initial Value= .05d0
- len :
real*8
- final length scale