Index: trunk/LMDZ.PLUTO/libf/phypluto/dyn1d/rcm1d.F
===================================================================
--- trunk/LMDZ.PLUTO/libf/phypluto/dyn1d/rcm1d.F	(revision 3814)
+++ trunk/LMDZ.PLUTO/libf/phypluto/dyn1d/rcm1d.F	(revision 3816)
@@ -14,6 +14,5 @@
       use comsoil_h, only: nsoilmx, layer, mlayer, inertiedat, volcapa
       use phyredem, only: physdem0,physdem1
-      use geometry_mod, only: init_geometry,
-     &                        init_geometry_cell_area_for_outputs
+      use geometry_mod, only: init_geometry
       use planete_mod, only: apoastr,periastr,year_day,peri_day,
      &         obliquit,z0,lmixmin,emin_turb,coefvis,coefir,
@@ -557,5 +556,4 @@
      &                   (/0.,0.,0.,0./),(/0.,0.,0.,0./),
      &                   cell_area,[1])
-      call init_geometry_cell_area_for_outputs(1,cell_area)
 ! Ehouarn: init_vertial_layers called later (because disvert not called yet)
 !      call init_vertical_layers(nlayer,preff,scaleheight,
Index: trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90
===================================================================
--- trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90	(revision 3814)
+++ trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90	(revision 3816)
@@ -2168,9 +2168,11 @@
       call write_output("emis","Emissivity","",emis)
       call write_output('phisfi','Surface Geopotential','m2s-2',phisfi)
-      if (grid_type == regular_lonlat) then
-          call write_output("area","Mesh area","m2", &
-                           cell_area_for_lonlat_outputs)
-        else ! unstructured grid (e.g. dynamico)
-          call write_output("area","Mesh area","m2",cell_area)
+      if (ngrid.ne.1) then
+         if (grid_type == regular_lonlat) then
+            call write_output("area","Mesh area","m2", &
+                              cell_area_for_lonlat_outputs)
+         else ! unstructured grid (e.g. dynamico)
+            call write_output("area","Mesh area","m2",cell_area)
+         endif
       endif