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

Timestamp:

05/22/12 18:06:53 (12 years ago)

Author:

ymipsl

Message:

some update

YM

File:

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

codes/icosagcm/trunk/src/metric.f90

@@ -194 +194 @@
    ENDDO
 
-  CALL dist_cart(vertex_glo(1,1,1)%xyz,vertex_glo(2,1,1)%xyz,d1) 
-  CALL dist_cart(vertex_glo(2,1,1)%xyz,vertex_glo(3,1,1)%xyz,d2) 
-  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 dist_cart(vertex_glo(1,1,1)%xyz,vertex_glo(2,1,1)%xyz,d1) 
+!  CALL dist_cart(vertex_glo(2,1,1)%xyz,vertex_glo(3,1,1)%xyz,d2) 
+!  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)
@@ -206 +206 @@
 !  PRINT *,"dist",vertex_glo(2,1,1)%xyz
 !  PRINT *,"dist",p1/sqrt(sum(p1**2))
-  CALL Centroide(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)
-  CALL Centroide(vertex_glo(1,2,1)%xyz,vertex_glo(1,1,1)%xyz,vertex_glo(2,1,1)%xyz,p1)
-  CALL dist_cart(vertex_glo(1,1,1)%xyz,p1,d1) 
-  CALL dist_cart(vertex_glo(2,1,1)%xyz,p1,d2) 
-  CALL dist_cart(vertex_glo(1,2,1)%xyz,p1,d3) 
+!  CALL Centroide(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)
+!  CALL Centroide(vertex_glo(1,2,1)%xyz,vertex_glo(1,1,1)%xyz,vertex_glo(2,1,1)%xyz,p1)
+!  CALL dist_cart(vertex_glo(1,1,1)%xyz,p1,d1) 
+!  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