Changeset 7447


Ignore:
Timestamp:
06/01/18 14:25:09 (4 months ago)
Author:
baker
Message:

Function NTRBL is default Integer type, not Integer*4

Location:
trunk/src/seismic_processing/rayloc_ew
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/seismic_processing/rayloc_ew/raylocator.f

    r1669 r7447  
    77        implicit none 
    88c 
     9        integer NTRBL 
     10c 
    911        character*8 phlst(2) 
    1012        character*1000 argv 
     
    1214        logical*4 prnt(3) 
    1315c 
    14         integer*4 ntrbl,ln,nphase,nexit 
     16        integer*4 ln,nphase,nexit 
    1517        integer*4 iargc,nargs 
    1618c       integer stdin, stdout 
     
    7981        if(argv.ne.'-i') go to 999 
    8082        call getarg(2,argv) 
    81         ln=ntrbl(argv) 
     83        ln=NTRBL(argv) 
    8284        tokenName=argv 
    8385 
     
    9698c 
    9799        call freeunit(mlu) 
    98         call tabin(mlu,dirname(1:ntrbl(dirname))//'/'//modnam) 
     100        call tabin(mlu,dirname(1:NTRBL(dirname))//'/'//modnam) 
    99101        call brnset(1,phlst,prnt) 
    100102c 
    101103        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', 
    103105        1 access='sequential',form='formatted',status='new') 
    104106c 
     
    106108c 
    107109        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', 
    109111        1       access='sequential',form='formatted',status='old') 
    110112        rewind(min) 
     
    121123c 
    122124        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', 
    124126        1       access='sequential',form='formatted',status='new') 
    125127        call output(nphase) 
  • trunk/src/seismic_processing/rayloc_ew/robust_util.f

    r4497 r7447  
    15151515c 
    15161516      implicit none 
     1517      integer NTRBL 
    15171518        real depmin,cn 
    15181519        parameter(depmin=1.,cn=5.) 
     
    15221523        character*20 prmmod 
    15231524        logical first 
    1524         integer*4 n,i,j,ntrbl,mprm,ln 
     1525        integer*4 n,i,j,mprm,ln 
    15251526c       integer*4 min,mout,merr 
    15261527      real*4 zs,tt1(maxtt),dtdd1(maxtt),dtdh1(maxtt),dddp1(maxtt), 
     
    16461647      else if(phase.eq.'pPKP'.or.phase.eq.'sPKP'.or.phase.eq.'SKP'.or. 
    16471648     1 phase.eq.'PKKP'.or.phase.eq.'P''P'''.or.phase.eq.'SKS') then 
    1648         ln=ntrbl(phase) 
     1649        ln=NTRBL(phase) 
    16491650        i=0 
    16501651        res=1e30 
     
    17651766c   Return the elevation correction. 
    17661767c 
    1767         implicit none 
     1768      implicit none 
     1769      integer NTRBL 
    17681770      real vp,vs,kmpd 
    17691771      parameter(vp=5.80,vs=3.46,kmpd=111.19493) 
    17701772        character*(*) phase 
    1771         integer*4 ntrbl,j 
     1773        integer*4 j 
    17721774        real*4 elev,dtdd 
    17731775c 
    17741776c   Figure out the arriving wave type. 
    1775         do j=ntrbl(phase),1,-1 
     1777        do j=NTRBL(phase),1,-1 
    17761778                if(phase(j:j).eq.'P') then 
    17771779c   This phase arrives as a P.  Do the correction using Vp. 
     
    17971799c 
    17981800      implicit none 
     1801      integer NTRBL 
    17991802      character*(*) modnam 
    18001803      character*20 flnm 
    1801       integer*4 nprm,i,j,k,ntrbl 
     1804      integer*4 nprm,i,j,k 
    18021805      include 'resprm.inc' 
    18031806      data flnm/'res_ac.prm'/ 
     
    18181821     >       status='old') 
    18191822      else 
    1820         open(nprm,file=dirname(1:ntrbl(dirname))//'/'//flnm, 
     1823        open(nprm,file=dirname(1:NTRBL(dirname))//'/'//flnm, 
    18211824     >       access='sequential',form='formatted',status='old') 
    18221825      endif 
     
    18411844 1    nbr=i-1 
    18421845      if(nbr.lt.mxbr) phnm(nbr+1)=' ' 
    1843 c     print *,'Getprm:  modnam nbr = ',modnam(1:ntrbl(modnam)),nbr 
     1846c     print *,'Getprm:  modnam nbr = ',modnam(1:NTRBL(modnam)),nbr 
    18441847      close(nprm) 
    18451848      return 
  • trunk/src/seismic_processing/rayloc_ew/robust_util_linux.f

    r3484 r7447  
    15151515c 
    15161516      implicit none 
    1517         real depmin,cn 
     1517      integer NTRBL 
     1518      real depmin,cn 
    15181519        parameter(depmin=1.,cn=5.) 
    15191520c   Miscellaneous constants. 
     
    15221523        character*20 prmmod 
    15231524        logical first 
    1524         integer*4 n,i,j,ntrbl,mprm,ln 
     1525        integer*4 n,i,j,mprm,ln 
    15251526c       integer*4 min,mout,merr 
    15261527      real*4 zs,tt1(maxtt),dtdd1(maxtt),dtdh1(maxtt),dddp1(maxtt), 
     
    16441645      else if(phase.eq.'pPKP'.or.phase.eq.'sPKP'.or.phase.eq.'SKP'.or. 
    16451646     1 phase.eq.'PKKP'.or.phase.eq.'P''P'''.or.phase.eq.'SKS') then 
    1646         ln=ntrbl(phase) 
     1647        ln=NTRBL(phase) 
    16471648        i=0 
    16481649        res=1e30 
     
    17631764c   Return the elevation correction. 
    17641765c 
    1765         implicit none 
     1766      implicit none 
     1767      integer NTRBL 
    17661768      real vp,vs,kmpd 
    17671769      parameter(vp=5.80,vs=3.46,kmpd=111.19493) 
    17681770        character*(*) phase 
    1769         integer*4 ntrbl,j 
     1771        integer*4 j 
    17701772        real*4 elev,dtdd 
    17711773c 
    17721774c   Figure out the arriving wave type. 
    1773         do j=ntrbl(phase),1,-1 
     1775        do j=NTRBL(phase),1,-1 
    17741776                if(phase(j:j).eq.'P') then 
    17751777c   This phase arrives as a P.  Do the correction using Vp. 
     
    17951797c 
    17961798      implicit none 
     1799      integer NTRBL 
    17971800      character*(*) modnam 
    17981801      character*20 flnm 
    1799       integer*4 nprm,i,j,k,ntrbl 
     1802      integer*4 nprm,i,j,k 
    18001803      include 'resprm.inc' 
    18011804      data flnm/'res_ac.prm'/ 
     
    18161819     >       status='old') 
    18171820      else 
    1818         open(nprm,file=dirname(1:ntrbl(dirname))//'/'//flnm, 
     1821        open(nprm,file=dirname(1:NTRBL(dirname))//'/'//flnm, 
    18191822     >       access='sequential',form='formatted',status='old') 
    18201823      endif 
     
    18391842 1    nbr=i-1 
    18401843      if(nbr.lt.mxbr) phnm(nbr+1)=' ' 
    1841 c     print *,'Getprm:  modnam nbr = ',modnam(1:ntrbl(modnam)),nbr 
     1844c     print *,'Getprm:  modnam nbr = ',modnam(1:NTRBL(modnam)),nbr 
    18421845      close(nprm) 
    18431846      return 
Note: See TracChangeset for help on using the changeset viewer.