Changeset 3078 for trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

Timestamp:

Oct 10, 2023, 11:41:27 AM (17 months ago)

Author:

emillour

Message:

Common physics/dynamics:
Add, for lon-lat outputs, a cell_area_for_lonlat_outputs()
field where polar mesh areas are adapted to match those
of the dynamics grid (where polar mesh points are replicated
along longitudes).

Mars PCM:
Enable outputing cell area with XIOS and make sure the polar grid point
is correctly handled (i.e. on lon-lat grid outputs polar mesh area must
be adjusted to account for the replicated polar meshes).

EM

File:

• trunk/LMDZ.MARS/libf/phymars/physiq_mod.F (modified) (3 diffs)

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -50 +50 @@
      &                     mlayer,layer ! soil mid layer depths
       use geometry_mod, only: longitude, latitude, cell_area,
-     &                        longitude_deg 
+     &                        cell_area_for_lonlat_outputs,longitude_deg
       use comgeomfi_h, only: sinlon, coslon, sinlat, coslat
       use surfdat_h, only: phisfi, albedodat, zmea, zstd, zsig, zgam,
@@ -110 +110 @@
       use wxios, only: wxios_context_init, xios_context_finalize
 #endif
-      USE mod_grid_phy_lmdz, ONLY: grid_type, unstructured
+      USE mod_grid_phy_lmdz, ONLY: grid_type, unstructured, 
+     &                             regular_lonlat
       use ioipsl_getin_p_mod, only: getin_p
       use comslope_mod, ONLY: nslope,def_slope,def_slope_mean,
@@ -3172 +3173 @@
          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)
+         endif
          call write_output("tsurf","Surface temperature","K",
      &                  tsurf(:,iflat))