--- trunk/libf/phylmd/Mobidic/press_coefoz.f90	2008/08/05 13:31:32	17
+++ trunk/Sources/phylmd/Mobidic/press_coefoz.f	2017/10/16 12:35:41	225
@@ -2,7 +2,7 @@
 
   implicit none
 
-  real, pointer, save:: plev(:)
+  real, allocatable, save:: plev(:)
   ! (pressure level of Mobidic input data, converted to Pa, in
   ! ascending order)
 
@@ -19,11 +19,11 @@
     ! 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 netcdf95, only: nf95_open, nf95_close, nf95_inq_varid, nf95_gw_var
     use netcdf, only: nf90_nowrite
 
     ! Variables local to the procedure:
-    integer ncid ! for NetCDF
+    integer ncid, varid ! for NetCDF
     integer n_plev ! number of pressure levels in the input data
     integer k
 
@@ -33,7 +33,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" and "regr_pr_int"
     ! require so:
     plev = plev * 100.