#include "swancpp.h"
! ALLOCATABLE DATA RELATED MODULES, file 2 of 3
!
! Contents of this file
!
! M_WCAP information for whitecapping
! OUTP_DATA information for output data
! M_SNL4 information for quadruplets
! M_BNDSPEC information for boundary conditions
! M_OBSTA information for obstacles
! M_GENARR contains a number of general arrays
! M_PARALL information for parallelisation with MPI
! M_DIFFR information for diffraction
!
#ifdef SWAN_MODEL
! For COAWST:
! M_COUPLING information for grid numbers
! M_CVMESH information for the cvmesh routine
! M_PARALL2 needed to have allocation of main vars.
! M_SREFINED used for refined grids
!
#endif
MODULE M_WCAP
!
!
! --|-----------------------------------------------------------|--
! | Delft University of Technology |
! | Faculty of Civil Engineering |
! | Environmental Fluid Mechanics Section |
! | P.O. Box 5048, 2600 GA Delft, The Netherlands |
! | |
! | Programmers: The SWAN team |
! --|-----------------------------------------------------------|--
!
!
! SWAN (Simulating WAves Nearshore); a third generation wave model
! Copyright (C) 1993-2017 Delft University of Technology
!
! This program 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 2 of
! the License, or (at your option) any later version.
!
! This program 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.
!
! A copy of the GNU General Public License is available at
! http://www.gnu.org/copyleft/gpl.html#SEC3
! or by writing to the Free Software Foundation, Inc.,
! 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
!
!
! 0. Authors
!
! 40.02: IJsbrand Haagsma
!
! 1. Updates
!
! Sep. 00: New Module
!
! 2. Purpose
!
! Create global variables used in whitecapping and integral parameter
! subroutines
!
! 3. Method
!
! MODULE construct
!
! 4. Modules used
!
! ---
!
IMPLICIT NONE
!
! 5. Argument variables
!
! ---
!
! 6. Parameter variables
!
! ---
!
! 7. Local variables
!
! ACTOT : Total action density per gridpoint
! EDRKTOT: Zeroth moment of energy / SQRT(wavenumber)
! EKTOT : Zeroth moment of energy * wavenumber
! ETOT1 : First moment of the energy density
! ETOT2 : Second moment of the energy density
! ETOT4 : Fourth moment of the energy density
! KM_WAM : Mean wavenumber according to the WAM-formulation
! KM01 : Mean wavenumber according to first and zeroth order moments
! SIGM_10: Mean frequency according to zeroth and first order moments
! SIGM01 : Mean frequency according to first and zeroth order moments
!
REAL, SAVE :: ACTOT
REAL, SAVE :: EDRKTOT
REAL, SAVE :: EKTOT
REAL, SAVE :: ETOT1
REAL, SAVE :: ETOT2
REAL, SAVE :: ETOT4
REAL, SAVE :: KM_WAM
REAL, SAVE :: KM01
REAL, SAVE :: SIGM_WAM
REAL, SAVE :: SIGM_10
REAL, SAVE :: SIGM01
!
!$OMP THREADPRIVATE(ACTOT, EDRKTOT, EKTOT, ETOT1, ETOT2, ETOT4,
!$OMP& KM_WAM, KM01, SIGM_WAM,
!$OMP& SIGM_10, SIGM01)
!
! SIGPOW : contains powers of relative frequencies
! second dimension indicates power of sigma
!
REAL, SAVE, ALLOCATABLE :: SIGPOW(:,:)
!
! 8. Subroutines and functions used
!
! ---
!
! 9. Subroutines and functions calling
!
! SSFILL :
! SINTGRL: Calculating integral paramters
! SWCAP : Calculating whitecapping source term
!
! 10. Error messages
!
! ---
!
! 11. Remarks
!
! ---
!
! 12. Structure
!
! ---
!
! 13. Source text
!
END MODULE M_WCAP
MODULE OUTP_DATA
!
!
! --|-----------------------------------------------------------|--
! | Delft University of Technology |
! | Faculty of Civil Engineering |
! | Environmental Fluid Mechanics Section |
! | P.O. Box 5048, 2600 GA Delft, The Netherlands |
! | |
! | Programmers: The SWAN team |
! --|-----------------------------------------------------------|--
!
!
! SWAN (Simulating WAves Nearshore); a third generation wave model
! Copyright (C) 1993-2017 Delft University of Technology
!
! This program 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 2 of
! the License, or (at your option) any later version.
!
! This program 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.
!
! A copy of the GNU General Public License is available at
! http://www.gnu.org/copyleft/gpl.html#SEC3
! or by writing to the Free Software Foundation, Inc.,
! 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
!
!
! 0. Authors
!
! 40.13: Nico Booij
! 40.31: Marcel Zijlema
!
! 1. Updates
!
! 40.13, July 01: New Module
! 40.13, Oct. 01: Longer filenames for output requests
! 40.31, Dec. 03: derive types OPSDAT, ORQDAT added
!
! 2. Purpose
!
! Contains data needed during generation of output
!
! 3. Method
!
! MODULE construct
!
! 4. Modules used
!
USE OCPCOMM2
!
IMPLICIT NONE
!
! 5. Argument variables
!
! ---
!
! 6. Parameter variables
!
! ---
!
! 7. Local variables
#ifdef SWAN_MODEL
INTEGER, PARAMETER :: MAX_OUTP_REQ = 25000 ! max. number of output requests
#else
INTEGER, PARAMETER :: MAX_OUTP_REQ = 250 ! max. number of output requests
#endif
CHARACTER (LEN=1) :: OUT_COMMENT = '%' ! comment sign for heading lines
LOGICAL :: LCOMPGRD
! formats for output:
CHARACTER (LEN=40) :: FLT_BLOCK = '(6E12.4)' ! floating point block
CHARACTER (LEN=40) :: FLT_TABLE = '(E11.4)' ! floating point table
CHARACTER (LEN=40) :: FIX_SPEC = '(200(1X,I4))' ! spectral output
! format for block output per process in case of collecting data
CHARACTER (LEN=40) :: FLT_BLKP = '(6E17.9)'
INTEGER :: FLD_TABLE = 12 ! field length for fixed-point table
INTEGER :: DEC_BLOCK = 4 ! number of decimals for fixed-point block
INTEGER :: DEC_SPEC = 4 ! number of decimals for spectral output
#ifdef SWAN_MODEL
! longer filenames for output requests
CHARACTER (LEN=LENFNM), POINTER :: OUTP_FILES(:)
CHARACTER (LEN=LENFNM), POINTER :: OUTP_FILES_G(:,:)
! filenames for output; index is output request sequence number
#else
! longer filenames for output requests
CHARACTER (LEN=LENFNM) :: OUTP_FILES(1:MAX_OUTP_REQ)
! filenames for output; index is output request sequence number
#endif
#ifdef SWAN_MODEL
INTEGER, POINTER, SAVE :: NREOQ != 0 ! actual number of requests saved
INTEGER, POINTER, SAVE :: NREOQ_G(:)
TYPE OPSDAT
CHARACTER (LEN=1) :: PSTYPE ! type (F, C, P, ...)
CHARACTER (LEN=8) :: PSNAME ! name of point set
INTEGER :: OPI(2) ! integer coefficients
REAL :: OPR(5) ! real coefficients
INTEGER :: MIP ! number of points
REAL, POINTER :: XP(:), YP(:), XQ(:), YQ(:) ! point coordinates
TYPE(OPSDAT), POINTER :: NEXTOPS
END TYPE OPSDAT
TYPE OPSDAT_Z
! TYPE(OPSDAT), POINTER :: CUOPS
TYPE(OPSDAT), POINTER :: FOPS
TYPE(OPSDAT), POINTER :: COPS
! TYPE(OPSDAT), POINTER :: ROPS
TYPE(OPSDAT), POINTER :: OPSTMP
END TYPE OPSDAT_Z
TYPE (OPSDAT_Z), ALLOCATABLE :: OPSDATZ_MOD(:)
LOGICAL, POINTER, SAVE :: LOPS ! = .FALSE.
LOGICAL, POINTER, SAVE :: LOPS_G(:)
TYPE ORQDAT
CHARACTER (LEN=4) :: RQTYPE ! type (BLK, TAB, SPC ...)
CHARACTER (LEN=8) :: PSNAME ! name of point set
INTEGER :: OQI(4) ! integer coefficients
REAL*8 :: OQR(2) ! real coefficients
INTEGER, POINTER :: IVTYP(:) ! type of output variable
REAL, POINTER :: FAC(:) ! multiplication factor of block output
TYPE(ORQDAT), POINTER :: NEXTORQ
END TYPE ORQDAT
! TYPE(ORQDAT), SAVE, TARGET :: FORQ
TYPE ORQDAT_Z
TYPE(ORQDAT), POINTER :: FORQ
! TYPE(ORQDAT), POINTER :: ORQTMP
END TYPE ORQDAT_Z
TYPE (ORQDAT_Z), ALLOCATABLE :: ORQDATZ_MOD(:)
LOGICAL, POINTER, SAVE :: LORQ
LOGICAL, POINTER, SAVE :: LORQ_G(:)
INTEGER, SAVE, POINTER :: IREC(:)
INTEGER, SAVE, POINTER :: IREC_G(:,:)
#else
INTEGER, SAVE :: NREOQ = 0 ! actual number of requests saved
TYPE OPSDAT
CHARACTER (LEN=1) :: PSTYPE ! type (F, C, P, ...)
CHARACTER (LEN=8) :: PSNAME ! name of point set
INTEGER :: OPI(2) ! integer coefficients
REAL :: OPR(5) ! real coefficients
INTEGER :: MIP ! number of points
REAL, POINTER :: XP(:), YP(:), XQ(:), YQ(:) ! point coordinates
TYPE(OPSDAT), POINTER :: NEXTOPS
END TYPE OPSDAT
TYPE(OPSDAT), SAVE, TARGET :: FOPS
TYPE(OPSDAT), SAVE, POINTER :: COPS
LOGICAL, SAVE :: LOPS = .FALSE.
TYPE ORQDAT
CHARACTER (LEN=4) :: RQTYPE ! type (BLK, TAB, SPC ...)
CHARACTER (LEN=8) :: PSNAME ! name of point set
INTEGER :: OQI(4) ! integer coefficients
REAL*8 :: OQR(2) ! real coefficients
INTEGER, POINTER :: IVTYP(:) ! type of output variable
REAL, POINTER :: FAC(:) ! multiplication factor of block output
TYPE(ORQDAT), POINTER :: NEXTORQ
END TYPE ORQDAT
TYPE(ORQDAT), SAVE, TARGET :: FORQ
#endif
!
! 8. Subroutines and functions used
!
! ---
!
! 9. Subroutines and functions calling
!
! SWREAD : reads data (command OUTPut OPTions)
! SWBLOK : produces block output
! SWTABP : produces table output
!
! 10. Error messages
!
! ---
!
! 11. Remarks
!
! ---
!
! 12. Structure
!
! ---
!
! 13. Source text
!
END MODULE OUTP_DATA
MODULE M_SNL4
!
!
! --|-----------------------------------------------------------|--
! | Delft University of Technology |
! | Faculty of Civil Engineering |
! | Environmental Fluid Mechanics Section |
! | P.O. Box 5048, 2600 GA Delft, The Netherlands |
! | |
! | Programmers: The SWAN team |
! --|-----------------------------------------------------------|--
!
!
! SWAN (Simulating WAves Nearshore); a third generation wave model
! Copyright (C) 1993-2017 Delft University of Technology
!
! This program 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 2 of
! the License, or (at your option) any later version.
!
! This program 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.
!
! A copy of the GNU General Public License is available at
! http://www.gnu.org/copyleft/gpl.html#SEC3
! or by writing to the Free Software Foundation, Inc.,
! 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
!
!
! 0. Authors
!
! 40.17: IJsbrand Haagsma
!
! 1. Updates
!
! Feb. 01: New Module
!
! 2. Purpose
!
! Create global variables used in quadruplet subroutines
!
! 3. Method
!
! MODULE construct
!
! 4. Modules used
!
! ---
!
IMPLICIT NONE
!
! 5. Argument variables
!
! ---
!
! 6. Parameter variables
!
! ---
!
! 7. Local variables
! MDIA : Number of quadruplets in the MDIA formulation
INTEGER, PUBLIC, SAVE :: MDIA = 1
! AF11 : Contains the scaling frequency for the DIA.
! CNL4_1: Contains the values for C1 in the MDIA formulation.
! CNL4_2: Contains the values for C2 in the MDIA formulation.
! LAMBDA: Contains the values for lambda in the MDIA formulation.
REAL, PUBLIC, SAVE, ALLOCATABLE :: AF11(:)
REAL, PUBLIC, SAVE, ALLOCATABLE :: CNL4_1(:)
REAL, PUBLIC, SAVE, ALLOCATABLE :: CNL4_2(:)
REAL, PUBLIC, SAVE, ALLOCATABLE :: LAMBDA(:)
! 8. Subroutines and functions used
!
! ---
!
! 9. Subroutines and functions calling
!
! ---
!
! 10. Error messages
!
! ---
!
! 11. Remarks
!
! ---
!
! 12. Structure
!
! ---
!
! 13. Source text
!
END MODULE M_SNL4
MODULE M_BNDSPEC
!
!
! --|-----------------------------------------------------------|--
! | Delft University of Technology |
! | Faculty of Civil Engineering and Geosciences |
! | Environmental Fluid Mechanics Section |
! | P.O. Box 5048, 2600 GA Delft, The Netherlands |
! | |
! | Programmer: Marcel Zijlema |
! --|-----------------------------------------------------------|--
!
!
! SWAN (Simulating WAves Nearshore); a third generation wave model
! Copyright (C) 1993-2017 Delft University of Technology
!
! This program 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 2 of
! the License, or (at your option) any later version.
!
! This program 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.
!
! A copy of the GNU General Public License is available at
! http://www.gnu.org/copyleft/gpl.html#SEC3
! or by writing to the Free Software Foundation, Inc.,
! 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
!
!
! 0. Authors
!
! 40.31: Marcel Zijlema
!
! 1. Updates
!
! Nov. 03: New Module
!
! 2. Purpose
!
! Contains data with respect to specification
! of boundary conditions
!
! 3. Method
!
! ---
!
! 4. Modules used
!
! ---
!
IMPLICIT NONE
!
! 5. Argument variables
!
! ---
!
! 6. Parameter variables
!
! ---
!
! 7. Local variables
!
! BFILED : data concerning boundary condition files
! BGP : array containing data w.r.t. boundary grid points
! BSPLOC : place in array BSPECS where to store interpolated spectra
! BSPDIR : spectral directions of input spectrum
! BSPFRQ : spectral frequencies of input spectrum
! CUBGP : current item in list of boundary grid points
! DSHAPE : indicates option for computation of directional distribution
! in the spectrum (boundary spectra etc.)
! =1: directional spread in degrees is given
! =2: power of COS is given
! FBNDFIL : first boundary condition file in list of files
! FBGP : first item in list of boundary grid points
! FBS : first item in list of boundary spectrum parameters
! FSHAPE : indicates option for computation of frequency distribution
! in the spectrum (boundary spectra etc.)
! =1: Pierson-Moskowitz
! =2: Jonswap
! =3: bin
! =4: Gaussian
! NBS : index of BSPECS
! NEXTBGP : pointer to next item in list of boundary grid points
! NEXTBS : pointer to next item in list of boundary spectrum parameters
! NEXTBSPC: pointer to next boundary condition file in list
! SPPARM : integral parameters used for computation of
! incident spectrum. Meaning:
! 1: significant wave height
! 2: wave period (peak or mean)
! 3: average wave direction
! 4: directional distribution coefficient
TYPE BSPCDAT
INTEGER :: BFILED(20)
INTEGER, POINTER :: BSPLOC(:)
REAL, POINTER :: BSPDIR(:), BSPFRQ(:)
TYPE(BSPCDAT), POINTER :: NEXTBSPC
END TYPE BSPCDAT
TYPE(BSPCDAT), SAVE, TARGET :: FBNDFIL
TYPE BSDAT
INTEGER :: NBS
INTEGER :: FSHAPE, DSHAPE
REAL :: SPPARM(4)
TYPE(BSDAT), POINTER :: NEXTBS
END TYPE BSDAT
TYPE(BSDAT), SAVE, TARGET :: FBS
#ifdef SWAN_MODEL
! LOGICAL, POINTER, SAVE :: LBGP ! = .FALSE.
LOGICAL, POINTER, SAVE :: LBGP_G(:)
!
! This is called a linked list.
!
TYPE BGPDAT
INTEGER :: BGP(6)
TYPE(BGPDAT), POINTER :: NEXTBGP
END TYPE BGPDAT
TYPE BGPDAT_Z
TYPE(BGPDAT), POINTER :: CURBGP
TYPE(BGPDAT), POINTER :: FBGP
TYPE(BGPDAT), POINTER :: CUBGP
TYPE(BGPDAT), POINTER :: BGPTMP
END TYPE BGPDAT_Z
TYPE (BGPDAT_Z), ALLOCATABLE :: BGPDATZ_MOD(:)
#else
TYPE BGPDAT
INTEGER :: BGP(6)
TYPE(BGPDAT), POINTER :: NEXTBGP
END TYPE BGPDAT
TYPE(BGPDAT), SAVE, TARGET :: FBGP
TYPE(BGPDAT), SAVE, POINTER :: CUBGP
#endif
! 8. Subroutines and functions used
!
! ---
!
! 9. Subroutines and functions calling
!
! ---
!
! 10. Error messages
!
! ---
!
! 11. Remarks
!
! ---
!
! 12. Structure
!
! ---
!
! 13. Source text
!
END MODULE M_BNDSPEC
MODULE M_OBSTA
!
!
! --|-----------------------------------------------------------|--
! | Delft University of Technology |
! | Faculty of Civil Engineering and Geosciences |
! | Environmental Fluid Mechanics Section |
! | P.O. Box 5048, 2600 GA Delft, The Netherlands |
! | |
! | Programmer: Marcel Zijlema |
! --|-----------------------------------------------------------|--
!
!
! SWAN (Simulating WAves Nearshore); a third generation wave model
! Copyright (C) 1993-2017 Delft University of Technology
!
! This program 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 2 of
! the License, or (at your option) any later version.
!
! This program 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.
!
! A copy of the GNU General Public License is available at
! http://www.gnu.org/copyleft/gpl.html#SEC3
! or by writing to the Free Software Foundation, Inc.,
! 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
!
!
! 0. Authors
!
! 40.31: Marcel Zijlema
!
! 1. Updates
!
! Nov. 03: New Module
!
! 2. Purpose
!
! Contains data with respect to obstacles
!
! 3. Method
!
! ---
!
! 4. Modules used
!
! ---
!
IMPLICIT NONE
!
! 5. Argument variables
!
! ---
!
! 6. Parameter variables
!
! ---
!
! 7. Local variables
!
! FOBSTAC : first obstacle in list of obstacles
! NCRPTS : number of corner points in obstacle
! NEXTOBST: pointer to next obstacle in list
! RFCOEF : reflection coefficients
! RFTYP1 : reflection type: standard (REFL)
! RFTYP2 : reflection type: diffusive/specular (RDIFF/RSPEC)
! RFTYP3 : reflection type: frequency-dependent (RFD)
! TRCF1D : frequency dependent transmission coefficients
! TRCF2D : frequency and direction dependent transmission coefficients
! TRCOEF : transmission coefficients
! TRTYPE : transmission type
! XCRP : x-coordinate of corner point
! YCRP : y-coordinate of corner point
TYPE OBSTDAT
INTEGER :: TRTYPE
REAL :: TRCOEF(3)
REAL, POINTER :: TRCF1D(:), TRCF2D(:,:)
INTEGER :: RFTYP1, RFTYP2, RFTYP3
REAL :: RFCOEF(6)
INTEGER :: NCRPTS
REAL, POINTER :: XCRP(:), YCRP(:)
TYPE(OBSTDAT), POINTER :: NEXTOBST
END TYPE OBSTDAT
TYPE(OBSTDAT), SAVE, TARGET :: FOBSTAC
! 8. Subroutines and functions used
!
! ---
!
! 9. Subroutines and functions calling
!
! ---
!
! 10. Error messages
!
! ---
!
! 11. Remarks
!
! ---
!
! 12. Structure
!
! ---
!
! 13. Source text
!
END MODULE M_OBSTA
MODULE M_GENARR
!
!
! --|-----------------------------------------------------------|--
! | Delft University of Technology |
! | Faculty of Civil Engineering and Geosciences |
! | Environmental Fluid Mechanics Section |
! | P.O. Box 5048, 2600 GA Delft, The Netherlands |
! | |
! | Programmer: Marcel Zijlema |
! --|-----------------------------------------------------------|--
!
!
! SWAN (Simulating WAves Nearshore); a third generation wave model
! Copyright (C) 1993-2017 Delft University of Technology
!
! This program 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 2 of
! the License, or (at your option) any later version.
!
! This program 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.
!
! A copy of the GNU General Public License is available at
! http://www.gnu.org/copyleft/gpl.html#SEC3
! or by writing to the Free Software Foundation, Inc.,
! 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
!
!
! 0. Authors
!
! 40.31: Marcel Zijlema
!
! 1. Updates
!
! Oct. 03: New Module
!
! 2. Purpose
!
! Create several allocatable arrays for SWAN computation
!
! 3. Method
!
! The following arrays will be created:
!
! KGRPNT, KGRBND
! XYTST
! AC2
! XCGRID, YCGRID
! SPCSIG, SPCDIR
! DEPTH , FRIC
! UXB , UYB
! WXI , WYI
! WLEVL , ASTDF
! NPLAF, TURBF
! MUDLF
!
! 4. Modules used
!
! ---
!
IMPLICIT NONE
!
! 5. Argument variables
!
! ---
!
! 6. Parameter variables
!
! ---
!
! 7. Local variables
!
! AC2 : Contains action density at present time step
! ASTDF : input field of air-sea temperature difference
! DEPTH : input field of depth
! FRIC : input field of friction
! KGRBND: array containing all boundary points
! (+ 2 extra zeros as area separator for all separated areas)
! KGRPNT: array containing indirect addresses for grid points
! LAYH : layer thickness for vegetation model
! MUDLF : input field containing fluid mud layer
! NPLAF : input field containing number of plants per square meter
! SPCDIR: (*,1); spectral directions (radians)
! (*,2); cosine of spectral directions
! (*,3); sine of spectral directions
! (*,4); cosine^2 of spectral directions
! (*,5); cosine*sine of spectral directions
! (*,6); sine^2 of spectral directions
! SPCSIG: Relative frequencies in computational domain in sigma-space
! TURBF : input field containing turbulent viscosity
! UXB : input field of contravariant U-velocity
! UYB : input field of contravariant V-velocity
! VEGDIL: vegetation diameter for each layer and grid point
! VEGDRL: drag coefficient for each layer and grid point
! VEGNSL: number of plants / m2 for each layer and grid point
! WLEVL : input field of water level
! WXI : input field of wind U-velocity (contravariant)
! WYI : input field of wind V-velocity (contravariant)
! XCGRID: Coordinates of computational grid in x-direction
! XYTST : Grid point indices of test points
! YCGRID: Coordinates of computational grid in y-direction
#ifdef SWAN_MODEL
TYPE T_M_GENARR_MOD
INTEGER, POINTER :: KGRPNT_G(:,:)
INTEGER, POINTER :: KGRBND_G(:)
INTEGER, POINTER :: XYTST_G(:)
REAL, POINTER :: AC2_G(:,:,:)
REAL, POINTER :: XCGRID_G(:,:)
REAL, POINTER :: YCGRID_G(:,:)
REAL, POINTER :: SPCSIG_G(:)
REAL, POINTER :: SPCDIR_G(:,:)
REAL, POINTER :: DEPTH_G(:)
REAL, POINTER :: FRIC_G(:)
REAL, POINTER :: UXB_G(:)
REAL, POINTER :: UYB_G(:)
REAL, POINTER :: WXI_G(:)
REAL, POINTER :: WYI_G(:)
REAL, POINTER :: WLEVL_G(:)
REAL, POINTER :: ASTDF_G(:)
REAL, POINTER :: NPLAF_G(:)
REAL, POINTER :: TURBF_G(:)
REAL, POINTER :: MUDLF_G(:)
# ifdef WRF_COUPLING
REAL, POINTER :: CosAngler_G(:)
REAL, POINTER :: SinAngler_G(:)
# endif
END TYPE T_M_GENARR_MOD
TYPE (T_M_GENARR_MOD), ALLOCATABLE :: M_GENARR_MOD(:)
INTEGER, POINTER :: KGRPNT(:,:)
INTEGER, POINTER :: KGRBND(:)
INTEGER, POINTER :: XYTST(:)
REAL , POINTER :: AC2(:,:,:)
REAL , POINTER :: XCGRID(:,:)
REAL , POINTER :: YCGRID(:,:)
REAL , POINTER :: SPCSIG(:)
REAL , POINTER :: SPCDIR(:,:)
REAL , POINTER :: DEPTH(:)
REAL , POINTER :: FRIC(:)
REAL , POINTER :: UXB(:)
REAL , POINTER :: UYB(:)
REAL , POINTER :: WXI(:)
REAL , POINTER :: WYI(:)
REAL , POINTER :: WLEVL(:)
REAL , POINTER :: ASTDF(:)
REAL , POINTER :: NPLAF(:)
REAL , POINTER :: TURBF(:)
REAL , POINTER :: MUDLF(:)
# ifdef WRF_COUPLING
REAL , POINTER :: CosAngler(:)
REAL , POINTER :: SinAngler(:)
# endif
#else
INTEGER, SAVE, ALLOCATABLE :: KGRPNT(:,:), KGRBND(:)
INTEGER, SAVE, ALLOCATABLE :: XYTST(:)
REAL , SAVE, ALLOCATABLE :: AC2(:,:,:)
REAL , SAVE, ALLOCATABLE :: XCGRID(:,:), YCGRID(:,:)
REAL , SAVE, ALLOCATABLE :: SPCSIG(:) , SPCDIR(:,:)
REAL , SAVE, ALLOCATABLE :: DEPTH(:) , FRIC(:)
REAL , SAVE, ALLOCATABLE :: UXB(:) , UYB(:)
REAL , SAVE, ALLOCATABLE :: WXI(:) , WYI(:)
REAL , SAVE, ALLOCATABLE :: WLEVL(:) , ASTDF(:)
REAL , SAVE, ALLOCATABLE :: MUDLF(:)
REAL , SAVE, ALLOCATABLE :: NPLAF(:) , TURBF(:)
#endif
REAL , SAVE, ALLOCATABLE :: LAYH(:), VEGDIL(:), VEGDRL(:),
& VEGNSL(:)
!ESMF!
!ESMF! added to save Sin exponential growth term for coupling
!ESMF LOGICAL, SAVE :: SAVE_SINBAC
!ESMF REAL , SAVE, ALLOCATABLE :: SINBAC(:,:,:)
! 8. Subroutines and functions used
!
! ---
!
! 9. Subroutines and functions calling
!
! ---
!
! 10. Error messages
!
! ---
!
! 11. Remarks
!
! ---
!
! 12. Structure
!
! ---
!
! 13. Source text
!
END MODULE M_GENARR
MODULE M_PARALL
!
!
! --|-----------------------------------------------------------|--
! | Delft University of Technology |
! | Faculty of Civil Engineering and Geosciences |
! | Environmental Fluid Mechanics Section |
! | P.O. Box 5048, 2600 GA Delft, The Netherlands |
! | |
! | Programmer: Marcel Zijlema |
! --|-----------------------------------------------------------|--
!
!
! SWAN (Simulating WAves Nearshore); a third generation wave model
! Copyright (C) 1993-2017 Delft University of Technology
!
! This program 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 2 of
! the License, or (at your option) any later version.
!
! This program 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.
!
! A copy of the GNU General Public License is available at
! http://www.gnu.org/copyleft/gpl.html#SEC3
! or by writing to the Free Software Foundation, Inc.,
! 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
!
!
! 0. Authors
!
! 40.31: Marcel Zijlema
! 40.41: Marcel Zijlema
!
! 1. Updates
!
! Dec. 03: New Module
! Jul. 04: introduction logicals
!
! 2. Purpose
!
! Contains data with respect to parallel process
! based on distributed-memory apprach using MPI
!
! 3. Method
!
! ---
!
! 4. Modules used
!
! ---
!
IMPLICIT NONE
!
! 5. Argument variables
!
! ---
!
! 6. Parameter variables
!
!JAC! IBLACK : integer used to colour subdomains 'black' for
!JAC! determining sequence of sweeps (=4,1,2,3)
!JAC! IGREEN : integer used to colour subdomains 'green' for
!JAC! determining sequence of sweeps (=3,4,1,2)
! IHALOX : width of halo area in x-direction
! IHALOY : width of halo area in y-direction
!JAC! IRED : integer used to colour subdomains 'red' for
!JAC! determining sequence of sweeps (=1,2,3,4)
!JAC! IYELOW : integer used to colour subdomains 'yellow' for
!JAC! determining sequence of sweeps (=2,3,4,1)
! MASTER : rank of master process
!
INTEGER MASTER
!JAC INTEGER IRED, IYELOW, IGREEN, IBLACK
INTEGER IHALOX, IHALOY
PARAMETER (MASTER = 1
!JAC & ,IRED = 1, IYELOW = 2,
!JAC & IGREEN = 3, IBLACK = 4,
!JAC & IHALOX = 1, IHALOY = 1
& ,IHALOX = 3, IHALOY = 3
& )
!
! 7. Local variables
!
! *** variables for parallel process with MPI:
!
! IAMMASTER: indicate whether this CPU is master or not
! INODE : rank of present node
! NPROC : number of nodes
! PARLL : flag to denote run as parallel (.TRUE.) or not (.FALSE.)
! SWINT : MPI datatype for integers
! SWMAX : MPI collective maximum operation
! SWMIN : MPI collective minimum operation
! SWREAL : MPI datatype for reals
! SWSUM : MPI collective summation
!
INTEGER INODE, NPROC
INTEGER SWINT, SWREAL
INTEGER SWMAX, SWMIN, SWSUM
LOGICAL IAMMASTER, PARLL
!
!COH! *** variables for coupling with COHERENS
!COH!
!COH! COMM : communicator for the wave subgroup
!COH!
!COH INTEGER COMM
!COH!
! *** information related to global domain and subdomains
!
!JAC! IBCOL : integer indicating the color of own subdomain
! IBLKAD : administration array for subdomain interfaces
! contents:
! pos. 1 number of neighbouring subdomains
! =m
! pos. 3*i-1 number of i-th neighbour
! pos. 3*i position of i-th neighbour with
! respect to present subdomain
! pos. 3*i+1 pointer of i-th neighbour in
! last part of this array
! pos. 3*m+2 number of overlapping unknowns
! on subdomain interface
! pos. 3*m+3 ... 3*m+2+n position of unknown in array
! to be sent to neighbour
! pos. 3*m+3+n ... 3*m+2*n+2 position of unknown in array
! to be received from neighbour
! IWEIG : weights to determine load per part
! KGRBGL : array containing all boundary points in global domain
! (+ 2 extra zeros as area separator for all separated areas)
! KGRPGL : indirect addressing for grid points in global domain
! =1: not active point
! >1: active point
! LENSPO : format length for spectral output
! LMXF : logical indicating whether first x-point of subdomain equals
! first x-point of global domain (=.TRUE.) or not (=.FALSE.)
! LMXL : logical indicating whether last x-point of subdomain equals
! last x-point of global domain (=.TRUE.) or not (=.FALSE.)
! LMYF : logical indicating whether first y-point of subdomain equals
! first y-point of global domain (=.TRUE.) or not (=.FALSE.)
! LMYL : logical indicating whether last y-point of subdomain equals
! last y-point of global domain (=.TRUE.) or not (=.FALSE.)
! MCGRDGL : number of wet grid points in global computational grid
! MXCGL : number of grid points in x-direction in global
! computational grid
! MXF : first index w.r.t. global grid in x-direction
! MXL : last index w.r.t. global grid in x-direction
! MYCGL : number of grid points in y-direction in global
! computational grid
! MYF : first index w.r.t. global grid in y-direction
! MYL : last index w.r.t. global grid in y-direction
! NBGGL : number of grid points for which boundary condition holds
! in global domain
! NGRBGL : number of boundary points in global domain
! XCLMAX : maximum x-coordinate in subdomain
! XCLMIN : minimum x-coordinate in subdomain
! YCLMAX : maximum y-coordinate in subdomain
! YCLMIN : minimum y-coordinate in subdomain
! XGRDGL : x-coordinate of computational grid in global domain
! YGRDGL : y-coordinate of computational grid in global domain
!
!JAC INTEGER IBCOL
#ifdef SWAN_MODEL
! Mine
INTEGER, SAVE :: WAV_COMM_WORLD
INTEGER :: Numspecl
INTEGER :: NBGRPT_last
INTEGER, SAVE, ALLOCATABLE :: BGRIDPG(:,:)
# ifdef NESTING
INTEGER, SAVE, ALLOCATABLE :: Numspec(:)
INTEGER, SAVE, ALLOCATABLE :: ac2size(:)
INTEGER, SAVE, ALLOCATABLE :: sstp(:)
INTEGER, SAVE, ALLOCATABLE :: snew(:)
# endif
!
INTEGER, POINTER :: MCGRDGL
INTEGER, SAVE, TARGET, ALLOCATABLE :: MCGRDGL_G(:)
INTEGER, POINTER :: NGRBGL
INTEGER, SAVE, TARGET, ALLOCATABLE :: NGRBGL_G(:)
INTEGER, POINTER :: MXCGL
INTEGER, SAVE, TARGET, ALLOCATABLE :: MXCGL_G(:)
INTEGER, POINTER :: MYCGL
INTEGER, SAVE, TARGET, ALLOCATABLE :: MYCGL_G(:)
INTEGER, POINTER :: MXF
INTEGER, SAVE, TARGET, ALLOCATABLE :: MXF_G(:)
INTEGER, POINTER :: MXL
INTEGER, SAVE, TARGET, ALLOCATABLE :: MXL_G(:)
INTEGER, POINTER :: MYF
INTEGER, SAVE, TARGET, ALLOCATABLE :: MYF_G(:)
INTEGER, POINTER :: MYL
INTEGER, SAVE, TARGET, ALLOCATABLE :: MYL_G(:)
INTEGER, POINTER :: NBGGL
INTEGER, SAVE, TARGET, ALLOCATABLE :: NBGGL_G(:)
TYPE T_PARALL_MOD
INTEGER, POINTER :: IBLKAD_G(:)
INTEGER, POINTER :: KGRPGL_G(:,:)
INTEGER, POINTER :: KGRBGL_G(:)
REAL, POINTER :: XGRDGL_G(:,:)
REAL, POINTER :: YGRDGL_G(:,:)
END TYPE T_PARALL_MOD
TYPE (T_PARALL_MOD), ALLOCATABLE :: PARALL_MOD(:)
INTEGER, POINTER :: IBLKAD(:)
INTEGER, POINTER :: KGRPGL(:,:)
INTEGER, POINTER :: KGRBGL(:)
REAL, POINTER :: XGRDGL(:,:)
REAL, POINTER :: YGRDGL(:,:)
REAL, POINTER :: XCLMAX
REAL, SAVE, TARGET, ALLOCATABLE :: XCLMAX_G(:)
REAL, POINTER :: YCLMAX
REAL, SAVE, TARGET, ALLOCATABLE :: YCLMAX_G(:)
REAL, POINTER :: XCLMIN
REAL, SAVE, TARGET, ALLOCATABLE :: XCLMIN_G(:)
REAL, POINTER :: YCLMIN
REAL, SAVE, TARGET, ALLOCATABLE :: YCLMIN_G(:)
LOGICAL, POINTER :: LMXF
LOGICAL, SAVE, TARGET, ALLOCATABLE :: LMXF_G(:)
LOGICAL, POINTER :: LMXL
LOGICAL, SAVE, TARGET, ALLOCATABLE :: LMXL_G(:)
LOGICAL, POINTER :: LMYF
LOGICAL, SAVE, TARGET, ALLOCATABLE :: LMYF_G(:)
LOGICAL, POINTER :: LMYL
LOGICAL, SAVE, TARGET, ALLOCATABLE :: LMYL_G(:)
!Orig
INTEGER :: LENSPO = 1000
INTEGER, SAVE, ALLOCATABLE :: IWEIG(:)
#else
INTEGER MCGRDGL, MXCGL, MYCGL
INTEGER MXF, MXL, MYF, MYL
INTEGER NGRBGL, NBGGL
REAL XCLMAX, XCLMIN, YCLMAX, YCLMIN
INTEGER :: LENSPO = 1000
LOGICAL LMXF, LMXL, LMYF, LMYL
INTEGER, SAVE, ALLOCATABLE :: IBLKAD(:)
INTEGER, SAVE, ALLOCATABLE :: IWEIG(:)
INTEGER, SAVE, ALLOCATABLE :: KGRPGL(:,:), KGRBGL(:)
REAL , SAVE, ALLOCATABLE :: XGRDGL(:,:), YGRDGL(:,:)
#endif
!
! 8. Subroutines and functions used
!
! ---
!
! 9. Subroutines and functions calling
!
! ---
!
! 10. Error messages
!
! ---
!
! 11. Remarks
!
! ---
!
! 12. Structure
!
! ---
!
! 13. Source text
!
END MODULE M_PARALL
#ifdef SWAN_MODEL
MODULE M_COUPLING
INTEGER, SAVE :: Gridnum
INTEGER, SAVE :: NUM_SGRIDS
integer, allocatable :: iics(:)
END MODULE M_COUPLING
# ifdef NESTING
MODULE M_SREFINED
TYPE T_REFINED_MOD
REAL, POINTER :: BOUNDHINDX_G(:)
REAL, POINTER :: BOUNDHINDY_G(:)
INTEGER, POINTER :: BOUNDHINDPF_G(:)
INTEGER, POINTER :: BOUNDHINDPR_G(:)
INTEGER, POINTER :: BOUNDHINDIX_G(:)
INTEGER, POINTER :: BOUNDHINDIY_G(:)
END TYPE T_REFINED_MOD
TYPE (T_REFINED_MOD), ALLOCATABLE :: REFINED_MOD(:)
REAL, POINTER :: BOUNDHINDX(:)
REAL, POINTER :: BOUNDHINDY(:)
INTEGER, POINTER :: BOUNDHINDPF(:)
INTEGER, POINTER :: BOUNDHINDPR(:)
INTEGER, POINTER :: BOUNDHINDIX(:)
INTEGER, POINTER :: BOUNDHINDIY(:)
END MODULE M_SREFINED
# endif
MODULE M_CVMESH
REAL, SAVE, POINTER :: XCSAVE_G(:)
REAL, SAVE, POINTER :: YCSAVE_G(:)
INTEGER, SAVE, POINTER :: MXITNR_G(:)
REAL, SAVE, POINTER :: XCSAVE
REAL, SAVE, POINTER :: YCSAVE
INTEGER, SAVE, POINTER :: MXITNR
END MODULE M_CVMESH
MODULE M_PARALL2
!
!
! --|-----------------------------------------------------------|--
! | Delft University of Technology |
! | Faculty of Civil Engineering and Geosciences |
! | Environmental Fluid Mechanics Section |
! | P.O. Box 5048, 2600 GA Delft, The Netherlands |
! | |
! | Programmer: M. Zijlema |
! --|-----------------------------------------------------------|--
!
!
! SWAN (Simulating WAves Nearshore); a third generation wave model
! Copyright (C) 2007 Delft University of Technology
!
! This program 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 2 of
! the License, or (at your option) any later version.
!
! This program 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.
!
! A copy of the GNU General Public License is available at
! http://www.gnu.org/copyleft/gpl.html#SEC3
! or by writing to the Free Software Foundation, Inc.,
! 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
!
!
! 0. Authors
!
! John C. Warner
!
! 1. Updates
!
! 2. Purpose
!
! Contains data with respect to parallel process
! based on distributed-memory apprach using MPI
!
! 3. Method
!
! ---
!
! 4. Modules used
!
! ---
!
IMPLICIT NONE
!
! 5. Argument variables
!
! ---
!
! 6. Parameter variables
!
!
#ifdef SWAN_MODEL
! INTEGER, SAVE :: Numgrids
TYPE T_PARALL2_MOD
REAL, POINTER :: BSPECS_G(:,:,:,:)
INTEGER, POINTER :: CROSS_G(:)
INTEGER, POINTER :: BGRIDP_G(:)
REAL, POINTER :: AC1_G(:,:,:)
REAL, POINTER :: COMPDA_G(:,:)
REAL, POINTER :: OURQT_G(:)
REAL, POINTER :: AC2_parent_G(:,:)
END TYPE T_PARALL2_MOD
TYPE (T_PARALL2_MOD), ALLOCATABLE :: PARALL2_MOD(:)
REAL, POINTER :: BSPECS(:,:,:,:)
INTEGER, POINTER :: CROSS(:)
INTEGER, POINTER :: BGRIDP(:)
REAL, POINTER :: AC1(:,:,:)
REAL, POINTER :: COMPDA(:,:)
REAL, POINTER :: OURQT(:)
REAL, POINTER :: AC2_parent(:,:)
#else
REAL, SAVE, ALLOCATABLE :: BSPECS(:,:,:,:)
INTEGER, SAVE, ALLOCATABLE :: CROSS(:)
INTEGER, SAVE, ALLOCATABLE :: BGRIDP(:)
REAL, SAVE, ALLOCATABLE :: AC1(:,:,:)
REAL, SAVE, ALLOCATABLE :: COMPDA(:,:)
REAL, SAVE, ALLOCATABLE :: OURQT(:)
#endif
END MODULE M_PARALL2
#endif
MODULE M_DIFFR
!
!
! --|-----------------------------------------------------------|--
! | Delft University of Technology |
! | Faculty of Civil Engineering |
! | Environmental Fluid Mechanics Section |
! | P.O. Box 5048, 2600 GA Delft, The Netherlands |
! | |
! | Programmers: The SWAN team |
! --|-----------------------------------------------------------|--
!
!
! SWAN (Simulating WAves Nearshore); a third generation wave model
! Copyright (C) 1993-2017 Delft University of Technology
!
! This program 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 2 of
! the License, or (at your option) any later version.
!
! This program 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.
!
! A copy of the GNU General Public License is available at
! http://www.gnu.org/copyleft/gpl.html#SEC3
! or by writing to the Free Software Foundation, Inc.,
! 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
!
!
! 0. Authors
!
! 40.21: Agnieszka Herman, Nico Booij
!
! 1. Updates
!
! Aug. 01: New Module
!
! 2. Purpose
!
! Contains global variables used in diffraction procedure
!
! 3. Method
!
! MODULE construct
!
! 4. Modules used
!
! ---
!
IMPLICIT NONE
!
! 5. Argument variables
!
! ---
!
! 6. Parameter variables
!
! ---
!
! 7. Local variables
!
! DIFPARAM: contains diffraction parameter (second derivative of Hs)
! DIFPARDX: derivative of DIFPARAM in x-direction
! DIFPARDY: derivative of DIFPARAM in y-direction
REAL, SAVE, ALLOCATABLE :: DIFPARAM(:)
REAL, SAVE, ALLOCATABLE :: DIFPARDX(:)
REAL, SAVE, ALLOCATABLE :: DIFPARDY(:)
! 8. Subroutines and functions used
!
! ---
!
! 9. Subroutines and functions calling
!
! DIFPAR : calculates the above arrays DIFPARAM, DIFPARDX, DIFPARDY
! SPROSD : computes propagation velocity in (x,y,theta) based on
! arrays DIFPARAM, DIFPARDX, DIFPARDY
!
! 10. Error messages
!
! ---
!
! 11. Remarks
!
! ---
!
! 12. Structure
!
! ---
!
! 13. Source text
!
END MODULE M_DIFFR
#ifdef SWAN_MODEL
MODULE M_MPI
INCLUDE 'mpif.h'
END MODULE M_MPI
#else
MODULE MPI
INCLUDE 'mpif.h'
END MODULE MPI
!PUN!/impi
!PUN!/impi MODULE MPI
!PUN!/impi INCLUDE 'mpif.h'
!PUN!/impi END MODULE MPI
#endif