Ignore:
Timestamp:
2014-02-04T13:11:35+01:00 (7 years ago)
Author:
trackstand2
Message:

Put jpkf into traldf_iso and traadv_tvd (inc nonosc)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_tvd.F90

    r3432 r4401  
    130130      END IF 
    131131      ! 
    132       zwi(:,:,:) = 0.e0 
     132      zwi(:,:,1:jpkf) = 0.e0 
    133133      ! 
    134134      !                                                          ! =========== 
     
    137137         ! 1. Bottom value : flux set to zero 
    138138         ! ---------------------------------- 
    139          zwx(:,:,jpk) = 0.e0    ;    zwz(:,:,jpk) = 0.e0 
    140          zwy(:,:,jpk) = 0.e0    ;    zwi(:,:,jpk) = 0.e0 
     139         ! ARPDBG: FINISS was using jpk correct here given that it can be  
     140         ! below the ocean floor?? 
     141         zwx(:,:,jpkf) = 0.e0    ;    zwz(:,:,jpkf) = 0.e0 
     142         zwy(:,:,jpkf) = 0.e0    ;    zwi(:,:,jpkf) = 0.e0 
    141143 
    142144         ! 2. upstream advection with initial mass fluxes & intermediate update 
     
    147149         DO jj = 1, jpjm1 
    148150            DO ji = 1, jpim1 
    149                DO jk = 1, jpkm1 
    150 #else 
    151          DO jk = 1, jpkm1 
     151               DO jk = 1, jpkfm1 
     152#else 
     153         DO jk = 1, jpkfm1 
    152154            DO jj = 1, jpjm1 
    153155               DO ji = 1, fs_jpim1   ! vector opt. 
     
    175177         DO jj = 1, jpj 
    176178            DO ji = 1, jpi 
    177                DO jk = 2, jpkm1 
    178 #else 
    179          DO jk = 2, jpkm1 
     179               DO jk = 2, jpkfm1 
     180#else 
     181         DO jk = 2, jpkfm1 
    180182            DO jj = 1, jpj 
    181183               DO ji = 1, jpi 
     
    194196         DO jj = 2, jpjm1 
    195197            DO ji = 2, jpim1 
    196                DO jk = 1, jpkm1 
     198               DO jk = 1, jpkfm1 
    197199                  z2dtt = p2dt(jk) 
    198200#else 
    199          DO jk = 1, jpkm1 
     201         DO jk = 1, jpkfm1 
    200202            z2dtt = p2dt(jk) 
    201203            DO jj = 2, jpjm1 
     
    223225         IF( l_trd )  THEN  
    224226            ! store intermediate advective trends 
    225             ztrdx(:,:,:) = zwx(:,:,:)   ;    ztrdy(:,:,:) = zwy(:,:,:)  ;   ztrdz(:,:,:) = zwz(:,:,:) 
     227            ztrdx(:,:,1:jpkf) = zwx(:,:,1:jpkf)   ;    ztrdy(:,:,1:jpkf) = zwy(:,:,1:jpkf)  ;   ztrdz(:,:,1:jpkf) = zwz(:,:,1:jpkf) 
    226228         END IF 
    227229         !                                 ! "Poleward" heat and salt transports (contribution of upstream fluxes) 
     
    238240         DO jj = 1, jpjm1 
    239241            DO ji = 1, jpim1 
    240                DO jk = 1, jpkm1 
    241 #else 
    242          DO jk = 1, jpkm1 
     242               DO jk = 1, jpkfm1 
     243#else 
     244         DO jk = 1, jpkfm1 
    243245            DO jj = 1, jpjm1 
    244246               DO ji = 1, fs_jpim1   ! vector opt. 
     
    257259            DO ji = 1, jpi 
    258260               zwz(ji,jj,1) = 0.e0   ! Surface value 
    259                DO jk = 2, jpkm1 
     261               DO jk = 2, jpkfm1 
    260262#else 
    261263         zwz(:,:,1) = 0.e0           ! Surface value 
    262264         ! 
    263          DO jk = 2, jpkm1            ! Interior value 
     265         DO jk = 2, jpkfm1            ! Interior value 
    264266            DO jj = 1, jpj 
    265267               DO ji = 1, jpi 
     
    289291         DO jj = 2, jpjm1 
    290292            DO ji = 2, jpim1 
    291                DO jk = 1, jpkm1 
    292 #else 
    293          DO jk = 1, jpkm1 
     293               DO jk = 1, jpkfm1 
     294#else 
     295         DO jk = 1, jpkfm1 
    294296            DO jj = 2, jpjm1 
    295297               DO ji = fs_2, fs_jpim1   ! vector opt.   
     
    390392      zbig  = 1.e+40_wp 
    391393      zrtrn = 1.e-15_wp 
    392       zbetup(:,:,jpk) = 0._wp   ;   zbetdo(:,:,jpk) = 0._wp 
     394      zbetup(:,:,jpkf) = 0._wp   ;   zbetdo(:,:,jpkf) = 0._wp 
    393395 
    394396 
     
    404406      DO jj = 2, jpjm1 
    405407         DO ji = 2, jpim1 
    406             DO jk = 1, jpkm1 
     408            DO jk = 1, jpkfm1 
    407409               ikm1 = MAX(jk-1,1) 
    408410               z2dtt = p2dt(jk) 
    409411#else 
    410       DO jk = 1, jpkm1 
     412      DO jk = 1, jpkfm1 
    411413         ikm1 = MAX(jk-1,1) 
    412414         z2dtt = p2dt(jk) 
     
    455457      DO jj = 2, jpjm1 
    456458         DO ji = 2, jpim1 
    457             DO jk = 1, jpkm1 
    458 #else 
    459       DO jk = 1, jpkm1 
     459            DO jk = 1, jpkfm1 
     460#else 
     461      DO jk = 1, jpkfm1 
    460462         DO jj = 2, jpjm1 
    461463            DO ji = fs_2, fs_jpim1   ! vector opt. 
Note: See TracChangeset for help on using the changeset viewer.