Bug fix in physics/GWD/ugwp_driver_v0.F to prevent using uninitialize…

…d Bnv2 array in line 350
climbfuji · Dec 4, 2024 · baff7cb · baff7cb
1 parent a3f4d93
commit baff7cb
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/physics/GWD/ugwp_driver_v0.F b/physics/GWD/ugwp_driver_v0.F
@@ -349,7 +349,9 @@ SUBROUTINE GWDPS_V0(IM,  km,    imx, do_tofd,
           BVF2 = grav2 * RDZ * (VTK(I,K+1)-VTK(I,K))
      &                       / (VTK(I,K+1)+VTK(I,K))
           bnv2(i,k+1) = max( BVF2, bnv2min )
-          RI_N(I,K+1) = Bnv2(i,k)/SHR2        ! Richardson number consistent with BNV2
+          ! https://github.com/NCAR/ccpp-physics/issues/1103
+          !RI_N(I,K+1) = Bnv2(i,k)/SHR2        ! Richardson number consistent with BNV2
+          RI_N(I,K+1) = Bnv2(i,max(k,2))/SHR2        ! Richardson number consistent with BNV2
 !
 ! add here computation for Ktur and OGW-dissipation fro VE-GFS
 !