- Timestamp:
- Feb 28, 2024, 12:26:25 PM (9 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.COMMON/libf/phy_common/geometry_mod.F90
r3241 r3242 36 36 !$OMP THREADPRIVATE(cell_area_for_lonlat_outputs) 37 37 38 REAL,SAVE,ALLOCATABLE :: real_area(:) ! real area of the cell (m2)39 ! multiplied by klon at poles40 !$OMP THREADPRIVATE(real_area)41 42 38 INTEGER,SAVE,ALLOCATABLE :: ind_cell_glo(:) ! global index of a local cell 43 39 !$OMP THREADPRIVATE(ind_cell_glo) … … 51 47 USE nrtype, ONLY : PI 52 48 IMPLICIT NONE 53 include "dimensions.h"54 49 INTEGER,INTENT(IN) :: klon ! number of columns for this MPI/OpenMP domain 55 50 REAL,INTENT(IN) :: longitude_(klon) … … 69 64 ALLOCATE(boundslat(klon,nvertex)) 70 65 ALLOCATE(cell_area(klon)) 71 ALLOCATE(real_area(klon))72 66 IF (PRESENT(ind_cell_glo_)) ALLOCATE(ind_cell_glo(klon)) 73 67 IF (PRESENT(dx_)) ALLOCATE(dx(klon)) … … 81 75 boundslat(:,:) = boundslat_(:,:) 82 76 cell_area(:) = cell_area_(:) 83 84 real_area(2:klon-1)=cell_area(2:klon-1)85 real_area(1)=cell_area(1)*iim86 real_area(klon)=cell_area(klon)*iim87 88 77 IF (PRESENT(ind_cell_glo_)) ind_cell_glo(:) = ind_cell_glo_(:) 89 78 IF (PRESENT(dx_)) dx(:) = dx_(:)
Note: See TracChangeset
for help on using the changeset viewer.