Changeset 970 for codes/icosagcm/trunk/src

Timestamp:

09/13/19 18:24:51 (5 years ago)

Author:

ymipsl

Message:

Bug fix in openMP for held&suarez physic.

YM

File:

• codes/icosagcm/trunk/src/initial/etat0_heldsz.f90 (modified) (3 diffs)

codes/icosagcm/trunk/src/initial/etat0_heldsz.f90

@@ -226 +226 @@
   SUBROUTINE compute_heldsz(ps,theta_eq,lat, theta_rhodz,u, theta) 
     USE theta2theta_rhodz_mod
+    USE omp_para
     REAL(rstd),INTENT(IN)    :: ps(iim*jjm) 
     REAL(rstd),INTENT(IN)    :: theta_eq(iim*jjm,llm) 
@@ -236 +237 @@
 
     CALL compute_theta_rhodz2theta(ps,theta_rhodz,theta,1)
-    DO l=1,llm
+    DO l=ll_begin,ll_end
        DO j=jj_begin-1,jj_end+1
           DO i=ii_begin-1,ii_end+1
@@ -247 +248 @@
     CALL compute_theta2theta_rhodz(ps,theta,theta_rhodz,1)
 
-    Do l=1,llm
+    Do l=ll_begin,ll_end
        u(:,l)=u(:,l)*(1.-itau_physics*dt*kfrict(l))
     END DO