c
write(*,*)'pressure in kPa',pr
c
write (*,*)trcn1(i)
end if
C HARD PACKED SNOW
if (idist.eq.6) then
if((ibase.eq.1).or.(ibase.eq.2).or.(ibase.eq.3))then
trcn1(i)=0.321*(pr)**0.97
end if
if(ibase.eq.5)then
trcn1(i)=0.321*(pr)**0.97
end if
end if
C SNOW ON ICE BASE
if((ibase.eq.4).and.(idist.eq.7))then
trcn1(i)=0.127*(pr)**1.06
else if((ibase.eq.4).and.(idist.eq.6))then
trcn1(i)=0.321*(pr)**0.97
end if
C ICE SURFACE
if ((ibase.eq.4).and.(hsnow.le.1.0))then
trcn1(i)=0.1*pr
end if
trcn1(i)=(trcn1(i)/6.8948)*area(i)
!CONVERT TO lbs
C NON-DRIVEN WHEEL
if(ip(i).eq.0) trcn1(i)=0.0
C SUM TRACTION FOR EACH ELEMENT
trsnow = trsnow + trcn1(i)*nwhls(i)
end do
return
end
C ******************************************************************
subroutine totalsnow(rsnow,trsnow,reshard,hsnow,sigma0,sumz)
common/m1/nwhls(10),rads(10),width(10),pload(10),area(10),aa(10)
common/m2/ip(10),pres1(10),tlength(10)
common/m3/xl2(10),cpffg(10)
common/m5/iveh,nunits,pmax,gclear,gvw
common/m6/def,sec
common/m7/ibase,rci,idist,slope
totalr = reshard + rsnow
totaldb = trsnow - totalr
32