Changeset c0884a8 in flexpart.git for src/obukhov.f90

Timestamp:

Jul 20, 2018, 2:40:39 PM (6 years ago)

Author:

pesei <petra seibert at univie ac at>

Branches:

univie

Children:

7ca2ef4

Parents:

f251e57

Message:

replace CTBTO code for checking type of GRIB

File:

• src/obukhov.f90 (modified) (4 diffs)

src/obukhov.f90

@@ -20 +20 @@
 !**********************************************************************
 
-real function obukhov(ps,tsurf,tdsurf,tlev,ustar,hf,akm,bkm,plev,metdata_format)
+real function obukhov(ps,tsurf,tdsurf,tlev,ustar,hf,akm,bkm,plev,id_centre)
 
   !********************************************************************
@@ -41 +41 @@
   !       if-then for meteo-type dependent code                       *
   !                                                                   *
+  !                                                                   *
+  !  Petra Seibert, 2018-06-26:                                       *
+  !    simplified version met data format detection                   *
+  !                                                                   *
   !********************************************************************
   !                                                                   *
@@ -54 +58 @@
   !     bkm     ECMWF vertical discretization parameter               *
   !     plev                                                          *
-  !     metdata_format format of metdata (ecmwf/gfs)                  *
+  !     id_centre format of metdata (ecmwf/gfs)                  *
   !                                                                   *
   !********************************************************************
 
   use par_mod
-  use class_gribfile
+  use check_gribfile_mod
 
   implicit none
 
-  integer :: metdata_format
+  integer :: id_centre
   real :: akm(nwzmax),bkm(nwzmax)
   real :: ps,tsurf,tdsurf,tlev,ustar,hf,e,ew,tv,rhoa,plev
@@ -72 +76 @@
   tv=tsurf*(1.+0.378*e/ps)               ! virtual temperature
   rhoa=ps/(r_air*tv)                      ! air density
-  if (metdata_format.eq.GRIBFILE_CENTRE_ECMWF) then
+  if (id_centre.eq.icg_id_ecmwf) then
   ak1=(akm(1)+akm(2))/2.
   bk1=(bkm(1)+bkm(2))/2.