Changeset 3790

Timestamp:

Dec 21, 2020, 4:43:24 PM (4 years ago)

Author:

oboucher

Message:

correct cell_area2D in the case of a rectangulat lat-lon grid for singularity at south and north poles
without this correction the calving flux sent to the coupler is not correct

File:

• LMDZ6/trunk/libf/phylmd/cpl_mod.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/cpl_mod.F90

@@ -242 +242 @@
     sum_error = sum_error + error
 
-
     CALL gather_omp(longitude_deg,rlon_mpi)
     CALL gather_omp(latitude_deg,rlat_mpi)
@@ -252 +251 @@
       CALL Grid1DTo2D_mpi(rlat_mpi,lat2D)
       CALL Grid1DTo2D_mpi(cell_area_mpi,cell_area2D)
+      !--the next line is required for lat-lon grid and should have no impact
+      !--for an unstructured grid for which nbp_lon=1
+      !--if north pole in process mpi then divide cell area of pole cell by number of replicates
+      IF (is_north_pole_dyn) cell_area2D(:,1)=cell_area2D(:,1)/FLOAT(nbp_lon)
+      !--if south pole in process mpi then divide cell area of pole cell by number of replicates
+      IF (is_south_pole_dyn) cell_area2D(:,jj_nb)=cell_area2D(:,jj_nb)/FLOAT(nbp_lon)
       mask_calving(:,:,:) = 0 
       WHERE ( lat2D >= 40) mask_calving(:,:,1) = 1
@@ -279 +284 @@
     ENDIF
     
-           
     IF (sum_error /= 0) THEN
        abort_message='Pb allocation variables couplees'
@@ -350 +354 @@
     ENDIF    ! is_sequential
     
-
 !*************************************************************************************
 ! compatibility test