Changeset 7753 for trunk/NEMOGCM/NEMO/TOP_SRC/trcini.F90
- Timestamp:
- 2017-03-03T12:46:59+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/TOP_SRC/trcini.F90
r7698 r7753 105 105 !! ** Purpose : passive tracers inventories at initialsation phase 106 106 !!---------------------------------------------------------------------- 107 INTEGER :: jk, jn , jj, ji! dummy loop indices107 INTEGER :: jk, jn ! dummy loop indices 108 108 CHARACTER (len=25) :: charout 109 109 !!---------------------------------------------------------------------- 110 110 ! ! masked grid volume 111 !$OMP PARALLEL112 !$OMP DO schedule(static) private(jk,jj,ji)113 111 DO jk = 1, jpk 114 DO jj = 1, jpj 115 DO ji = 1, jpi 116 cvol(ji,jj,jk) = e1e2t(ji,jj) * e3t_n(ji,jj,jk) * tmask(ji,jj,jk) 117 END DO 118 END DO 112 cvol(:,:,jk) = e1e2t(:,:) * e3t_n(:,:,jk) * tmask(:,:,jk) 119 113 END DO 120 !121 !$OMP DO schedule(static) private(jn)122 DO jn = 1, jptra123 trai(jn) = 0._wp ! initial content of all tracers124 END DO125 !$OMP END PARALLEL126 114 ! ! total volume of the ocean 127 115 areatot = glob_sum( cvol(:,:,:) ) 128 116 ! 117 trai(:) = 0._wp ! initial content of all tracers 129 118 DO jn = 1, jptra 130 119 trai(jn) = trai(jn) + glob_sum( trn(:,:,:,jn) * cvol(:,:,:) ) … … 231 220 USE trcdta ! initialisation from files 232 221 ! 233 INTEGER :: jn, jl , jk, jj, ji! dummy loop indices222 INTEGER :: jn, jl ! dummy loop indices 234 223 !!---------------------------------------------------------------------- 235 224 ! … … 265 254 ENDIF 266 255 ! 267 !$OMP PARALLEL DO schedule(static) private(jn,jk,jj,ji) 268 DO jn = 1, jptra 269 DO jk = 1, jpk 270 DO jj = 1, jpj 271 DO ji = 1, jpi 272 trb(ji,jj,jk,jn) = trn(ji,jj,jk,jn) 273 END DO 274 END DO 275 END DO 276 END DO 256 trb(:,:,:,:) = trn(:,:,:,:) 277 257 ! 278 258 ENDIF 279 259 280 !$OMP PARALLEL DO schedule(static) private(jn,jk,jj,ji) 281 DO jn = 1, jptra 282 DO jk = 1, jpk 283 DO jj = 1, jpj 284 DO ji = 1, jpi 285 tra(ji,jj,jk,jn) = 0._wp 286 END DO 287 END DO 288 END DO 289 END DO 260 tra(:,:,:,:) = 0._wp 290 261 ! ! Partial top/bottom cell: GRADh(trn) 291 262 END SUBROUTINE trc_ini_state
Note: See TracChangeset
for help on using the changeset viewer.