phylmd/Mobidic/press_coefoz.f90

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	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