Update legacy code in GMP (#180)

Oden-EAG · Mar 1, 2024 · efef205 · efef205
1 parent d2b3b25
commit efef205
Show file tree

Hide file tree

Showing 54 changed files with 3,626 additions and 3,628 deletions.
diff --git a/trunk/src/GMP/base/GMP_kernels.F90 b/trunk/src/GMP/base/GMP_kernels.F90
@@ -797,15 +797,15 @@ subroutine rectaB_back(No,T,Norient, X,Dxdt)
 !----------------------------------------------------------------------
 !       2nd TERM OF COMPARISON: xw, dxwdse                            !
 !  .....evaluate local edge parameterization of master triangle       !
-        call edge_param(QUAD,ie,se, sw,dswdse)                      !
+        call edge_param(QUAD,ie,se, sw,dswdse)                        !
 !  .....evaluate rectangle parameterization                           !
         call recta(No,sw, xw,dxwds)                                   !
 !  .....compute derivative wrt edge parameter                         !
         dxwdse(1:3) = dxwds(1:3,1)*dswdse(1) + dxwds(1:3,2)*dswdse(2) !
 !  .....evaluate the edge bubble                                      !
         xw(1:3) = xw(1:3)                           &                 !
                 - (xv(1:3,iv1)*(1.d0 - se) + xv(1:3,iv2)*se)          !
-        dxwdse(1:3) = dxwdse(1:3)                            &         !
+        dxwdse(1:3) = dxwdse(1:3)                            &        !
                 - (xv(1:3,iv2) -  xv(1:3,iv1))                        !
 !----------------------------------------------------------------------
 !       COMPARE                                                       !
@@ -814,17 +814,17 @@ subroutine rectaB_back(No,T,Norient, X,Dxdt)
           smax = max(smax,abs(xw(ivar) - xc(ivar)))                   !
           dmax = max(dmax,abs(dxcdse(ivar) - dxwdse(ivar)))           !
         enddo                                                         !
-        if ((smax.gt.GEOM_TOL).or.(dmax.gt.GEOM_TOL)) then      !
+        if ((smax.gt.GEOM_TOL).or.(dmax.gt.GEOM_TOL)) then            !
           write(*,7002) No,ie,smax,dmax                               !
  7002     format(' rectaB: No,ie,smax,dmax = ',i5,i2,2e12.5)          !
-!c          call pause                                                  !
-!cc          if (iprint .eq. 1) call my_tests                            !
-!cc          iprint = 1                                                  !
-!cc          write(*,9000) No                                            !
- 9000     format(' trianB: rerunning rountine for triangle No = ',I6, &!
-      ' with printing flag on.')                                      !
-!cc          call pause                                                  !
-!cc          go to 10                                                    !
+!!          call pause                                                !
+!!!          if (iprint .eq. 1) call my_tests                         !
+!!!          iprint = 1                                               !
+!!!          write(*,9000) No                                         !
+ 9000     format(' trianB: rerunning rountine for triangle No = ', &  !
+                           I6,' with printing flag on.')              !
+!!!          call pause                                               !
+!!!          go to 10                                                 !
         endif                                                         !
 !======================================================================
 !

diff --git a/trunk/src/GMP/base/block_hexa.F90 b/trunk/src/GMP/base/block_hexa.F90
@@ -296,7 +296,7 @@ end subroutine hexa_TraHex
 !
 !-----------------------------------------------------------------------
 !
-      subroutine hexa_CylHex(No,Eta, X,Dxdeta)
+   subroutine hexa_CylHex(No,Eta, X,Dxdeta)
 !
       use control
       use GMP          , only : HEXAS,POINTS,NDIM
@@ -394,4 +394,4 @@ subroutine hexa_CylHex(No,Eta, X,Dxdeta)
       endif
 !
 !
-      end subroutine hexa_CylHex
+   end subroutine hexa_CylHex
diff --git a/trunk/src/GMP/base/block_prism.F90 b/trunk/src/GMP/base/block_prism.F90
@@ -72,7 +72,7 @@ subroutine prism(No,Eta, X,Dxdeta)
       endselect
 !
 !
-      end subroutine prism
+   end subroutine prism
 !
 !
 !