Merge pull request #125 from OpenSEMBA/dev

Merge dev with feature/sgbc changes into main

Author: Alberto-o

src_main_pub/healer.F90

@@ -222,7 +222,7 @@ SUBROUTINE CreateVolumeMM (layoutnumber, Mtag, tags, numertag, MMiEx, MMiEy, MMi
             DO i = punto%XI, puntoPlus1%XE
                medio = MMiHx (i, j, k)
 !               IF (medio /= 0) THEN   !ojo esto estaba antes de 031016 y daba maxima prioridad al medio 0 PEC. Ahora puedo tener medios con mas prioridad!!! !?!? cambio agresivo 031016!!!
-                  IF (med(indicemedio)%Priority >= med(medio)%Priority) THEN
+                  IF (med(indicemedio)%Priority > med(medio)%Priority) THEN
                      MMiHx (i, j, k) = indicemedio
                      Mtag(i,j,k)=64*numertag 
                      tags%face%x(i,j,k) = 64*numertag
@@ -238,7 +238,7 @@ SUBROUTINE CreateVolumeMM (layoutnumber, Mtag, tags, numertag, MMiEx, MMiEy, MMi
             DO i = punto%XI, punto%XE
                medio = MMiHy (i, j, k)
 !               IF (medio /= 0) THEN   !ojo esto estaba antes de 031016 y daba maxima prioridad al medio 0 PEC. Ahora puedo tener medios con mas prioridad!!! !?!? cambio agresivo 031016!!!
-                  IF (med(indicemedio)%Priority >= med(medio)%Priority) THEN
+                  IF (med(indicemedio)%Priority > med(medio)%Priority) THEN
                      MMiHy (i, j, k) = indicemedio
                      Mtag(i,j,k)=64*numertag 
                      tags%face%y(i,j,k) = 64*numertag
@@ -254,7 +254,7 @@ SUBROUTINE CreateVolumeMM (layoutnumber, Mtag, tags, numertag, MMiEx, MMiEy, MMi
             DO i = punto%XI, punto%XE
                medio = MMiHz (i, j, k)
 !               IF (medio /= 0) THEN   !ojo esto estaba antes de 031016 y daba maxima prioridad al medio 0 PEC. Ahora puedo tener medios con mas prioridad!!! !?!? cambio agresivo 031016!!!
-                  IF (med(indicemedio)%Priority >= med(medio)%Priority) THEN
+                  IF (med(indicemedio)%Priority > med(medio)%Priority) THEN
                      MMiHz (i, j, k) = indicemedio
                      Mtag(i,j,k)=64*numertag 
                      tags%face%z(i,j,k) = 64*numertag
@@ -361,12 +361,12 @@ SUBROUTINE CreateSurfaceMM (layoutnumber, Mtag, tags, numertag, MMiEx, MMiEy, MM
                      CALL AddToShared (iEz, i, j, k, indicemedio, medio, Eshared)
                   END IF
                END DO
-            END DO
+            END DO  
             DO j = punto%YI, punto%YE
                DO k = punto%ZI, punto%ZE
                   medio = MMiHx (i, j, k)
 !                  IF (medio /= 0) THEN   !ojo esto estaba antes de 031016 y daba maxima prioridad al medio 0 PEC. Ahora puedo tener medios con mas prioridad!!! !?!? cambio agresivo 031016!!!
-                     IF (med(indicemedio)%Priority >= med(medio)%Priority) then
+                     IF (med(indicemedio)%Priority > med(medio)%Priority) then
                          MMiHx (i, j, k) = indicemedio; 
                          Mtag(i,j,k)=64*numertag ! if (.true..or.(Mtag(i,j,k)==0).or.(int(Mtag(i,j,k)/64) == numertag)) Mtag(i,j,k) = IBSET(64*numertag,3);
                          tags%face%x(i,j,k) = 64*numertag
@@ -408,7 +408,7 @@ SUBROUTINE CreateSurfaceMM (layoutnumber, Mtag, tags, numertag, MMiEx, MMiEy, MM
                DO k = punto%ZI, punto%ZE
                   medio = MMiHy (i, j, k)
 !                  IF (medio /= 0) THEN   !ojo esto estaba antes de 031016 y daba maxima prioridad al medio 0 PEC. Ahora puedo tener medios con mas prioridad!!! !?!? cambio agresivo 031016!!!
-                     IF (med(indicemedio)%Priority >= med(medio)%Priority) then
+                     IF (med(indicemedio)%Priority > med(medio)%Priority) then
                          MMiHy (i, j, k) = indicemedio; 
                          Mtag(i,j,k)=64*numertag ! if (.true..or.(Mtag(i,j,k)==0).or.(int(Mtag(i,j,k)/64) == numertag)) Mtag(i,j,k) = IBSET(64*numertag,4);;
                          tags%face%y(i,j,k) = 64*numertag
@@ -449,7 +449,7 @@ SUBROUTINE CreateSurfaceMM (layoutnumber, Mtag, tags, numertag, MMiEx, MMiEy, MM
                DO j = punto%YI, punto%YE
                   medio = MMiHz (i, j, k)
 !                  IF (medio /= 0) THEN   !ojo esto estaba antes de 031016 y daba maxima prioridad al medio 0 PEC. Ahora puedo tener medios con mas prioridad!!! !?!? cambio agresivo 031016!!!
-                     IF (med(indicemedio)%Priority >= med(medio)%Priority) then
+                     IF (med(indicemedio)%Priority > med(medio)%Priority) then
                          MMiHz (i, j, k) = indicemedio; 
                          Mtag(i,j,k)=64*numertag ! if (.true..or.(Mtag(i,j,k)==0).or.(int(Mtag(i,j,k)/64) == numertag)) Mtag(i,j,k) = IBSET(64*numertag,5);
                          tags%face%z(i,j,k) = 64*numertag