- Timestamp:
- 2011-11-18T22:41:44+01:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_qck.F90
r3116 r3159 26 26 USE diaptr ! poleward transport diagnostics 27 27 USE trc_oce ! share passive tracers/Ocean variables 28 USE wrk_nemo_2 ! Memory Allocation 29 USE timing ! Timing 28 30 29 31 IMPLICIT NONE … … 91 93 !!---------------------------------------------------------------------- 92 94 95 ! 96 IF( nn_timing == 1 ) CALL timing_start('tra_adv_qck') 97 ! 93 98 IF( kt == kit000 ) THEN 94 99 IF(lwp) WRITE(numout,*) … … 108 113 CALL tra_adv_cen2_k( kt, cdtype, pwn, ptn, pta, kjpt ) 109 114 ! 115 IF( nn_timing == 1 ) CALL timing_stop('tra_adv_qck') 116 ! 110 117 END SUBROUTINE tra_adv_qck 111 118 … … 116 123 !! 117 124 !!---------------------------------------------------------------------- 118 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released119 125 USE oce , ONLY: zwx => ua ! ua used as workspace 120 USE wrk_nemo, ONLY: zfu => wrk_3d_11 , zfc => wrk_3d_12, zfd => wrk_3d_13 ! 3D workspace121 126 ! 122 127 INTEGER , INTENT(in ) :: kt ! ocean time-step index … … 130 135 INTEGER :: ji, jj, jk, jn ! dummy loop indices 131 136 REAL(wp) :: ztra, zbtr, zdir, zdx, zdt, zmsk ! local scalars 137 REAL(wp), POINTER, DIMENSION(:,:,:) :: zfu, zfc, zfd 132 138 !---------------------------------------------------------------------- 133 139 ! 134 IF( wrk_in_use(3, 11,12,13) ) THEN 135 CALL ctl_stop('tra_adv_qck_i: requested workspace arrays unavailable') ; RETURN 136 ENDIF 140 CALL wrk_alloc( jpi, jpj, jpk, zfu, zfc, zfd ) 137 141 ! ! =========== 138 142 DO jn = 1, kjpt ! tracer loop … … 229 233 END DO 230 234 ! 231 IF( wrk_not_released(3, 11,12,13) ) CALL ctl_stop('tra_adv_qck_i: failed to release workspace arrays')235 CALL wrk_dealloc( jpi, jpj, jpk, zfu, zfc, zfd ) 232 236 ! 233 237 END SUBROUTINE tra_adv_qck_i … … 239 243 !! 240 244 !!---------------------------------------------------------------------- 241 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released242 245 USE oce , ONLY: zwy => ua ! ua used as workspace 243 USE wrk_nemo, ONLY: zfu => wrk_3d_11 , zfc => wrk_3d_12, zfd => wrk_3d_13 ! 3D workspace244 246 ! 245 247 INTEGER , INTENT(in ) :: kt ! ocean time-step index … … 253 255 INTEGER :: ji, jj, jk, jn ! dummy loop indices 254 256 REAL(wp) :: ztra, zbtr, zdir, zdx, zdt, zmsk ! local scalars 257 REAL(wp), POINTER, DIMENSION(:,:,:) :: zfu, zfc, zfd 255 258 !---------------------------------------------------------------------- 256 259 ! 257 IF(wrk_in_use(3, 11,12,13))THEN 258 CALL ctl_stop('tra_adv_qck_j: ERROR: requested workspace arrays unavailable') 259 RETURN 260 END IF 260 CALL wrk_alloc( jpi, jpj, jpk, zfu, zfc, zfd ) 261 ! 261 262 ! ! =========== 262 263 DO jn = 1, kjpt ! tracer loop … … 360 361 END DO 361 362 ! 362 IF( wrk_not_released(3, 11,12,13) ) CALL ctl_stop('tra_adv_qck_j: failed to release workspace arrays')363 CALL wrk_dealloc( jpi, jpj, jpk, zfu, zfc, zfd ) 363 364 ! 364 365 END SUBROUTINE tra_adv_qck_j … … 437 438 REAL(wp) :: zc, zcurv, zfho ! - - 438 439 !---------------------------------------------------------------------- 439 440 ! 441 IF( nn_timing == 1 ) CALL timing_start('quickest') 442 ! 440 443 DO jk = 1, jpkm1 441 444 DO jj = 1, jpj … … 468 471 END DO 469 472 ! 473 IF( nn_timing == 1 ) CALL timing_stop('quickest') 474 ! 470 475 END SUBROUTINE quickest 471 476
Note: See TracChangeset
for help on using the changeset viewer.