- Timestamp:
- 2016-11-21T11:23:19+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/r5936_INGV1_WAVE-coupling/NEMOGCM/NEMO/OPA_SRC/step.F90
r7168 r7279 134 134 IF( lk_zdfgls ) CALL zdf_gls( kstp ) ! GLS closure scheme for Kz 135 135 IF( ln_zdfqiao ) THEN 136 CALL zdf_qiao(kstp ) ! Qiao vertical mixing 137 DO jk = 1, jpkm1 138 DO jj = 1, jpj 139 DO ji = 1, jpi 140 avmu(ji,jj,jk) = (avmu(ji,jj,jk) + QBvu(ji,jj,jk)) * umask(ji,jj,jk) 141 avmv(ji,jj,jk) = (avmv(ji,jj,jk) + QBvv(ji,jj,jk)) * vmask(ji,jj,jk) 142 avt( ji,jj,jk) = (avt( ji,jj,jk) + QBv(ji,jj,jk)) * tmask(ji,jj,jk) 136 !Activated Qiao enhanced turbulence but neither ln_wave or ln_sdw are activated 137 IF ( .NOT.( ln_wave .AND. ln_sdw ) ) THEN 138 CALL ctl_stop ( 'Ask for wave Qiao enhanced turbulence but ln_wave and ln_sdw have to be activated') 139 ELSE 140 CALL zdf_qiao(kstp ) ! Qiao vertical mixing 141 DO jk = 1, jpkm1 142 DO jj = 1, jpj 143 DO ji = 1, jpi 144 avmu(ji,jj,jk) = (avmu(ji,jj,jk) + QBvu(ji,jj,jk)) * umask(ji,jj,jk) 145 avmv(ji,jj,jk) = (avmv(ji,jj,jk) + QBvv(ji,jj,jk)) * vmask(ji,jj,jk) 146 avt( ji,jj,jk) = (avt( ji,jj,jk) + QBv(ji,jj,jk)) * tmask(ji,jj,jk) 147 END DO 148 END DO 143 149 END DO 144 END DO 145 END DO 150 ENDIF 146 151 ENDIF 147 152 ! … … 231 236 CALL dyn_adv ( kstp ) ! advection (vector or flux form) 232 237 CALL dyn_vor ( kstp ) ! vorticity term including Coriolis 233 IF( ln_stcor ) CALL dyn_stcor ( kstp ) ! Stokes-Coriolis forcing 238 IF( ln_wave .AND. ln_sdw .AND. ln_stcor ) & 239 & CALL dyn_stcor ( kstp ) ! Stokes-Coriolis forcing 234 240 CALL dyn_ldf ( kstp ) ! lateral mixing 235 241 CALL dyn_hpg ( kstp ) ! horizontal gradient of Hydrostatic pressure
Note: See TracChangeset
for help on using the changeset viewer.