From 5c58ba4225ad147158658da89eca666c6c2d7e4f Mon Sep 17 00:00:00 2001 From: Magnus Ulimoen Date: Mon, 16 Oct 2023 16:36:03 +0200 Subject: [PATCH] Fix computation of vgrav tables --- src/common/vgravtables.f90 | 20 ++++++++++++++++++-- 1 file changed, 18 insertions(+), 2 deletions(-) diff --git a/src/common/vgravtables.f90 b/src/common/vgravtables.f90 index 2ee33bcc..b1d28cd0 100644 --- a/src/common/vgravtables.f90 +++ b/src/common/vgravtables.f90 @@ -77,13 +77,13 @@ subroutine vgravtables select case(drydep_scheme) case (DRYDEP_SCHEME_EMEP,DRYDEP_SCHEME_ZHANG,DRYDEP_SCHEME_EMERSON) ! expected kg/m3 - rho_part = def_comp(m)%densitygcm3 / 1000.0 + rho_part = def_comp(m)%densitygcm3 * 1e3 ! expected m diam_part = 2 * def_comp(m)%radiusmym / 1e6 do ip=1,numpresvg ! Expecting pascal - p = max(1.0, pbasevg + ip*pincrvg) / 100.0 + p = max(1.0, pbasevg + ip*pincrvg) * 100.0 do it=1,numtempvg t = tbasevg + it*tincrvg roa = p / (real(t, kind=real64) * R) @@ -113,6 +113,22 @@ subroutine vgravtables end do end do end select + + block + use snapdebug, only: iulog, idebug + if (idebug == 1) then + write(iulog,*) "VGRAV PRINTOUT START COMPONENT ", n + write(iulog,*) "p t vg" + do it=1,size(vgtable,1) + do ip=1,size(vgtable,2) + t = tbasevg + it*tincrvg + p = pbasevg + ip*pincrvg + write(iulog,'(F7.1 F7.1,F14.6)') p, t, vgtable(it,ip,n) + end do + end do + write(iulog,*) "VGRAV PRINTOUT END" + endif + end block end do do_comp end subroutine vgravtables