Reads the parameters unique to this module.
32 type(ocean_grid_type),
intent(in) :: G
33 type(verticalGrid_type),
intent(in) :: GV
34 type(unit_scale_type),
intent(in) :: US
35 type(param_file_type),
intent(in) :: param_file
36 real,
intent(out) :: S_ref
37 real,
intent(out) :: dSdz
38 real,
intent(out) :: delta_S
39 real,
intent(out) :: dSdx
40 real,
intent(out) :: T_ref
41 real,
intent(out) :: dTdz
42 real,
intent(out) :: delta_T
43 real,
intent(out) :: dTdx
44 real,
intent(out) :: L_zone
45 logical,
optional,
intent(in) :: just_read_params
50 just_read = .false. ;
if (
present(just_read_params)) just_read = just_read_params
53 call log_version(param_file, mdl, version,
'Initialization of an analytic baroclinic zone')
54 call openparameterblock(param_file,
'BCZIC')
55 call get_param(param_file, mdl,
"S_REF", s_ref,
'Reference salinity', units=
'ppt', &
56 default=35., do_not_log=just_read)
57 call get_param(param_file, mdl,
"DSDZ", dsdz,
'Salinity stratification', &
58 units=
'ppt/m', default=0.0, scale=us%Z_to_m, do_not_log=just_read)
59 call get_param(param_file, mdl,
"DELTA_S",delta_s,
'Salinity difference across baroclinic zone', &
60 units=
'ppt', default=0.0, do_not_log=just_read)
61 call get_param(param_file, mdl,
"DSDX",dsdx,
'Meridional salinity difference', &
62 units=
'ppt/'//trim(g%x_axis_units), default=0.0, do_not_log=just_read)
63 call get_param(param_file, mdl,
"T_REF",t_ref,
'Reference temperature',units=
'C', &
64 default=10., do_not_log=just_read)
65 call get_param(param_file, mdl,
"DTDZ", dtdz,
'Temperature stratification', &
66 units=
'C/m', default=0.0, scale=us%Z_to_m, do_not_log=just_read)
67 call get_param(param_file, mdl,
"DELTA_T",delta_t,
'Temperature difference across baroclinic zone', &
68 units=
'C', default=0.0, do_not_log=just_read)
69 call get_param(param_file, mdl,
"DTDX",dtdx,
'Meridional temperature difference', &
70 units=
'C/'//trim(g%x_axis_units), default=0.0, do_not_log=just_read)
71 call get_param(param_file, mdl,
"L_ZONE",l_zone,
'Width of baroclinic zone', &
72 units=g%x_axis_units, default=0.5*g%len_lat, do_not_log=just_read)
73 call closeparameterblock(param_file)