Commit 2363b24c authored by Sabine's avatar Sabine
Browse files

added hmix interpolation option to com_mod

parent cea125d8
......@@ -249,7 +249,7 @@ subroutine advance(itime,nrelpoint,ldt,up,vp,wp, &
if (ngrid.le.0) then
do k=1,2
mind=memind(k) ! eso: compatibility with 3-field version
if (interpolhmix)
if (interpolhmix) then
h1(k)=p1*hmix(ix ,jy ,1,mind) &
+ p2*hmix(ixp,jy ,1,mind) &
+ p3*hmix(ix ,jyp,1,mind) &
......
......@@ -752,7 +752,7 @@ module com_mod
logical :: lroot=.true. ! true if serial version, or if MPI .and. root process
logical :: usekernel=.false. ! true if the output kernel shall be switched on
logical :: interpolhmxi=.false. ! true if the hmix shall be interpolated
logical :: interpolhmix=.false. ! true if the hmix shall be interpolated
contains
......
......@@ -131,9 +131,6 @@ subroutine timemanager
loutstart=loutnext-loutaver/2
loutend=loutnext+loutaver/2
! open(127,file=path(2)(1:length(2))//'depostat.dat'
! + ,form='unformatted')
!write (*,*) 'writing deposition statistics depostat.dat!'
!**********************************************************************
! Loop over the whole modelling period in time steps of mintime seconds
......@@ -146,6 +143,8 @@ subroutine timemanager
! print*, 'Initialized lifetime'
!CGZ-lifetime: set lifetime to 0
if (interpolhmix) write(*,*) 'Interpolating hmix'
write(*,46) float(itime)/3600,itime,numpart
if (verbosity.gt.0) then
......
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