~asc/fluidity/iceshelfcavity

« back to all changes in this revision

Viewing changes to legacy_reservoir_prototype/tests/CompModel_1d_2Phases_3Components_LVI_NonEquilibrium_DG/InputComp2P3C_EOSCoeff.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:
1
 
program main
2
 
  implicit none
3
 
  integer :: ndim1, ndim2
4
 
  real, dimension( :, : ), allocatable :: field
5
 
  integer :: idim1, idim2
6
 
 
7
 
  open( 5, file = 'filedim', status = 'unknown' )
8
 
  read( 5, * ) ndim1, ndim2
9
 
  close( 5 )
10
 
 
11
 
  allocate( field( ndim1, ndim2 ))
12
 
 
13
 
  field = 0
14
 
  field( 1 : ndim1, 1 ) = 1.
15
 
 
16
 
  do idim1 = 1, ndim1
17
 
     do idim2 = 1, ndim2
18
 
        print*, field( idim1, idim2 )
19
 
     end do
20
 
  end do
21
 
 
22
 
end program main