      SUBROUTINE VDPL(QTTL,SY,SZ,D1,X,DX,DXA,XS,XCH,DP,DPA,SDEPL,INT)
      COMMON VDM(41),UT,VDM1(23),FP,VDM2(8),Z0,VDM3(13),
     +IPO,I2MC,IMA,VDM4(11),IM,IR,IL,VDM5(3),K33,K42
      DIMENSION FKMS(6),Z0T(12),DDD(10),DDL(10),DDBR(10),
     +ARD(10),PT(6,3),ITT(12),Y(10)
      DATA PT/.05,.05,.1,.1,.15,.15,.05,.1,.15,.2,.25,.3,.1,
     +.15,.2,.25,.3,.35/
      DATA ITT/3,2,3,1,3,2,3,2,2,3,3,0/
      DATA FKMS/.9,.8,.6,.4,.2,.05/
      DATA Z0T/100.,.03,100.,.005,100.,10.,100.,100.,100.,
     +100.,100.,100./
      GOTO(101,102,103,104)INT
 101  ITI=ITT(IL)
      IF(K33.EQ.0.AND.IL.LT.12.AND.IM.LT.7) FP=PT(IM,ITI)
      IF (K42.EQ.0.AND.IL.LT.12) Z0=Z0T(IL)
      UZ=UT*.005**FP
      USTR=FP*FKMS(IM)*UZ
      DVF2=USTR/UZ
      DVF1=DVF2*USTR
      BINV=0.06*(19600.*USTR*Z0)**.45
      DEPV=DVF1/(1.+DVF2*BINV)
      WRITE(*,50)
 50   FORMAT(/5X,'FP',7X,'FKMS',7X,'UT',7X,'USTR',
     +6X,'BINV',6X,'DEPV',6X,'QTTL')
      WRITE(*,60) FP,FKMS(IM),UT,USTR,BINV,DEPV,QTTL
 60   FORMAT(1X,1P7E10.3)
      DEPV60=DEPV*60.
      RETURN
C     ENTRY VDPL1(SY,SZ,D1)
 102  DO 10 I=1,10
 10   ARD(I)=0.
      DDMX=QTTL/(188.496*SY*SZ*UT)
      DDD(1)=D1/10.
      DDDL=(ALOG(DDMX)-ALOG(DDD(1)))/9.
      DDL(1)=ALOG(DDD(1))
      DO 20 I=2,10
      DDL(I)=DDL(I-1)+DDDL
      DDD(I)=EXP(DDL(I))
 20   DDBR(I-1)=(DDD(I-1)+DDD(I))/2.
      IF (IPO.NE.4) RETURN
      IF (IMA.GT.0) GO TO 35
      WRITE(*,30) (DDD(I),I=1,10)
 30   FORMAT(/6X,'DOSAGE CONTOURS',10F5.0)
      WRITE(*,31)
 31   FORMAT(/7X,'X',7X,'DP',4X,'CONTOUR HALF-WIDTHS')
      RETURN
 35   WRITE(*,37) (DDD(I),I=1,10)
 37   FORMAT(/7X,'CONC. CONTOURS',10F5.0)
      WRITE(*,38)
 38   FORMAT(/7X,'X',7X,'CP',4X,'CONTOUR HALF-WIDTHS')
      RETURN
C     ENTRY VDPL2(X,DX,DXA,XS,XCH,DP,DPA,SDEPL)
 103  DP=DP*SDEPL
      DPA=DPA*SDEPL
      IF (DPA.LE.0) RETURN
      DPLA=ALOG(DPA)
      DO 70 IY=1,10
      ARG=DPLA-DDL(IY)
      IF (ARG.LT.0.) GO TO 80
      Y(IY)=1.41421*SY*SQRT(ARG)
      IF (X.NE.XCH) ARD(IY)=ARD(IY)+Y(IY)*DXA
      IF (X.EQ.XCH) ARD(IY)=ARD(IY)-Y(IY)*(XS+DX-XCH)
 70   CONTINUE
      IY=11
 80   QD=0.
      IF (IPO.GT.0) WRITE(*,'(F5.2)') SDEPL
      IY=IY-1
      DO 90 I=1,9
 90   QD=QD+DDBR(I)*(ARD(I)-ARD(I+1))*DEPV60
      SDEPL=(QTTL-QD)/QTTL
      IF (SDEPL.LT.0.) SDEPL=1.E-20
      IF (IPO.EQ.4) WRITE(*,100) X,DPA,(Y(I),I=1,IY)
 100  FORMAT(1X,F10.0,1PE10.3,0P10F5.0)
      RETURN
C     ENTRY VDPL3
 104  WRITE(*,*) '   C OR D     AREA'
      WRITE(*,'(1X,1P2E10.3)') (DDD(I),ARD(I),I=10,1,-1)
      RETURN
      END
