- Timestamp:
- 2016-11-01T15:26:15+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_r6711_SIMPLIF_6_aerobulk/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfbfr.F90
r6140 r7163 99 99 ! 100 100 IF( nn_timing == 1 ) CALL timing_start('zdf_bfr') 101 !102 IF( kt == nit000 .AND. lwp ) THEN103 WRITE(numout,*)104 WRITE(numout,*) 'zdf_bfr : Set bottom friction coefficient (non-linear case)'105 WRITE(numout,*) '~~~~~~~~'106 ENDIF107 101 ! 108 102 IF( nn_bfr == 2 ) THEN ! quadratic bottom friction only … … 259 253 IF(lwp) WRITE(numout,*) 260 254 IF(lwp) WRITE(numout,*) 'zdf_bfr_init : momentum bottom friction' 261 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~ ~'255 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 262 256 IF(lwp) WRITE(numout,*) ' Namelist nam_bfr : set bottom friction parameters' 263 257 ! … … 266 260 CASE( 0 ) 267 261 IF(lwp) WRITE(numout,*) ' free-slip ' 268 bfrua(:,:) = 0. e0269 bfrva(:,:) = 0. e0270 tfrua(:,:) = 0. e0271 tfrva(:,:) = 0. e0262 bfrua(:,:) = 0._wp 263 bfrva(:,:) = 0._wp 264 tfrua(:,:) = 0._wp 265 tfrva(:,:) = 0._wp 272 266 ! 273 267 CASE( 1 ) … … 321 315 IF(lwp) WRITE(numout,*) ' log formulation ln_bfr2d = ', ln_loglayer 322 316 IF(lwp) WRITE(numout,*) ' bottom roughness rn_bfrz0 [m] = ', rn_bfrz0 323 IF( rn_bfrz0 <=0.e0) THEN317 IF( rn_bfrz0 <= 0._wp ) THEN 324 318 WRITE(ctmp1,*) ' bottom roughness must be strictly positive' 325 319 CALL ctl_stop( ctmp1 ) … … 336 330 IF(lwp) WRITE(numout,*) ' log formulation ln_tfr2d = ', ln_loglayer 337 331 IF(lwp) WRITE(numout,*) ' top roughness rn_tfrz0 [m] = ', rn_tfrz0 338 IF( rn_tfrz0 <=0.e0) THEN332 IF( rn_tfrz0 <= 0._wp ) THEN 339 333 WRITE(ctmp1,*) ' top roughness must be strictly positive' 340 334 CALL ctl_stop( ctmp1 )
Note: See TracChangeset
for help on using the changeset viewer.