ProfileWithDoubleGaussianExponentialPlusPycnocline
Build a mixed Gaussian-exponential pycnocline profile used by the built-in benchmarks.
Developer documentation: this item describes internal implementation details.
Declaration
[rhoFunc,N2Func,zIn] = ProfileWithDoubleGaussianExponentialPlusPycnocline(rho_0,D,delta_p,z_p,L_s,L_d,z_T,L_deep,N0,Nq,Np)
Parameters
rho_0reference surface densityDwater-column depthdelta_ppycnocline thickness scalez_ppycnocline center depthL_sshallow Gaussian scaleL_ddeep Gaussian scalez_Ttransition depth to the exponential tailL_deepdeep exponential scaleN0surface buoyancy frequencyNqmixed-layer buoyancy-frequency scaleNppycnocline buoyancy-frequency scale
Returns
rhoFuncdensity function handleN2Funcbuoyancy-frequency function handlezIndepth domain for the constructed profile