Changeset 7447
 Timestamp:
 06/01/18 14:25:09 (4 months ago)
 Location:
 trunk/src/seismic_processing/rayloc_ew
 Files:

 3 edited
Legend:
 Unmodified
 Added
 Removed

trunk/src/seismic_processing/rayloc_ew/raylocator.f
r1669 r7447 7 7 implicit none 8 8 c 9 integer NTRBL 10 c 9 11 character*8 phlst(2) 10 12 character*1000 argv … … 12 14 logical*4 prnt(3) 13 15 c 14 integer*4 ntrbl,ln,nphase,nexit16 integer*4 ln,nphase,nexit 15 17 integer*4 iargc,nargs 16 18 c integer stdin, stdout … … 79 81 if(argv.ne.'i') go to 999 80 82 call getarg(2,argv) 81 ln= ntrbl(argv)83 ln=NTRBL(argv) 82 84 tokenName=argv 83 85 … … 96 98 c 97 99 call freeunit(mlu) 98 call tabin(mlu,dirname(1: ntrbl(dirname))//'/'//modnam)100 call tabin(mlu,dirname(1:NTRBL(dirname))//'/'//modnam) 99 101 call brnset(1,phlst,prnt) 100 102 c 101 103 call freeunit(merr) 102 open(merr,file=dirname(1: ntrbl(dirname))//'/'//'RayLocError'//tokenName(1:ln)//'.txt',104 open(merr,file=dirname(1:NTRBL(dirname))//'/'//'RayLocError'//tokenName(1:ln)//'.txt', 103 105 1 access='sequential',form='formatted',status='new') 104 106 c … … 106 108 c 107 109 call freeunit(min) 108 open(min,file=dirname(1: ntrbl(dirname))//'/'//'RayLocInput'//tokenName(1:ln)//'.txt',110 open(min,file=dirname(1:NTRBL(dirname))//'/'//'RayLocInput'//tokenName(1:ln)//'.txt', 109 111 1 access='sequential',form='formatted',status='old') 110 112 rewind(min) … … 121 123 c 122 124 call freeunit(mout) 123 open(mout,file=dirname(1: ntrbl(dirname))//'/'//'RayLocOutput'//tokenName(1:ln)//'.txt',125 open(mout,file=dirname(1:NTRBL(dirname))//'/'//'RayLocOutput'//tokenName(1:ln)//'.txt', 124 126 1 access='sequential',form='formatted',status='new') 125 127 call output(nphase) 
trunk/src/seismic_processing/rayloc_ew/robust_util.f
r4497 r7447 1515 1515 c 1516 1516 implicit none 1517 integer NTRBL 1517 1518 real depmin,cn 1518 1519 parameter(depmin=1.,cn=5.) … … 1522 1523 character*20 prmmod 1523 1524 logical first 1524 integer*4 n,i,j, ntrbl,mprm,ln1525 integer*4 n,i,j,mprm,ln 1525 1526 c integer*4 min,mout,merr 1526 1527 real*4 zs,tt1(maxtt),dtdd1(maxtt),dtdh1(maxtt),dddp1(maxtt), … … 1646 1647 else if(phase.eq.'pPKP'.or.phase.eq.'sPKP'.or.phase.eq.'SKP'.or. 1647 1648 1 phase.eq.'PKKP'.or.phase.eq.'P''P'''.or.phase.eq.'SKS') then 1648 ln= ntrbl(phase)1649 ln=NTRBL(phase) 1649 1650 i=0 1650 1651 res=1e30 … … 1765 1766 c Return the elevation correction. 1766 1767 c 1767 implicit none 1768 implicit none 1769 integer NTRBL 1768 1770 real vp,vs,kmpd 1769 1771 parameter(vp=5.80,vs=3.46,kmpd=111.19493) 1770 1772 character*(*) phase 1771 integer*4 ntrbl,j1773 integer*4 j 1772 1774 real*4 elev,dtdd 1773 1775 c 1774 1776 c Figure out the arriving wave type. 1775 do j= ntrbl(phase),1,11777 do j=NTRBL(phase),1,1 1776 1778 if(phase(j:j).eq.'P') then 1777 1779 c This phase arrives as a P. Do the correction using Vp. … … 1797 1799 c 1798 1800 implicit none 1801 integer NTRBL 1799 1802 character*(*) modnam 1800 1803 character*20 flnm 1801 integer*4 nprm,i,j,k ,ntrbl1804 integer*4 nprm,i,j,k 1802 1805 include 'resprm.inc' 1803 1806 data flnm/'res_ac.prm'/ … … 1818 1821 > status='old') 1819 1822 else 1820 open(nprm,file=dirname(1: ntrbl(dirname))//'/'//flnm,1823 open(nprm,file=dirname(1:NTRBL(dirname))//'/'//flnm, 1821 1824 > access='sequential',form='formatted',status='old') 1822 1825 endif … … 1841 1844 1 nbr=i1 1842 1845 if(nbr.lt.mxbr) phnm(nbr+1)=' ' 1843 c print *,'Getprm: modnam nbr = ',modnam(1: ntrbl(modnam)),nbr1846 c print *,'Getprm: modnam nbr = ',modnam(1:NTRBL(modnam)),nbr 1844 1847 close(nprm) 1845 1848 return 
trunk/src/seismic_processing/rayloc_ew/robust_util_linux.f
r3484 r7447 1515 1515 c 1516 1516 implicit none 1517 real depmin,cn 1517 integer NTRBL 1518 real depmin,cn 1518 1519 parameter(depmin=1.,cn=5.) 1519 1520 c Miscellaneous constants. … … 1522 1523 character*20 prmmod 1523 1524 logical first 1524 integer*4 n,i,j, ntrbl,mprm,ln1525 integer*4 n,i,j,mprm,ln 1525 1526 c integer*4 min,mout,merr 1526 1527 real*4 zs,tt1(maxtt),dtdd1(maxtt),dtdh1(maxtt),dddp1(maxtt), … … 1644 1645 else if(phase.eq.'pPKP'.or.phase.eq.'sPKP'.or.phase.eq.'SKP'.or. 1645 1646 1 phase.eq.'PKKP'.or.phase.eq.'P''P'''.or.phase.eq.'SKS') then 1646 ln= ntrbl(phase)1647 ln=NTRBL(phase) 1647 1648 i=0 1648 1649 res=1e30 … … 1763 1764 c Return the elevation correction. 1764 1765 c 1765 implicit none 1766 implicit none 1767 integer NTRBL 1766 1768 real vp,vs,kmpd 1767 1769 parameter(vp=5.80,vs=3.46,kmpd=111.19493) 1768 1770 character*(*) phase 1769 integer*4 ntrbl,j1771 integer*4 j 1770 1772 real*4 elev,dtdd 1771 1773 c 1772 1774 c Figure out the arriving wave type. 1773 do j= ntrbl(phase),1,11775 do j=NTRBL(phase),1,1 1774 1776 if(phase(j:j).eq.'P') then 1775 1777 c This phase arrives as a P. Do the correction using Vp. … … 1795 1797 c 1796 1798 implicit none 1799 integer NTRBL 1797 1800 character*(*) modnam 1798 1801 character*20 flnm 1799 integer*4 nprm,i,j,k ,ntrbl1802 integer*4 nprm,i,j,k 1800 1803 include 'resprm.inc' 1801 1804 data flnm/'res_ac.prm'/ … … 1816 1819 > status='old') 1817 1820 else 1818 open(nprm,file=dirname(1: ntrbl(dirname))//'/'//flnm,1821 open(nprm,file=dirname(1:NTRBL(dirname))//'/'//flnm, 1819 1822 > access='sequential',form='formatted',status='old') 1820 1823 endif … … 1839 1842 1 nbr=i1 1840 1843 if(nbr.lt.mxbr) phnm(nbr+1)=' ' 1841 c print *,'Getprm: modnam nbr = ',modnam(1: ntrbl(modnam)),nbr1844 c print *,'Getprm: modnam nbr = ',modnam(1:NTRBL(modnam)),nbr 1842 1845 close(nprm) 1843 1846 return
Note: See TracChangeset
for help on using the changeset viewer.