Changeset 6140 for trunk/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfric.F90
- Timestamp:
- 2015-12-21T12:35:23+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfric.F90
r5836 r6140 55 55 56 56 !! * Substitutions 57 # include " domzgr_substitute.h90"57 # include "vectopt_loop_substitute.h90" 58 58 !!---------------------------------------------------------------------- 59 59 !! NEMO/OPA 4.0 , NEMO Consortium (2011) … … 133 133 ! ----------------- 134 134 DO jj = 2, jpjm1 135 DO ji = 2,jpim1136 zcoef = 0.5 / fse3w(ji,jj,jk)135 DO ji = fs_2, fs_jpim1 136 zcoef = 0.5 / e3w_n(ji,jj,jk) 137 137 ! ! shear of horizontal velocity 138 138 zdku = zcoef * ( ub(ji-1,jj,jk-1) + ub(ji,jj,jk-1) & … … 151 151 z05alp = 0.5_wp * rn_alp 152 152 DO jj = 1, jpjm1 ! Eddy viscosity coefficients (avm) 153 DO ji = 1, jpim1153 DO ji = 1, fs_jpim1 154 154 avmu(ji,jj,jk) = umask(ji,jj,jk) * rn_avmri / ( 1. + z05alp*( zwx(ji+1,jj)+zwx(ji,jj) ) )**nn_ric 155 155 avmv(ji,jj,jk) = vmask(ji,jj,jk) * rn_avmri / ( 1. + z05alp*( zwx(ji,jj+1)+zwx(ji,jj) ) )**nn_ric … … 157 157 END DO 158 158 DO jj = 2, jpjm1 ! Eddy diffusivity coefficients (avt) 159 DO ji = 2,jpim1159 DO ji = fs_2, fs_jpim1 160 160 avt(ji,jj,jk) = tmric(ji,jj,jk) / ( 1._wp + rn_alp * zwx(ji,jj) ) & 161 161 & * ( avmu(ji,jj,jk) + avmu(ji-1,jj,jk) & … … 176 176 ! ------------------------------------------------------- 177 177 zflageos = ( 0.5 + SIGN( 0.5, nn_eos - 1. ) ) * rau0 178 DO jj = 1, jpj179 DO ji = 1, jpi178 DO jj = 2, jpjm1 179 DO ji = fs_2, fs_jpim1 180 180 zrhos = rhop(ji,jj,1) + zflageos * ( 1. - tmask(ji,jj,1) ) 181 181 zustar = SQRT( taum(ji,jj) / ( zrhos + rsmall ) ) … … 189 189 ! are always equal to the namelist values rn_wtmix/rn_wvmix 190 190 ! ------------------------------------------------------- 191 DO jj = 1, jpj192 DO ji = 1, jpi191 DO jj = 2, jpjm1 192 DO ji = fs_2, fs_jpim1 193 193 avmv(ji,jj,1) = MAX( avmv(ji,jj,1), rn_wvmix ) 194 194 avmu(ji,jj,1) = MAX( avmu(ji,jj,1), rn_wvmix ) … … 200 200 ! ------------------------------------------------------- 201 201 DO jk = 2, jpkm1 202 DO jj = 1, jpj203 DO ji = 1, jpi202 DO jj = 2, jpjm1 203 DO ji = fs_2, fs_jpim1 204 204 IF( fsdept(ji,jj,jk) < ekm_dep(ji,jj) ) THEN 205 205 avmv(ji,jj,jk) = MAX( avmv(ji,jj,jk), rn_wvmix ) … … 212 212 213 213 DO jk = 1, jpkm1 214 DO jj = 1, jpj215 DO ji = 1, jpi214 DO jj = 2, jpjm1 215 DO ji = fs_2, fs_jpim1 216 216 avmv(ji,jj,jk) = avmv(ji,jj,jk) * vmask(ji,jj,jk) 217 217 avmu(ji,jj,jk) = avmu(ji,jj,jk) * umask(ji,jj,jk)
Note: See TracChangeset
for help on using the changeset viewer.