view release on metacpan or  search on metacpan

lib/App/SeismicUnixGui/fortran/src/wrimod2.for  view on Meta::CPAN

	SUBROUTINE WRIMOD2(NL,VT,VB,DZ,VST,VSB,RHOT,RHOB)
	DIMENSION VT(*),VB(*),DZ(*),VST(*),VSB(*),RHOT(*),RHOB(*)
C
C WRITES VELOCITY MODEL ON THE SCREEN 
C
	ICH = 6
	WRITE(ICH,50)
	WRITE(ICH,100)
	WRITE(ICH,50)
C
	A1=0.
	A2=0.
	DO 20 I=1,NL
	A1=A1+DZ(I)
C COMPUTES 2-WAY TRAVELTIME ASSUMING EITHER CONSTANT VELOCITY OR
C CONSTANT VELOCITY GRADIENT LAYERS.
C
	IF(VT(I).EQ.VB(I)) THEN
		A2=A2+2.*DZ(I)/VT(I)
	ELSE
		A2=A2+2.*DZ(I)*ALOG(VB(I)/VT(I))/(VB(I)-VT(I))
	ENDIF
C
  	WRITE(ICH,200) I,VT(I),VB(I),DZ(I),A1,A2,VST(I),VSB(I),
     +  RHOT(I),RHOB(I)
20	CONTINUE
	WRITE(ICH,50)
C
50	FORMAT('')
100	FORMAT(3X,'#',1X,'VTOP(km/s)','VBOTT',X,'DZ',1X,'(km)',1X,'ZBOTT',
     +  3X,'TWTT(s)',3X,'VSTOP',2X,'VSBOTT',4X,'RHOT',4X,'RHOB')
200	FORMAT(I4,11F8.5)
C
	RETURN
        END