1 SUBROUTINE snfrac_gfs(SNEQV,IVEG,SNCOVR)
16 integer,
intent(in) :: iveg
17 REAL,
intent(in) :: sneqv
18 REAL,
intent(out) :: sncovr
19 REAL salp,snup(13),rsnow
23 DATA snup /0.080, 0.080, 0.080, 0.080, 0.080, 0.080, &
24 & 0.040, 0.040, 0.040, 0.040, 0.025, 0.040, &
30 IF (sneqv < snup(iveg))
THEN
31 rsnow = sneqv/snup(iveg)
32 sncovr = 1. - (exp(-salp*rsnow) - rsnow*exp(-salp))
36 sncovr = max(0.,min(sncovr,1.))