Changeset 7698 for trunk/NEMOGCM/NEMO/OPA_SRC/DOM/domhgr.F90

Timestamp:

2017-02-18T10:02:03+01:00 (7 years ago)

Author:

mocavero

Message:

update trunk with OpenMP parallelization

File:

• trunk/NEMOGCM/NEMO/OPA_SRC/DOM/domhgr.F90 (modified) (3 diffs)

trunk/NEMOGCM/NEMO/OPA_SRC/DOM/domhgr.F90

@@ -40 +40 @@
    !!----------------------------------------------------------------------
    !! NEMO/OPA 3.7 , NEMO Consortium (2016)
-   !! $Id$ 
+   !! $Id$
    !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
@@ -117 +117 @@
       IF( iff == 0 ) THEN                 ! Coriolis parameter has not been defined 
          IF(lwp) WRITE(numout,*) '          Coriolis parameter calculated on the sphere from gphif & gphit'
-         ff_f(:,:) = 2. * omega * SIN( rad * gphif(:,:) )     ! compute it on the sphere at f-point
-         ff_t(:,:) = 2. * omega * SIN( rad * gphit(:,:) )     !    -        -       -    at t-point
+!$OMP PARALLEL DO schedule(static) private(jj, ji)
+         DO jj = 1, jpj
+            DO ji = 1, jpi
+               ff_f(ji,jj) = 2. * omega * SIN( rad * gphif(ji,jj) )     ! compute it on the sphere at f-point
+               ff_t(ji,jj) = 2. * omega * SIN( rad * gphit(ji,jj) )     !    -        -       -    at t-point
+            END DO
+         END DO
       ELSE
          IF( ln_read_cfg ) THEN
@@ -130 +135 @@
       !                             !==  associated horizontal metrics  ==!
       !
-      r1_e1t(:,:) = 1._wp / e1t(:,:)   ;   r1_e2t (:,:) = 1._wp / e2t(:,:)
-      r1_e1u(:,:) = 1._wp / e1u(:,:)   ;   r1_e2u (:,:) = 1._wp / e2u(:,:)
-      r1_e1v(:,:) = 1._wp / e1v(:,:)   ;   r1_e2v (:,:) = 1._wp / e2v(:,:)
-      r1_e1f(:,:) = 1._wp / e1f(:,:)   ;   r1_e2f (:,:) = 1._wp / e2f(:,:)
-      !
-      e1e2t (:,:) = e1t(:,:) * e2t(:,:)   ;   r1_e1e2t(:,:) = 1._wp / e1e2t(:,:)
-      e1e2f (:,:) = e1f(:,:) * e2f(:,:)   ;   r1_e1e2f(:,:) = 1._wp / e1e2f(:,:)
+!$OMP PARALLEL DO schedule(static) private(jj, ji)
+      DO jj = 1, jpj
+         DO ji = 1, jpi
+            r1_e1t(ji,jj) = 1._wp / e1t(ji,jj)   ;   r1_e2t (ji,jj) = 1._wp / e2t(ji,jj)
+            r1_e1u(ji,jj) = 1._wp / e1u(ji,jj)   ;   r1_e2u (ji,jj) = 1._wp / e2u(ji,jj)
+            r1_e1v(ji,jj) = 1._wp / e1v(ji,jj)   ;   r1_e2v (ji,jj) = 1._wp / e2v(ji,jj)
+            r1_e1f(ji,jj) = 1._wp / e1f(ji,jj)   ;   r1_e2f (ji,jj) = 1._wp / e2f(ji,jj)
+            !
+            e1e2t (ji,jj) = e1t(ji,jj) * e2t(ji,jj)   ;   r1_e1e2t(ji,jj) = 1._wp / e1e2t(ji,jj)
+            e1e2f (ji,jj) = e1f(ji,jj) * e2f(ji,jj)   ;   r1_e1e2f(ji,jj) = 1._wp / e1e2f(ji,jj)
+         END DO
+      END DO
       IF( ie1e2u_v == 0 ) THEN               ! u- & v-surfaces have not been defined
          IF(lwp) WRITE(numout,*) '          u- & v-surfaces calculated as e1 e2 product'
-         e1e2u (:,:) = e1u(:,:) * e2u(:,:)         ! compute them
-         e1e2v (:,:) = e1v(:,:) * e2v(:,:) 
+!$OMP PARALLEL DO schedule(static) private(jj, ji)
+         DO jj = 1, jpj
+            DO ji = 1, jpi
+               e1e2u (ji,jj) = e1u(ji,jj) * e2u(ji,jj)         ! compute them
+               e1e2v (ji,jj) = e1v(ji,jj) * e2v(ji,jj) 
+            END DO
+         END DO
       ELSE
          IF(lwp) WRITE(numout,*) '          u- & v-surfaces have been read in "mesh_mask" file:'
          IF(lwp) WRITE(numout,*) '                     grid size reduction in strait(s) is used'
       ENDIF
-      r1_e1e2u(:,:) = 1._wp / e1e2u(:,:)     ! compute their invert in any cases
-      r1_e1e2v(:,:) = 1._wp / e1e2v(:,:)
-      !   
-      e2_e1u(:,:) = e2u(:,:) / e1u(:,:)
-      e1_e2v(:,:) = e1v(:,:) / e2v(:,:)
+!$OMP PARALLEL DO schedule(static) private(jj, ji)
+      DO jj = 1, jpj
+         DO ji = 1, jpi
+            r1_e1e2u(ji,jj) = 1._wp / e1e2u(ji,jj)     ! compute their invert in any cases
+            r1_e1e2v(ji,jj) = 1._wp / e1e2v(ji,jj)
+            !   
+            e2_e1u(ji,jj) = e2u(ji,jj) / e1u(ji,jj)
+            e1_e2v(ji,jj) = e1v(ji,jj) / e2v(ji,jj)
+         END DO
+      END DO
       !
       !