C Output:
C MRMULT
Motion resistance multiplier
C
MRMULT = 1.0
C
THAWCM = DTHAW / XCM2IN
C
The following 3 lines handle cases for which THAWCM < 2.5 cm
IF( THAWCM .LE. 2.5 )THEN
MRMULT = 1.0
ELSE IF( (THAWCM .GT. 2.5) .AND. (THAWCM .LT. 4.0)) THEN
MRMULT = 2.883 * THAWCM - 6.056
ELSE IF( (THAWCM .GE. 4.0) .AND. (THAWCM .LT. 8.0)) THEN
MRMULT = -0.22*THAWCM*THAWCM + 3.54*THAWCM - 5.24
ELSE IF( (THAWCM .GE. 8.) .AND. (THAWCM .LT. 12.0)) THEN
MRMULT = 0.225*THAWCM + 7.2167
ELSE
MRMULT = 10.0
END IF
IF( NVEH(I).EQ. TRACK )
& CALL ERROR( `MRMULT',
& `Multiplier may not be valid for TRACKED elements' )
IF( .NOT.SVALID(KUSCS) ) THEN
WRITE(ERRMSG,9901)KUSCS,ZUSCS(KUSCS)(1:LENSTR(ZUSCS(KUSCS))),
&
NVEH(I),ZVTYPE
CALL ERROR( `MRMULT', ERRMSG )
END IF
C
RETURN
END
52
Previous Page
