mass_balance__snow.F90

Go to the documentation of this file.

module mass_balance__snow
 
  use iso_c_binding, only        : c_short, c_int, c_float, c_double
  implicit none
 
contains
 
  elemental subroutine calculate_snow_mass_balance( snow_storage,         &
                                                    potential_snowmelt,   &
                                                    snowmelt,             &
                                                    net_snowfall )
 
    real (c_float), intent(inout)      :: snow_storage
    real (c_float), intent(inout)      :: snowmelt
    real (c_float), intent(in)         :: potential_snowmelt
    real (c_float), intent(in)         :: net_snowfall
 
 
    snow_storage = max(0.0_c_float, snow_storage + net_snowfall )
 
    if( snow_storage > potential_snowmelt ) then
 
      snowmelt = potential_snowmelt
      snow_storage = snow_storage - snowmelt
 
    else   ! not enough snowcover to satisfy the amount that *could* melt
 
      snowmelt = snow_storage
      snow_storage = 0.0_c_float
 
    end if
 
  end subroutine calculate_snow_mass_balance
 
end module mass_balance__snow