[cig-commits] commit: fix the dimension of interfaces in paraview export.

Mercurial hg at geodynamics.org
Thu Jul 5 15:58:13 PDT 2012


changeset:   135:e5de65e3c5dc
tag:         tip
user:        Sylvain Barbot <sbarbot at caltech.edu>
date:        Thu Jul 05 15:58:09 2012 -0700
files:       src/input.f90
description:
fix the dimension of interfaces in paraview export.


diff -r 4b4b129d1b47 -r e5de65e3c5dc src/input.f90
--- a/src/input.f90	Wed Jul 04 02:39:05 2012 -0700
+++ b/src/input.f90	Thu Jul 05 15:58:09 2012 -0700
@@ -519,7 +519,7 @@ CONTAINS
 
           rffilename=trim(in%wdir)//"/linearlayer-"//digit//".vtp"
           CALL exportvtk_rectangle(0.d0,0.d0,in%linearlayer(k)%z, &
-                                   DBLE(in%sx1)*in%dx1,DBLE(in%sx2)*in%dx2, &
+                                   DBLE(in%sx2)*in%dx2,DBLE(in%sx1)*in%dx1, &
                                    0._8,1.5708d0,rffilename)
 #endif
        END DO
@@ -633,7 +633,7 @@ CONTAINS
           ! export the viscous layer in VTK format
           rffilename=trim(in%wdir)//"/nonlinearlayer-"//digit//".vtp"
           CALL exportvtk_rectangle(0.d0,0.d0,in%nonlinearlayer(k)%z, &
-                                   DBLE(in%sx1)*in%dx1,DBLE(in%sx2)*in%dx2, &
+                                   DBLE(in%sx2)*in%dx2,DBLE(in%sx1)*in%dx1, &
                                    0._8,1.57d0,rffilename)
 #endif
        END DO



More information about the CIG-COMMITS mailing list