!********************************************************************** ! Copyright 1998,1999,2000,2001,2002,2005,2007,2008,2009,2010 * ! Andreas Stohl, Petra Seibert, A. Frank, Gerhard Wotawa, * ! Caroline Forster, Sabine Eckhardt, John Burkhart, Harald Sodemann * ! * ! This file is part of FLEXPART. * ! * ! FLEXPART is free software: you can redistribute it and/or modify * ! it under the terms of the GNU General Public License as published by* ! the Free Software Foundation, either version 3 of the License, or * ! (at your option) any later version. * ! * ! FLEXPART is distributed in the hope that it will be useful, * ! but WITHOUT ANY WARRANTY; without even the implied warranty of * ! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * ! GNU General Public License for more details. * ! * ! You should have received a copy of the GNU General Public License * ! along with FLEXPART. If not, see . * !********************************************************************** integer function detectformat() !***************************************************************************** ! * ! This routine reads the 1st file with windfields to determine * ! the format. * ! * ! Authors: M. Harustak * ! * ! 6 May 2015 * ! * ! Unified ECMWF and GFS builds * ! Marian Harustak, 12.5.2017 * ! - Added routine to FP10 Flexpart distribution * !***************************************************************************** ! * ! Variables: * ! fname file name of file to check * ! * !***************************************************************************** use par_mod use com_mod use class_gribfile implicit none character(len=255) :: filename character(len=255) :: wfname1(maxwf) integer :: metdata_format ! If no file is available if ( maxwf.le.0 ) then print*,'No wind file available' detectformat = GRIBFILE_CENTRE_UNKNOWN return endif ! construct filename filename = path(3)(1:length(3)) // trim(wfname(1)) ! get format detectformat = gribfile_centre(TRIM(filename)) end