← Back to branch summary

~ubuntu-branches/ubuntu/wily/ifrit/wily

« back to all changes in this revision

Viewing changes to docs/writeIFrIT.f

  • Committer: Bazaar Package Importer
  • Author(s): Mark Hymers
  • Date: 2009-09-13 14:53:24 UTC
  • mfrom: (1.1.11 upstream) (2.1.4 sid)
  • Revision ID: james.westby@ubuntu.com-20090913145324-084ivqvf29arkc1b
Tags: 3.3.2-1
* New upstream release.
* Bump Standards-Version to 3.8.3.  No changes needed.

Show diffs side-by-side

added added

removed removed

Lines of Context:
docs/writeIFrIT.f
1
 
C
2
 
C  Write to text uniform scalars data file
3
 
C
4
 
      subroutine WriteIFrITUniformScalarsTxtFile(n1,n2,n3,var1,var2,var3,
5
 
     .     filename)
6
 
      integer n1, n2, n3  ! Size of the computational mesh in 3 directions
7
 
      real*4 var1(n1,n2,n3)
8
 
      real*4 var2(n1,n2,n3) ! Three scalar variables
9
 
      real*4 var3(n1,n2,n3)  
10
 
      character*(*) filename  ! Name of the file
11
 
      open(unit=1, file=filename)
12
 
      write(1,*) n1, n2, n3
13
 
      do k=1,n3
14
 
         do j=1,n2
15
 
            do i=1,n1
16
 
               write(1,*) var1(i,j,k), var2(i,j,k), var3(i,j,k)
17
 
            enddo
18
 
         enddo
19
 
      enddo
20
 
      close(1)
21
 
      return
22
 
      end
23
 
C
24
 
C  Write to binary uniform scalars data file
25
 
C
26
 
      subroutine WriteIFrITUniformScalarsBinFile(n1,n2,n3,var1,var2,var3,
27
 
     .     filename)
28
 
      integer n1, n2, n3  ! Size of the computational mesh in 3 directions
29
 
      real*4 var1(n1,n2,n3)
30
 
      real*4 var2(n1,n2,n3) ! Three scalar variables
31
 
      real*4 var3(n1,n2,n3)  
32
 
      character*(*) filename  ! Name of the file
33
 
      open(unit=1, file=filename, form='unformatted')
34
 
      write(1) n1, n2, n3
35
 
      write(1) (((var1(i,j,k),i=1,n1),j=1,n2),k=1,n3)
36
 
      write(1) (((var2(i,j,k),i=1,n1),j=1,n2),k=1,n3)
37
 
      write(1) (((var3(i,j,k),i=1,n1),j=1,n2),k=1,n3)
38
 
      close(1)
39
 
      return
40
 
      end
41
 
C
42
 
C  Write to text uniform vectors data file
43
 
C
44
 
      subroutine WriteIFrITUniformVectorsTxtFile(n1,n2,n3,vect,filename)
45
 
      integer n1, n2, n3  ! Size of the computational mesh in 3 directions
46
 
      real*4 vect(3,n1,n2,n3)  ! Vector field
47
 
      character*(*) filename  ! Name of the file
48
 
      open(unit=1, file=filename)
49
 
      write(1,*) n1, n2, n3
50
 
      do k=1,n3
51
 
         do j=1,n2
52
 
            do i=1,n1
53
 
               write(1,*) vect(1,i,j,k), vect(2,i,j,k), vect(3,i,j,k)
54
 
            enddo
55
 
         enddo
56
 
      enddo
57
 
      close(1)
58
 
      return
59
 
      end
60
 
C
61
 
C  Write to binary uniform vectors data file
62
 
C
63
 
      subroutine WriteIFrITUniformVectorsBinFile(n1,n2,n3,vect,filename)
64
 
      integer n1, n2, n3  ! Size of the computational mesh in 3 directions
65
 
      real*4 vect(3,n1,n2,n3)  ! Vector field
66
 
      character*(*) filename  ! Name of the file
67
 
      open(unit=1, file=filename, form='unformatted')
68
 
      write(1) n1, n2, n3
69
 
      write(1) (((vect(1,i,j,k),i=1,n1),j=1,n2),k=1,n3)
70
 
      write(1) (((vect(2,i,j,k),i=1,n1),j=1,n2),k=1,n3)
71
 
      write(1) (((vect(3,i,j,k),i=1,n1),j=1,n2),k=1,n3)
72
 
      close(1)
73
 
      return
74
 
      end
75
 
C
76
 
C  Write to text uniform tensors data file
77
 
C
78
 
      subroutine WriteIFrITUniformTensorsTxtFile(n1,n2,n3,tens,filename)
79
 
      integer n1, n2, n3  ! Size of the computational mesh in 3 directions
80
 
      real*4 tens(6,n1,n2,n3)  ! Tensor field
81
 
      character*(*) filename  ! Name of the file
82
 
      open(unit=1, file=filename)
83
 
      write(1,*) n1, n2, n3
84
 
      do k=1,n3
85
 
         do j=1,n2
86
 
            do i=1,n1
87
 
               write(1,*) tens(1,i,j,k), tens(2,i,j,k), tens(3,i,j,k),
88
 
     .              tens(4,i,j,k), tens(5,i,j,k), tens(6,i,j,k)
89
 
            enddo
90
 
         enddo
91
 
      enddo
92
 
      close(1)
93
 
      return
94
 
      end
95
 
C
96
 
C  Write to binary uniform tensors data file
97
 
C
98
 
      subroutine WriteIFrITUniformTensorsBinFile(n1,n2,n3,tens,filename)
99
 
      integer n1, n2, n3  ! Size of the computational mesh in 3 directions
100
 
      real*4 tens(6,n1,n2,n3)  ! Tensor field
101
 
      character*(*) filename  ! Name of the file
102
 
      open(unit=1, file=filename, form='unformatted')
103
 
      write(1) n1, n2, n3
104
 
      write(1) (((tens(1,i,j,k),i=1,n1),j=1,n2),k=1,n3)
105
 
      write(1) (((tens(2,i,j,k),i=1,n1),j=1,n2),k=1,n3)
106
 
      write(1) (((tens(3,i,j,k),i=1,n1),j=1,n2),k=1,n3)
107
 
      write(1) (((tens(4,i,j,k),i=1,n1),j=1,n2),k=1,n3)
108
 
      write(1) (((tens(5,i,j,k),i=1,n1),j=1,n2),k=1,n3)
109
 
      write(1) (((tens(6,i,j,k),i=1,n1),j=1,n2),k=1,n3)
110
 
      close(1)
111
 
      return
112
 
      end
113
 
C
114
 
C  Write to text basic particles data file
115
 
C
116
 
      subroutine WriteIFrITBasicParticlesTxtFile(n,xl,yl,zl,xh,yh,zh,
117
 
     .     x,y,z,attr1,attr2,attr3,filename)
118
 
      integer n  ! Number of particles
119
 
      real*4 xl, yl, zl, xh, yh, zh  ! Bounding box 
120
 
      real*4 x(n), y(n), z(n)  ! Particle positions (can be real*8)
121
 
      real*4 attr1(n), attr2(n), attr3(n)  ! Particle attributes
122
 
      character*(*) filename  ! Name of the file
123
 
      open(unit=1, file=filename)
124
 
      write(1,*) n
125
 
      write(1,*) xl, yl, zl, xh, yh, zh
126
 
      do i=1,n
127
 
         write(1,*) x(i), y(i), z(i), attr1(i), attr2(i), attr3(i)
128
 
      enddo
129
 
      close(1)
130
 
      return
131
 
      end
132
 
C
133
 
C  Write to binary basic particles data file
134
 
C
135
 
      subroutine WriteIFrITBasicParticlesBinFile(n,xl,yl,zl,xh,yh,zh,
136
 
     .     x,y,z,attr1,attr2,attr3,filename)
137
 
      integer n  ! Number of particles
138
 
      real*4 xl, yl, zl, xh, yh, zh  ! Bounding box
139
 
      real*4 x(n), y(n), z(n)  ! Particle positions (can be real*8)
140
 
      real*4 attr1(n), attr2(n), attr3(n)  ! Particle attributes
141
 
      character*(*) filename  ! Name of the file
142
 
      open(unit=1, file=filename, form='unformatted')
143
 
      write(1) n
144
 
      write(1) xl, yl, zl, xh, yh, zh
145
 
      write(1) (x(i),i=1,n)
146
 
      write(1) (y(i),i=1,n)
147
 
      write(1) (z(i),i=1,n)
148
 
      write(1) (attr1(i),i=1,n)
149
 
      write(1) (attr2(i),i=1,n)
150
 
      write(1) (attr3(i),i=1,n)
151
 
      close(1)
152
 
      return
153
 
      end
 
1
C
 
2
C  Write to text uniform scalars data file
 
3
C
 
4
      subroutine WriteIFrITUniformScalarsTxtFile(n1,n2,n3,var1,var2,var3,
 
5
     .     filename)
 
6
      integer n1, n2, n3  ! Size of the computational mesh in 3 directions
 
7
      real*4 var1(n1,n2,n3)
 
8
      real*4 var2(n1,n2,n3) ! Three scalar variables
 
9
      real*4 var3(n1,n2,n3)  
 
10
      character*(*) filename  ! Name of the file
 
11
      open(unit=1, file=filename)
 
12
      write(1,*) n1, n2, n3
 
13
      do k=1,n3
 
14
         do j=1,n2
 
15
            do i=1,n1
 
16
               write(1,*) var1(i,j,k), var2(i,j,k), var3(i,j,k)
 
17
            enddo
 
18
         enddo
 
19
      enddo
 
20
      close(1)
 
21
      return
 
22
      end
 
23
C
 
24
C  Write to binary uniform scalars data file
 
25
C
 
26
      subroutine WriteIFrITUniformScalarsBinFile(n1,n2,n3,var1,var2,var3,
 
27
     .     filename)
 
28
      integer n1, n2, n3  ! Size of the computational mesh in 3 directions
 
29
      real*4 var1(n1,n2,n3)
 
30
      real*4 var2(n1,n2,n3) ! Three scalar variables
 
31
      real*4 var3(n1,n2,n3)  
 
32
      character*(*) filename  ! Name of the file
 
33
      open(unit=1, file=filename, form='unformatted')
 
34
      write(1) n1, n2, n3
 
35
      write(1) (((var1(i,j,k),i=1,n1),j=1,n2),k=1,n3)
 
36
      write(1) (((var2(i,j,k),i=1,n1),j=1,n2),k=1,n3)
 
37
      write(1) (((var3(i,j,k),i=1,n1),j=1,n2),k=1,n3)
 
38
      close(1)
 
39
      return
 
40
      end
 
41
C
 
42
C  Write to text uniform vectors data file
 
43
C
 
44
      subroutine WriteIFrITUniformVectorsTxtFile(n1,n2,n3,vect,filename)
 
45
      integer n1, n2, n3  ! Size of the computational mesh in 3 directions
 
46
      real*4 vect(3,n1,n2,n3)  ! Vector field
 
47
      character*(*) filename  ! Name of the file
 
48
      open(unit=1, file=filename)
 
49
      write(1,*) n1, n2, n3
 
50
      do k=1,n3
 
51
         do j=1,n2
 
52
            do i=1,n1
 
53
               write(1,*) vect(1,i,j,k), vect(2,i,j,k), vect(3,i,j,k)
 
54
            enddo
 
55
         enddo
 
56
      enddo
 
57
      close(1)
 
58
      return
 
59
      end
 
60
C
 
61
C  Write to binary uniform vectors data file
 
62
C
 
63
      subroutine WriteIFrITUniformVectorsBinFile(n1,n2,n3,vect,filename)
 
64
      integer n1, n2, n3  ! Size of the computational mesh in 3 directions
 
65
      real*4 vect(3,n1,n2,n3)  ! Vector field
 
66
      character*(*) filename  ! Name of the file
 
67
      open(unit=1, file=filename, form='unformatted')
 
68
      write(1) n1, n2, n3
 
69
      write(1) (((vect(1,i,j,k),i=1,n1),j=1,n2),k=1,n3)
 
70
      write(1) (((vect(2,i,j,k),i=1,n1),j=1,n2),k=1,n3)
 
71
      write(1) (((vect(3,i,j,k),i=1,n1),j=1,n2),k=1,n3)
 
72
      close(1)
 
73
      return
 
74
      end
 
75
C
 
76
C  Write to text uniform tensors data file
 
77
C
 
78
      subroutine WriteIFrITUniformTensorsTxtFile(n1,n2,n3,tens,filename)
 
79
      integer n1, n2, n3  ! Size of the computational mesh in 3 directions
 
80
      real*4 tens(6,n1,n2,n3)  ! Tensor field
 
81
      character*(*) filename  ! Name of the file
 
82
      open(unit=1, file=filename)
 
83
      write(1,*) n1, n2, n3
 
84
      do k=1,n3
 
85
         do j=1,n2
 
86
            do i=1,n1
 
87
               write(1,*) tens(1,i,j,k), tens(2,i,j,k), tens(3,i,j,k),
 
88
     .              tens(4,i,j,k), tens(5,i,j,k), tens(6,i,j,k)
 
89
            enddo
 
90
         enddo
 
91
      enddo
 
92
      close(1)
 
93
      return
 
94
      end
 
95
C
 
96
C  Write to binary uniform tensors data file
 
97
C
 
98
      subroutine WriteIFrITUniformTensorsBinFile(n1,n2,n3,tens,filename)
 
99
      integer n1, n2, n3  ! Size of the computational mesh in 3 directions
 
100
      real*4 tens(6,n1,n2,n3)  ! Tensor field
 
101
      character*(*) filename  ! Name of the file
 
102
      open(unit=1, file=filename, form='unformatted')
 
103
      write(1) n1, n2, n3
 
104
      write(1) (((tens(1,i,j,k),i=1,n1),j=1,n2),k=1,n3)
 
105
      write(1) (((tens(2,i,j,k),i=1,n1),j=1,n2),k=1,n3)
 
106
      write(1) (((tens(3,i,j,k),i=1,n1),j=1,n2),k=1,n3)
 
107
      write(1) (((tens(4,i,j,k),i=1,n1),j=1,n2),k=1,n3)
 
108
      write(1) (((tens(5,i,j,k),i=1,n1),j=1,n2),k=1,n3)
 
109
      write(1) (((tens(6,i,j,k),i=1,n1),j=1,n2),k=1,n3)
 
110
      close(1)
 
111
      return
 
112
      end
 
113
C
 
114
C  Write to text basic particles data file
 
115
C
 
116
      subroutine WriteIFrITBasicParticlesTxtFile(n,xl,yl,zl,xh,yh,zh,
 
117
     .     x,y,z,attr1,attr2,attr3,filename)
 
118
      integer n  ! Number of particles
 
119
      real*4 xl, yl, zl, xh, yh, zh  ! Bounding box 
 
120
      real*4 x(n), y(n), z(n)  ! Particle positions (can be real*8)
 
121
      real*4 attr1(n), attr2(n), attr3(n)  ! Particle attributes
 
122
      character*(*) filename  ! Name of the file
 
123
      open(unit=1, file=filename)
 
124
      write(1,*) n
 
125
      write(1,*) xl, yl, zl, xh, yh, zh
 
126
      do i=1,n
 
127
         write(1,*) x(i), y(i), z(i), attr1(i), attr2(i), attr3(i)
 
128
      enddo
 
129
      close(1)
 
130
      return
 
131
      end
 
132
C
 
133
C  Write to binary basic particles data file
 
134
C
 
135
      subroutine WriteIFrITBasicParticlesBinFile(n,xl,yl,zl,xh,yh,zh,
 
136
     .     x,y,z,attr1,attr2,attr3,filename)
 
137
      integer n  ! Number of particles
 
138
      real*4 xl, yl, zl, xh, yh, zh  ! Bounding box
 
139
      real*4 x(n), y(n), z(n)  ! Particle positions (can be real*8)
 
140
      real*4 attr1(n), attr2(n), attr3(n)  ! Particle attributes
 
141
      character*(*) filename  ! Name of the file
 
142
      open(unit=1, file=filename, form='unformatted')
 
143
      write(1) n
 
144
      write(1) xl, yl, zl, xh, yh, zh
 
145
      write(1) (x(i),i=1,n)
 
146
      write(1) (y(i),i=1,n)
 
147
      write(1) (z(i),i=1,n)
 
148
      write(1) (attr1(i),i=1,n)
 
149
      write(1) (attr2(i),i=1,n)
 
150
      write(1) (attr3(i),i=1,n)
 
151
      close(1)
 
152
      return
 
153
      end