bulk_snow_module.f90

Go to the documentation of this file.

module bulk_snow_module

 implicit none

 private

 public calc_density

contains


 subroutine calc_density(lensfc, landmask, swe, snd, rank, density)

       implicit none

       integer, intent(in) :: lensfc, rank
       integer, intent(in) :: landmask(lensfc)
       real, intent(in)    :: swe(lensfc), snd(lensfc)
       real, intent(out)   :: density(lensfc)

       real :: dens_mean
       integer :: n

        ! density = swe/snd
        do n =1,lensfc
                if (snd(n) > 0.001 ) then
                        density(n) = swe(n)/snd(n)
                else 
                        density(n)=0.1 
                endif 
        enddo

        where (density < 0.0001) density = 0.1

        ! calculate mean density of snow over land
        if (count(landmask==2) > 0) then
                ! mean density over snow-covered land
                dens_mean = sum(density, mask = (landmask==2 )) &
                         / count(landmask==2)
                print *, 'mean density on rank ', rank,': ', dens_mean
        else
                dens_mean = 0.1  ! default value if have no snow in tile
                print *, 'no snow on rank ', rank, ' using default density ', dens_mean
        endif

        ! for grid cells with no valid density, fill in the average snodens
        where( swe <= 0.001 ) density = dens_mean

 end subroutine calc_density

end module bulk_snow_module