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

Timestamp:

Feb 27, 2024, 11:40:35 PM (10 months ago)

Author:

afalco

Message:

Pluto PCM:
outputs for CH4, CO, etc
AF

File:

• trunk/LMDZ.COMMON/libf/phy_common/geometry_mod.F90 (modified) (4 diffs)

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

@@ -36 +36 @@
 !$OMP THREADPRIVATE(cell_area_for_lonlat_outputs)
 
+  REAL,SAVE,ALLOCATABLE :: real_area(:) ! real area of the cell (m2)
+                                        ! multiplied by klon at poles
+  !$OMP THREADPRIVATE(real_area)
+
   INTEGER,SAVE,ALLOCATABLE :: ind_cell_glo(:)      ! global index of a local cell
 !$OMP THREADPRIVATE(ind_cell_glo)
@@ -47 +51 @@
   USE nrtype, ONLY : PI
   IMPLICIT NONE
+  include "dimensions.h"
     INTEGER,INTENT(IN) :: klon ! number of columns for this MPI/OpenMP domain
     REAL,INTENT(IN) :: longitude_(klon)
@@ -64 +69 @@
     ALLOCATE(boundslat(klon,nvertex))
     ALLOCATE(cell_area(klon))
+    ALLOCATE(real_area(klon))
     IF (PRESENT(ind_cell_glo_)) ALLOCATE(ind_cell_glo(klon))
     IF (PRESENT(dx_)) ALLOCATE(dx(klon))
@@ -75 +81 @@
     boundslat(:,:) = boundslat_(:,:)
     cell_area(:) = cell_area_(:)
+
+    real_area(2:klon-1)=cell_area(2:klon-1)
+    real_area(1)=cell_area(1)*iim
+    real_area(klon)=cell_area(klon)*iim
+
     IF (PRESENT(ind_cell_glo_)) ind_cell_glo(:) = ind_cell_glo_(:)
     IF (PRESENT(dx_)) dx(:) = dx_(:)