diff --git a/src/getfields.f90 b/src/getfields.f90 index fc4061484314c614ce844f37bdcd66372c3f3489..bfc24d25a8b029f958827da19ddc9070806189dc 100644 --- a/src/getfields.f90 +++ b/src/getfields.f90 @@ -137,6 +137,10 @@ subroutine getfields(itime,nstop) end do 40 indmin=indj + if (WETBKDEP) then + call writeprecip(itime,memind(1)) + endif + else ! No wind fields, which can be used, are currently in memory @@ -168,6 +172,10 @@ subroutine getfields(itime,nstop) end do 60 indmin=indj + if (WETBKDEP) then + call writeprecip(itime,memind(1)) + endif + endif lwindinterv=abs(memtime(2)-memtime(1)) diff --git a/src/makefile b/src/makefile index 35c81e469a047a197b8bba091e1ffbe4a37d50e1..7c5d15fb5afd48e8df388cd582a8d16fe6157d46 100644 --- a/src/makefile +++ b/src/makefile @@ -143,6 +143,7 @@ OBJECTS_GFS = \ OBJECTS = \ advance.o initialize.o \ writeheader.o writeheader_txt.o \ +writeprecip.o \ writeheader_surf.o assignland.o\ part0.o gethourlyOH.o\ caldate.o partdep.o \ @@ -431,6 +432,7 @@ wetdepo.o: com_mod.o par_mod.o point_mod.o wetdepokernel.o: com_mod.o par_mod.o unc_mod.o wetdepokernel_nest.o: com_mod.o par_mod.o unc_mod.o writeheader.o: com_mod.o outg_mod.o par_mod.o point_mod.o +writeprecip.o: com_mod.o par_mod.o point_mod.o writeheader_nest.o: com_mod.o outg_mod.o par_mod.o point_mod.o writeheader_nest_surf.o: com_mod.o outg_mod.o par_mod.o point_mod.o writeheader_surf.o: com_mod.o outg_mod.o par_mod.o point_mod.o diff --git a/src/par_mod.f90 b/src/par_mod.f90 index 5f45767882d3e7061d5efb4b6d9547042a9a4df3..0f1a645b6d8f16421898d1b21d2e00f2ffca6634 100644 --- a/src/par_mod.f90 +++ b/src/par_mod.f90 @@ -251,7 +251,7 @@ module par_mod integer,parameter :: unitspecies=1, unitoutrecept=91, unitoutreceptppt=92 integer,parameter :: unitlsm=1, unitsurfdata=1, unitland=1, unitwesely=1 integer,parameter :: unitOH=1 - integer,parameter :: unitdates=94, unitheader=90,unitheader_txt=100, unitshortpart=95 + integer,parameter :: unitdates=94, unitheader=90,unitheader_txt=100, unitshortpart=95, unitprecip=101 integer,parameter :: unitboundcond=89 integer,parameter :: unittmp=101 diff --git a/src/timemanager.f90 b/src/timemanager.f90 index 21f6ca3b934154ec5be23886136af234355fd518..7ae91dbe0401a95b6f40e80b0657934093994abc 100644 --- a/src/timemanager.f90 +++ b/src/timemanager.f90 @@ -147,6 +147,8 @@ subroutine timemanager !CGZ-lifetime: set lifetime to 0 if (.not.usekernel) write(*,*) 'Not using the kernel' + if (turboff) write(*,*) 'Turbulence switched off' + write(*,46) float(itime)/3600,itime,numpart if (verbosity.gt.0) then