Changeset 12495 for NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/src/OCE/TRD/trdglo.F90
- Timestamp:
- 2020-03-02T09:10:34+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/src/OCE/TRD/trdglo.F90
r12377 r12495 76 76 INTEGER :: ji, jj, jk ! dummy loop indices 77 77 INTEGER :: ikbu, ikbv ! local integers 78 REAL(wp):: zvm, zvt, zvs, z1_2r au0 ! local scalars78 REAL(wp):: zvm, zvt, zvs, z1_2rho0 ! local scalars 79 79 REAL(wp), DIMENSION(jpi,jpj) :: ztswu, ztswv, z2dx, z2dy ! 2D workspace 80 80 !!---------------------------------------------------------------------- … … 125 125 ! 126 126 IF( ktrd == jpdyn_zdf ) THEN ! zdf trend: compute separately the surface forcing trend 127 z1_2r au0 = 0.5_wp / rau0127 z1_2rho0 = 0.5_wp / rho0 128 128 DO_2D_10_10 129 129 zvt = ( utau_b(ji,jj) + utau(ji,jj) ) * tmask_i(ji+1,jj) * tmask_i(ji,jj) * umask(ji,jj,jk) & 130 & * z1_2r au0 * e1e2u(ji,jj)130 & * z1_2rho0 * e1e2u(ji,jj) 131 131 zvs = ( vtau_b(ji,jj) + vtau(ji,jj) ) * tmask_i(ji,jj+1) * tmask_i(ji,jj) * vmask(ji,jj,jk) & 132 & * z1_2r au0 * e1e2v(ji,jj)132 & * z1_2rho0 * e1e2v(ji,jj) 133 133 umo(jpdyn_tau) = umo(jpdyn_tau) + zvt 134 134 vmo(jpdyn_tau) = vmo(jpdyn_tau) + zvs … … 141 141 ! ! 142 142 ! IF( ln_drgimp ) THEN ! implicit drag case: compute separately the bottom friction 143 ! z1_2r au0 = 0.5_wp / rau0143 ! z1_2rho0 = 0.5_wp / rho0 144 144 ! DO jj = 1, jpjm1 145 145 ! DO ji = 1, jpim1 … … 203 203 CALL eos( ts(:,:,:,:,Kmm), rhd, rhop ) ! now potential density 204 204 205 zcof = 0.5_wp / r au0 ! Density flux at w-point205 zcof = 0.5_wp / rho0 ! Density flux at w-point 206 206 zkz(:,:,1) = 0._wp 207 207 DO jk = 2, jpk … … 209 209 END DO 210 210 211 zcof = 0.5_wp / r au0 ! Density flux at u and v-points211 zcof = 0.5_wp / rho0 ! Density flux at u and v-points 212 212 DO_3D_10_10( 1, jpkm1 ) 213 213 zkx(ji,jj,jk) = zcof * e2u(ji,jj) * e3u(ji,jj,jk,Kmm) * uu(ji,jj,jk,Kmm) * ( rhop(ji,jj,jk) + rhop(ji+1,jj,jk) ) … … 347 347 9546 FORMAT(' 0 < horizontal diffusion : ', e20.13) 348 348 9547 FORMAT(' 0 < vertical diffusion : ', e20.13) 349 9548 FORMAT(' pressure gradient u2 = - 1/r au0 u.dz(rhop) : ', e20.13, ' u.dz(rhop) =', e20.13)349 9548 FORMAT(' pressure gradient u2 = - 1/rho0 u.dz(rhop) : ', e20.13, ' u.dz(rhop) =', e20.13) 350 350 ! 351 351 ! Save potential to kinetic energy conversion for next time step
Note: See TracChangeset
for help on using the changeset viewer.