Changeset 3816

Timestamp:

Jun 25, 2025, 11:10:57 AM (4 days 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

Files:

• LMDZ.GENERIC/libf/phystd/physiq_mod.F90 (modified) (3 diffs)
• LMDZ.PLUTO/libf/phypluto/dyn1d/rcm1d.F (modified) (2 diffs)
• LMDZ.PLUTO/libf/phypluto/physiq_mod.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

@@ -31 +31 @@
       use comsaison_h, only: mu0, fract, dist_star, declin, right_ascen
       use comsoil_h, only: nsoilmx, layer, mlayer, inertiedat
-      use geometry_mod, only: latitude, longitude, cell_area
+      use geometry_mod, only: latitude, longitude, cell_area, &
+                          cell_area_for_lonlat_outputs
       USE comgeomfi_h, only: totarea, totarea_planet
       USE tracer_h, only: noms, mmol, radius, rho_q, qext, &
@@ -55 +56 @@
       use planetwide_mod, only: planetwide_minval,planetwide_maxval,planetwide_sumval
       use mod_phys_lmdz_para, only : is_master
+      USE mod_grid_phy_lmdz, ONLY: regular_lonlat, grid_type, unstructured
       use planete_mod, only: apoastr, periastr, year_day, peri_day, &
                             obliquit, nres, z0
@@ -2546 +2548 @@
       call write_output("p","Pressure","Pa",pplay)
       call write_output("omega","omega","Pa/s",omega)
-
+      call write_output('phisfi','Surface Geopotential','m2s-2',phisfi)
+      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
+      
 !     Subsurface temperatures
 !        call writediagsoil(ngrid,"tsurf","Surface temperature","K",2,tsurf)

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