C
Hard surface
IBASE = 1
C ELSE IF( DTHAW .LE. 0.0 )THEN
ELSE IF( DFREEZ .GT. 0.0 .AND. DTHAW .LE. 0.0 )THEN
C
Frozen soil surface
IBASE = 2
ELSE
C
Unfrozen soil
IBASE = 3
END IF
C
SIGMA0 = SIGMA
C
C ... Find maximum contact pressure
PMAX = 0.0
DO 10 I=1, NAMBLY
IF( NVEH(I) .EQ. WHEEL )THEN
AREA1 = GCA(I,JPSI)
XNWORT = NWHL(I)
ELSE
AREA1 = GCA(I,JPSI) / 2.0
XNWORT = 2.0
END IF
C
Weight on single traction element
WGHT1 = WGHT(I) / XNWORT
C
Maximum contact pressure of all traction elements
PMAX = AMAX1( PMAX, WGHT1/AREA1 )
C
10 CONTINUE
C
C ... Compute final density sigmaf
IF( PMAX .GT. 101.5 ) THEN
SIGMAF = 0.65
ELSE IF( PMAX .GT. 50.8 ) THEN
SIGMAF = 0.60
ELSE IF( PMAX .GT. 30.5 ) THEN
SIGMAF = 0.55
ELSE
SIGMAF = 0.5
END IF
C
FLAG = 1.0
SUMZ = 0.0
C
C Initial density and final density information
ZMAX = AMAX1( 0.0, DSNOW * (1. - (SIGMA0/SIGMAF)) )
IF(KIV(3).GT.YES)WRITE(LUN1,601)PMAX,SIGMA0,SIGMAF,ZMAX
601 FORMAT(` Max elmnt cntct pres PMAX =',F8.4/
&
` Inp snow density
SIGMA0 =',F8.4/
&
` Final snow density SIGMAF =',F8.4/
&
` Max avail sinkage
ZMAX =',F8.4)
C
DO 20 I = 1, NAMBLY
43
Previous Page
