Index: trunk/LMDZ.PLUTO/libf/phypluto/testconservmass.F90
===================================================================
--- trunk/LMDZ.PLUTO/libf/phypluto/testconservmass.F90	(revision 3573)
+++ trunk/LMDZ.PLUTO/libf/phypluto/testconservmass.F90	(revision 3577)
@@ -4,5 +4,5 @@
       use comcstfi_mod, only: g
       use geometry_mod, only: cell_area
-      USE comgeomfi_h, only: totarea
+      USE comgeomfi_h, only: totarea,totarea_planet
       implicit none
 
@@ -49,5 +49,5 @@
       pstot=0.0
 
-      do ig = 2, ngrid-1
+      do ig = 1, ngrid
 
          ! mass atm kg/m2
@@ -62,16 +62,8 @@
 
       enddo
-      atmmasstot= atmmasstot + ps(1)/g*cell_area(1)*iim
-      atmmasstot= atmmasstot + ps(ngrid)/g*cell_area(ngrid)*iim
-
-      surfmasstot= surfmasstot + surfmass(1)*cell_area(1)*iim
-      surfmasstot= surfmasstot + surfmass(ngrid)*cell_area(ngrid)*iim
-
-      pstot= pstot + ps(1)*cell_area(1)*iim
-      pstot= pstot + ps(ngrid)*cell_area(ngrid)*iim
                
-      atmmasstot=  atmmasstot/(totarea+(iim-1)*(cell_area(1)+cell_area(ngrid)))
-      pstot=  pstot/(totarea+(iim-1)*(cell_area(1)+cell_area(ngrid)))
-      surfmasstot = surfmasstot/(totarea+(iim-1)*(cell_area(1)+cell_area(ngrid)))
+      atmmasstot=  atmmasstot/totarea_planet
+      pstot=  pstot/totarea_planet
+      surfmasstot = surfmasstot/totarea_planet
 
       totmass=surfmasstot+atmmasstot