Changeset 882 for codes/icosagcm/trunk/src/sphere

Timestamp:

06/11/19 14:59:17 (5 years ago)

Author:

ymipsl

Message:

Metric is now write in start.nc/restart.nc
Metric can be read at restart if read_metric=y.

YM

File:

• codes/icosagcm/trunk/src/sphere/geometry.f90 (modified) (1 diff)

codes/icosagcm/trunk/src/sphere/geometry.f90

@@ -295 +295 @@
     CALL update_circumcenters
 
-    DO ind=1,ndomain
-      IF (.NOT. assigned_domain(ind)  .OR. .NOT. is_omp_level_master ) CYCLE
-      d=>domain(ind)
-      CALL swap_dimensions(ind)
-      CALL swap_geometry(ind)
-      DO j=jj_begin,jj_end
-        DO i=ii_begin,ii_end
-          n=(j-1)*iim+i
-          DO k=0,5
-            x1(:) = xyz_v(n+z_pos(k+1),:)
-            x2(:) = d%vertex(:,k,i,j) 
-            IF (norm(x1-x2)>1e-10) THEN
-              PRINT*,"vertex diff ",ind,i,j,k
-              PRINT*,x1
-              PRINT*,x2
-            ENDIF
-          ENDDO
-        ENDDO
-      ENDDO
-    ENDDO
+!ym you get discrepency if you try to read metric from file, vertex are not saved in file
+!    DO ind=1,ndomain
+!      IF (.NOT. assigned_domain(ind)  .OR. .NOT. is_omp_level_master ) CYCLE
+!      d=>domain(ind)
+!      CALL swap_dimensions(ind)
+!      CALL swap_geometry(ind)
+!      DO j=jj_begin,jj_end
+!        DO i=ii_begin,ii_end
+!          n=(j-1)*iim+i
+!          DO k=0,5
+!            x1(:) = xyz_v(n+z_pos(k+1),:)
+!            x2(:) = d%vertex(:,k,i,j) 
+!            IF (norm(x1-x2)>1e-10) THEN
+!              PRINT*,"vertex diff ",ind,i,j,k
+!              PRINT*,x1
+!              PRINT*,x2
+!            ENDIF
+!          ENDDO
+!        ENDDO
+!      ENDDO
+!   ENDDO