Changeset 925 for codes/icosagcm/devel/src
- Timestamp:
- 06/21/19 00:11:29 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/src/dissip/dissip_gcm.f90
r924 r925 376 376 END SELECT 377 377 END SELECT 378 378 379 IF(is_master) PRINT *, 'vertical profile of horizontal dissipation : ', zvert(:) 380 379 381 DO l=1,llm 380 382 tau_graddiv(l) = tau_graddiv(l)/zvert(l) … … 407 409 itau_dissip=MAX(1,itau_dissip) 408 410 dtdissip=itau_dissip*dt 411 409 412 IF (is_master) THEN 410 413 PRINT *,"init_dissip: rayleigh_tau",rayleigh_tau, "mintau ",mintau … … 412 415 ENDIF 413 416 417 IF (dtdissip>2.*mintau) THEN 418 IF(is_master) PRINT *, 'The CFL condition for dissipation dtdissip<2*mintau is violated : dtdissip, mintau ', dtdissip, mintau 419 STOP 420 END IF 421 414 422 END SUBROUTINE dissip_timescale 415 423
Note: See TracChangeset
for help on using the changeset viewer.