Changes in src/obukhov.f90 [6ecb30a:e200b7a] in flexpart.git

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)
 
   !********************************************************************
@@ -27 +27 @@
   !                       Date:   1994-06-27                          *
   !                                                                   *
-  !     This program calculates Obukhov scale height from surface     *
-  !     meteorological data and sensible heat flux.                   *
+  !  Update: A. Stohl, 2000-09-25, avoid division by zero by          *
+  !  setting ustar to minimum value                                   *
   !                                                                   *
   !********************************************************************
   !                                                                   *
-  !  Update: A. Stohl, 2000-09-25, avoid division by zero by          *
-  !  setting ustar to minimum value                                   *
-  !  CHANGE: 17/11/2005 Caroline Forster NCEP GFS version             *
-  !                                                                   *
-  !   Unified ECMWF and GFS builds                                    *
-  !   Marian Harustak, 12.5.2017                                      *
-  !     - Merged obukhov and obukhov_gfs into one routine using       *
-  !       if-then for meteo-type dependent code                       *
+  !     This program calculates Obukhov scale height from surface     *
+  !     meteorological data and sensible heat flux.                   *
   !                                                                   *
   !********************************************************************
@@ -53 +47 @@
   !     akm     ECMWF vertical discretization parameter               *
   !     bkm     ECMWF vertical discretization parameter               *
-  !     plev                                                          *
-  !     metdata_format format of metdata (ecmwf/gfs)                  *
   !                                                                   *
   !********************************************************************
 
   use par_mod
-  use class_gribfile
 
   implicit none
 
-  integer :: metdata_format
   real :: akm(nwzmax),bkm(nwzmax)
   real :: ps,tsurf,tdsurf,tlev,ustar,hf,e,ew,tv,rhoa,plev
@@ -72 +62 @@
   tv=tsurf*(1.+0.378*e/ps)               ! virtual temperature
   rhoa=ps/(r_air*tv)                      ! air density
-  if (metdata_format.eq.GRIBFILE_CENTRE_ECMWF) then
   ak1=(akm(1)+akm(2))/2.
   bk1=(bkm(1)+bkm(2))/2.
   plev=ak1+bk1*ps                        ! Pressure level 1
-  end if
   theta=tlev*(100000./plev)**(r_air/cpa) ! potential temperature
   if (ustar.le.0.) ustar=1.e-8