Ignore:
Timestamp:
2015-11-24T18:34:11+01:00 (5 years ago)
Author:
acc
Message:

Branch dev_r5836_NOC3_vvl_by_default. Changes to ensure successful SETTE testing. Affects tranxt.F90, trcnxt.F90 and traqsr.F90 only. Reduces dangerous dependencies on fields set in haloes but introduces additional lbc_lnks instead.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90

    r5883 r5916  
    5050   USE agrif_opa_interp 
    5151#endif 
     52#  include "vectopt_loop_substitute.h90" 
    5253 
    5354   IMPLICIT NONE 
     
    9293      INTEGER, INTENT(in) ::   kt    ! ocean time-step index 
    9394      !! 
    94       INTEGER  ::   jk, jn    ! dummy loop indices 
    95       REAL(wp) ::   zfact     ! local scalars 
     95      INTEGER  ::   ji, jj, jk, jn   ! dummy loop indices 
     96      REAL(wp) ::   zfact            ! local scalars 
    9697      REAL(wp), POINTER, DIMENSION(:,:,:) ::  ztrdt, ztrds 
    9798      !!---------------------------------------------------------------------- 
     
    142143            END DO 
    143144         END DO 
     145         ! 
    144146      ELSE                                            ! Leap-Frog + Asselin filter time stepping 
    145147         ! 
     
    148150           &                                                                sbc_tsc, sbc_tsc_b, jpts )  ! non-linear free surface 
    149151         ENDIF 
     152         ! 
     153         DO jn = 1, jpts 
     154            CALL lbc_lnk( tsb(:,:,:,jn), 'T', 1._wp )  
     155            CALL lbc_lnk( tsn(:,:,:,jn), 'T', 1._wp ) 
     156            CALL lbc_lnk( tsa(:,:,:,jn), 'T', 1._wp ) 
     157         END DO 
    150158      ENDIF      
    151159      ! 
     
    220228         ! 
    221229         DO jk = 1, jpkm1 
    222             DO jj = 1, jpj 
    223                DO ji = 1, jpi 
     230            DO jj = 2, jpjm1 
     231               DO ji = fs_2, fs_jpim1 
    224232                  ztn = ptn(ji,jj,jk,jn)                                     
    225233                  ztd = pta(ji,jj,jk,jn) - 2. * ztn + ptb(ji,jj,jk,jn)      ! time laplacian on tracers 
     
    305313            zfact1 = atfp * p2dt(jk) 
    306314            zfact2 = zfact1 / rau0 
    307             DO jj = 1, jpj 
    308                DO ji = 1, jpi 
     315            DO jj = 2, jpjm1 
     316               DO ji = fs_2, fs_jpim1 
    309317                  ze3t_b = e3t_b(ji,jj,jk) 
    310318                  ze3t_n = e3t_n(ji,jj,jk) 
Note: See TracChangeset for help on using the changeset viewer.