diff --git a/mo_grotz.f90 b/mo_grotz.f90 index 996b73f62917b5a77398267fe3507f0b5c58cc3c..21edfcb581561f65cf22fbfde10541a0c9f02b76 100644 --- a/mo_grotz.f90 +++ b/mo_grotz.f90 @@ -227,8 +227,8 @@ CONTAINS CALL snow_precip (m_snow, H_abs_snow, thick_snow, dt, liquid_precip, T2m) ELSE IF(MAX(liquid_precip, solid_precip)>0.0_wp .AND. N_Active==1) THEN !Precip in water layer CALL snow_precip_0 (H_abs(1), S_abs(1), m(1), T(1), dt, liquid_precip, T2m) - ELSE IF(MIN(liquid_precip, solid_precip)<0.0_wp .AND. thick_snow>0._wp .AND. N_Active>1) THEN !Erosion of snow layer - CALL snow_erosion (m_snow, H_abs_snow, thick_snow, dt, liquid_precip, T2m) + ELSE IF(MIN(liquid_precip, solid_precip)<0.0_wp .AND. thick_snow>0._wp .AND. N_Active>1) THEN !Erosion of snow layer + CALL snow_erosion (m_snow, H_abs_snow, thick_snow, dt, liquid_precip, T2m) END IF ELSE IF (precip_flag==0) THEN @@ -334,8 +334,8 @@ CONTAINS grav_drain = grav_drain / time_out end if - grav_heat_flux_down = grav_heat_flux_down / time_out - grav_heat_flux_up = grav_heat_flux_up / time_out + grav_heat_flux_down = grav_heat_flux_down / time_out + grav_heat_flux_up = grav_heat_flux_up / time_out !Calling standard output CALL output(Nlayer, T, psi_s, psi_l, thick, S_bu, ray, format_T, format_psi, & @@ -385,8 +385,8 @@ CONTAINS grav_drain = 0.0_wp grav_salt = 0.0_wp grav_temp = 0.0_wp - grav_heat_flux_down = 0._wp - grav_heat_flux_up = 0._wp + grav_heat_flux_down = 0._wp + grav_heat_flux_up = 0._wp melt_thick_output(:) = 0._wp n_time_out = 0