Changes in src/readwind_gfs.f90 [db91eb7:6ecb30a] in flexpart.git

File:

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

src/readwind_gfs.f90

@@ -1 @@
-
+!**********************************************************************
 ! Copyright 1998,1999,2000,2001,2002,2005,2007,2008,2009,2010         *
 ! Andreas Stohl, Petra Seibert, A. Frank, Gerhard Wotawa,             *
@@ -108 +108 @@
   character(len=24) :: gribErrorMsg = 'Error reading grib file'
   character(len=20) :: gribFunction = 'readwind_gfs'
-  character(len=20) :: shortname
 
 
@@ -121 +120 @@
 
   !HSO
-  call grib_open_file(ifile,path(3)(1:length(3)) &
+5   call grib_open_file(ifile,path(3)(1:length(3)) &
          //trim(wfname(indj)),'r',iret)
   if (iret.ne.GRIB_SUCCESS) then
@@ -163 +162 @@
 
   !read the grib2 identifiers
-  call grib_get_string(igrib,'shortName',shortname,iret)
-
   call grib_get_int(igrib,'discipline',discipl,iret)
 !  call grib_check(iret,gribFunction,gribErrorMsg)
@@ -176 +173 @@
        valSurf,iret)
 !  call grib_check(iret,gribFunction,gribErrorMsg)
-  
-!  write(*,*) 'Field: ',ifield,parCat,parNum,typSurf,shortname
+
   !convert to grib1 identifiers
   isec1(6)=-1
@@ -218 +214 @@
     isec1(7)=105         ! indicatorOfTypeOfLevel
     isec1(8)=10
-  elseif ((parCat.eq.1).and.(parNum.eq.22).and.(typSurf.eq.100)) then ! CLWMR Cloud Mixing Ratio [kg/kg]:
-    isec1(6)=153         ! indicatorOfParameter
-    isec1(7)=100         ! indicatorOfTypeOfLevel
-    isec1(8)=valSurf/100 ! level, convert to hPa
   elseif ((parCat.eq.3).and.(parNum.eq.1).and.(typSurf.eq.101)) then ! SLP
     isec1(6)=2           ! indicatorOfParameter
@@ -556 +548 @@
         endif
       endif
-! SEC & IP 12/2018 read GFS clouds
-      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.
-      endif
-
 
     end do
@@ -690 +675 @@
     call shift_field(vvh,nxfield,ny,nuvzmax,nuvz,1,1)
     call shift_field(wwh,nxfield,ny,nwzmax,nwz,1,1)
-! IP & SEC adding GFS Clouds 20181205
-    call shift_field(clwch,nxfield,ny,nuvzmax,nuvz,2,n)
   endif