Fix in unit conversion in computing srflag for GFDL-MP

lisa-bengtsson · Mar 29, 2019 · 32d2b1d · 32d2b1d
1 parent 384ee3d
commit 32d2b1d
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/gfsphysics/GFS_layer/GFS_physics_driver.F90 b/gfsphysics/GFS_layer/GFS_physics_driver.F90
@@ -3775,9 +3775,9 @@ subroutine GFS_physics_driver                                 &
 !          if ((snow0(i,1)+ice0(i,1)+graupel0(i,1)+csnow) > 0.0) then
 !            Sfcprop%srflag(i) = 1.                   ! clu: set srflag to 'snow' (i.e. 1)
 !          endif
-           total_precip = snow0(i,1)+ice0(i,1)+graupel0(i,1)+rain0(i,1)+Diag%rainc(i)
-           if (total_precip*tem > rainmin) then
-             Sfcprop%srflag(i) = (snow0(i,1)+ice0(i,1)+graupel0(i,1)+csnow)/total_precip
+           total_precip = (snow0(i,1)+ice0(i,1)+graupel0(i,1)+rain0(i,1))*tem+Diag%rainc(i)
+           if (total_precip > rainmin) then
+             Sfcprop%srflag(i) = ((snow0(i,1)+ice0(i,1)+graupel0(i,1))*tem+csnow)/total_precip
            endif
         enddo
       elseif( .not. Model%cal_pre) then