Changeset 5916


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.

Location:
branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO
Files:
3 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) 
  • branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/OPA_SRC/TRA/traqsr.F90

    r5883 r5916  
    3232   USE restart        ! ocean restart 
    3333   USE lib_mpp        ! MPP library 
     34   USE lbclnk         ! ocean lateral boundary conditions (or mpp link) 
    3435   USE wrk_nemo       ! Memory Allocation 
    3536   USE timing         ! Timing 
     
    255256            END DO 
    256257         END DO 
     258         ! Update haloes since lim_thd needs fraqsr_1lev to be defined everywhere 
     259         CALL lbc_lnk( fraqsr_1lev(:,:), 'T', 1._wp ) 
    257260      ENDIF 
    258261      ! 
  • branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/TOP_SRC/TRP/trcnxt.F90

    r5883 r5916  
    133133           &                                                                   sbc_trc, sbc_trc_b, jptra )  ! non-linear ssh 
    134134         ENDIF 
     135         ! 
     136         DO jn = 1, jptra 
     137            CALL lbc_lnk( trb(:,:,:,jn), 'T', 1._wp )  
     138            CALL lbc_lnk( trn(:,:,:,jn), 'T', 1._wp ) 
     139            CALL lbc_lnk( tra(:,:,:,jn), 'T', 1._wp ) 
     140         END DO 
    135141      ENDIF 
    136142      ! 
Note: See TracChangeset for help on using the changeset viewer.