Commit 7873bf79 authored by Sabine's avatar Sabine

added isec1 for CLWR in gfs, pointsource, smaller range according to Chrstian Maurer

parent 437c545f
......@@ -133,6 +133,7 @@ subroutine readwind_gfs(indj,n,uuh,vvh,wwh)
numpv=0
numpw=0
numprh=0
numpclwch=0
ifield=0
10 ifield=ifield+1
!
......@@ -556,17 +557,17 @@ subroutine readwind_gfs(indj,n,uuh,vvh,wwh)
endif
endif
! SEC & IP 12/2018 read GFS clouds
if(isec1(6).eq.153) then !! CLWCR Cloud liquid water content [kg/kg]
if((isec1(6).eq.153).and.(isec1(7).eq.100)) then !! CLWCR Cloud liquid water content [kg/kg]
if((i.eq.0).and.(j.eq.0)) then
do ii=1,nuvz
if ((isec1(8)*100.0).eq.akz(ii)) numpu=ii
if ((isec1(8)*100.0).eq.akz(ii)) numpclwch=ii
end do
endif
help=zsec4(nxfield*(ny-j-1)+i+1)
if(i.le.i180) then
clwch(i179+i,j,numpu,n)=help
clwch(i179+i,j,numpclwch,n)=help
else
clwch(i-i181,j,numpu,n)=help
clwch(i-i181,j,numpclwch,n)=help
endif
readclouds=.true.
sumclouds=.true.
......
......@@ -113,7 +113,8 @@ subroutine releaseparticles(itime)
!*****************************************************************************
average_timecorrect=0.
do k=1,nspec
if (zpoint1(i).gt.0.5) then ! point source
if(abs(xpoint2(i)-xpoint1(i)).lt.1.E-4.and.abs(ypoint2(i)-ypoint1(i)).lt.1.E-4) then
! if (zpoint1(i).gt.0.5) then ! point source
timecorrect(k)=point_hour(k,nhour)*point_dow(k,ndayofweek)
else ! area source
timecorrect(k)=area_hour(k,nhour)*area_dow(k,ndayofweek)
......
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