- Timestamp:
- 2017-09-27T16:29:24+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfgls.F90
r8215 r8568 159 159 !!-------------------------------------------------------------------- 160 160 ! 161 IF( nn_timing == 1) CALL timing_start('zdf_gls')161 IF( ln_timing ) CALL timing_start('zdf_gls') 162 162 ! 163 163 ! Preliminary computing … … 822 822 ENDIF 823 823 ! 824 IF( nn_timing == 1) CALL timing_stop('zdf_gls')824 IF( ln_timing ) CALL timing_stop('zdf_gls') 825 825 ! 826 826 END SUBROUTINE zdf_gls … … 852 852 !!---------------------------------------------------------- 853 853 ! 854 IF( nn_timing == 1 )CALL timing_start('zdf_gls_init')854 IF( ln_timing ) CALL timing_start('zdf_gls_init') 855 855 ! 856 856 REWIND( numnam_ref ) ! Namelist namzdf_gls in reference namelist : Vertical eddy diffivity and viscosity using gls turbulent closure scheme … … 1077 1077 rl_sf = vkarmn 1078 1078 ELSE 1079 rl_sf = rc0 * SQRT(rc0/rcm_sf) * SQRT( ( (1._wp + 4._wp*rmm + 8._wp*rmm**2_wp)*rsc_tke & 1080 & + 12._wp * rsc_psi0*rpsi2 - (1._wp + 4._wp*rmm) & 1081 & *SQRT(rsc_tke*(rsc_tke & 1082 & + 24._wp*rsc_psi0*rpsi2)) ) & 1083 & /(12._wp*rnn**2.) & 1084 & ) 1079 rl_sf = rc0 * SQRT(rc0/rcm_sf) * SQRT( ( (1._wp + 4._wp*rmm + 8._wp*rmm**2_wp) * rsc_tke & 1080 & + 12._wp*rsc_psi0*rpsi2 - (1._wp + 4._wp*rmm) & 1081 & *SQRT(rsc_tke*(rsc_tke & 1082 & + 24._wp*rsc_psi0*rpsi2)) ) & 1083 & /(12._wp*rnn**2.) ) 1085 1084 ENDIF 1086 1085 … … 1130 1129 CALL gls_rst( nit000, 'READ' ) ! (en, avt_k, avm_k, hmxl_n) 1131 1130 ! 1132 IF( nn_timing == 1 )CALL timing_stop('zdf_gls_init')1131 IF( ln_timing ) CALL timing_stop('zdf_gls_init') 1133 1132 ! 1134 1133 END SUBROUTINE zdf_gls_init
Note: See TracChangeset
for help on using the changeset viewer.