C
SIGMA0
Input snow specific gravity
C
SIGMAF
Final snow specific gravity
C
SLIP
Vehicle wheel slip
C
SOILZ
Sinkage in soil of current traction element (in)
C
SUMZ
Total vehicle sinkage
C
TFOW
Temporary traction coefficient
C
TRCN1
Maximum traction for i-th tire (lb)
C
XL2
Snow contact length for i-th suspension element (in)
C
XNWORT
Number of wheels or tracks on suspension assembly
C
XWID1
Width of one traction element (in)
C
WGHT1
Weight on single traction element (wheel or track) (lb)
C
X2KPA
Conversion factor from psi to KPA
C
ZMAX
Maximum sinkage in snow (in)
C
IMPLICIT NONE
INCLUDE `nrmmdefs.inc'
INCLUDE `nrmmdrvc.inc'
INCLUDE `nrmmcntc.inc'
INCLUDE `nrmmvehc.inc'
INCLUDE `nrmmvppc.inc'
INCLUDE `nrmmscnc.inc'
INCLUDE `nrmmterc.inc'
INCLUDE `nrmmtppc.inc'
C
INTEGER I
REAL ACOS
REAL AA, AREA1, DEFSEC, FLAG
REAL PMAX, PN, PRES1, RADSI, RRES1, RTOW
REAL SIGMA0, SIGMAF, SLIP, SOILZ, SUMZ
REAL TFOW, TRCN1, WGHT1, X2KPA, XL2
REAL XNWORT, XWID1, ZMAX
C
REAL RHARD
EXTERNAL RHARD
C
C Conversion factor for psi to KPA
PARAMETER ( X2KPA = XM2IN * XM2IN / 1000.0 / XN2LB )
C
C Diagnostic printout of input data:
IF( KIV(3) .GT. YES )THEN
CALL DIAGIN(`!0***Routine CRLSNOW(IV3) `// `CRREL snow model')
CALL DIAG( `CL DIAW DSNOW DTHAW GCA' )
CALL DIAG( `ITSURF ITUT NAMBLY NVEH NWHL' )
CALL DIAG( `SECTW SIGMA TRAKLN TRAKWD WGHT' )
CALL DIAG( `WT WTE' )
END IF
C
CRLSNOW = .TRUE.
C
C ... Set Snow model flag, IBASE
IF( ITSURF.EQ.3 )THEN
C
Ice
IBASE = 4
ELSE IF( ITUT.GE.11.AND.ITUT.LE.13 )THEN
42
Previous Page
