Changeset 3816 for trunk/LMDZ.PLUTO/libf

Timestamp:

Jun 25, 2025, 11:10:57 AM (3 weeks ago)

Author:

afalco

Message:

Generic: enable writing cell_area with XIOSin physiq_mod, only for 3D runs.
Pluto: only write variable in 3D runs (do not calculate cell_area for 1D runs).
AF

Location:

trunk/LMDZ.PLUTO/libf/phypluto

Files:

• dyn1d/rcm1d.F (modified) (2 diffs)
• physiq_mod.F90 (modified) (1 diff)

trunk/LMDZ.PLUTO/libf/phypluto/dyn1d/rcm1d.F

@@ -14 +14 @@
       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 +556 @@
      &                   (/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,

trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90

@@ -2168 +2168 @@
       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