Changeset 3816 for trunk/LMDZ.GENERIC/libf

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

File:

• trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90 (modified) (3 diffs)

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)