write(*,*)'are you using wheeled (0) or a tracked vehicle (1)'
read(*,*)ivehs
write(*,*)'input the thaw depth (cm)'
read(*,*)zthaw
trmlt = 1.
if((iwet.eq.3).and.((isoil.lt.13).or.(isoil.gt.15)))then
if(zthaw.le.2.5)then
trmlt = trmlt
else if((zthaw.gt.2.5).and.(zthaw.lt.15.))then
trmlt = trmlt*(2.379*(1/zthaw**2.)+0.619)
else if(zthaw.ge.15.)then
trmlt = 0.63*trmlt
end if
else if(iwet.eq.2) then
trmlt = 1.
else if(iwet.eq.1)then
write(*,*)'*********************************************************'
write(*,*)'for dry soils, the effect of thawing on traction is a'
write(*,*)'function of applied stress; these functions have not
write(*,*)yet been developed so a multiplier of 1.0 is used'
write(*,*)'*********************************************************'
end if
if((isoil.ge.13).and.(isoil.le.15))then
trmlt = 1.
write(*,*)'********************************************************'
write(*,*)'A multiplier of 1.0 is used for water, rock, or pavement'
write(*,*)'********************************************************'
end if
if(((isoil.ge.6).and.(isoil.le.12)).or.(isoil.eq.20))then
write(*,*)'*******************************************************'
write(*,*)'calculations are based on sandy soils'
write(*,*)'*******************************************************'
end if
if(ivehs.eq.1)then
write(*,*)'*******************************************'
write(*,*)'calculations are based on wheeled vehicles `
write(*,*)'*******************************************'
end if
return
end
C ************************************************************************
subroutine mrtmult(zthaw,svalid,isoil,iwet,rmult,ivehs)
c this subroutine calculates the motion resistance multipler
rmult = 1.
c
if(iwet.eq.3)then
if((iwet.eq.3).and.((isoil.lt.13).or.(isoil.gt.15)))then
if(zthaw.le.2.5)then
rmult = 1.
else if((zthaw.gt.2.5).and.(zthaw.lt.4.))then
34