- Timestamp:
- 2011-12-21T10:11:11+01:00 (13 years ago)
- Location:
- branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/OPA_SRC/ZDF
- Files:
-
- 2 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') -
branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90
r3229 r3231 192 192 !! (= Kz dz[Ub] * dz[Un] ) 193 193 !! --------------------------------------------------------------------- 194 USE oce , ONLY: zdiag => ua ! (ua,va) used as workspace195 USE oce , ONLY: tsa ! (tsa) used as workspace196 !!197 194 INTEGER :: ji, jj, jk ! dummy loop arguments 198 195 !!bfr INTEGER :: ikbu, ikbv, ikbum1, ikbvm1 ! temporary scalar … … 207 204 REAL(wp) :: zzd_up, zzd_lw ! - - 208 205 !!bfr REAL(wp) :: zebot ! - - 209 REAL(wp), POINTER, DIMENSION(:,:,:) :: zd_up, zd_lw210 206 INTEGER , POINTER, DIMENSION(:,: ) :: imlc 211 207 REAL(wp), POINTER, DIMENSION(:,: ) :: zhlc 212 REAL(wp), POINTER, DIMENSION(:,:,:) :: zpelc 208 REAL(wp), POINTER, DIMENSION(:,:,:) :: zpelc, zdiag, zd_up, zd_lw 213 209 !!-------------------------------------------------------------------- 214 210 ! … … 217 213 CALL wrk_alloc( jpi,jpj, imlc ) ! integer 218 214 CALL wrk_alloc( jpi,jpj, zhlc ) 219 CALL wrk_alloc( jpi,jpj,jpk, zpelc ) 220 ! 221 zd_up => tsa(:,:,:,1) 222 zd_lw => tsa(:,:,:,2) 223 215 CALL wrk_alloc( jpi,jpj,jpk, zpelc, zdiag, zd_up, zd_lw ) 216 ! 224 217 zbbrau = rn_ebb / rau0 ! Local constant initialisation 225 218 zfact1 = -.5_wp * rdt … … 437 430 CALL wrk_dealloc( jpi,jpj, imlc ) ! integer 438 431 CALL wrk_dealloc( jpi,jpj, zhlc ) 439 CALL wrk_dealloc( jpi,jpj,jpk, zpelc )432 CALL wrk_dealloc( jpi,jpj,jpk, zpelc, zdiag, zd_up, zd_lw ) 440 433 ! 441 434 IF( nn_timing == 1 ) CALL timing_stop('tke_tke') … … 479 472 !! - avmu, avmv : now vertical eddy viscosity at uw- and vw-points 480 473 !!---------------------------------------------------------------------- 481 USE oce, ONLY: zmpdl => ua ! ua used as workspace482 USE oce, ONLY: tsa ! use tsa as workspace483 !!484 474 INTEGER :: ji, jj, jk ! dummy loop indices 485 475 REAL(wp) :: zrn2, zraug, zcoef, zav ! local scalars 486 476 REAL(wp) :: zdku, zpdlr, zri, zsqen ! - - 487 477 REAL(wp) :: zdkv, zemxl, zemlm, zemlp ! - - 488 REAL(wp), POINTER, DIMENSION(:,:,:) :: zm xlm, zmxld478 REAL(wp), POINTER, DIMENSION(:,:,:) :: zmpdl, zmxlm, zmxld 489 479 !!-------------------------------------------------------------------- 490 480 ! 491 481 IF( nn_timing == 1 ) CALL timing_start('tke_avn') 492 ! 493 zmxlm => tsa(:,:,:,1) 494 zmxld => tsa(:,:,:,2) 482 483 CALL wrk_alloc( jpi,jpj,jpk, zmpdl, zmxlm, zmxld ) 495 484 496 485 ! !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< … … 664 653 & tab3d_2=avmv, clinfo2= ' v: ', mask2=vmask, ovlap=1, kdim=jpk ) 665 654 ENDIF 655 ! 656 CALL wrk_dealloc( jpi,jpj,jpk, zmpdl, zmxlm, zmxld ) 666 657 ! 667 658 IF( nn_timing == 1 ) CALL timing_stop('tke_avn')
Note: See TracChangeset
for help on using the changeset viewer.