Changeset 889

Timestamp:

06/13/19 10:12:33 (5 years ago)

Author:

ymipsl

Message:

OpenMP fix when using "phys_smooth_tendency=y"

YM

File:

• codes/icosagcm/trunk/src/physics/physics_external.F90 (modified) (6 diffs)

codes/icosagcm/trunk/src/physics/physics_external.F90

@@ -34 +34 @@
 
     CALL initialize_external_physics
+!$OMP PARALLEL    
     CALL allocate_field(f_theta_rhodz0, field_t, type_real, llm, nqdyn, name='theta_rhodz0')
     CALL allocate_field(f_u0,field_u,type_real,llm,name='u0')
@@ -46 +47 @@
     phys_smooth_tendency=.FALSE.
     CALL getin("phys_smooth_tendency",phys_smooth_tendency)
-    
+!$OMP END PARALLEL    
         
   END SUBROUTINE init_physics
@@ -88 +89 @@
 
     IF (phys_smooth_tendency) THEN
-    
+
       IF (MOD(it_,itau_physics)==1) THEN
         DO ind=1, ndomain
@@ -103 +104 @@
           rhodz=f_rhodz(ind)
        
-          theta_rhodz0(:,:,1)=theta_rhodz(:,:,1)
-          u0=u
-          q0=q
+          theta_rhodz0(:,ll_begin:ll_end,1)=theta_rhodz(:,ll_begin:ll_end,1)
+          u0(:,ll_begin:ll_end)=u(:,ll_begin:ll_end)
+          q0(:,ll_begin:ll_end,:)=q(:,ll_begin:ll_end,:)
           CALL compute_rhodz(.TRUE., ps, rhodz)
         ENDDO
@@ -129 +130 @@
           rhodz=f_rhodz(ind)
        
-          dtheta_rhodz(:,:,1)=(theta_rhodz(:,:,1)-theta_rhodz0(:,:,1))/itau_physics
-          du=(u-u0)/itau_physics
-   
+          dtheta_rhodz(:,ll_begin:ll_end,1)=(theta_rhodz(:,ll_begin:ll_end,1)-theta_rhodz0(:,ll_begin:ll_end,1))/itau_physics
+          
+          du(:,ll_begin:ll_end)=(u(:,ll_begin:ll_end)-u0(:,ll_begin:ll_end))/itau_physics
+         
           DO iq=1, nqtot
-            dq(:,:,iq)=((q(:,:,iq)-q0(:,:,iq))/itau_physics)*rhodz(:,:)
+            dq(:,ll_begin:ll_end,iq)=((q(:,ll_begin:ll_end,iq)-q0(:,ll_begin:ll_end,iq))/itau_physics)*rhodz(:,ll_begin:ll_end)
           ENDDO
          
-          theta_rhodz(:,:,1)=theta_rhodz0(:,:,1)
-          u=u0
-          q=q0
+          theta_rhodz(:,ll_begin:ll_end,1)=theta_rhodz0(:,ll_begin:ll_end,1)
+          u(:,ll_begin:ll_end)=u0(:,ll_begin:ll_end)
+          q(:,ll_begin:ll_end,:)=q0(:,ll_begin:ll_end,:)
         ENDDO
      ENDIF
@@ -156 +158 @@
        ps=f_ps(ind)
 
-       u=u+du
-       theta_rhodz=theta_rhodz+dtheta_rhodz
+       u(:,ll_begin:ll_end)=u(:,ll_begin:ll_end)+du(:,ll_begin:ll_end)
+       theta_rhodz(:,ll_begin:ll_end,1)=theta_rhodz(:,ll_begin:ll_end,1)+dtheta_rhodz(:,ll_begin:ll_end,1)
        CALL compute_rhodz(.TRUE., ps, rhodz)
        DO iq=1, nqtot
-         q(:,:,iq)=q(:,:,iq)+dq(:,:,iq)/rhodz(:,:)
+         q(:,ll_begin:ll_end,iq)=q(:,ll_begin:ll_end,iq)+dq(:,ll_begin:ll_end,iq)/rhodz(:,ll_begin:ll_end)
        ENDDO
      ENDDO
+!$OMP BARRIER
    
    ELSE