call thaw(trmlt,rmult,zsnow)
else
call freeze(isoil)
end if
end if
write (*,*)'enter y to run another condition'
read(*,300)ans
300
format (a1)
if ((ans.eq."Y") .or. (ans.eq."y")) then
goto 1000
end if
c
call exit(0)
end
C ########################################################################
subroutine vehdata
C THIS SUBROUTINE READS DATA FROM THE VEHICLE DATA FILE
common/m1/nwhls(10),rads(10),width(10),pload(10),area(10),aa(10)
common/m2/ip(10),pres1(10),tlength(10)
common/m5/iveh,nunits,pmax,gclear,gvw
common/m6/def,sec
character*12 mfile
write(*,*)'enter the vehicle data file name'
read(*,10) mfile
open(21,status='old',file=mfile,form='formatted')
C READ IN VEHICLE DATA
read(21,100)title
read(21,110)iveh
!WHEELED=0, TRACKED=1
if(iveh.eq.0) then
C
READ WHEEL VEHICLE DATA - ENGLISH UNITS
read(21,115)nunits
!NO. AXLES,
read(21,220)gclear
!GROUND CLEARANCE
read(21,110)(nwhls(i),i=1,nunits)
!NUMBER OF WHEELS PER AXLE
read(21,110)(ip(i),i=1,nunits)
!AXLE, O=NO POWER ,1=POWERED
read(21,120)(rads(i),i=1,nunits)
!RADIUS OF EACH WHEEL, IN.
read(21,120)(width(i),i=1,nunits)
!WIDTH OF EACH WHEEL, IN.
read(21,120)(pload(i),i=1,nunits)
!LOAD ON EACH WHEEL, LBS.
read(21,120)(area(i),i=1,nunits)
!TIRE CONTACT AREA, SQ. IN.
read(21,*)sec,def
!UNLOADED SECTION HEIGHT, LOADED
TIRE DEFLECTION
24
Previous Page
