Skip to content

Commit

Permalink
Bugfix for tau
Browse files Browse the repository at this point in the history
With damping option ‘JFV’, the stratospheric damping rate now actually
gets set to the JFV equilibrium. Before this bug fix, the JFV tau_stat
has computed, but not used in setting tdamp.
  • Loading branch information
mjucker committed Sep 15, 2015
1 parent 3db4c4a commit 38816fd
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions fms_riga_hs_jucker/src/atmos_param/hs_forcing/hs_forcing.f90
Original file line number Diff line number Diff line change
Expand Up @@ -798,7 +798,7 @@ subroutine newtonian_damping ( Time, lat, ps, p_full, p_half, t, tdt, teq, tau,
teq = teq_strat
end where
end if
!! same for tamping rate
!! same for damping rate
if ( trim(equilibrium_tau_option) == 'JFV' ) then
tau_strat = 0.
do k=1,size(t,3)
Expand Down Expand Up @@ -838,7 +838,7 @@ subroutine newtonian_damping ( Time, lat, ps, p_full, p_half, t, tdt, teq, tau,
tau_strat = 1./tau_strat
endif
if ( trim(equilibrium_tau_option) == 'strat_file' .or. &
&trim(equilibrium_tau_option) == 'strat_file' ) then
&trim(equilibrium_tau_option) == 'JFV' ) then
! merge with HS troposphere
Pp = 0.
where ( p_full .le. p_hs .and. p_full .ge. p_bd )
Expand Down

0 comments on commit 38816fd

Please sign in to comment.