common/m5/iveh,nunits,pmax,gclear,gvw
if (iveh.eq.1) then
C CALCULATE VALUES FOR TRACKED VEHICLES
do i=1,nunits
nwhls(i)=2
!ALL TRACKED VEHICLES HAVE 2 TRACKS PER AXLE
ip(i)=1
!ALL TRACKS ARE DRIVEN
area(i)=width(i)*tlength(i)
pload(i)=gvw/nunits/nwhls(i)
pres1(i)=pload(i)/area(i)
c
write(*,*)'pressure',pres1(i)
end do
else
C CALCULATE GROUND CONTACT PRESSURES FOR WHEELED VEHICLES
do i=1,nunits
pres1(i)=pload(i)/area(i)
c
write(*,*)'pressure1',pres1(i)
end do
end if
C
CALCULATE THE MAXIMUM PRESSURE
pmax=pres1(1)
do i=1,nunits
if (pres1(i).gt.pmax) pmax=pres1(i)
end do
return
end
C **********************************************************************
subroutine density(hsnow,sigma0,sigmaf,zmax)
C THIS SUBROUTINE CALCULATES THE FINAL DENSITY AND SINKAGE
common/m5/iveh,nunits,pmax,gclear,gvw
C COMPUTE FINAL DENSITY DATA
sigmaf=.5
if(pmax.gt.30.5) sigmaf=.55
if(pmax.gt.50.8) sigmaf=.60
if(pmax.gt.101.5) sigmaf=.65
C COMPUTE SINKAGE
zmax= hsnow * (1. - (sigma0/sigmaf))
if(zmax.lt.0.)zmax=0.
return
end
C *******************************************************************
subroutine rterrain(hsnow,sigma0,sigmaf,zmax,sumz,rsnow)
27
Previous Page
