--- /tmp/RAM_NRL/ram1.5.f 2000-09-13 00:00:00.000000000 -0700 +++ ram.f 2011-10-19 17:11:34.000000000 -0700 @@ -3,6 +3,8 @@ c ****************************************************************** c ***** Range-dependent Acoustic Model, Version 1.5, 13-Sep-00 ***** c ****************************************************************** +c +c \$Id: \$ c c This code was developed by Michael D. Collins at the Naval c Research Laboratory in Washington, DC. It solves range-dependent @@ -48,7 +50,7 @@ c c mr=bathymetry points, mz=depth grid, mp=pade terms. c - parameter (mr=100,mz=8000,mp=10) + parameter (mr=5000,mz=100000,mp=16) dimension rb(mr),zb(mr),cw(mz),cb(mz),rhob(mz),attn(mz),alpw(mz), > alpb(mz),f1(mz),f2(mz),f3(mz),ksq(mz),ksqw(mz),ksqb(mz),u(mz), > v(mz),tlg(mz),r1(mz,mp),r2(mz,mp),r3(mz,mp),s1(mz,mp), @@ -118,7 +120,8 @@ dir=ri-float(ir) k0=omega/c0 nz=zmax/dz-0.5 - nzplt=zmplt/dz-0.5 +c last depth in tl.grid was NOT being output (jcp) + nzplt=zmplt/dz z=zb(1) iz=1.0+z/dz iz=max(2,iz) @@ -150,7 +153,10 @@ do 5 i=ndz,nzplt,ndz lz=lz+1 5 continue - write(3)lz +c write(3)lz +c header containing a more detailed description of the problem (jcp) + write(3) freq,zs,zr,rmax,dr,ndr,zmax,dz,ndz,zmplt,c0,np,ns,rs,lz + c c The initial profiles and starting field. c @@ -460,14 +466,14 @@ complex ur,u(mz) real f3(mz),tlg(mz) c - ur=(1.0-dir)*f3(ir)*u(ir)+dir*f3(ir+1)*u(ir+1) - tl=-20.0*alog10(cabs(ur)+eps)+10.0*alog10(r+eps) - write(2,*)r,tl -c mdr=mdr+1 if(mdr.eq.ndr)then mdr=0 c + ur=(1.0-dir)*f3(ir)*u(ir)+dir*f3(ir+1)*u(ir+1) + tl=-20.0*alog10(cabs(ur)+eps)+10.0*alog10(r+eps) + write(2,*)r,tl +c j=0 do 1 i=ndz,nzplt,ndz ur=u(i)*f3(i)