Changeset 10874 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/ZDF/zdfphy.F90
- Timestamp:
- 2019-04-15T15:57:37+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/ZDF/zdfphy.F90
r10829 r10874 218 218 219 219 220 SUBROUTINE zdf_phy( kt , ktlev1, ktlev2, kt2lev1, kt2lev2)220 SUBROUTINE zdf_phy( kt ) 221 221 !!---------------------------------------------------------------------- 222 222 !! *** ROUTINE zdf_phy *** … … 231 231 !!---------------------------------------------------------------------- 232 232 INTEGER, INTENT(in) :: kt ! ocean time-step index 233 INTEGER, INTENT(in) :: ktlev1, ktlev2 ! time level indices for 3-time-level source terms234 INTEGER, INTENT(in) :: kt2lev1, kt2lev2 ! time level indices for 2-time-level source terms235 233 ! 236 234 INTEGER :: ji, jj, jk ! dummy loop indice … … 256 254 ! 257 255 IF( l_zdfsh2 ) & !* shear production at w-points (energy conserving form) 258 CALL zdf_sh2( u u(:,:,:,ktlev1), vv(:,:,:,ktlev1), uu(:,:,:,ktlev2), vv(:,:,:,ktlev2), avm_k, & ! <<== in256 CALL zdf_sh2( ub, vb, un, vn, avm_k, & ! <<== in 259 257 & zsh2 ) ! ==>> out : shear production 260 258 ! 261 259 SELECT CASE ( nzdf_phy ) !* Vertical eddy viscosity and diffusivity coefficients at w-points 262 CASE( np_RIC ) ; CALL zdf_ric( kt, gdept (:,:,:,kt2lev2), zsh2, avm_k, avt_k ) ! Richardson number dependent Kz263 CASE( np_TKE ) ; CALL zdf_tke( kt , ktlev1, ktlev2, kt2lev1, kt2lev2, zsh2, avm_k, avt_k ) ! TKE closure scheme for Kz260 CASE( np_RIC ) ; CALL zdf_ric( kt, gdept_n, zsh2, avm_k, avt_k ) ! Richardson number dependent Kz 261 CASE( np_TKE ) ; CALL zdf_tke( kt , zsh2, avm_k, avt_k ) ! TKE closure scheme for Kz 264 262 CASE( np_GLS ) ; CALL zdf_gls( kt , zsh2, avm_k, avt_k ) ! GLS closure scheme for Kz 265 263 CASE( np_OSM ) ; CALL zdf_osm( kt , avm_k, avt_k ) ! OSMOSIS closure scheme for Kz … … 320 318 IF( ln_zdfgls ) CALL gls_rst( kt, 'WRITE' ) 321 319 IF( ln_zdfric ) CALL ric_rst( kt, 'WRITE' ) 322 ! NB. OSMOSIS restart (osm_rst) will be called in step.F90 after w whas been updated320 ! NB. OSMOSIS restart (osm_rst) will be called in step.F90 after wn has been updated 323 321 ENDIF 324 322 !
Note: See TracChangeset
for help on using the changeset viewer.