Index: /trunk/LMDZ.MARS/libf/phymars/planetwide_mod.F90
===================================================================
--- /trunk/LMDZ.MARS/libf/phymars/planetwide_mod.F90	(revision 1462)
+++ /trunk/LMDZ.MARS/libf/phymars/planetwide_mod.F90	(revision 1463)
@@ -20,4 +20,9 @@
 end interface
 
+interface planetwide_sumval ! sum() , over the entire planet
+  module procedure planetwide_sumval_i1, planetwide_sumval_i2, &
+                   planetwide_sumval_r1, planetwide_sumval_r2
+end interface
+
 contains
 
@@ -200,4 +205,94 @@
   end subroutine planetwide_minval_r2
 
-
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
+  subroutine planetwide_sumval_i1(values,values_sum)
+  use dimphy, only: klon
+  use mod_grid_phy_lmdz, only : klon_glo
+  implicit none
+  integer,intent(in) :: values(:) ! local grid (klon)
+  integer,intent(out) :: values_sum
+#ifdef CPP_PARA
+  integer :: values_glo(klon_glo) ! global grid
+  
+  ! gather field on master:
+  call gather(values,values_glo)
+  ! calculate sum value
+  if (is_master) then
+    values_sum=SUM(values_glo(:))
+  endif
+  ! broadcast information to all cores
+  call bcast(values_sum)
+#else
+  values_sum=SUM(values(:))
+#endif
+  end subroutine planetwide_sumval_i1
+
+  subroutine planetwide_sumval_i2(values,values_sum)
+  use dimphy, only: klon, klev
+  use mod_grid_phy_lmdz, only : klon_glo
+  implicit none
+  integer,intent(in) :: values(:,:) ! local grid (klon,klev)
+  integer,intent(out) :: values_sum
+#ifdef CPP_PARA
+  integer :: values_glo(klon_glo,klev) ! global grid
+  
+  ! gather field on master:
+  call gather(values,values_glo)
+  ! calculate sum value
+  if (is_master) then
+    values_sum=SUM(values_glo)
+  endif
+  ! broadcast information to all cores
+  call bcast(values_sum)
+#else
+  values_sum=SUM(values)
+#endif
+  end subroutine planetwide_sumval_i2
+
+  subroutine planetwide_sumval_r1(values,values_sum)
+  use dimphy, only: klon
+  use mod_grid_phy_lmdz, only : klon_glo
+  implicit none
+  real,intent(in) :: values(:) ! local grid (klon)
+  real,intent(out) :: values_sum
+#ifdef CPP_PARA
+  real :: values_glo(klon_glo) ! global grid
+  
+  ! gather field on master:
+  call gather(values,values_glo)
+  ! calculate sum value
+  if (is_master) then
+    values_sum=SUM(values_glo)
+  endif
+  ! broadcast information to all cores
+  call bcast(values_sum)
+#else
+  values_sum=SUM(values)
+#endif
+  end subroutine planetwide_sumval_r1
+
+  subroutine planetwide_sumval_r2(values,values_sum)
+  use dimphy, only: klon, klev
+  use mod_grid_phy_lmdz, only : klon_glo
+  implicit none
+  real,intent(in) :: values(:,:) ! local grid (klon,klev)
+  real,intent(out) :: values_sum
+#ifdef CPP_PARA
+  real :: values_glo(klon_glo,klev) ! global grid
+  
+  ! gather field on master:
+  call gather(values,values_glo)
+  ! calculate sum value
+  if (is_master) then
+    values_sum=SUM(values_glo)
+  endif
+  ! broadcast information to all cores
+  call bcast(values_sum)
+#else
+  values_sum=SUM(values)
+#endif
+  end subroutine planetwide_sumval_r2
+  
+  
 end module planetwide_mod