mom_coriolisadv mom_coriolisadv::coriolisadv_cs integer, parameter integer, parameter mom_coriolisadv::sadourny75_energy sadourny75_energy = 1 Enumeration values for Coriolis_Scheme. coradcalc coriolisadv_init integer, parameter integer, parameter mom_coriolisadv::arakawa_hsu90 arakawa_hsu90 = 2 Enumeration values for Coriolis_Scheme. coradcalc coriolisadv_init integer, parameter integer, parameter mom_coriolisadv::robust_enstro robust_enstro = 3 Enumeration values for Coriolis_Scheme. coradcalc coriolisadv_init integer, parameter integer, parameter mom_coriolisadv::sadourny75_enstro sadourny75_enstro = 4 Enumeration values for Coriolis_Scheme. coradcalc coriolisadv_init integer, parameter integer, parameter mom_coriolisadv::arakawa_lamb81 arakawa_lamb81 = 5 Enumeration values for Coriolis_Scheme. coradcalc coriolisadv_init integer, parameter integer, parameter mom_coriolisadv::al_blend al_blend = 6 Enumeration values for Coriolis_Scheme. coradcalc coriolisadv_init character *(20), parameter character*(20), parameter mom_coriolisadv::sadourny75_energy_string sadourny75_energy_string = "SADOURNY75_ENERGY" Enumeration values for Coriolis_Scheme. coriolisadv_init character *(20), parameter character*(20), parameter mom_coriolisadv::arakawa_hsu_string arakawa_hsu_string = "ARAKAWA_HSU90" Enumeration values for Coriolis_Scheme. coriolisadv_init character *(20), parameter character*(20), parameter mom_coriolisadv::robust_enstro_string robust_enstro_string = "ROBUST_ENSTRO" Enumeration values for Coriolis_Scheme. coriolisadv_init character *(20), parameter character*(20), parameter mom_coriolisadv::sadourny75_enstro_string sadourny75_enstro_string = "SADOURNY75_ENSTRO" Enumeration values for Coriolis_Scheme. coriolisadv_init character *(20), parameter character*(20), parameter mom_coriolisadv::arakawa_lamb_string arakawa_lamb_string = "ARAKAWA_LAMB81" Enumeration values for Coriolis_Scheme. coriolisadv_init character *(20), parameter character*(20), parameter mom_coriolisadv::al_blend_string al_blend_string = "ARAKAWA_LAMB_BLEND" Enumeration values for Coriolis_Scheme. coriolisadv_init integer, parameter integer, parameter mom_coriolisadv::ke_arakawa ke_arakawa = 10 Enumeration values for KE_Scheme. coriolisadv_init gradke integer, parameter integer, parameter mom_coriolisadv::ke_simple_gudonov ke_simple_gudonov = 11 Enumeration values for KE_Scheme. coriolisadv_init gradke integer, parameter integer, parameter mom_coriolisadv::ke_gudonov ke_gudonov = 12 Enumeration values for KE_Scheme. coriolisadv_init gradke character *(20), parameter character*(20), parameter mom_coriolisadv::ke_arakawa_string ke_arakawa_string = "KE_ARAKAWA" Enumeration values for KE_Scheme. coriolisadv_init character *(20), parameter character*(20), parameter mom_coriolisadv::ke_simple_gudonov_string ke_simple_gudonov_string = "KE_SIMPLE_GUDONOV" Enumeration values for KE_Scheme. coriolisadv_init character *(20), parameter character*(20), parameter mom_coriolisadv::ke_gudonov_string ke_gudonov_string = "KE_GUDONOV" Enumeration values for KE_Scheme. coriolisadv_init integer, parameter integer, parameter mom_coriolisadv::pv_adv_centered pv_adv_centered = 21 Enumeration values for PV_Adv_Scheme. coradcalc coriolisadv_init integer, parameter integer, parameter mom_coriolisadv::pv_adv_upwind1 pv_adv_upwind1 = 22 Enumeration values for PV_Adv_Scheme. coradcalc coriolisadv_init character *(20), parameter character*(20), parameter mom_coriolisadv::pv_adv_centered_string pv_adv_centered_string = "PV_ADV_CENTERED" Enumeration values for PV_Adv_Scheme. coriolisadv_init character *(20), parameter character*(20), parameter mom_coriolisadv::pv_adv_upwind1_string pv_adv_upwind1_string = "PV_ADV_UPWIND1" Enumeration values for PV_Adv_Scheme. coriolisadv_init subroutine, public subroutine, public mom_coriolisadv::coradcalc (u, v, h, uh, vh, CAu, CAv, OBC, AD, G, GV, US, CS) coradcalc u u v v h h uh uh vh vh CAu CAu CAv CAv OBC OBC AD AD G G GV GV US US CS CS Calculates the Coriolis and momentum advection contributions to the acceleration. g Ocen grid structure gv Vertical grid structure u Zonal velocity [L T-1 ~> m s-1] v Meridional velocity [L T-1 ~> m s-1] h Layer thickness [H ~> m or kg m-2] uh Zonal transport u*h*dy [H L2 T-1 ~> m3 s-1 or kg s-1] vh Meridional transport v*h*dx [H L2 T-1 ~> m3 s-1 or kg s-1] cau Zonal acceleration due to Coriolis and momentum advection [L T-2 ~> m s-2]. cav Meridional acceleration due to Coriolis and momentum advection [L T-2 ~> m s-2]. obc Open boundary control structure ad Storage for acceleration diagnostics us A dimensional unit scaling type cs Control structure for MOM_CoriolisAdv al_blend arakawa_hsu90 arakawa_lamb81 gradke mom_error_handler::mom_error mom_open_boundary::obc_direction_n pv_adv_centered pv_adv_upwind1 robust_enstro sadourny75_energy sadourny75_enstro mom_dynamics_split_rk2::step_mom_dyn_split_rk2 mom_dynamics_unsplit::step_mom_dyn_unsplit mom_dynamics_unsplit_rk2::step_mom_dyn_unsplit_rk2 subroutine subroutine mom_coriolisadv::gradke (u, v, h, KE, KEx, KEy, k, OBC, G, US, CS) gradke u u v v h h KE KE KEx KEx KEy KEy k k OBC OBC G G US US CS CS Calculates the acceleration due to the gradient of kinetic energy. g Ocen grid structure u Zonal velocity [L T-1 ~> m s-1] v Meridional velocity [L T-1 ~> m s-1] h Layer thickness [H ~> m or kg m-2] ke Kinetic energy per unit mass [L2 T-2 ~> m2 s-2] kex Zonal acceleration due to kinetic energy gradient [L T-2 ~> m s-2] key Meridional acceleration due to kinetic energy gradient [L T-2 ~> m s-2] k Layer number to calculate for obc Open boundary control structure us A dimensional unit scaling type cs Control structure for MOM_CoriolisAdv ke_arakawa ke_gudonov ke_simple_gudonov coradcalc subroutine, public subroutine, public mom_coriolisadv::coriolisadv_init (Time, G, GV, US, param_file, diag, AD, CS) coriolisadv_init Time Time G G GV GV US US param_file param_file diag diag AD AD CS CS Initializes the control structure for coriolisadv_cs. time Current model time g Ocean grid structure gv Vertical grid structure us A dimensional unit scaling type param_file Runtime parameter handles diag Diagnostics control structure ad Strorage for acceleration diagnostics cs Control structure fro MOM_CoriolisAdv al_blend al_blend_string arakawa_hsu90 arakawa_hsu_string arakawa_lamb81 arakawa_lamb_string ke_arakawa ke_arakawa_string ke_gudonov ke_gudonov_string ke_simple_gudonov ke_simple_gudonov_string mom_error_handler::mom_error mom_error_handler::mom_mesg pv_adv_centered pv_adv_centered_string pv_adv_upwind1 pv_adv_upwind1_string mom_diag_mediator::register_diag_field robust_enstro robust_enstro_string sadourny75_energy sadourny75_energy_string sadourny75_enstro sadourny75_enstro_string mom_string_functions::uppercase mom_dynamics_split_rk2::initialize_dyn_split_rk2 mom_dynamics_unsplit::initialize_dyn_unsplit mom_dynamics_unsplit_rk2::initialize_dyn_unsplit_rk2 subroutine, public subroutine, public mom_coriolisadv::coriolisadv_end (CS) coriolisadv_end CS CS Destructor for coriolisadv_cs. cs Control structure fro MOM_CoriolisAdv Accelerations due to the Coriolis force and momentum advection. This file contains the subroutine that calculates the time derivatives of the velocities due to Coriolis acceleration and momentum advection. This subroutine uses either a vorticity advection scheme from Arakawa and Hsu, Mon. Wea. Rev. 1990, or Sadourny's (JAS 1975) energy conserving scheme. Both have been modified to use general orthogonal coordinates as described in Arakawa and Lamb, Mon. Wea. Rev. 1981. Both schemes are second order accurate, and allow for vanishingly small layer thicknesses. The Arakawa and Hsu scheme globally conserves both total energy and potential enstrophy in the limit of nondivergent flow. Sadourny's energy conserving scheme conserves energy if the flow is nondivergent or centered difference thickness fluxes are used. A small fragment of the grid is shown below: j+1 x ^ x ^ x At x: q, CoriolisBu j+1 > o > o > At ^: v, CAv, vh j x ^ x ^ x At >: u, CAu, uh, a, b, c, d j > o > o > At o: h, KE j-1 x ^ x ^ x i-1 i i+1 At x & ^: i i+1 At > & o: The boundaries always run through q grid points (x).