--- trunk/libf/phylmd/Mobidic/regr_pr_o3.f90 2008/10/15 16:19:57 20 +++ trunk/libf/phylmd/Mobidic/regr_pr_o3.f90 2011/07/01 15:00:48 47 @@ -25,11 +25,11 @@ use conf_gcm_m, only: dayref use dimens_m, only: iim, jjm, llm use netcdf95, only: nf95_open, nf95_close, nf95_inq_varid, handle_err, & - nf95_get_coord + nf95_gw_var use netcdf, only: nf90_nowrite, nf90_get_var - use numer_rec, only: assert + use nr_util, only: assert use grid_change, only: dyn_phy - use regr1_step_av_m, only: regr1_step_av + use l_util, only: regr1_step_av use pressure_var, only: p3d real, intent(out):: o3_mob_regr(:, :, :) ! (iim + 1, jjm + 1, llm) @@ -61,7 +61,8 @@ call nf95_open("coefoz_LMDZ.nc", nf90_nowrite, ncid) - call nf95_get_coord(ncid, "plev", plev) + call nf95_inq_varid(ncid, "plev", varid) + call nf95_gw_var(ncid, varid, plev) ! Convert from hPa to Pa because "regr_pr_av" requires so: plev = plev * 100. n_plev = size(plev)