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

[willic3 at geodynamics.org]

Author: willic3
Date: 2006-06-20 19:10:19 -0700 (Tue, 20 Jun 2006)
New Revision: 3830

Modified:
   short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/write_ucd_header.F
Log:
I left off the units in the field descriptions, which caused ParaView
to barf.  This has been fixed.  It should already be correct for
binary UCD output.


Modified: short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/write_ucd_header.F
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/write_ucd_header.F	2006-06-21 01:42:06 UTC (rev 3829)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/write_ucd_header.F	2006-06-21 02:10:19 UTC (rev 3830)
@@ -57,11 +57,18 @@
       character stress*6,strain*4,strate*9
       data stress,strain,strate/"Pascal","None","1/seconds"/
 c
+c...  external functions
+c
+      integer nchar
+      external nchar
+c
 c...  local variables
 c
       integer i,ibyte,indl,indu
+      integer il
       integer iout(3*nstatesmax)
       character nlabels*1024,nunits*1024
+      character slabel*100
 c
 c...  included variable definitions
 c
@@ -75,7 +82,15 @@
         write(kucd,"(100i5)") nstatestot,(iout(i),i=1,nstatestot)
 c
         do i=1,nstatestot
-          write(kucd,"(a11)") labels(istatoutc(i))
+          il=nchar(labels(istatoutc(i)))
+          if(istatoutc(i).le.nstatesmax) then
+            slabel=labels(istatoutc(i))(1:il)//","//stress
+          else if(istatoutc(i).le.2*nstatesmax) then
+            slabel=labels(istatoutc(i))(1:il)//","//strain
+          else if(istatoutc(i).le.3*nstatesmax) then
+            slabel=labels(istatoutc(i))(1:il)//","//strate
+          end if
+          write(kucd,"(a100)") slabel
         end do
       else if(iucd.eq.itwo) then
         indl=ione