Commit af2898ae authored by Espen Sollum's avatar Espen Sollum
Browse files

Removed some debugging statements

parent e3790867
......@@ -41,7 +41,7 @@ subroutine ohreaction(itime,ltsample,loutnext)
real :: ohrate,temp,pres
real, parameter :: smallnum = tiny(0.0) ! smallest number that can be handled
real, parameter :: g2mol = 1.e-3*6.0221E23/17.00 ! conversion factor for OH to molec/cm3
real :: w1, maxdiff
real :: w1!, maxdiff
real(kind=dp) :: jul
! Default value (for old OH fields) TODO not needed
......@@ -76,7 +76,7 @@ subroutine ohreaction(itime,ltsample,loutnext)
! ESO DEBUG
idm=0
maxdiff=0.0
! maxdiff=0.0
do jpart=1,numpart
......@@ -152,9 +152,9 @@ subroutine ohreaction(itime,ltsample,loutnext)
if (ohfields_version == 2) then
oh_average=OH_hourly(OHx,OHy,OHz,ohmemind(1)) + &
& w1*(OH_hourly(OHx,OHy,OHz,ohmemind(2))-OH_hourly(OHx,OHy,OHz,ohmemind(1)))
if (abs(OH_hourly(OHx,OHy,OHz,ohmemind(2))-OH_hourly(OHx,OHy,OHz,ohmemind(1))) > maxdiff) then
maxdiff=abs(OH_hourly(OHx,OHy,OHz,ohmemind(2))-OH_hourly(OHx,OHy,OHz,ohmemind(1)))
end if
! if (abs(OH_hourly(OHx,OHy,OHz,ohmemind(2))-OH_hourly(OHx,OHy,OHz,ohmemind(1))) > maxdiff) then
! maxdiff=abs(OH_hourly(OHx,OHy,OHz,ohmemind(2))-OH_hourly(OHx,OHy,OHz,ohmemind(1)))
! end if
else
oh_average=OH_hourly(OHx,OHy,OHz,1)+&
&(OH_hourly(OHx,OHy,OHz,2)-OH_hourly(OHx,OHy,OHz,1))*&
......@@ -167,16 +167,8 @@ subroutine ohreaction(itime,ltsample,loutnext)
!**********************************************************
temp=tt(ix,jy,indz,n)
if (ohfields_version == 2) then
! Convert kg/kg to molec/cm**3
!*****************************
! pres=presOH3D(OHx,OHy,OHz)
! oh_average=oh_average*pres/101325.*298./temp*2.46E19*28.96/17.00
! oh_average=oh_average*rho(ix,jy,indz,n)*g2mol! /1.e3/17.00*6.0221E23
oh_fact=rho(ix,jy,indz,n)*g2mol
end if
if (ohfields_version == 2) oh_fact=rho(ix,jy,indz,n)*g2mol ! Convert kg/kg to molec/cm**3
do k=1,nspec
if (ohcconst(k).gt.0.) then
......@@ -200,7 +192,7 @@ subroutine ohreaction(itime,ltsample,loutnext)
endif ! oh_average.gt.smallnum
end do !continue loop over all particles
write(*,*) 'maxdiff', maxdiff
! write(*,*) 'maxdiff', maxdiff
end subroutine ohreaction
......
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