Changeset 131 for codes/icosagcm/trunk/src/metric.f90

Timestamp:

02/08/13 12:09:35 (11 years ago)

Author:

ymipsl

Message:

Some operations must be only done by the mpi master task.

YM

File:

• codes/icosagcm/trunk/src/metric.f90 (modified) (4 diffs)

codes/icosagcm/trunk/src/metric.f90

@@ -124 +124 @@
   CALL dist_cart(vertex_glo(3,1,1)%xyz,vertex_glo(4,1,1)%xyz,d3) 
   CALL div_arc(vertex_glo(1,1,1)%xyz,vertex_glo(3,1,1)%xyz,0.5,p1)
-!  CALL div_arc(vertex_glo(2,1,1)%xyz,vertex_glo(1,3,1)%xyz,1./3,p1)
-!  CALL div_arc(vertex_glo(1,2,1)%xyz,vertex_glo(3,1,1)%xyz,1./3,p2)
-!  CALL div_arc(vertex_glo(2,2,1)%xyz,vertex_glo(1,1,1)%xyz,1./3,p3)
-!  PRINT *, "dist",d1
-!  PRINT *, "dist",d2
-!  PRINT *, "dist",d3
-!  PRINT *,"dist",vertex_glo(2,1,1)%xyz
-!  PRINT *,"dist",p1/sqrt(sum(p1**2))
+
   CALL circumcenter(vertex_glo(1,1,1)%xyz,vertex_glo(2,1,1)%xyz,vertex_glo(1,2,1)%xyz,p1)
 !  CALL Centroide(vertex_glo(1,2,1)%xyz,vertex_glo(2,1,1)%xyz,vertex_glo(1,1,1)%xyz,p1)
@@ -138 +131 @@
   CALL dist_cart(vertex_glo(2,1,1)%xyz,p1,d2) 
   CALL dist_cart(vertex_glo(1,2,1)%xyz,p1,d3) 
-!  PRINT *, "dist",d1
-!  PRINT *, "dist",d2
-!  PRINT *, "dist",d3
 
   END SUBROUTINE compute_face
@@ -337 +327 @@
          
         ind=vertex_glo(i,j,nf)%ind
-        IF (ind==0) THEN
-          PRINT *,"ind=0",i,j,nf
-        ENDIF
         delta=MOD(vertex_glo(i,j,nf)%delta+neighbour+6,6)
         ind2=cell_glo(ind)%neighbour(delta)
@@ -758 +745 @@
     cell_glo(ind)%neighbour(3)=cell_glo(ind)%neighbour(2)
       
-
-    ind=vertex_glo(1,jjm_glo,1)%ind
-    DO i=0,5
-      ind2=cell_glo(ind)%neighbour(i)
-      IF (ind2>0) PRINT *,"neighbour",i,cell_glo(ind2)%assign_face,cell_glo(ind2)%assign_i,cell_glo(ind2)%assign_j
-    ENDDO
-    
     
  !! assignation des delta