rmult = rmult*(2.883*zthaw - 6.056)
else if((zthaw.ge.4.).and.(zthaw.le.8.))then
rmult = rmult*(-0.22*zthaw**2. + 3.54*zthaw - 5.24)
else if((zthaw.gt.8.).and.(zthaw.lt.12.))then
rmult = rmult*(0.225*zthaw + 7.2167)
else if(zthaw.ge.12)then
rmult = rmult*10.0
end if
end if
return
end
C ****************************************************************************
subroutine soil(isoil,iwet)
c this subroutine asks the USCS soil type and wetness
write(*,*)'choose an USCS soil type code:'
write(*,*)' 1 = SW'
write(*,*)' 2 = SP'
write(*,*)' 3 = SM'
write(*,*)' 4 = SC'
write(*,*)' 5 = SMSC'
write(*,*)' 6 = CL'
write(*,*)' 7 = ML'
write(*,*)' 8 = CLML'
write(*,*)' 9 = CH'
write(*,*)'10 = MH'
write(*,*)'11 = OL'
write(*,*)'12 = OH'
write(*,*)'13 = WATER'
write(*,*)'14 = PAVEMENT'
write(*,*)'15 = ROCK'
write(*,*)'16 = GW'
write(*,*)'17 = GP'
write(*,*)'18 = GM'
write(*,*)'19 = GC'
write(*,*)'20 = PT'
read(*,*)isoil
write(*,*)'input the soil wetness'
write(*,*)' 1 = dry'
write(*,*)' 2 = moist'
write(*,*)' 3 = wet (water content is above liquid limit)'
read(*,*)iwet
return
end
C **************************************************************************
subroutine freeze
write(*,*)'FROZEN GROUND MENU'
write(*,*)'1 = input the frost depth'
write(*,*)'2 = input the vehicle class'
35
Previous Page
