Changeset 7873bf7 in flexpart.git

Timestamp:

Apr 19, 2019, 10:40:41 AM (5 years ago)

Author:

Sabine <sabine.eckhardt@…>

Branches:

master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug

Children:

5d74ed9

Parents:

437c545

Message:

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

Location:

src

Files:

• readwind_gfs.f90 (modified) (2 diffs)
• releaseparticles.f90 (modified) (1 diff)

src/readwind_gfs.f90

@@ -134 +134 @@
   numpw=0
   numprh=0
+  numpclwch=0
   ifield=0
 10   ifield=ifield+1
@@ -557 +558 @@
       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
-        else
-          clwch(i-i181,j,numpu,n)=help
+          clwch(i179+i,j,numpclwch,n)=help
+        else
+          clwch(i-i181,j,numpclwch,n)=help
         endif
         readclouds=.true.

src/releaseparticles.f90

@@ -114 +114 @@
       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