trunk/libf/test_disvert.f90

program test_disvert

  ! This is a program in Fortran 95.
  ! This program writes the vertical coordinates at level interfaces.

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

  implicit none

  integer unit, l

  REAL s(llm+1) ! hybrid sigma-pressure coordinate
  ! "s(l)" : à l'interface des couches "l" et "l-1"

  real p(llm + 1) ! pressure at level, in Pa
  real z(llm) ! pressure-altitude (km)

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

  print *, "llm = ", llm
  pa = 5e4
  call disvert
  s = ap / pa + bp
  p = ap + bp * preff
  z = - 7. * log(ap(:llm) / preff + bp(:llm))

  call new_unit(unit)
  open(unit, file="test_disvert.csv", status="replace", action="write")
  ! Title line:
  write(unit, fmt=*) '"ap (hPa)" "bp" "s" "p (hPa)" "z (km)"'
  do l = 1, llm
     write(unit, fmt=*) ap(l) / 100., bp(l), s(l), p(l) / 100., z(l)
  end do
  close(unit)
  print *, 'The file "test_disvert.csv" has been created.'

  open(unit, file="delta_z.csv", status="replace", action="write")
  ! Title line:
  write(unit, fmt=*) '"z middle (km)" "delta z (km)"'
  do l = 1, llm - 1
     write(unit, fmt=*) (z(l) + z(l+1)) / 2, z(l+1) - z(l)
  end do
  close(unit)
  print *, 'The file "delta_z.csv" has been created.'

end program test_disvert
1	program test_disvert
2
3	! This is a program in Fortran 95.
4	! This program writes the vertical coordinates at level interfaces.
5
6	use dimens_m, only: llm
7	use comvert, only: pa, disvert, ap, bp, preff
8	use new_unit_m, only: new_unit
9
10	implicit none
11
12	integer unit, l
13
14	REAL s(llm+1) ! hybrid sigma-pressure coordinate
15	! "s(l)" : à l'interface des couches "l" et "l-1"
16
17	real p(llm + 1) ! pressure at level, in Pa
18	real z(llm) ! pressure-altitude (km)
19
20	!---------------------
21
22	print *, "llm = ", llm
23	pa = 5e4
24	call disvert
25	s = ap / pa + bp
26	p = ap + bp * preff
27	z = - 7. * log(ap(:llm) / preff + bp(:llm))
28
29	call new_unit(unit)
30	open(unit, file="test_disvert.csv", status="replace", action="write")
31	! Title line:
32	write(unit, fmt=*) '"ap (hPa)" "bp" "s" "p (hPa)" "z (km)"'
33	do l = 1, llm
34	write(unit, fmt=*) ap(l) / 100., bp(l), s(l), p(l) / 100., z(l)
35	end do
36	close(unit)
37	print *, 'The file "test_disvert.csv" has been created.'
38
39	open(unit, file="delta_z.csv", status="replace", action="write")
40	! Title line:
41	write(unit, fmt=*) '"z middle (km)" "delta z (km)"'
42	do l = 1, llm - 1
43	write(unit, fmt=*) (z(l) + z(l+1)) / 2, z(l+1) - z(l)
44	end do
45	close(unit)
46	print *, 'The file "delta_z.csv" has been created.'
47
48	end program test_disvert