Changeset 2549 for trunk/LMDZ.COMMON/libf/dyn3dpar

Timestamp:

Jul 9, 2021, 5:18:17 PM (3 years ago)

Author:

emillour

Message:

Common dynamics:
Fix a weird issue that only arises in mixed MPI/OpenMP mode with
deallocation (via deallocate_buffer called by WaitRequest?) if
only the halo of surface pressure is registered. Adding the registration
of a companion 3D field seems to solve the problem.
EM

File:

• trunk/LMDZ.COMMON/libf/dyn3dpar/leapfrog_p.F (modified) (1 diff)

trunk/LMDZ.COMMON/libf/dyn3dpar/leapfrog_p.F

@@ -1303 +1303 @@
         ! sponge layer
         if (callsponge) then
-          call Register_Hallo(ps,ip1jm,1,1,1,1,1,Request_Dissip)
+! NB: Call to Register_Hallo is in principle only required for ps
+!     but it seems that then WaitRequest crashes (only in mixed MPI/OpenMP)
+!     unless a 3D field has also been registered
+          call Register_Hallo(teta,ip1jmp1,llm,0,1,1,0,Request_Dissip)
+          call Register_Hallo(ps,ip1jmp1,1,0,1,1,0,Request_Dissip)
           call SendRequest(Request_Dissip)
 c$OMP BARRIER