rossby_front_2d_initialization character(len=40) character(len=40) rossby_front_2d_initialization::mdl mdl = "Rossby_front_2d_initialization" This module's name. rossby_front_initialize_temperature_salinity rossby_front_initialize_thickness rossby_front_initialize_velocity real, parameter real, parameter rossby_front_2d_initialization::frontfractionalwidth frontfractionalwidth = 0.5 Width of front as fraction of domain [nondim]. dtdy ypseudo real, parameter real, parameter rossby_front_2d_initialization::hmlmin hmlmin = 0.25 Shallowest ML as fractional depth of ocean [nondim]. dtdy hml real, parameter real, parameter rossby_front_2d_initialization::hmlmax hmlmax = 0.75 Deepest ML as fractional depth of ocean [nondim]. dtdy hml subroutine, public subroutine, public rossby_front_2d_initialization::rossby_front_initialize_thickness (h, G, GV, US, param_file, just_read_params) rossby_front_initialize_thickness h h G G GV GV US US param_file param_file just_read_params just_read_params Initialization of thicknesses in 2D Rossby front test. g Grid structure gv Vertical grid structure us A dimensional unit scaling type h The thickness that is being initialized [H ~> m or kg m-2] param_file A structure indicating the open file to parse for model parameter values. just_read_params If present and true, this call will only read parameters without changing h. hml mdl mom_error_handler::mom_error mom_error_handler::mom_mesg regrid_consts::regridding_rho regrid_consts::regridding_sigma subroutine, public subroutine, public rossby_front_2d_initialization::rossby_front_initialize_temperature_salinity (T, S, h, G, GV, param_file, eqn_of_state, just_read_params) rossby_front_initialize_temperature_salinity T T S S h h G G GV GV param_file param_file eqn_of_state eqn_of_state just_read_params just_read_params Initialization of temperature and salinity in the Rossby front test. g Grid structure gv The ocean's vertical grid structure. t Potential temperature [degC] s Salinity [ppt] h Thickness [H ~> m or kg m-2] param_file Parameter file handle eqn_of_state Equation of state structure just_read_params If present and true, this call will only read parameters without changing T & S. hml mdl subroutine, public subroutine, public rossby_front_2d_initialization::rossby_front_initialize_velocity (u, v, h, G, GV, US, param_file, just_read_params) rossby_front_initialize_velocity u u v v h h G G GV GV US US param_file param_file just_read_params just_read_params Initialization of u and v in the Rossby front test. g Grid structure gv Vertical grid structure u i-component of velocity [L T-1 ~> m s-1] v j-component of velocity [L T-1 ~> m s-1] h Thickness [H ~> m or kg m-2] us A dimensional unit scaling type param_file A structure indicating the open file to parse for model parameter values. just_read_params If present and true, this call will only read parameters without setting u & v. dtdy hml mdl mom_state_initialization::mom_initialize_state real function real function rossby_front_2d_initialization::ypseudo (G, lat) ypseudo G G lat lat Pseudo coordinate across domain used by Hml() and dTdy() returns a coordinate from -PI/2 .. PI/2 squashed towards the center of the domain. g Grid structure lat Latitude frontfractionalwidth dtdy hml real function real function rossby_front_2d_initialization::hml (G, lat) hml G G lat lat Analytic prescription of mixed layer depth in 2d Rossby front test, in the same units as Gmax_depth. g Grid structure lat Latitude hmlmax hmlmin ypseudo rossby_front_initialize_temperature_salinity rossby_front_initialize_thickness rossby_front_initialize_velocity real function real function rossby_front_2d_initialization::dtdy (G, dT, lat) dtdy G G dT dT lat lat Analytic prescription of mixed layer temperature gradient in 2d Rossby front test. g Grid structure dt Top to bottom temperature difference lat Latitude frontfractionalwidth hmlmax hmlmin ypseudo rossby_front_initialize_velocity Initial conditions for the 2D Rossby front test. Consistent with a linear equation of state, the system has a constant stratification below the mixed layer, stratified in temperature only. Isotherms are flat below the mixed layer and vertical within. Salinity is constant. The mixed layer has a half sine form so that there are no mixed layer or temperature gradients at the side walls.Below the mixed layer the potential temperature, $\theta(z)$, is given by \[ \theta(z) = \theta_0 - \Delta \theta \left( z + h_{ML} \right) \] where $ \theta_0 $ and $ \Delta \theta $ are external model parameters.The depth of the mixed layer, $H_{ML}$ is \[ h_{ML}(y) = h_{min} + \left( h_{max} - h_{min} \right) \cos{\pi y/L} \]. The temperature in mixed layer is given by the reference temperature at $z=h_{ML}$ so that \begin{eqnarray} \theta(y,z) = \theta_0 - \Delta \theta \left( z + h_{ML} \right) & \forall & z < h_{ML}(y) T.B.D. \end{eqnarray}