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 3231 for branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfgls.F90 – NEMO

Ignore:
Timestamp:
2011-12-21T10:11:11+01:00 (12 years ago)
Author:
smasson
Message:

dev_NEMO_MERGE_2011: supress TARGET attribute for tsa and use work arrays

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfgls.F90

    r3229 r3231  
    133133      !!              coefficients using the GLS turbulent closure scheme. 
    134134      !!---------------------------------------------------------------------- 
    135       USE oce     , ONLY: tsa                 ! use tsa as workspace 
    136       USE oce     , ONLY: z_elem_a  =>   ua   ! use ua as workspace 
    137       USE oce     , ONLY: z_elem_b  =>   va   ! use va as workspace 
    138       ! 
    139135      INTEGER, INTENT(in) ::   kt ! ocean time step 
    140136      INTEGER  ::   ji, jj, jk, ibot, ibotm1, dir  ! dummy loop arguments 
     
    152148      REAL(wp), POINTER, DIMENSION(:,:,:) ::   eps         ! dissipation rate 
    153149      REAL(wp), POINTER, DIMENSION(:,:,:) ::   zwall_psi   ! Wall function use in the wb case (ln_sigpsi.AND.ln_crban=T) 
    154       REAL(wp), POINTER, DIMENSION(:,:,:) ::   z_elem_c, psi 
     150      REAL(wp), POINTER, DIMENSION(:,:,:) ::   z_elem_a, z_elem_b, z_elem_c, psi 
    155151      !!-------------------------------------------------------------------- 
    156152      ! 
     
    158154      ! 
    159155      CALL wrk_alloc( jpi,jpj, zdep, zflxs, zhsro ) 
    160       CALL wrk_alloc( jpi,jpj,jpk, eb, mxlb, shear, eps, zwall_psi ) 
    161       ! 
    162       z_elem_c  => tsa(:,:,:,1) 
    163       psi       => tsa(:,:,:,2) 
     156      CALL wrk_alloc( jpi,jpj,jpk, eb, mxlb, shear, eps, zwall_psi, z_elem_a, z_elem_b, z_elem_c, psi ) 
    164157 
    165158      ! Preliminary computing 
     
    890883      ! 
    891884      CALL wrk_dealloc( jpi,jpj, zdep, zflxs, zhsro ) 
    892       CALL wrk_dealloc( jpi,jpj,jpk, eb, mxlb, shear, eps, zwall_psi ) 
     885      CALL wrk_dealloc( jpi,jpj,jpk, eb, mxlb, shear, eps, zwall_psi, z_elem_a, z_elem_b, z_elem_c, psi ) 
    893886      ! 
    894887      IF( nn_timing == 1 )  CALL timing_stop('zdf_gls') 
Note: See TracChangeset for help on using the changeset viewer.