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 3147 for branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/TRP/trcadv.F90 – NEMO

Ignore:
Timestamp:
2011-11-17T17:28:03+01:00 (13 years ago)
Author:
cetlod
Message:

branch dev_NEMO_MERGE_2011:New dynamical allocation for TOP

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/TOP_SRC/TRP/trcadv.F90

    r3124 r3147  
    2626   USE ldftra_oce      ! lateral diffusion coefficient on tracers 
    2727   USE prtctl_trc      ! Print control 
     28   USE wrk_nemo_2      ! Memory allocation 
    2829 
    2930   IMPLICIT NONE 
     
    6768      !! ** Method  : - Update the tracer with the advection term following nadv 
    6869      !!---------------------------------------------------------------------- 
    69       USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released 
    70       USE wrk_nemo, ONLY: zun => wrk_3d_4, zvn => wrk_3d_5, zwn => wrk_3d_6   ! effective velocity 
    7170      !! 
    7271      INTEGER, INTENT(in) ::   kt   ! ocean time-step index 
     
    7473      INTEGER ::   jk  
    7574      CHARACTER (len=22) ::   charout 
    76       !!---------------------------------------------------------------------- 
    77       ! 
    78       IF( wrk_in_use(3, 4,5,6) ) THEN 
    79          CALL ctl_stop('trc_adv : requested workspace arrays unavailable')   ;   RETURN 
    80       ENDIF 
     75      REAL(wp), POINTER, DIMENSION(:,:,:) :: zun, zvn, zwn  ! effective velocity 
     76      !!---------------------------------------------------------------------- 
     77      ! 
     78      CALL wrk_alloc( jpi, jpj, jpk, zun, zvn, zwn ) 
     79      ! 
    8180 
    8281      IF( kt == nittrc000 )   CALL trc_adv_ctl          ! initialisation & control of options 
     
    141140      END IF 
    142141      ! 
    143       IF( wrk_not_released(3, 4,5,6) ) CALL ctl_stop('trc_adv : failed to release workspace arrays.') 
     142      CALL wrk_dealloc( jpi, jpj, jpk, zun, zvn, zwn ) 
    144143      ! 
    145144   END SUBROUTINE trc_adv 
Note: See TracChangeset for help on using the changeset viewer.