--- trunk/phylmd/Orography/start_init_orog_m.f	2014/03/05 14:57:53	82
+++ trunk/Sources/phylmd/Orography/start_init_orog.f	2018/01/05 14:45:45	247
@@ -17,7 +17,7 @@
        zthe_2d, zpic_2d, zval_2d)
 
     use conf_dat2d_m, only: conf_dat2d
-    use comgeom, only: rlatu, rlonv
+    use dynetat0_m, only: rlatu, rlonv
     use grid_noro_m, only: grid_noro
     use indicesol, only: epsfra
     use netcdf, only: nf90_nowrite
@@ -53,9 +53,9 @@
     INTEGER iml_rel
     INTEGER jml_rel
     INTEGER ncid, varid
-    REAL, pointer:: relief(:, :)
+    REAL, ALLOCATABLE:: relief(:, :) ! in m
     REAL, ALLOCATABLE:: lon_rad(:), lat_rad(:)
-    REAL, pointer:: lon_ini(:), lat_ini(:)
+    REAL, ALLOCATABLE:: lon_ini(:), lat_ini(:)
 
     !-----------------------------------
 
@@ -92,13 +92,10 @@
 
     CALL conf_dat2d(lon_ini, lat_ini, lon_rad, lat_rad, relief , &
          interbar=.FALSE.)
-    deallocate(lon_ini, lat_ini) ! pointers
 
     print *, 'Compute all the parameters needed for the gravity wave drag code'
-
     CALL grid_noro(lon_rad, lat_rad, relief, rlonv, rlatu, phis, zmea_2d, &
          zstd_2d, zsig_2d, zgam_2d, zthe_2d, zpic_2d, zval_2d, mask)
-    deallocate(relief) ! pointer
     phis(:, :) = phis(:, :) * 9.81
 
     mask(2:, 1) = mask(1, 1) ! north pole