Changeset db91eb7 in flexpart.git for src/readwind_gfs.f90

Timestamp:

Dec 5, 2018, 11:52:16 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:

e9e0f06, 3d7eebf

Parents:

79e0349

Message:

reading clouds in NCEP

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
 
 
@@ -120 +121 @@
 
   !HSO
-5   call grib_open_file(ifile,path(3)(1:length(3)) &
+  call grib_open_file(ifile,path(3)(1:length(3)) &
          //trim(wfname(indj)),'r',iret)
   if (iret.ne.GRIB_SUCCESS) then
@@ -162 +163 @@
 
   !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)
@@ -173 +176 @@
        valSurf,iret)
 !  call grib_check(iret,gribFunction,gribErrorMsg)
-
+  
+!  write(*,*) 'Field: ',ifield,parCat,parNum,typSurf,shortname
   !convert to grib1 identifiers
   isec1(6)=-1
@@ -214 +218 @@
     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
@@ -548 +556 @@
         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
@@ -675 +690 @@
     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