phylmd/Mobidic/press_coefoz.f

module press_coefoz_m

  implicit none

  real, pointer, save:: plev(:)
  ! (pressure level of Mobidic input data, converted to Pa, in
  ! ascending order)

  real, allocatable, save:: press_in_edg(:)
  ! (edges of pressure intervals for Mobidic input data, in Pa, in
  ! ascending order)

contains

  subroutine press_coefoz

    ! This procedure is called once per run.
    ! It reads the pressure levels from "coefoz_LMDZ.nc".
    ! We assume that, in "coefoz_LMDZ.nc", the pressure levels are in hPa
    ! and strictly increasing.

    use netcdf95, only: nf95_open, nf95_close, nf95_get_coord
    use netcdf, only: nf90_nowrite

    ! Variables local to the procedure:
    integer ncid ! for NetCDF
    integer n_plev ! number of pressure levels in the input data
    integer k

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

    print *, "Call sequence information: press_coefoz"

    call nf95_open("coefoz_LMDZ.nc", nf90_nowrite, ncid)

    call nf95_get_coord(ncid, "plev", plev)
    ! Convert from hPa to Pa because "regr_pr_av" and "regr_pr_int"
    ! require so:
    plev = plev * 100.
    n_plev = size(plev)

    call nf95_close(ncid)

    ! Compute edges of pressure intervals:
    allocate(press_in_edg(n_plev + 1))
    press_in_edg(1) = 0.
    ! We choose edges halfway in logarithm:
    forall (k = 2:n_plev) press_in_edg(k) = sqrt(plev(k - 1) * plev(k))
    press_in_edg(n_plev + 1) = huge(0.)
    ! (infinity, but any value guaranteed to be greater than the
    ! surface pressure would do)

  end subroutine press_coefoz

end module press_coefoz_m
1	guez	17	module press_coefoz_m
2
3			implicit none
4
5			real, pointer, save:: plev(:)
6			! (pressure level of Mobidic input data, converted to Pa, in
7			! ascending order)
8
9			real, allocatable, save:: press_in_edg(:)
10			! (edges of pressure intervals for Mobidic input data, in Pa, in
11			! ascending order)
12
13			contains
14
15			subroutine press_coefoz
16
17			! This procedure is called once per run.
18			! It reads the pressure levels from "coefoz_LMDZ.nc".
19			! We assume that, in "coefoz_LMDZ.nc", the pressure levels are in hPa
20			! and strictly increasing.
21
22			use netcdf95, only: nf95_open, nf95_close, nf95_get_coord
23			use netcdf, only: nf90_nowrite
24
25			! Variables local to the procedure:
26			integer ncid ! for NetCDF
27			integer n_plev ! number of pressure levels in the input data
28			integer k
29
30			!---------------------------------------
31
32			print *, "Call sequence information: press_coefoz"
33
34			call nf95_open("coefoz_LMDZ.nc", nf90_nowrite, ncid)
35
36			call nf95_get_coord(ncid, "plev", plev)
37			! Convert from hPa to Pa because "regr_pr_av" and "regr_pr_int"
38			! require so:
39			plev = plev * 100.
40			n_plev = size(plev)
41
42			call nf95_close(ncid)
43
44			! Compute edges of pressure intervals:
45			allocate(press_in_edg(n_plev + 1))
46			press_in_edg(1) = 0.
47			! We choose edges halfway in logarithm:
48			forall (k = 2:n_plev) press_in_edg(k) = sqrt(plev(k - 1) * plev(k))
49			press_in_edg(n_plev + 1) = huge(0.)
50			! (infinity, but any value guaranteed to be greater than the
51			! surface pressure would do)
52
53			end subroutine press_coefoz
54
55			end module press_coefoz_m