Changeset 3294 for trunk/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf_imp.F90
- Timestamp:
- 2012-01-28T17:44:18+01:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf_imp.F90
r2715 r3294 35 35 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 36 36 USE lib_mpp ! MPP library 37 USE wrk_nemo ! Memory Allocation 38 USE timing ! Timing 37 39 38 40 IMPLICIT NONE … … 55 57 CONTAINS 56 58 57 SUBROUTINE tra_zdf_imp( kt, cdtype, p2dt, ptb, pta, kjpt )59 SUBROUTINE tra_zdf_imp( kt, kit000, cdtype, p2dt, ptb, pta, kjpt ) 58 60 !!---------------------------------------------------------------------- 59 61 !! *** ROUTINE tra_zdf_imp *** … … 74 76 !! ** Action : - pta becomes the after tracer 75 77 !!--------------------------------------------------------------------- 76 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released77 78 USE oce , ONLY: zwd => ua , zws => va ! (ua,va) used as 3D workspace 78 USE wrk_nemo, ONLY: zwi => wrk_3d_6 , zwt => wrk_3d_7 ! 3D workspace79 79 ! 80 80 INTEGER , INTENT(in ) :: kt ! ocean time-step index 81 INTEGER , INTENT(in ) :: kit000 ! first time step index 81 82 CHARACTER(len=3) , INTENT(in ) :: cdtype ! =TRA or TRC (tracer indicator) 82 83 INTEGER , INTENT(in ) :: kjpt ! number of tracers … … 87 88 INTEGER :: ji, jj, jk, jn ! dummy loop indices 88 89 REAL(wp) :: zrhs, ze3tb, ze3tn, ze3ta ! local scalars 90 REAL(wp), POINTER, DIMENSION(:,:,:) :: zwi, zwt 89 91 !!--------------------------------------------------------------------- 90 91 IF( wrk_in_use(3, 6,7) ) THEN92 CALL ctl_stop('tra_zdf_imp : requested workspace arrays unavailable.') ; RETURN93 ENDIF94 95 IF( kt == nit000 ) THEN92 ! 93 IF( nn_timing == 1 ) CALL timing_start('tra_zdf_imp') 94 ! 95 CALL wrk_alloc( jpi, jpj, jpk, zwi, zwt ) 96 ! 97 IF( kt == kit000 ) THEN 96 98 IF(lwp)WRITE(numout,*) 97 99 IF(lwp)WRITE(numout,*) 'tra_zdf_imp : implicit vertical mixing on ', cdtype … … 228 230 ! ! ================= ! 229 231 ! 230 IF( wrk_not_released(3, 6,7) ) CALL ctl_stop('tra_zdf_imp: failed to release workspace arrays') 232 CALL wrk_dealloc( jpi, jpj, jpk, zwi, zwt ) 233 ! 234 IF( nn_timing == 1 ) CALL timing_stop('tra_zdf_imp') 231 235 ! 232 236 END SUBROUTINE tra_zdf_imp
Note: See TracChangeset
for help on using the changeset viewer.