Changeset 985 for LMDZ4/trunk/libf/dyn3dpar/gradiv2_p.F

Timestamp:

Jul 30, 2008, 5:50:03 PM (16 years ago)

Author:

Laurent Fairhead

Message:

Mise a jour de dyn3dpar par rapport a dyn3d, inclusion OpenMP et filtre FFT YM
LF

File:

• LMDZ4/trunk/libf/dyn3dpar/gradiv2_p.F (modified) (6 diffs)

LMDZ4/trunk/libf/dyn3dpar/gradiv2_p.F

@@ -16 +16 @@
       USE times
       USE Write_field_p
+      USE mod_hallo
       IMPLICIT NONE
 c
@@ -33 +34 @@
 c
       REAL,SAVE :: div(ip1jmp1,llm)
+      REAL      :: tmp_div2(ip1jmp1,llm)
       REAL signe, nugrads
       INTEGER l,ij,iter,ld
       INTEGER :: ijb,ije,jjb,jje
+      Type(Request)  :: request_dissip
       
 c    ........................................................
@@ -63 +66 @@
 
 c$OMP BARRIER
-c$OMP MASTER      
-      call suspend_timer(timer_dissip)
-      call exchange_Hallo(gdy,ip1jm,llm,1,0)
-      call resume_timer(timer_dissip)
-c$OMP END MASTER
+       call Register_Hallo(gdy,ip1jm,llm,1,0,0,1,Request_dissip)
+       call SendRequest(Request_dissip)
+c$OMP BARRIER
+       call WaitRequest(Request_dissip)
 c$OMP BARRIER
 c
@@ -81 +83 @@
       IF( ld.GT.1 )   THEN
 c$OMP BARRIER
-c$OMP MASTER        
-        call suspend_timer(timer_dissip)
-        call exchange_Hallo(div,ip1jmp1,llm,1,1)
-        call resume_timer(timer_dissip)
-c$OMP END MASTER        
+       call Register_Hallo(div,ip1jmp1,llm,1,1,1,1,Request_dissip)
+       call SendRequest(Request_dissip)
+c$OMP BARRIER
+       call WaitRequest(Request_dissip)
 c$OMP BARRIER
         CALL laplacien_p ( klevel, div,  div     )
@@ -94 +95 @@
         DO iter = 1, ld -2
 c$OMP BARRIER
-c$OMP MASTER
-         call suspend_timer(timer_dissip)
-         call exchange_Hallo(div,ip1jmp1,llm,1,1)
-         call resume_timer(timer_dissip)
-c$OMP END MASTER
+       call Register_Hallo(div,ip1jmp1,llm,1,1,1,1,Request_dissip)
+       call SendRequest(Request_dissip)
 c$OMP BARRIER
+       call WaitRequest(Request_dissip)
+
+c$OMP BARRIER
+
          CALL laplacien_gam ( klevel,cuvscvgam1,cvuscugam1,unsair_gam1,
      *                       unsapolnga1, unsapolsga1,  div, div       )
@@ -112 +114 @@
 c       call exchange_Hallo(div,ip1jmp1,llm,0,1)
 c$OMP BARRIER
-c$OMP MASTER        
-        call suspend_timer(timer_dissip)
-        call exchange_Hallo(div,ip1jmp1,llm,1,1)
-        call resume_timer(timer_dissip)
-c$OMP END MASTER
+       call Register_Hallo(div,ip1jmp1,llm,1,1,1,1,Request_dissip)
+       call SendRequest(Request_dissip)
 c$OMP BARRIER
-c       call write_field3d_p('div4',reshape(div,(/iip1,jjp1,llm/)))
+       call WaitRequest(Request_dissip)
+
+c$OMP BARRIER
+
+
        CALL  grad_p  ( klevel,  div,   gdx,  gdy             )