Changeset 4451 for branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90

Timestamp:

2014-02-04T13:16:21+01:00 (10 years ago)

Author:

trackstand2

Message:

Add use of mbkmax to tra_nxt and traswp

File:

• branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90 (modified) (6 diffs)

branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90

@@ -48 +48 @@
    USE agrif_opa_interp
 #endif
+   USE timing
 
    IMPLICIT NONE
@@ -108 +109 @@
       !!----------------------------------------------------------------------
 
+      CALL timing_start('tra_nxt')
+
       IF( kt == nit000 ) THEN
          IF(lwp) WRITE(numout,*)
@@ -155 +158 @@
             DO jj = 1, jpj
                DO ji = 1, jpi
-                  DO jk = 1, jpkm1
+                  DO jk = 1, mbkmax(ji,jj)-1 ! jpkm1
                      tsn(ji,jj,jk,jn) = tsa(ji,jj,jk,jn)    
                   END DO
@@ -185 +188 @@
          DO jj = 1, jpj
             DO ji = 1, jpi
-               DO jk = 1, jpkm1
+               DO jk = 1, mbkmax(ji,jj)-1 ! jpkm1
                   zfact = 1.e0 / r2dtra(jk)             
                   ztrdt(ji,jj,jk) = ( tsb(ji,jj,jk,jp_tem) - ztrdt(ji,jj,jk) ) * zfact
@@ -207 +210 @@
       IF(ln_ctl)   CALL prt_ctl( tab3d_1=tsn(:,:,:,jp_tem), clinfo1=' nxt  - Tn: ', mask1=tmask,   &
          &                       tab3d_2=tsn(:,:,:,jp_sal), clinfo2=       ' Sn: ', mask2=tmask )
+      !
+      CALL timing_stop('tra_nxt','section')
       !
    END SUBROUTINE tra_nxt
@@ -371 +376 @@
          DO jj = 1, jpj
             DO ji = 1, jpi
-               DO jk = 1, jpkm1
+               DO jk = 1, mbkmax(ji,jj)-1 ! jpkm1
                   !! DCSE_NEMO: could try promoting these scalars to vectors
                   zfact1 = atfp * rdttra(jk)