Changeset 595 for codes/icosagcm/devel/src/diagnostics

Timestamp:

10/19/17 10:32:52 (7 years ago)

Author:

dubos

Message:

devel : fix halo issues with computation of energy fluxes

Location:

codes/icosagcm/devel/src/diagnostics

Files:

• diagflux.F90 (modified) (1 diff)
• observable.f90 (modified) (1 diff)

codes/icosagcm/devel/src/diagnostics/diagflux.F90

@@ -150 +150 @@
     INTEGER :: ij, l, ij_omp_begin_ext, ij_omp_end_ext
     Rd = kappa*cpp
-    CALL distrib_level(ij_end_ext-ij_begin_ext+1,ij_omp_begin_ext,ij_omp_end_ext)
-    ij_omp_begin_ext = ij_omp_begin_ext+ij_begin_ext-1
-    ij_omp_end_ext = ij_omp_end_ext+ij_begin_ext-1
+    ! even if loops are of the _ext variant, we still need halo exchanges before reconstructing fluxes at cell centers
+    ! => loop over interior region
+    CALL distrib_level(ij_end-ij_begin+1,ij_omp_begin_ext,ij_omp_end_ext)
+    ij_omp_begin_ext = ij_omp_begin_ext+ij_begin-1
+    ij_omp_end_ext = ij_omp_end_ext+ij_begin-1
 #include "../kernels/energy_fluxes.k90"
   END SUBROUTINE compute_diagflux_energy

codes/icosagcm/devel/src/diagnostics/observable.f90

@@ -167 +167 @@
           CALL output_field("massflux_lon",f_buf_ulon)
           CALL output_field("massflux_lat",f_buf_ulat)
+
+          CALL transfert_request(f_epotfluxt,req_e1_vect) 
           CALL flux_centered_lonlat(1./(itau_out*dt) , f_epotfluxt, f_buf_ulon, f_buf_ulat)
           CALL output_field("epot_t", f_epot)
           CALL output_field("epotflux_lon",f_buf_ulon)
           CALL output_field("epotflux_lat",f_buf_ulat)
+
+          CALL transfert_request(f_ekinfluxt,req_e1_vect) 
           CALL flux_centered_lonlat(1./(itau_out*dt) , f_ekinfluxt, f_buf_ulon, f_buf_ulat)
           CALL output_field("ekin_t", f_ekin)
           CALL output_field("ekinflux_lon",f_buf_ulon)
           CALL output_field("ekinflux_lat",f_buf_ulat)
+
+          CALL transfert_request(f_enthalpyfluxt,req_e1_vect) 
           CALL flux_centered_lonlat(1./(itau_out*dt) , f_enthalpyfluxt, f_buf_ulon, f_buf_ulat)
           CALL output_field("enthalpy_t", f_enthalpy)
           CALL output_field("enthalpyflux_lon",f_buf_ulon)
           CALL output_field("enthalpyflux_lat",f_buf_ulat)
+
           CALL qflux_centered_lonlat(1./(itau_out*dt) , f_qfluxt, f_qfluxt_lon, f_qfluxt_lat)
           CALL output_field("qmass_t", f_qmasst)