Changeset 295 for codes/icosagcm/trunk/src/wind.f90

Timestamp:

10/31/14 14:52:01 (10 years ago)

Author:

ymipsl

Message:

Merging OpenMP parallisme mode : by subdomain and on vertical level.
This feature is actually experimental but may be retro-compatible with the last method based only on subdomain

YM

File:

• codes/icosagcm/trunk/src/wind.f90 (modified) (20 diffs)

codes/icosagcm/trunk/src/wind.f90

@@ -47 +47 @@
   SUBROUTINE compute_wind_centered(ue,ucenter)
   USE icosa
-  
+  USE omp_para
   IMPLICIT NONE
   REAL(rstd) :: ue(3*iim*jjm,llm)
@@ -53 +53 @@
   INTEGER :: i,j,ij,l    
  
-    DO l=1,llm
+    DO l=ll_begin,ll_end
       DO j=jj_begin,jj_end
         DO i=ii_begin,ii_end
@@ -73 +73 @@
  SUBROUTINE compute_wind_on_edge(ue,uedge)
   USE icosa
+  USE omp_para
     
   IMPLICIT NONE
@@ -83 +84 @@
     CALL compute_tangential_compound(ue,ut)
   
-    DO l=1,llm
+    DO l=ll_begin,ll_end
       DO j=jj_begin,jj_end
         DO i=ii_begin,ii_end
@@ -100 +101 @@
  SUBROUTINE compute_tangential_compound(ue,ut)
   USE icosa  
+  USE omp_para
   IMPLICIT NONE
   REAL(rstd) :: ue(3*iim*jjm,llm)
@@ -105 +107 @@
   INTEGER :: i,j,l,ij
     
-  DO l=1,llm
+  DO l=ll_begin,ll_end
     DO j=jj_begin,jj_end
       DO i=ii_begin,ii_end
@@ -155 +157 @@
  SUBROUTINE compute_wind_lonlat_compound(u, ulon, ulat)
   USE icosa  
+  USE omp_para
     
   IMPLICIT NONE
@@ -164 +167 @@
     
   
-    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
@@ -184 +187 @@
   SUBROUTINE compute_wind_from_lonlat_compound(ulon, ulat, u)
   USE icosa  
+  USE omp_para
     
   IMPLICIT NONE
@@ -192 +196 @@
   INTEGER :: i,j,ij,l     
   
-    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
@@ -207 +211 @@
   SUBROUTINE compute_wind_centered_from_lonlat_compound(ulon, ulat, u)
   USE icosa  
+  USE omp_para
     
   IMPLICIT NONE
@@ -214 +219 @@
 
   INTEGER :: i,j,ij,l     
-  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
@@ -248 +253 @@
   SUBROUTINE compute_wind_perp_from_lonlat_compound(ulon, ulat, up)
   USE icosa  
+  USE omp_para
     
   IMPLICIT NONE
@@ -259 +265 @@
    CALL compute_wind_from_lonlat_compound(ulon, ulat, u)
 
-    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
@@ -297 +303 @@
  SUBROUTINE compute_wind_centered_lonlat_compound(uc, ulon, ulat)
   USE icosa  
+  USE omp_para
     
   IMPLICIT NONE
@@ -306 +313 @@
     
   
-    DO l=1,llm
+    DO l=ll_begin,ll_end
       DO j=jj_begin,jj_end
         DO i=ii_begin,ii_end
@@ -320 +327 @@
  SUBROUTINE compute_wind_centered_from_wind_lonlat_centered(ulon, ulat,uc)
   USE icosa  
+  USE omp_para
     
   IMPLICIT NONE
@@ -329 +337 @@
     
   
-    DO l=1,llm
+    DO l=ll_begin,ll_end
       DO j=jj_begin,jj_end
         DO i=ii_begin,ii_end
@@ -344 +352 @@
  SUBROUTINE compute_wind_perp_from_wind_centered(uc,un)
   USE icosa  
+  USE omp_para
     
   IMPLICIT NONE
@@ -352 +361 @@
     
   
-    DO l=1,llm
+    DO l=ll_begin,ll_end
       DO j=jj_begin,jj_end
         DO i=ii_begin,ii_end