New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 9045 for branches/2017/dev_merge_2017/NEMOGCM/NEMO – NEMO

Ignore:
Timestamp:
2017-12-14T13:03:31+01:00 (6 years ago)
Author:
timgraham
Message:

Add exlicit loops

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_ts.F90

    r9043 r9045  
    523523      ! Note that the "unclipped" bottom friction parameter is used even with explicit drag 
    524524      IF( ln_wd_il ) THEN 
    525         zu_frc(:,:) = zu_frc(:,:) + MAX(r1_hu_n(:,:) * zCdU_u(:,:),-1._wp / rdtbt) * zwx(:,:) *  wdrampu(:,:) 
    526         zv_frc(:,:) = zv_frc(:,:) + MAX(r1_hv_n(:,:) * zCdU_v(:,:),-1._wp / rdtbt) * zwy(:,:) *  wdrampv(:,:) 
     525         DO jj = 2, jpjm1 
     526            DO ji = fs_2, fs_jpim1   ! vector opt. 
     527               zu_frc(ji,jj) = zu_frc(ji,jj) + MAX(r1_hu_n(ji,jj) * zCdU_u(ji,jj),-1._wp / rdtbt) * zwx(ji,jj) *  wdrampu(ji,jj) 
     528               zv_frc(ji,jj) = zv_frc(ji,jj) + MAX(r1_hv_n(ji,jj) * zCdU_v(ji,jj),-1._wp / rdtbt) * zwy(ji,jj) *  wdrampv(ji,jj) 
     529            END DO 
     530         END DO 
    527531      ELSE 
    528         zu_frc(:,:) = zu_frc(:,:) + r1_hu_n(:,:) * zCdU_u(:,:) * zwx(:,:) 
    529         zv_frc(:,:) = zv_frc(:,:) + r1_hv_n(:,:) * zCdU_v(:,:) * zwy(:,:) 
     532         DO jj = 2, jpjm1 
     533            DO ji = fs_2, fs_jpim1   ! vector opt. 
     534               zu_frc(ji,jj) = zu_frc(ji,jj) + r1_hu_n(ji,jj) * zCdU_u(ji,jj) * zwx(ji,jj) 
     535               zv_frc(ji,jj) = zv_frc(ji,jj) + r1_hv_n(ji,jj) * zCdU_v(ji,jj) * zwy(ji,jj) 
     536            END DO 
     537         END DO 
    530538      END IF 
    531539      ! 
     
    10361044!jth do implicitly instead 
    10371045         IF ( .NOT. ll_wd ) THEN ! Revert to explicit for bit comparison tests in non wad runs 
    1038             zu_trd(:,:) = zu_trd(:,:) + zCdU_u(:,:) * un_e(:,:) * hur_e(:,:) 
    1039             zv_trd(:,:) = zv_trd(:,:) + zCdU_v(:,:) * vn_e(:,:) * hvr_e(:,:) 
     1046            DO jj = 2, jpjm1 
     1047               DO ji = fs_2, fs_jpim1   ! vector opt. 
     1048                  zu_trd(ji,jj) = zu_trd(ji,jj) + zCdU_u(ji,jj) * un_e(ji,jj) * hur_e(ji,jj) 
     1049                  zv_trd(ji,jj) = zv_trd(ji,jj) + zCdU_v(ji,jj) * vn_e(ji,jj) * hvr_e(ji,jj) 
     1050               END DO 
     1051            END DO 
    10401052         ENDIF  
    10411053         ! 
Note: See TracChangeset for help on using the changeset viewer.