19 |
! We assume that, in "coefoz_LMDZ.nc", the pressure levels are in hPa |
! We assume that, in "coefoz_LMDZ.nc", the pressure levels are in hPa |
20 |
! and strictly increasing. |
! and strictly increasing. |
21 |
|
|
22 |
use netcdf95, only: nf95_open, nf95_close, nf95_get_coord |
use netcdf95, only: nf95_open, nf95_close, nf95_inq_varid, nf95_gw_var |
23 |
use netcdf, only: nf90_nowrite |
use netcdf, only: nf90_nowrite |
24 |
|
|
25 |
! Variables local to the procedure: |
! Variables local to the procedure: |
26 |
integer ncid ! for NetCDF |
integer ncid, varid ! for NetCDF |
27 |
integer n_plev ! number of pressure levels in the input data |
integer n_plev ! number of pressure levels in the input data |
28 |
integer k |
integer k |
29 |
|
|
33 |
|
|
34 |
call nf95_open("coefoz_LMDZ.nc", nf90_nowrite, ncid) |
call nf95_open("coefoz_LMDZ.nc", nf90_nowrite, ncid) |
35 |
|
|
36 |
call nf95_get_coord(ncid, "plev", plev) |
call nf95_inq_varid(ncid, "plev", varid) |
37 |
|
call nf95_gw_var(ncid, varid, plev) |
38 |
! Convert from hPa to Pa because "regr_pr_av" and "regr_pr_int" |
! Convert from hPa to Pa because "regr_pr_av" and "regr_pr_int" |
39 |
! require so: |
! require so: |
40 |
plev = plev * 100. |
plev = plev * 100. |