Index: LMDZ6/trunk/libf/dyn3d/iniacademic.F90
===================================================================
--- LMDZ6/trunk/libf/dyn3d/iniacademic.F90	(revision 4374)
+++ LMDZ6/trunk/libf/dyn3d/iniacademic.F90	(revision 4376)
@@ -18,5 +18,5 @@
   USE comconst_mod, ONLY: cpp, kappa, g, daysec, dtvr, pi, im, jm
   USE logic_mod, ONLY: iflag_phys, read_start
-  USE comvert_mod, ONLY: ap, bp, preff, presnivs, pressure_exner
+  USE comvert_mod, ONLY: ap, bp, preff, pa, presnivs, pressure_exner
   USE temps_mod, ONLY: annee_ref, day_ini, day_ref
   USE ener_mod, ONLY: etot0,ptot0,ztot0,stot0,ang0
@@ -114,4 +114,14 @@
   endif
 
+
+  if ( ( .not. read_start ) .and. ( iflag_phys > 2 ) ) then
+        ! 101080 : specific value for CMIP5 aqua/terra planets
+        ! "Specify the initial dry mass to be equivalent to
+        !  a global mean surface pressure (101325 minus 245) Pa."
+        preff=101080.
+        call getin('psurf',preff)
+        pa=preff/2.
+  endif
+
   CALL iniconst
   CALL inigeom
@@ -120,14 +130,5 @@
   ! Initialize pressure and mass field if read_start=.false.
   IF (.NOT. read_start) THEN
-     ! surface pressure
-     if (iflag_phys>2) then
-        ! specific value for CMIP5 aqua/terra planets
-        ! "Specify the initial dry mass to be equivalent to
-        !  a global mean surface pressure (101325 minus 245) Pa."
-        ps(:)=101080.  
-     else
-        ! use reference surface pressure
-        ps(:)=preff
-     endif
+     ps(:)=preff
 
      ! ground geopotential
Index: LMDZ6/trunk/libf/dyn3dmem/iniacademic_loc.F90
===================================================================
--- LMDZ6/trunk/libf/dyn3dmem/iniacademic_loc.F90	(revision 4374)
+++ LMDZ6/trunk/libf/dyn3dmem/iniacademic_loc.F90	(revision 4376)
@@ -19,5 +19,5 @@
   USE comconst_mod, ONLY: cpp, kappa, g, daysec, dtvr, pi, im, jm
   USE logic_mod, ONLY: iflag_phys, read_start
-  USE comvert_mod, ONLY: ap, bp, preff, presnivs, pressure_exner
+  USE comvert_mod, ONLY: ap, bp, preff, pa, presnivs, pressure_exner
   USE temps_mod, ONLY: annee_ref, day_ini, day_ref
   USE ener_mod, ONLY: etot0,ptot0,ztot0,stot0,ang0
@@ -118,4 +118,13 @@
   endif
 
+  if ( ( .not. read_start ) .and. (iflag_phys>2) ) then
+     ! 101080 : specific value for CMIP5 aqua/terra planets
+     ! "Specify the initial dry mass to be equivalent to
+     !  a global mean surface pressure (101325 minus 245) Pa."
+     preff=101080.
+     call getin('psurf',preff)
+     pa=preff/2.
+  endif
+
   CALL iniconst
   CALL inigeom
@@ -133,13 +142,5 @@
 
      ! surface pressure
-     if (iflag_phys>2) then
-        ! specific value for CMIP5 aqua/terra planets
-        ! "Specify the initial dry mass to be equivalent to
-        !  a global mean surface pressure (101325 minus 245) Pa."
-        ps_glo(:)=101080.  
-     else
-        ! use reference surface pressure
-        ps_glo(:)=preff
-     endif
+     ps_glo(:)=preff
 
      ! ground geopotential