Additional renames
authorloizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 23 Feb 2011 15:19:40 +0000 (15:19 +0000)
committerloizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 23 Feb 2011 15:19:40 +0000 (15:19 +0000)
TAmpt/AMPT/linana.f
TAmpt/AMPT/zpc.f

index d0bc107..142017c 100644 (file)
@@ -322,12 +322,12 @@ c     boost the reference frame up by beta (pznew=gam(pz+beta e)):
          write(6,*) '1',dbex,dbey,dbez,beta2,gam
       endif
 c
-      call lorenz(de1,dpx,dpy,dpz,-dbex,-dbey,-dbez)
+      call lorenza(de1,dpx,dpy,dpz,-dbex,-dbey,-dbez)
         ptwo(1,1)=sngl(pxnew)
         ptwo(1,2)=sngl(pynew)
         ptwo(1,3)=sngl(pznew)
         ptwo(1,4)=sngl(enenew)
-      call lorenz(de2,-dpx,-dpy,-dpz,-dbex,-dbey,-dbez)
+      call lorenza(de2,-dpx,-dpy,-dpz,-dbex,-dbey,-dbez)
         ptwo(2,1)=sngl(pxnew)
         ptwo(2,2)=sngl(pynew)
         ptwo(2,3)=sngl(pznew)
@@ -376,7 +376,7 @@ cc      SAVE /NOPREC/
       COMMON/HPARNT/HIPR1(100),IHPR2(50),HINT1(100),IHNT2(50)
 cc      SAVE /HPARNT/
 c     7/20/01: use double precision
-c     otherwise sometimes beta>1 and gamma diverge in lorenz():
+c     otherwise sometimes beta>1 and gamma diverge in lorenza():
       COMMON/SOFT/PXSGS(MAXSTR,3),PYSGS(MAXSTR,3),PZSGS(MAXSTR,3),
      &     PESGS(MAXSTR,3),PMSGS(MAXSTR,3),GXSGS(MAXSTR,3),
      &     GYSGS(MAXSTR,3),GZSGS(MAXSTR,3),FTSGS(MAXSTR,3),
@@ -997,7 +997,7 @@ c
                write(6,*) '2',bex,bey,bez,beta2,gam
             endif
 c
-            call lorenz(ftavg0,gxavg0,gyavg0,gzavg0,-bex,-bey,-bez)
+            call lorenza(ftavg0,gxavg0,gyavg0,gzavg0,-bex,-bey,-bez)
               GXAR(inatt)=sngl(pxnew)
               GYAR(inatt)=sngl(pynew)
               GZAR(inatt)=sngl(pznew)
@@ -1365,12 +1365,12 @@ c     boost the reference frame down by beta to get to the pair rest frame:
                write(6,*) '4b',bex,bey,bez,beta2,gam
             endif
 c
-            call lorenz(ftp(j),gxp(j),gyp(j),gzp(j),bex,bey,bez)
+            call lorenza(ftp(j),gxp(j),gyp(j),gzp(j),bex,bey,bez)
             gxp0(j)=pxnew
             gyp0(j)=pynew
             gzp0(j)=pznew
             ftp0(j)=enenew
-            call lorenz(pep(j),pxp(j),pyp(j),pzp(j),bex,bey,bez)
+            call lorenza(pep(j),pxp(j),pyp(j),pzp(j),bex,bey,bez)
             pxp0(j)=pxnew
             pyp0(j)=pynew
             pzp0(j)=pznew
@@ -1406,12 +1406,12 @@ c     for 3-body kinematics, used for baryons formation:
          endif
 c     boost the reference frame down by beta to get to the 3-parton rest frame:
          do 1002 j=1,3
-            call lorenz(ftp(j),gxp(j),gyp(j),gzp(j),bex,bey,bez)
+            call lorenza(ftp(j),gxp(j),gyp(j),gzp(j),bex,bey,bez)
             gxp0(j)=pxnew
             gyp0(j)=pynew
             gzp0(j)=pznew
             ftp0(j)=enenew
-            call lorenz(pep(j),pxp(j),pyp(j),pzp(j),bex,bey,bez)
+            call lorenza(pep(j),pxp(j),pyp(j),pzp(j),bex,bey,bez)
             pxp0(j)=pxnew
             pyp0(j)=pynew
             pzp0(j)=pznew
index cffdd43..b980d69 100644 (file)
@@ -445,7 +445,7 @@ cc        external ran1
         do 1001 i = mul + 1, mul + incmul
            ityp0(i) = 21
            xmass0(i) = xmp
-           call energy(e, temp)
+           call energya(e, temp)
            call momntm(px, py, pz, e)
 c     7/20/01:
 c           e = sqrt(e ** 2 + xmp ** 2)
@@ -455,7 +455,7 @@ c           e = sqrt(e ** 2 + xmp ** 2)
               bex = 0d0
               bey = 0d0
               bez = -tanh(eta(i))
-              call lorenz(e, px, py, pz, bex, bey, bez)
+              call lorenza(e, px, py, pz, bex, bey, bez)
               px0(i) = pxnew
               py0(i) = pynew
               pz0(i) = pznew
@@ -564,7 +564,7 @@ cc      SAVE /rndm3/
         return
         end
         
-        subroutine energy(e, temp)
+        subroutine energya(e, temp)
 
 c       to generate the magnitude of the momentum e,
 c       knowing the temperature of the local thermal distribution temp
@@ -640,7 +640,7 @@ cc      SAVE /prec1/
 cc      SAVE /lora/
         SAVE   
 
-        external lorenz
+        external lorenza
 
         bex = 0d0 
         bey = 0d0
@@ -652,7 +652,7 @@ c       save data for many runs of the same initial condition
            py1 = gy0(i)
            pz1 = gz0(i)
            e1 = ft0(i)
-           call lorenz(e1, px1, py1, pz1, bex, bey, bez)
+           call lorenza(e1, px1, py1, pz1, bex, bey, bez)
            gx0(i) = pxnew
            gy0(i) = pynew
            gz0(i) = pznew
@@ -661,7 +661,7 @@ c       save data for many runs of the same initial condition
            py1 = py0(i)
            pz1 = pz0(i)
            e1 = e0(i)
-           call lorenz(e1, px1, py1, pz1, bex, bey, bez)
+           call lorenza(e1, px1, py1, pz1, bex, bey, bez)
            px0(i) = pxnew
            py0(i) = pynew
            pz0(i) = pznew
@@ -2021,7 +2021,7 @@ ctransend
         bey = (py1 + py2) / (e1 + e2)
         bez = (pz1 + pz2) / (e1 + e2)
 
-        call lorenz(e1, px1, py1, pz1, bex, bey, bez)
+        call lorenza(e1, px1, py1, pz1, bex, bey, bez)
 cc      SAVE pxnew, ..., values for later use.
         px1 = pxnew
         py1 = pynew
@@ -2036,13 +2036,13 @@ cc      SAVE pxnew, ..., values for later use.
 c       we boost to the cm frame, get rotation axis, and rotate 1 particle 
 c       momentum
 
-        call lorenz(t1, x1, y1, z1, bex, bey, bez)
+        call lorenza(t1, x1, y1, z1, bex, bey, bez)
 
         x1 = pxnew
         y1 = pynew
         z1 = pznew
 
-        call lorenz(t2, x2, y2, z2, bex, bey, bez)
+        call lorenza(t2, x2, y2, z2, bex, bey, bez)
 
         x2 = pxnew
         y2 = pynew
@@ -2067,12 +2067,12 @@ c        e2 = e1
         e2 = dsqrt(px2**2+py2**2+pz2**2+xmass(jscat)**2)
 
 c       boost the 2 particle 4 momentum back to lab frame
-        call lorenz(e1, px1, py1, pz1, -bex, -bey, -bez)
+        call lorenza(e1, px1, py1, pz1, -bex, -bey, -bez)
         px(iscat) = pxnew
         py(iscat) = pynew
         pz(iscat) = pznew
         e(iscat) = enenew
-        call lorenz(e2, px2, py2, pz2, -bex, -bey, -bez)        
+        call lorenza(e2, px2, py2, pz2, -bex, -bey, -bez)        
         px(jscat) = pxnew
         py(jscat) = pynew
         pz(jscat) = pznew
@@ -6429,7 +6429,7 @@ cms  &                        em(i, 4, ian) / vol / ntotal
 
 ******************************************************************************
 
-      subroutine lorenz(energy, px, py, pz, bex, bey, bez)
+      subroutine lorenza(energy, px, py, pz, bex, bey, bez)
 
 c     add in a cut for beta2 to prevent gam to be nan (infinity)