- Timestamp:
- 2011-12-21T10:11:11+01:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfgls.F90
r3229 r3231 133 133 !! coefficients using the GLS turbulent closure scheme. 134 134 !!---------------------------------------------------------------------- 135 USE oce , ONLY: tsa ! use tsa as workspace136 USE oce , ONLY: z_elem_a => ua ! use ua as workspace137 USE oce , ONLY: z_elem_b => va ! use va as workspace138 !139 135 INTEGER, INTENT(in) :: kt ! ocean time step 140 136 INTEGER :: ji, jj, jk, ibot, ibotm1, dir ! dummy loop arguments … … 152 148 REAL(wp), POINTER, DIMENSION(:,:,:) :: eps ! dissipation rate 153 149 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, psi150 REAL(wp), POINTER, DIMENSION(:,:,:) :: z_elem_a, z_elem_b, z_elem_c, psi 155 151 !!-------------------------------------------------------------------- 156 152 ! … … 158 154 ! 159 155 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 ) 164 157 165 158 ! Preliminary computing … … 890 883 ! 891 884 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 ) 893 886 ! 894 887 IF( nn_timing == 1 ) CALL timing_stop('zdf_gls')
Note: See TracChangeset
for help on using the changeset viewer.