[cig-commits] r5485 - short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d

willic3 at geodynamics.org willic3 at geodynamics.org
Wed Dec 6 13:25:42 PST 2006


Author: willic3
Date: 2006-12-06 13:25:42 -0800 (Wed, 06 Dec 2006)
New Revision: 5485

Modified:
   short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/write_ucd_mesh.F
Log:
Changed integer output format to allow larger values for ASCII UCD files.



Modified: short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/write_ucd_mesh.F
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/write_ucd_mesh.F	2006-12-06 20:59:49 UTC (rev 5484)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/write_ucd_mesh.F	2006-12-06 21:25:42 UTC (rev 5485)
@@ -145,12 +145,12 @@
 c
 c...  header info
 c
-        write(kucd,"(5i7)") numnp,numelv,nnattr,neattr,nmattr
+        write(kucd,"(5i9)") numnp,numelv,nnattr,neattr,nmattr
 c
 c...  write nodal coordinates
 c
         do i=1,numnp
-          write(kucd,"(i7,3(2x,1pe15.8))") i,(x(j,i),j=1,nsd)
+          write(kucd,"(i9,3(2x,1pe15.8))") i,(x(j,i),j=1,nsd)
         end do
 c
 c...  write element connectivities
@@ -164,7 +164,7 @@
             do j=1,nen
               itmp(j)=ien(j,ielg)
             end do
-            write(kucd,"(2i7,2x,a4,20i7)") ielg,matid,
+            write(kucd,"(2i9,2x,a4,20i9)") ielg,matid,
      &       eltype(indtype),(itmp(inducd(j,indtype)),j=1,nen)
           end do
         end do
@@ -174,7 +174,7 @@
 c
         iopt=itwo
         call open_ucd(kucd,iprestress,nstep,ucdroot,iopt,iucd)
-        write(kucd,"(5i7)") ngpts,numelg,nngattr,negattr,nmattr
+        write(kucd,"(5i9)") ngpts,numelg,nngattr,negattr,nmattr
 c
 c...  write Gauss coordinates
 c
@@ -191,7 +191,7 @@
               xg(2)=xg(2)+xl(2,j)*sh(4,j,l)
               xg(3)=xg(3)+xl(3,j)*sh(4,j,l)
             end do
-            write(kucd,"(i7,3(2x,1pe15.8))") igpt,(xg(j),j=1,nsd)
+            write(kucd,"(i9,3(2x,1pe15.8))") igpt,(xg(j),j=1,nsd)
           end do
         end do
         close(kucd)



More information about the cig-commits mailing list