diff --git a/README.md b/README.md index 5765206560..148dcf2a76 100644 --- a/README.md +++ b/README.md @@ -1,4 +1,4 @@ -MPAS-v8.0.1 +MPAS-v8.0.2 ==== The Model for Prediction Across Scales (MPAS) is a collaborative project for diff --git a/src/core_atmosphere/Registry.xml b/src/core_atmosphere/Registry.xml index 4e50bea204..42630440c5 100644 --- a/src/core_atmosphere/Registry.xml +++ b/src/core_atmosphere/Registry.xml @@ -1,5 +1,5 @@ - + diff --git a/src/core_atmosphere/physics/mpas_atmphys_driver_cloudiness.F b/src/core_atmosphere/physics/mpas_atmphys_driver_cloudiness.F index 9e444ae2f7..93e833b328 100644 --- a/src/core_atmosphere/physics/mpas_atmphys_driver_cloudiness.F +++ b/src/core_atmosphere/physics/mpas_atmphys_driver_cloudiness.F @@ -51,7 +51,10 @@ module mpas_atmphys_driver_cloudiness ! Laura D. Fowler (laura@ucar.edu) / 2016-07-05. ! * since we removed the local variable radt_cld_scheme from mpas_atmphys_vars.F, now defines radt_cld_scheme ! as a pointer to config_radt_cld_scheme. -! Laura D. Fowler (laura@ucar.edu) / 2917-02-16. +! Laura D. Fowler (laura@ucar.edu) / 2017-02-16. +! * this is a bug fix. dx_p is converted from meters to kilometers prior to calling the thompson parameterization +! of the cloud fraction. +! Laura D. Fowler (laura@ucar.edu) / 2024-03-23. contains @@ -119,6 +122,8 @@ subroutine cloudiness_from_MPAS(configs,mesh,diag_physics,sfc_input,its,ite) do j = jts,jte do i = its,ite dx_p(i,j) = len_disp / meshDensity(i)**0.25 + !conversion of dx_p from meters to kilometers. + dx_p(i,j) = dx_p(i,j)*0.001 xland_p(i,j) = xland(i) enddo diff --git a/src/core_atmosphere/physics/physics_wrf/module_mp_thompson.F b/src/core_atmosphere/physics/physics_wrf/module_mp_thompson.F index b5b9189b95..9abf279048 100644 --- a/src/core_atmosphere/physics/physics_wrf/module_mp_thompson.F +++ b/src/core_atmosphere/physics/physics_wrf/module_mp_thompson.F @@ -1593,6 +1593,8 @@ subroutine mp_thompson (qv1d, qc1d, qi1d, qr1d, qs1d, qg1d, ni1d, & smod(k) = 0. smoe(k) = 0. smof(k) = 0. + mvd_r(k) = 0. + mvd_c(k) = 0. enddo !+---+-----------------------------------------------------------------+ diff --git a/src/core_init_atmosphere/Registry.xml b/src/core_init_atmosphere/Registry.xml index e40f5c3722..9f5abc147e 100644 --- a/src/core_init_atmosphere/Registry.xml +++ b/src/core_init_atmosphere/Registry.xml @@ -1,5 +1,5 @@ - + diff --git a/src/core_landice/Registry.xml b/src/core_landice/Registry.xml index 42e2a99254..91db32ee9a 100644 --- a/src/core_landice/Registry.xml +++ b/src/core_landice/Registry.xml @@ -1,5 +1,5 @@ - + diff --git a/src/core_ocean/Registry.xml b/src/core_ocean/Registry.xml index c084373916..a5f42e4d4a 100644 --- a/src/core_ocean/Registry.xml +++ b/src/core_ocean/Registry.xml @@ -1,5 +1,5 @@ - + - + - + diff --git a/src/core_test/Registry.xml b/src/core_test/Registry.xml index f2ffbbf18d..102cce6de3 100644 --- a/src/core_test/Registry.xml +++ b/src/core_test/Registry.xml @@ -1,5 +1,5 @@ - + diff --git a/src/operators/mpas_geometry_utils.F b/src/operators/mpas_geometry_utils.F index ba9d49522a..1a50a02e60 100644 --- a/src/operators/mpas_geometry_utils.F +++ b/src/operators/mpas_geometry_utils.F @@ -42,9 +42,10 @@ real (kind=RKIND) function mpas_sphere_angle(ax, ay, az, bx, by, bz, cx, cy, cz) real (kind=RKIND) :: s ! Semiperimeter of the triangle real (kind=RKIND) :: sin_angle - a = acos(max(min(bx*cx + by*cy + bz*cz,1.0_RKIND),-1.0_RKIND)) ! Eqn. (3) - b = acos(max(min(ax*cx + ay*cy + az*cz,1.0_RKIND),-1.0_RKIND)) ! Eqn. (2) - c = acos(max(min(ax*bx + ay*by + az*bz,1.0_RKIND),-1.0_RKIND)) ! Eqn. (1) + + a = mpas_arc_length(bx, by, bz, cx, cy, cz) + b = mpas_arc_length(ax, ay, az, cx, cy, cz) + c = mpas_arc_length(ax, ay, az, bx, by, bz) ABx = bx - ax ABy = by - ay