IF( (facz.GE.dd).AND.(facz.LT.dd20) ) THEN
plumek(i) = one - min(one,(facz-dd)/(dd20-dd))
tempN = facz*tempN !raise to the Npower+1
plumek(i) = one - min(one,two/dd20*facz-tempN/tempN20)
& three/dd20*facz - three/(dd20*dd20)*facz*facz
_RL facz, dd, dd20
facz = abs(fact*plumek(i))
IF (SPDepth(i).GT.facz) THEN
& (exp(facz/dd)-one)/
facz = fact*swdk(i)
IF ( facz .LT. -200. _d 0 ) THEN
swdk(i) = rfac(jwtype) * exp( facz/a1(jwtype) )
& + (1. _d 0 - rfac(jwtype)) * exp( facz/a2(jwtype) )
_RL facz