Output updraft value of convective cloud condensate

hafs-community · Jul 16, 2024 · e2107d2 · e2107d2
1 parent c700e1b
commit e2107d2
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 2 deletions.
diff --git a/physics/CONV/SAMF/samfdeepcnv.f b/physics/CONV/SAMF/samfdeepcnv.f
@@ -3430,7 +3430,10 @@ subroutine samfdeepcnv_run (im,km,first_time_step,restart,        &
         do i = 1, im
           if (cnvflg(i) .and. rn(i) > 0.) then
             if (k >= kbcon(i) .and. k < ktcon(i)) then
-              cnvw(i,k) = cnvwt(i,k) * xmb(i) * dt2
+               cnvw(i,k) = cnvwt(i,k) * xmb(i) * dt2
+               tem=max(sigmaout(i,k),0.)
+               tem1=min(tem,1.0)
+               cnvw(i,k)=cnvw(i,k)*tem1
             endif
           endif
         enddo

diff --git a/physics/CONV/SAMF/samfshalcnv.f b/physics/CONV/SAMF/samfshalcnv.f
@@ -2405,7 +2405,10 @@ subroutine samfshalcnv_run(im,km,itc,ntc,cliq,cp,cvap,            &
           if (cnvflg(i)) then
             if (k >= kbcon(i) .and. k < ktcon(i)) then
               cnvw(i,k) = cnvwt(i,k) * xmb(i) * dt2
-            endif
+              tem=max(sigmaout(i,k),0.)
+              tem1=min(tem,1.0)
+              cnvw(i,k)=cnvw(i,k)*tem1
+           endif
           endif
         enddo
       enddo