Changeset 3078 for trunk/LMDZ.COMMON/libf/phy_common

Timestamp:

Oct 10, 2023, 11:41:27 AM (15 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.COMMON/libf/phy_common/geometry_mod.F90 (modified) (2 diffs)

trunk/LMDZ.COMMON/libf/phy_common/geometry_mod.F90

@@ -27 +27 @@
 !$OMP THREADPRIVATE(dy)
 
-  REAL,SAVE,ALLOCATABLE :: cell_area(:)      ! area of the cell
+  REAL,SAVE,ALLOCATABLE :: cell_area(:) ! area of the cell (m2)
 !$OMP THREADPRIVATE(cell_area)
+
+  REAL,SAVE,ALLOCATABLE :: cell_area_for_lonlat_outputs(:) 
+                           ! for lon-lat outputs only: area of the cell
+                           ! but with polar values as in dyn grid
+                           ! (with replicated polar grid point)
+!$OMP THREADPRIVATE(cell_area_for_lonlat_outputs)
 
   INTEGER,SAVE,ALLOCATABLE :: ind_cell_glo(:)      ! global index of a local cell
@@ -76 +82 @@
 
 
+  SUBROUTINE init_geometry_cell_area_for_outputs(klon, &
+                                                 cell_area_for_lonlat_outputs_)
+  IMPLICIT NONE
+  INTEGER,INTENT(IN) :: klon ! number of columns for this MPI/OpenMP domain
+  REAL,INTENT(IN) :: cell_area_for_lonlat_outputs_(klon) ! tweaked lon-lat mesh
+                     ! cell areas where polar values are as on dyn lon-lat grid
+                     ! for outputs
+  
+  ALLOCATE(cell_area_for_lonlat_outputs(klon))
+  cell_area_for_lonlat_outputs(1:klon)=cell_area_for_lonlat_outputs_(1:klon)
+  
+  END SUBROUTINE init_geometry_cell_area_for_outputs
+
 END MODULE geometry_mod