IcePenetSWFluxFrac (I,J) = 0. _d 0
IcePenetSWFluxFrac(I,J) = ZERO
IcePenetSWFluxFrac(I,J) =
& (1.0 - IcePenetSWFluxFrac(I,J))
& IcePenetSWFluxFrac(I,J)
^ IcePenetSWFluxFrac(I,J), ALB(I,J)
_RL IcePenetSWFluxFrac