Program and driver, and data modules roughly split, next step will be…

… cleaning and debugging to get it compiled

src/programs/CMakeLists.txt

@@ -42,7 +42,10 @@ foreach( prec sp dp )
   endif()
   if( HAVE_${prec} )
     ecbuild_add_executable(TARGET  ectrans-benchmark-uni-${prec}
-                           SOURCES ectrans-benchmark-program.F90 ectrans-benchmark-driver.F90
+                           SOURCES 
+			   ectrans-benchmark-program.F90 
+			   ectrans-benchmark-driver.F90
+			   ectrans-benchmark-data.F90
                            LIBS
                              fiat
                              parkind_${prec}

src/programs/ectrans-benchmark-data_mod.f90

@@ -0,0 +1,171 @@
+! (C) Copyright 2014- ECMWF.
+!
+! This software is licensed under the terms of the Apache Licence Version 2.0
+! which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
+! In applying this licence, ECMWF does not waive the privileges and immunities
+! granted to it by virtue of its status as an intergovernmental organisation
+! nor does it submit to any jurisdiction.
+!
+
+module transform_driver_data_mod
+
+use parkind1, only: jpim, jprb, jprd
+
+implicit none
+
+
+integer(kind=jpim) :: istack, getstackusage
+real(kind=jprb), dimension(1) :: zmaxerr(5), zerr(5)
+real(kind=jprb) :: zmaxerrg
+
+! Output unit numbers
+integer(kind=jpim), parameter :: nerr     = 0 ! Unit number for STDERR
+integer(kind=jpim), parameter :: nout     = 6 ! Unit number for STDOUT
+integer(kind=jpim), parameter :: noutdump = 7 ! Unit number for field output
+
+! Default parameters
+integer(kind=jpim) :: nsmax   = 79  ! Spectral truncation
+integer(kind=jpim) :: iters   = 10  ! Number of iterations for transform test
+integer(kind=jpim) :: nfld    = 1   ! Number of scalar fields 
+integer(kind=jpim) :: nlev    = 1   ! Number of vertical levels
+
+integer(kind=jpim) :: nflevg
+integer(kind=jpim) :: ndgl ! Number of latitudes
+integer(kind=jpim) :: nspec2
+integer(kind=jpim) :: ngptot
+integer(kind=jpim) :: ngptotg
+integer(kind=jpim) :: ifld
+integer(kind=jpim) :: jroc
+integer(kind=jpim) :: jb
+integer(kind=jpim) :: nspec2g
+integer(kind=jpim) :: i
+integer(kind=jpim) :: ja
+integer(kind=jpim) :: ib
+integer(kind=jpim) :: jprtrv
+
+integer(kind=jpim), allocatable :: nloen(:), nprcids(:)
+integer(kind=jpim) :: myproc, jj
+integer :: jstep
+
+real(kind=jprd) :: ztinit, ztloop, timef, ztstepmax, ztstepmin, ztstepavg, ztstepmed
+real(kind=jprd) :: ztstepmax1, ztstepmin1, ztstepavg1, ztstepmed1
+real(kind=jprd) :: ztstepmax2, ztstepmin2, ztstepavg2, ztstepmed2
+real(kind=jprd), allocatable :: ztstep(:), ztstep1(:), ztstep2(:)
+
+real(kind=jprb), allocatable :: znormsp(:), znormsp1(:), znormdiv(:), znormdiv1(:)
+real(kind=jprb), allocatable :: znormvor(:), znormvor1(:), znormt(:), znormt1(:)
+real(kind=jprd) :: zaveave(0:jpmaxstat)
+
+! Grid-point space data structures
+real(kind=jprb), allocatable, target :: zgmv   (:,:,:,:) ! Multilevel fields at t and t-dt
+real(kind=jprb), allocatable, target :: zgmvs  (:,:,:)   ! Single level fields at t and t-dt
+real(kind=jprb), pointer :: zgp3a (:,:,:,:) ! Multilevel fields at t and t-dt
+real(kind=jprb), pointer :: zgpuv   (:,:,:,:) ! Multilevel fields at t and t-dt
+real(kind=jprb), pointer :: zgp2 (:,:,:) ! Single level fields at t and t-dt
+
+! Spectral space data structures
+real(kind=jprb), allocatable, target :: sp3d(:,:,:)
+real(kind=jprb), pointer :: zspvor(:,:) => null()
+real(kind=jprb), pointer :: zspdiv(:,:) => null()
+real(kind=jprb), pointer :: zspsc3a(:,:,:) => null()
+real(kind=jprb), allocatable :: zspsc2(:,:)
+
+logical :: lstack = .false. ! Output stack info
+logical :: luserpnm = .false.
+logical :: lkeeprpnm = .false.
+logical :: luseflt = .false. ! Use fast legendre transforms
+logical :: ltrace_stats = .false.
+logical :: lstats_omp = .false.
+logical :: lstats_comms = .false.
+logical :: lstats_mpl = .false.
+logical :: lstats = .true. ! gstats statistics
+logical :: lbarrier_stats = .false.
+logical :: lbarrier_stats2 = .false.
+logical :: ldetailed_stats = .false.
+logical :: lstats_alloc = .false.
+logical :: lsyncstats = .false.
+logical :: lstatscpu = .false.
+logical :: lstats_mem = .false.
+logical :: lxml_stats = .false.
+logical :: lfftw = .true. ! Use FFTW for Fourier transforms
+logical :: lvordiv = .false.
+logical :: lscders = .false.
+logical :: luvders = .false.
+logical :: lprint_norms = .false. ! Calculate and print spectral norms
+logical :: lmeminfo = .false. ! Show information from FIAT routine ec_meminfo at the end
+
+integer(kind=jpim) :: nstats_mem = 0
+integer(kind=jpim) :: ntrace_stats = 0
+integer(kind=jpim) :: nprnt_stats = 1
+
+! The multiplier of the machine epsilon used as a tolerance for correctness checking
+! ncheck = 0 (the default) means that correctness checking is disabled
+integer(kind=jpim) :: ncheck = 0
+
+logical :: lmpoff = .false. ! Message passing switch
+
+! Verbosity level (0 or 1)
+integer :: verbosity = 0
+
+real(kind=jprb) :: zra = 6371229._jprb
+
+integer(kind=jpim) :: nmax_resol = 37 ! Max number of resolutions
+integer(kind=jpim) :: npromatr = 0 ! nproma for trans lib
+integer(kind=jpim) :: ncombflen = 1800000 ! Size of comm buffer
+
+integer(kind=jpim) :: nproc ! Number of procs
+integer(kind=jpim) :: nthread
+integer(kind=jpim) :: nprgpns ! Grid-point decomp
+integer(kind=jpim) :: nprgpew ! Grid-point decomp
+integer(kind=jpim) :: nprtrv = 0 ! Spectral decomp
+integer(kind=jpim) :: nprtrw = 0 ! Spectral decomp
+integer(kind=jpim) :: nspecresmin = 80 ! Minimum spectral resolution, for controlling nprtrw
+integer(kind=jpim) :: mysetv
+integer(kind=jpim) :: mysetw
+integer(kind=jpim) :: mp_type = 2 ! Message passing type
+integer(kind=jpim) :: mbx_size = 150000000 ! Mailbox size
+
+integer(kind=jpim) :: nflevl
+
+! sumpini
+integer(kind=jpim) :: isqr
+logical :: lsync_trans = .true. ! Activate barrier sync
+logical :: leq_regions = .true. ! Eq regions flag
+
+
+integer(kind=jpim) :: nproma = 0
+integer(kind=jpim) :: ngpblks
+! locals
+integer(kind=jpim) :: iprtrv
+integer(kind=jpim) :: iprtrw
+integer(kind=jpim) :: iprused, ilevpp, irest, ilev, jlev
+
+integer(kind=jpim) :: ndimgmv  = 0 ! Third dim. of gmv "(nproma,nflevg,ndimgmv,ngpblks)"
+integer(kind=jpim) :: ndimgmvs = 0 ! Second dim. gmvs "(nproma,ndimgmvs,ngpblks)"
+
+integer(kind=jpim) :: jbegin_uv = 0
+integer(kind=jpim) :: jend_uv   = 0
+integer(kind=jpim) :: jbegin_sc = 0
+integer(kind=jpim) :: jend_sc   = 0
+integer(kind=jpim) :: jbegin_scder_NS = 0
+integer(kind=jpim) :: jend_scder_NS = 0
+integer(kind=jpim) :: jbegin_scder_EW = 0
+integer(kind=jpim) :: jend_scder_EW = 0
+integer(kind=jpim) :: jbegin_uder_EW = 0
+integer(kind=jpim) :: jend_uder_EW = 0
+integer(kind=jpim) :: jbegin_vder_EW = 0
+integer(kind=jpim) :: jend_vder_EW = 0
+
+logical :: ldump_values = .false.
+
+integer, external :: ec_mpirank
+logical :: luse_mpi = .true.
+
+character(len=16) :: cgrid = ''
+
+integer(kind=jpim) :: ierr
+
+
+end module transform_driver_data_mod
+
+!===================================================================================================