Commit a76d954a authored by Sabine's avatar Sabine
Browse files

bugfix: wetscav was used as array and real at the same time

parent db375cc3
...@@ -106,7 +106,7 @@ subroutine timemanager ...@@ -106,7 +106,7 @@ subroutine timemanager
integer :: loutnext,loutstart,loutend integer :: loutnext,loutstart,loutend
integer :: ix,jy,ldeltat,itage,nage,idummy integer :: ix,jy,ldeltat,itage,nage,idummy
integer :: i_nan=0,ii_nan,total_nan_intl=0 !added by mc to check instability in CBL scheme integer :: i_nan=0,ii_nan,total_nan_intl=0 !added by mc to check instability in CBL scheme
real :: outnum,weight,prob_rec(maxspec),prob(maxspec),decfact,wetscav(maxspec) real :: outnum,weight,prob_rec(maxspec),prob(maxspec),decfact,wetscav
! real :: uap(maxpart),ucp(maxpart),uzp(maxpart) ! real :: uap(maxpart),ucp(maxpart),uzp(maxpart)
! real :: us(maxpart),vs(maxpart),ws(maxpart) ! real :: us(maxpart),vs(maxpart),ws(maxpart)
! integer(kind=2) :: cbt(maxpart) ! integer(kind=2) :: cbt(maxpart)
...@@ -570,8 +570,8 @@ subroutine timemanager ...@@ -570,8 +570,8 @@ subroutine timemanager
do ks=1,nspec do ks=1,nspec
if ((xscav_frac1(j,ks).lt.0)) then if ((xscav_frac1(j,ks).lt.0)) then
call get_wetscav(itime,lsynctime,loutnext,j,ks,grfraction,idummy,idummy,wetscav) call get_wetscav(itime,lsynctime,loutnext,j,ks,grfraction,idummy,idummy,wetscav)
if (wetscav(ks).gt.0) then if (wetscav.gt.0) then
xscav_frac1(j,ks)=wetscav(ks)* & xscav_frac1(j,ks)=wetscav* &
(zpoint2(npoint(j))-zpoint1(npoint(j)))*grfraction(1) (zpoint2(npoint(j))-zpoint1(npoint(j)))*grfraction(1)
else else
xmass1(j,ks)=0. xmass1(j,ks)=0.
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment