CAD in Civil Engineering - Baustatik-Info-Server

Page 46 Computer Languages for Engineering - SS 13131132 write(io,*,iostat=ioerr) n1,n2 ! write the dimensions133 if (ioerr /= 0) then134 write(*,*) ’*** Error: writing not possible’135 iwritemat = -1 ! exit, if io error136 return137 endif138139 ! over the rows140 do i=1,n1141 write(io,*) (m(i,j),j=1,n2)142 enddo143144 iwritemat = 0 ! 0 return: everything ok145 end function iwritemat146147 ! Read the dimension of a matrix from a file148 integer function ireadmatdim(io,n1,n2)149150 integer::io ! io channel151 integer::n1,n2 ! dimensions of the matrix152 integer::ioerr ! error flag153154 read(io,*,iostat=ioerr) n1,n2155 if (ioerr /= 0) then ! if io-error, perhaps a wrong format156 write(*,*) ’*** Error: wrong file format’157 ireadmatdim = -1158 return159 endif160161 ireadmatdim = 0162 end function ireadmatdim163164 ! read matrix data from a file165 !166 integer function ireadmat(io,m,n1,n2)167 implicit none168169 integer::io ! io channel number170 integer::n1 ! number of rows171 integer::n2 ! number of columns172 real(8), dimension(n1,n2):: m ! matrix data173 integer::ioerr, i, j174175 ! over the rows176 do i=1,n1177 read(io,*,iostat=ioerr) (m(i,j),j=1,n2)178 if (ioerr /= 0) then ! important to check the read-io179 write(*,*)’*** Error: format’180 ireadmat = -1181 return182 endif183 enddo184185 ireadmat = 0186 end function ireadmatE. Baeck