- Timestamp:
- 2015-12-02T17:12:45+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5936_INGV1_WAVE/NEMOGCM/NEMO/OPA_SRC/step.F90
r5930 r5983 26 26 !! 3.6 ! 2012-07 (J. Simeon, G. Madec. C. Ethe) Online coarsening of outputs 27 27 !! 3.6 ! 2014-04 (F. Roquet, G. Madec) New equations of state 28 !! 3.6 ! 2014-10 (E. Clementi, P. Oddo) Add Qiao vertical mixing in case of waves 28 29 !! 3.7 ! 2014-10 (G. Madec) LDF simplication 29 30 !! - ! 2014-12 (G. Madec) remove KPP scheme … … 75 76 !! -8- Outputs and diagnostics 76 77 !!---------------------------------------------------------------------- 77 INTEGER :: j k! dummy loop indice78 INTEGER :: ji,jj,jk ! dummy loop indice 78 79 INTEGER :: indic ! error indicator if < 0 79 80 INTEGER :: kcall ! optional integer argument (dom_vvl_sf_nxt) … … 132 133 IF( lk_zdftke ) CALL zdf_tke( kstp ) ! TKE closure scheme for Kz 133 134 IF( lk_zdfgls ) CALL zdf_gls( kstp ) ! GLS closure scheme for Kz 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) 143 END DO 144 END DO 145 END DO 146 ENDIF 147 ! 134 148 IF( lk_zdfcst ) THEN ! Constant Kz (reset avt, avm[uv] to the background value) 135 149 avt (:,:,:) = rn_avt0 * wmask (:,:,:) … … 218 232 CALL dyn_vor ( kstp ) ! vorticity term including Coriolis 219 233 CALL dyn_ldf ( kstp ) ! lateral mixing 234 IF( ln_stcor ) CALL dyn_stcor ( kstp ) ! Stokes-Coriolis forcing 220 235 CALL dyn_hpg ( kstp ) ! horizontal gradient of Hydrostatic pressure 221 236 CALL dyn_spg ( kstp ) ! surface pressure gradient
Note: See TracChangeset
for help on using the changeset viewer.