42 real(RP),
private :: OCEAN_PHY_ROUGHNESS_miller92_CM0 = 1.0e-3_rp
43 real(RP),
private :: OCEAN_PHY_ROUGHNESS_miller92_Z0MI = 0.0e-0_rp
44 real(RP),
private :: OCEAN_PHY_ROUGHNESS_miller92_Z0MR = 1.8e-2_rp
45 real(RP),
private :: OCEAN_PHY_ROUGHNESS_miller92_Z0MS = 1.1e-1_rp
46 real(RP),
private :: OCEAN_PHY_ROUGHNESS_miller92_Z0HI = 1.4e-5_rp
47 real(RP),
private :: OCEAN_PHY_ROUGHNESS_miller92_Z0HR = 0.0e-0_rp
48 real(RP),
private :: OCEAN_PHY_ROUGHNESS_miller92_Z0HS = 4.0e-1_rp
49 real(RP),
private :: OCEAN_PHY_ROUGHNESS_miller92_Z0EI = 1.3e-4_rp
50 real(RP),
private :: OCEAN_PHY_ROUGHNESS_miller92_Z0ER = 0.0e-0_rp
51 real(RP),
private :: OCEAN_PHY_ROUGHNESS_miller92_Z0ES = 6.2e-1_rp
63 namelist / param_ocean_phy_roughness_miller92 / &
64 ocean_phy_roughness_miller92_cm0, &
65 ocean_phy_roughness_miller92_z0mi, &
66 ocean_phy_roughness_miller92_z0mr, &
67 ocean_phy_roughness_miller92_z0ms, &
68 ocean_phy_roughness_miller92_z0hi, &
69 ocean_phy_roughness_miller92_z0hr, &
70 ocean_phy_roughness_miller92_z0hs, &
71 ocean_phy_roughness_miller92_z0ei, &
72 ocean_phy_roughness_miller92_z0er, &
73 ocean_phy_roughness_miller92_z0es
83 read(
io_fid_conf,nml=param_ocean_phy_roughness_miller92,iostat=ierr)
85 log_info(
"OCEAN_PHY_ROUGHNESS_miller92_setup",*)
'Not found namelist. Default used.'
86 elseif( ierr > 0 )
then
87 log_error(
"OCEAN_PHY_ROUGHNESS_miller92_setup",*)
'Not appropriate names in namelist PARAM_OCEAN_PHY_ROUGHNESS_MILLER92. Check!'
90 log_nml(param_ocean_phy_roughness_miller92)
111 integer,
intent(in) :: oia, ois, oie
112 integer,
intent(in) :: oja, ojs, oje
113 real(rp),
intent(in) :: uabs(oia,oja)
114 real(rp),
intent(out) :: z0m (oia,oja)
115 real(rp),
intent(out) :: z0h (oia,oja)
116 real(rp),
intent(out) :: z0e (oia,oja)
128 z0m(i,j) = max( ocean_phy_roughness_miller92_z0mi &
129 + ocean_phy_roughness_miller92_z0mr / grav * ustar * ustar &
132 z0h(i,j) = max( ocean_phy_roughness_miller92_z0hi &
133 + ocean_phy_roughness_miller92_z0hr / grav * ustar * ustar &
136 z0e(i,j) = max( ocean_phy_roughness_miller92_z0ei &
137 + ocean_phy_roughness_miller92_z0er / grav * ustar * ustar &