Changeset 706
- Timestamp:
- 07/17/18 17:24:56 (6 years ago)
- Location:
- codes/icosagcm/trunk/src
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/dissip/dissip_gcm.f90
r648 r706 453 453 454 454 IF(mintau>0) THEN 455 itau_dissip=INT(mintau/dt) 456 dtdissip=itau_dissip*dt 455 IF (itau_dissip==0) THEN 456 IF (is_master) PRINT *,"init_dissip: Automatic computation of itau_dissip..." 457 itau_dissip=INT(mintau/dt) 458 ENDIF 457 459 ELSE 458 IF (is_master) PRINT *," No dissipation time set, setting itau_dissip to 1000000000"460 IF (is_master) PRINT *,"init_dissip: No dissipation time set, setting itau_dissip to 1000000000" 459 461 itau_dissip=100000000 460 462 END IF 461 463 itau_dissip=MAX(1,itau_dissip) 462 IF (is_master) PRINT *,"rayleigh_tau",rayleigh_tau, "mintau ",mintau, & 463 "itau_dissip",itau_dissip," dtdissip ",dtdissip 464 dtdissip=itau_dissip*dt 465 IF (is_master) THEN 466 PRINT *,"init_dissip: rayleigh_tau",rayleigh_tau, "mintau ",mintau 467 PRINT *,"init_dissip: itau_dissip",itau_dissip," dtdissip ",dtdissip 468 ENDIF 464 469 465 470 END SUBROUTINE init_dissip … … 557 562 CALL trace_end("dissip") 558 563 564 CALL write_dissip_tendencies 559 565 !$OMP BARRIER 560 566 … … 588 594 END SUBROUTINE relax 589 595 596 SUBROUTINE write_dissip_tendencies 597 USE observable_mod, ONLY : f_buf_ulon, f_buf_ulat 598 USE wind_mod 599 USE output_field_mod 600 601 CALL transfert_request(f_due_diss1,req_e1_vect) 602 CALL un2ulonlat(f_due_diss1, f_buf_ulon, f_buf_ulat, (1./(tau_graddiv(1)))) 603 CALL output_field("dulon_diss1",f_buf_ulon) 604 CALL output_field("dulat_diss1",f_buf_ulat) 605 ! 606 CALL transfert_request(f_due_diss2,req_e1_vect) 607 CALL un2ulonlat(f_due_diss2, f_buf_ulon, f_buf_ulat, (1./(tau_graddiv(1)))) 608 CALL output_field("dulon_diss2",f_buf_ulon) 609 CALL output_field("dulat_diss2",f_buf_ulat) 610 END SUBROUTINE write_dissip_tendencies 611 590 612 END SUBROUTINE dissip 591 613 -
codes/icosagcm/trunk/src/time/time.f90
r548 r706 75 75 CALL getin('itau_adv',itau_adv) 76 76 77 itau_dissip= 177 itau_dissip=0 ! set to zero which implies itau_dissip will be automatically computed (see init_dissip) 78 78 CALL getin('itau_dissip',itau_dissip) 79 79
Note: See TracChangeset
for help on using the changeset viewer.