Changeset 1394 for LMDZ4/branches

Timestamp:

May 27, 2010, 12:07:16 PM (15 years ago)

Author:

Ehouarn Millour

Message:

OpenMP bug fix for Newtonian case.

EM

File:

• LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/leapfrog_p.F (modified) (2 diffs)

LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/leapfrog_p.F

@@ -974 +974 @@
        ijb=ij_begin
        ije=ij_end
-       teta(ijb:ije,:)=teta(ijb:ije,:)
-     s  -iphysiq*dtvr*(teta(ijb:ije,:)-tetarappel(ijb:ije,:))/taurappel
+!$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 
+       do l=1,llm
+       teta(ijb:ije,l)=teta(ijb:ije,l)
+     &  -iphysiq*dtvr*(teta(ijb:ije,l)-tetarappel(ijb:ije,l))/taurappel
+       enddo
+!$OMP END DO
 
        call Register_Hallo(ucov,ip1jmp1,llm,0,1,1,0,Request_Physic)
@@ -983 +987 @@
        call WaitRequest(Request_Physic)     
 c$OMP BARRIER
+!$OMP MASTER
        call friction_p(ucov,vcov,iphysiq*dtvr)
+!$OMP END MASTER
+!$OMP BARRIER
       ENDIF ! of IF(iflag_phys.EQ.2)