trunk/libf/test_disvert.f90

program test_disvert

  use dimens_m, only: llm
  use comvert, only: pa, disvert, ap, bp, preff
  use comconst, only: initialize
  use new_unit_m, only: new_unit

  implicit none

  integer unit, l
  REAL sigma(llm+1) ! "sigma(l)" : à l'interface des couches "l" et "l-1"
  real p(llm + 1) ! pressure at level, in Pa

  !---------------------

  print *, "llm = ", llm
  call initialize
  pa = 5e4
  call disvert
  sigma = ap / pa + bp
  p = ap + bp * preff

  unit = new_unit()
  open(unit, file="test_disvert.csv", status="replace", action="write")
  write(unit, fmt=*) '"ap (Pa)" "bp" "sigma" "p (Pa)"' ! title line
  do l = 1, llm
     write(unit, fmt=*) ap(l), bp(l), sigma(l), p(l)
  end do
  close(unit)
  print *, 'The file "test_disvert.csv" has been created.'

end program test_disvert
1	guez	3	program test_disvert
2
3			use dimens_m, only: llm
4			use comvert, only: pa, disvert, ap, bp, preff
5			use comconst, only: initialize
6	guez	7	use new_unit_m, only: new_unit
7	guez	3
8			implicit none
9
10			integer unit, l
11			REAL sigma(llm+1) ! "sigma(l)" : à l'interface des couches "l" et "l-1"
12			real p(llm + 1) ! pressure at level, in Pa
13
14			!---------------------
15
16			print *, "llm = ", llm
17			call initialize
18			pa = 5e4
19			call disvert
20			sigma = ap / pa + bp
21			p = ap + bp * preff
22
23			unit = new_unit()
24			open(unit, file="test_disvert.csv", status="replace", action="write")
25			write(unit, fmt=*) '"ap (Pa)" "bp" "sigma" "p (Pa)"' ! title line
26			do l = 1, llm
27			write(unit, fmt=*) ap(l), bp(l), sigma(l), p(l)
28			end do
29			close(unit)
30			print *, 'The file "test_disvert.csv" has been created.'
31
32			end program test_disvert