~asc/fluidity/iceshelfcavity

« back to all changes in this revision

Viewing changes to legacy_reservoir_prototype/tests/BL_test_10elements_1d/input_satura_BL1dg_func.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
 
 
2
 
program input_satura_funct
3
 
  implicit none
4
 
  ! Allocating  SATURA( CV_PHA_NONODS )
5
 
  integer, parameter :: cv_nloc = 3, totele = 50, nphase = 2
6
 
  integer, parameter :: cv_nonods = cv_nloc * totele, cv_pha_nonods = cv_nonods * nphase
7
 
  real, dimension( cv_pha_nonods ) :: satura
8
 
  ! Local
9
 
  integer :: i
10
 
 
11
 
 ! print*, cv_pha_nonods
12
 
  do i = 1, cv_pha_nonods
13
 
     if( i <= cv_nonods ) then
14
 
        satura( i ) = 0.
15
 
     else
16
 
        satura( i ) = 1.
17
 
     endif
18
 
     print*, satura( i )
19
 
  end do
20
 
 
21
 
  return
22
 
end program input_satura_funct
23