- Timestamp:
- 2017-12-19T09:26:25+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsink.F90
r7753 r9124 69 69 !!--------------------------------------------------------------------- 70 70 ! 71 IF( nn_timing == 1 )CALL timing_start('p4z_sink')71 IF( ln_timing ) CALL timing_start('p4z_sink') 72 72 73 73 … … 261 261 ENDIF 262 262 ! 263 IF( nn_timing == 1 )CALL timing_stop('p4z_sink')263 IF( ln_timing ) CALL timing_stop('p4z_sink') 264 264 ! 265 265 END SUBROUTINE p4z_sink 266 266 267 267 268 SUBROUTINE p4z_sink_init 268 269 !!---------------------------------------------------------------------- … … 270 271 !!---------------------------------------------------------------------- 271 272 INTEGER :: jk 272 273 !!---------------------------------------------------------------------- 274 ! 273 275 ik100 = 10 ! last level where depth less than 100 m 274 276 DO jk = jpkm1, 1, -1 … … 282 284 ! 283 285 END SUBROUTINE p4z_sink_init 286 284 287 285 288 SUBROUTINE p4z_sink2( pwsink, psinkflx, jp_tra, kiter ) … … 294 297 !! transport term, i.e. div(u*tra). 295 298 !!--------------------------------------------------------------------- 296 !297 299 INTEGER , INTENT(in ) :: jp_tra ! tracer index index 298 300 INTEGER , INTENT(in ) :: kiter ! number of iterations for time-splitting 299 301 REAL(wp), INTENT(in ), DIMENSION(jpi,jpj,jpk) :: pwsink ! sinking speed 300 302 REAL(wp), INTENT(inout), DIMENSION(jpi,jpj,jpk) :: psinkflx ! sinking fluxe 301 ! !303 ! 302 304 INTEGER :: ji, jj, jk, jn 303 305 REAL(wp) :: zigma,zew,zign, zflx, zstep … … 305 307 !!--------------------------------------------------------------------- 306 308 ! 307 IF( nn_timing == 1 )CALL timing_start('p4z_sink2')309 IF( ln_timing ) CALL timing_start('p4z_sink2') 308 310 ! 309 311 ! Allocate temporary workspace … … 391 393 CALL wrk_dealloc( jpi, jpj, jpk, ztraz, zakz, zwsink2, ztrb ) 392 394 ! 393 IF( nn_timing == 1) CALL timing_stop('p4z_sink2')395 IF( ln_timing ) CALL timing_stop('p4z_sink2') 394 396 ! 395 397 END SUBROUTINE p4z_sink2 … … 401 403 !!---------------------------------------------------------------------- 402 404 INTEGER :: ierr(3) 403 405 !!---------------------------------------------------------------------- 406 ! 404 407 ierr(:) = 0 405 408 !
Note: See TracChangeset
for help on using the changeset viewer.