Changeset 4444
- Timestamp:
- 2014-02-04T13:15:42+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_tvd.F90
r4409 r4444 130 130 END IF 131 131 ! 132 zwi(:,:,1:jpkf) = 0.e0 132 #if defined key_z_first 133 DO jj = 1, jpj, 1 134 DO ji = 1, jpi, 1 135 zwi(ji,jj,1:mbkmax(ji,jj)) = 0.e0_wp 136 END DO 137 END DO 138 #else 139 zwi(:,:,1:jpk) = 0.e0_wp 140 #endif 133 141 ! 134 142 ! ! =========== … … 139 147 ! ARPDBG: FINISS was using jpk correct here given that it can be 140 148 ! below the ocean floor?? 149 #if defined key_z_first 150 DO jj = 1, jpj, 1 151 DO ji = 1, jpi, 1 152 zwx(ji,jj,mbkmax(ji,jj)) = 0.0_wp 153 zwy(ji,jj,mbkmax(ji,jj)) = 0.0_wp 154 zwz(ji,jj,mbkmax(ji,jj)) = 0.0_wp 155 zwi(ji,jj,mbkmax(ji,jj)) = 0.0_wp 156 END DO 157 END DO 158 #else 141 159 zwx(:,:,jpkf) = 0.e0 ; zwz(:,:,jpkf) = 0.e0 142 160 zwy(:,:,jpkf) = 0.e0 ; zwi(:,:,jpkf) = 0.e0 143 161 #endif 144 162 ! 2. upstream advection with initial mass fluxes & intermediate update 145 163 ! -------------------------------------------------------------------- … … 149 167 DO jj = 1, jpjm1 150 168 DO ji = 1, jpim1 151 DO jk = 1, jpkfm1169 DO jk = 1, mbkmax(ji,jj)-1 ! jpkfm1 152 170 #else 153 171 DO jk = 1, jpkfm1 … … 177 195 DO jj = 1, jpj 178 196 DO ji = 1, jpi 179 DO jk = 2, jpkfm1197 DO jk = 2, mbkmax(ji,jj)-1 ! jpkfm1 180 198 #else 181 199 DO jk = 2, jpkfm1 … … 196 214 DO jj = 2, jpjm1 197 215 DO ji = 2, jpim1 198 DO jk = 1, jpkfm1216 DO jk = 1, mbkmax(ji,jj)-1 ! jpkfm1 199 217 z2dtt = p2dt(jk) 200 218 #else … … 240 258 DO jj = 1, jpjm1 241 259 DO ji = 1, jpim1 242 DO jk = 1, jpkfm1260 DO jk = 1, mbkmax(ji,jj)-1 ! jpkfm1 243 261 #else 244 262 DO jk = 1, jpkfm1 … … 259 277 DO ji = 1, jpi 260 278 zwz(ji,jj,1) = 0.e0 ! Surface value 261 DO jk = 2, jpkfm1279 DO jk = 2, mbkmax(ji,jj)-1 ! jpkfm1 262 280 #else 263 281 zwz(:,:,1) = 0.e0 ! Surface value … … 291 309 DO jj = 2, jpjm1 292 310 DO ji = 2, jpim1 293 DO jk = 1, jpkfm1311 DO jk = 1, mbkmax(ji,jj)-1 ! jpkfm1 294 312 #else 295 313 DO jk = 1, jpkfm1
Note: See TracChangeset
for help on using the changeset viewer.