← Back to branch summary

~asc/fluidity/iceshelfcavity

« back to all changes in this revision

Viewing changes to legacy_reservoir_prototype/tests/CompModel_1d_2Phases_3Components_NonEquilibrium/InputComp2P3C_permability.f90

  • Committer: Adam Candy
  • Date: 2011-10-12 20:39:04 UTC
  • Revision ID: adam.candy@imperial.ac.uk-20111012203904-9gwqzs4z798zcw0l
Automatic sync to launchpad

Show diffs side-by-side

added added

removed removed

Lines of Context:
legacy_reservoir_prototype/tests/CompModel_1d_2Phases_3Components_NonEquilibrium/InputComp2P3C_permability.f90
1
 
program main
2
 
  implicit none
3
 
  integer :: ndim1, ndim2, ndim3
4
 
  real, dimension( :, :, : ), allocatable :: field
5
 
  integer :: idim1, idim2, idim3
6
 
 
7
 
  open( 5, file = 'filedim', status = 'unknown' )
8
 
  read( 5, * ) ndim1, ndim2, ndim3
9
 
  close( 5 )
10
 
 
11
 
  allocate( field( ndim1, ndim2, ndim3 ))
12
 
 
13
 
  field = 0
14
 
  field( 1 : ndim1, 1, 1 ) = 1.
15
 
 
16
 
  do idim1 = 1, ndim1
17
 
     do idim2 = 1, ndim2
18
 
        do idim3 = 1, ndim3
19
 
           print*, field( idim1, idim2, idim3 )
20
 
        end do
21
 
     end do
22
 
  end do
23
 
 
24
 
 
25
 
end program main