[cig-commits] r5493 - short/3D/lithomop/trunk/lithomop3d/liblithomop3d

[willic3 at geodynamics.org]

Author: willic3
Date: 2006-12-07 07:20:49 -0800 (Thu, 07 Dec 2006)
New Revision: 5493

Modified:
   short/3D/lithomop/trunk/lithomop3d/liblithomop3d/write_ucd_mesh.F
Log:
Changed default integer output format to prevent overflows.



Modified: short/3D/lithomop/trunk/lithomop3d/liblithomop3d/write_ucd_mesh.F
===================================================================
--- short/3D/lithomop/trunk/lithomop3d/liblithomop3d/write_ucd_mesh.F	2006-12-07 11:56:23 UTC (rev 5492)
+++ short/3D/lithomop/trunk/lithomop3d/liblithomop3d/write_ucd_mesh.F	2006-12-07 15:20:49 UTC (rev 5493)
@@ -144,12 +144,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
@@ -163,7 +163,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
@@ -173,7 +173,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
@@ -190,7 +190,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)