Changes in src/readwind_gfs.f90 [5d74ed9:db91eb7] in flexpart.git

File:

• src/readwind_gfs.f90 (modified) (3 diffs)

src/readwind_gfs.f90

@@ -83 +83 @@
 
   ! NCEP
-  integer :: numpt,numpu,numpv,numpw,numprh,numpclwch
+  integer :: numpt,numpu,numpv,numpw,numprh
   real :: help, temp, ew
   real :: elev
@@ -134 +134 @@
   numpw=0
   numprh=0
-  numpclwch=0
   ifield=0
 10   ifield=ifield+1
@@ -558 +557 @@
       endif
 ! SEC & IP 12/2018 read GFS clouds
-      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)) numpclwch=ii
-            end do
-        endif
-        help=zsec4(nxfield*(ny-j-1)+i+1)
-        if(i.le.i180) then
-          clwch(i179+i,j,numpclwch,n)=help
-        else
-          clwch(i-i181,j,numpclwch,n)=help
-        endif
+      if(isec1(6).eq.153) then  !! CLWCR  Cloud liquid water content [kg/kg]
+        clwch(i,j,nlev_ec-k+2,n)=zsec4(nxfield*(ny-j-1)+i+1)
         readclouds=.true.
         sumclouds=.true.
-!        readclouds=.false.
-!       sumclouds=.false.
       endif