New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 3294 for trunk/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf_imp.F90 – NEMO

Ignore:
Timestamp:
2012-01-28T17:44:18+01:00 (12 years ago)
Author:
rblod
Message:

Merge of 3.4beta into the trunk

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf_imp.F90

    r2715 r3294  
    3535   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
    3636   USE lib_mpp         ! MPP library 
     37   USE wrk_nemo        ! Memory Allocation 
     38   USE timing          ! Timing 
    3739 
    3840   IMPLICIT NONE 
     
    5557CONTAINS 
    5658  
    57    SUBROUTINE tra_zdf_imp( kt, cdtype, p2dt, ptb, pta, kjpt )  
     59   SUBROUTINE tra_zdf_imp( kt, kit000, cdtype, p2dt, ptb, pta, kjpt )  
    5860      !!---------------------------------------------------------------------- 
    5961      !!                  ***  ROUTINE tra_zdf_imp  *** 
     
    7476      !! ** Action  : - pta  becomes the after tracer 
    7577      !!--------------------------------------------------------------------- 
    76       USE wrk_nemo, ONLY:   wrk_in_use, wrk_not_released 
    7778      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 workspace  
    7979      ! 
    8080      INTEGER                              , INTENT(in   ) ::   kt       ! ocean time-step index 
     81      INTEGER                              , INTENT(in   ) ::   kit000          ! first time step index 
    8182      CHARACTER(len=3)                     , INTENT(in   ) ::   cdtype   ! =TRA or TRC (tracer indicator) 
    8283      INTEGER                              , INTENT(in   ) ::   kjpt     ! number of tracers 
     
    8788      INTEGER  ::  ji, jj, jk, jn   ! dummy loop indices 
    8889      REAL(wp) ::  zrhs, ze3tb, ze3tn, ze3ta   ! local scalars 
     90      REAL(wp), POINTER, DIMENSION(:,:,:) ::  zwi, zwt 
    8991      !!--------------------------------------------------------------------- 
    90  
    91       IF( wrk_in_use(3, 6,7) ) THEN 
    92          CALL ctl_stop('tra_zdf_imp : requested workspace arrays unavailable.')   ;   RETURN 
    93       ENDIF 
    94  
    95       IF( kt == nit000 )  THEN 
     92      ! 
     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 
    9698         IF(lwp)WRITE(numout,*) 
    9799         IF(lwp)WRITE(numout,*) 'tra_zdf_imp : implicit vertical mixing on ', cdtype 
     
    228230      !                                               ! ================= ! 
    229231      ! 
    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') 
    231235      ! 
    232236   END SUBROUTINE tra_zdf_imp 
Note: See TracChangeset for help on using the changeset viewer.