C SLIP DUE TO TRACTION
if (totalr.le.trsnow) then
slipt = .2*totalr / trsnow
slipw = 1./(slipt + 1.) !RATIO OF VEHICLE SPEED TO WHEEL SPEED
else
slipw=0.0
end if
C WRITE TOTAL VEHICLE TRACTION AND RESISTANCE
write(*,210)hsnow,sigma0,slope
write(*,*)'ratio of vehicle speed to wheel speed: `,slipw
write(*,*)'
units are lbs:'
write(*,260)sumz,trsnow,totalr,totaldb
write(*,*)'
units are N:'
write(*,260)sumz,trsnow*4.4482,totalr*4.4482,totaldb*4.4482
write (*,*)'drawbar coefficient = ;(totaldb/gvw)
200
format(1x,a12)
210
format(/1x,"snow depth, cm.= ",f5.2,2x,"density, g/cc= ",f5.4,
& " slope= ",f6.2)
260
format(1x,"rut,cm.= ",f4.1,2x,"tr =",f8.1,2x,"rr= ",f8.1,
& 2x,"db= ",f8.1)
return
end
C ***************************************************************
subroutine thaw(trmlt,rmult,zsnow)
C
if(layer.eq.1.and.zsnow.lt.1.0/xcm2in)then
do i = 1, nunits
c
trmltn = (dowpb(i) + rtowpb(i))*trmlt
c
rtowpb(i) = rtowpb(i)*rmult
c
rtowt(i) = rtowt(i)*rmult
c
dowpb(i) = trmltn - rtowpb(i)
end do
C
else
C
write(*,*)'routine thaw is not implemented for:'
C
if(layer.ne.1.)then
C
write(*,*)'critical depth > 6 in.'
C
else
C
write(*,*)'snow cover > 1 cm'
C
end if
C
end if
write(*,*)'traction multiplier ',trmlt
write(*,*)'motion resistance multiplier',rmult
C
return
end
C********************************************************************
subroutine trmult(isoil,iwet,trmlt,zthaw,ivehs)
c this subroutine calculates the traction multiplier
33
Previous Page
