Commit d1a87077 authored by Sabine's avatar Sabine
Browse files

When running in wet deposition mode a textfile with precipitation data is written out

parent 58882983
...@@ -137,6 +137,10 @@ subroutine getfields(itime,nstop) ...@@ -137,6 +137,10 @@ subroutine getfields(itime,nstop)
end do end do
40 indmin=indj 40 indmin=indj
if (WETBKDEP) then
call writeprecip(itime,memind(1))
endif
else else
! No wind fields, which can be used, are currently in memory ! No wind fields, which can be used, are currently in memory
...@@ -168,6 +172,10 @@ subroutine getfields(itime,nstop) ...@@ -168,6 +172,10 @@ subroutine getfields(itime,nstop)
end do end do
60 indmin=indj 60 indmin=indj
if (WETBKDEP) then
call writeprecip(itime,memind(1))
endif
endif endif
lwindinterv=abs(memtime(2)-memtime(1)) lwindinterv=abs(memtime(2)-memtime(1))
......
...@@ -143,6 +143,7 @@ OBJECTS_GFS = \ ...@@ -143,6 +143,7 @@ OBJECTS_GFS = \
OBJECTS = \ OBJECTS = \
advance.o initialize.o \ advance.o initialize.o \
writeheader.o writeheader_txt.o \ writeheader.o writeheader_txt.o \
writeprecip.o \
writeheader_surf.o assignland.o\ writeheader_surf.o assignland.o\
part0.o gethourlyOH.o\ part0.o gethourlyOH.o\
caldate.o partdep.o \ caldate.o partdep.o \
...@@ -431,6 +432,7 @@ wetdepo.o: com_mod.o par_mod.o point_mod.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.o: com_mod.o par_mod.o unc_mod.o
wetdepokernel_nest.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 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.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_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 writeheader_surf.o: com_mod.o outg_mod.o par_mod.o point_mod.o
......
...@@ -251,7 +251,7 @@ module par_mod ...@@ -251,7 +251,7 @@ module par_mod
integer,parameter :: unitspecies=1, unitoutrecept=91, unitoutreceptppt=92 integer,parameter :: unitspecies=1, unitoutrecept=91, unitoutreceptppt=92
integer,parameter :: unitlsm=1, unitsurfdata=1, unitland=1, unitwesely=1 integer,parameter :: unitlsm=1, unitsurfdata=1, unitland=1, unitwesely=1
integer,parameter :: unitOH=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 :: unitboundcond=89
integer,parameter :: unittmp=101 integer,parameter :: unittmp=101
......
...@@ -147,6 +147,8 @@ subroutine timemanager ...@@ -147,6 +147,8 @@ subroutine timemanager
!CGZ-lifetime: set lifetime to 0 !CGZ-lifetime: set lifetime to 0
if (.not.usekernel) write(*,*) 'Not using the kernel' if (.not.usekernel) write(*,*) 'Not using the kernel'
if (turboff) write(*,*) 'Turbulence switched off'
write(*,46) float(itime)/3600,itime,numpart write(*,46) float(itime)/3600,itime,numpart
if (verbosity.gt.0) then 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