source: palm/trunk/SOURCE/gust_mod.f90 @ 2895

!> @file gust_mod.f90
!------------------------------------------------------------------------------!
! This file is part of PALM.
!
! PALM 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.
!
! PALM 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
! PALM. If not, see <http://www.gnu.org/licenses/>.
!
! Copyright 1997-2018 Leibniz Universitaet Hannover
!------------------------------------------------------------------------------!
!
! Current revisions:
! -----------------
! 
! 
! Former revisions:
! -----------------
!
! Initial interface definition
!
! 
!
! Description:
! ------------
!> Gust model.
!>
!> @todo This is just a dummy module. The actual module ist not released yet.
!------------------------------------------------------------------------------!
 MODULE gust_mod

    USE indices,                                                               &
        ONLY:  nxl, nxlg, nxr, nxrg, nys, nysg, nyn, nyng, nzb, nzt

    USE kinds

    IMPLICIT NONE

    LOGICAL  ::  gust_module_enabled = .FALSE.       !< switch, if the entire module is used at all

    SAVE

    PRIVATE

!
!-- Public functions
    PUBLIC &
       gust_parin, &
       gust_check_parameters, &
       gust_check_data_output_pr, &
       gust_check_data_output, &
       gust_init_arrays, &
       gust_init, &
       gust_define_netcdf_grid, &
       gust_header, &
       gust_actions, &
       gust_swap_timelevel, &
       gust_3d_data_averaging, &
       gust_data_output_2d, &
       gust_data_output_3d, &
       gust_statistics, &
       gust_read_restart_data, &
       gust_write_restart_data
!
!-- Public parameters, constants and initial values
    PUBLIC &
       gust_module_enabled


    INTERFACE gust_parin
       MODULE PROCEDURE gust_parin
    END INTERFACE gust_parin

    INTERFACE gust_check_parameters
       MODULE PROCEDURE gust_check_parameters
    END INTERFACE gust_check_parameters

    INTERFACE gust_check_data_output_pr
       MODULE PROCEDURE gust_check_data_output_pr
    END INTERFACE gust_check_data_output_pr

    INTERFACE gust_check_data_output
       MODULE PROCEDURE gust_check_data_output
    END INTERFACE gust_check_data_output

    INTERFACE gust_init_arrays
       MODULE PROCEDURE gust_init_arrays
    END INTERFACE gust_init_arrays

    INTERFACE gust_init
       MODULE PROCEDURE gust_init
    END INTERFACE gust_init

    INTERFACE gust_define_netcdf_grid
       MODULE PROCEDURE gust_define_netcdf_grid
    END INTERFACE gust_define_netcdf_grid

    INTERFACE gust_header
       MODULE PROCEDURE gust_header
    END INTERFACE gust_header

    INTERFACE gust_actions
       MODULE PROCEDURE gust_actions
       MODULE PROCEDURE gust_actions_ij
    END INTERFACE gust_actions

    INTERFACE gust_swap_timelevel
       MODULE PROCEDURE gust_swap_timelevel
    END INTERFACE gust_swap_timelevel

    INTERFACE gust_3d_data_averaging
       MODULE PROCEDURE gust_3d_data_averaging
    END INTERFACE gust_3d_data_averaging

    INTERFACE gust_data_output_2d
       MODULE PROCEDURE gust_data_output_2d
    END INTERFACE gust_data_output_2d

    INTERFACE gust_data_output_3d
       MODULE PROCEDURE gust_data_output_3d
    END INTERFACE gust_data_output_3d

    INTERFACE gust_statistics
       MODULE PROCEDURE gust_statistics
    END INTERFACE gust_statistics

    INTERFACE gust_read_restart_data
       MODULE PROCEDURE gust_read_restart_data
    END INTERFACE gust_read_restart_data

    INTERFACE gust_write_restart_data
       MODULE PROCEDURE gust_write_restart_data
    END INTERFACE gust_write_restart_data

 CONTAINS


!------------------------------------------------------------------------------!
! Description:
! ------------
!> Parin for &gust_par for gust module
!------------------------------------------------------------------------------!
    SUBROUTINE gust_parin


       IMPLICIT NONE

       CHARACTER (LEN=80)  ::  line  !< dummy string that contains the current line of the parameter file

       NAMELIST /gust_par/  &
          gust_module_enabled

       line = ' '
!
!--    Try to find gust module package
       REWIND ( 11 )
       line = ' '
       DO   WHILE ( INDEX( line, '&gust_par' ) == 0 )
          READ ( 11, '(A)', END=10 )  line
       ENDDO
       BACKSPACE ( 11 )
!
!--    Read user-defined namelist
       READ ( 11, gust_par )
!
!--    Set flag that indicates that the gust module is switched on
       gust_module_enabled = .TRUE.

10     CONTINUE


    END SUBROUTINE gust_parin


!------------------------------------------------------------------------------!
! Description:
! ------------
!> Check parameters routine for gust module
!------------------------------------------------------------------------------!
    SUBROUTINE gust_check_parameters


       IMPLICIT NONE


    END SUBROUTINE gust_check_parameters


!------------------------------------------------------------------------------!
! Description:
! ------------
!> Check data output of profiles for gust module
!------------------------------------------------------------------------------!
    SUBROUTINE gust_check_data_output_pr( variable, var_count, unit, dopr_unit )


       IMPLICIT NONE

       CHARACTER (LEN=*) ::  unit      !<
       CHARACTER (LEN=*) ::  variable  !<
       CHARACTER (LEN=*) ::  dopr_unit !< local value of dopr_unit

       INTEGER(iwp) ::  var_count      !<


    END SUBROUTINE gust_check_data_output_pr

!------------------------------------------------------------------------------!
! Description:
! ------------
!> Check data output for gust module
!------------------------------------------------------------------------------!
    SUBROUTINE gust_check_data_output( var, unit )


       IMPLICIT NONE

       CHARACTER (LEN=*) ::  unit  !<
       CHARACTER (LEN=*) ::  var   !<


    END SUBROUTINE gust_check_data_output


!------------------------------------------------------------------------------!
! Description:
! ------------
!> Allocate gust module arrays and define pointers
!------------------------------------------------------------------------------!
    SUBROUTINE gust_init_arrays


       IMPLICIT NONE


    END SUBROUTINE gust_init_arrays


!------------------------------------------------------------------------------!
! Description:
! ------------
!> Initialization of the gust module
!------------------------------------------------------------------------------!
    SUBROUTINE gust_init( dots_label, dots_unit, dots_num, dots_max )


       IMPLICIT NONE

       INTEGER(iwp) ::  dots_num
       INTEGER(iwp) ::  dots_max
       CHARACTER (LEN=13), DIMENSION(dots_max) :: dots_unit
       CHARACTER (LEN=13), DIMENSION(dots_max) :: dots_label


    END SUBROUTINE gust_init


!------------------------------------------------------------------------------!
!
! Description:
! ------------
!> Subroutine defining appropriate grid for netcdf variables.
!> It is called out from subroutine netcdf.
!------------------------------------------------------------------------------!
    SUBROUTINE gust_define_netcdf_grid( var, found, grid_x, grid_y, grid_z )


       IMPLICIT NONE

       CHARACTER (LEN=*), INTENT(IN)  ::  var         !<
       LOGICAL, INTENT(IN)           ::  found       !<
       CHARACTER (LEN=*), INTENT(IN) ::  grid_x      !<
       CHARACTER (LEN=*), INTENT(IN) ::  grid_y      !<
       CHARACTER (LEN=*), INTENT(IN) ::  grid_z      !<


    END SUBROUTINE gust_define_netcdf_grid


!------------------------------------------------------------------------------!
! Description:
! ------------
!> Header output for gust module
!------------------------------------------------------------------------------!
    SUBROUTINE gust_header ( io )


       IMPLICIT NONE

       INTEGER(iwp), INTENT(IN) ::  io  !< Unit of the output file


    END SUBROUTINE gust_header


!------------------------------------------------------------------------------!
! Description:
! ------------
!> Call for all grid points
!------------------------------------------------------------------------------!
    SUBROUTINE gust_actions( location )


       IMPLICIT NONE

       CHARACTER (LEN=*) ::  location !<


    END SUBROUTINE gust_actions


!------------------------------------------------------------------------------!
! Description:
! ------------
!> Call for grid point i,j
!------------------------------------------------------------------------------!
    SUBROUTINE gust_actions_ij( i, j, location )


       IMPLICIT NONE

       CHARACTER (LEN=*) ::  location

       INTEGER(iwp) ::  i
       INTEGER(iwp) ::  j


    END SUBROUTINE gust_actions_ij


!------------------------------------------------------------------------------!
! Description:
! ------------
!> Swapping of timelevels
!------------------------------------------------------------------------------!
    SUBROUTINE gust_swap_timelevel ( mod_count )


       IMPLICIT NONE

       INTEGER, INTENT(IN) :: mod_count


    END SUBROUTINE gust_swap_timelevel


!------------------------------------------------------------------------------!
!
! Description:
! ------------
!> Subroutine for averaging 3D data
!------------------------------------------------------------------------------!
    SUBROUTINE gust_3d_data_averaging( mode, variable )


       IMPLICIT NONE

       CHARACTER (LEN=*) ::  mode    !<
       CHARACTER (LEN=*) :: variable !<


    END SUBROUTINE gust_3d_data_averaging

!------------------------------------------------------------------------------!
!
! Description:
! ------------
!> Subroutine defining 2D output variables
!------------------------------------------------------------------------------!
    SUBROUTINE gust_data_output_2d( av, variable, found, grid, local_pf,       &
                                    two_d, nzb_do, nzt_do )


       IMPLICIT NONE

       CHARACTER (LEN=*) ::  grid     !<
       CHARACTER (LEN=*) ::  variable !<

       INTEGER(iwp) ::  av     !< flag to control data output of instantaneous or time-averaged data
       INTEGER(iwp) ::  nzb_do !< lower limit of the data output (usually 0)
       INTEGER(iwp) ::  nzt_do !< vertical upper limit of the data output (usually nz_do3d)

       LOGICAL      ::  found !<
       LOGICAL      ::  two_d !< flag parameter that indicates 2D variables (horizontal cross sections)

       REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb:nzt+1) ::  local_pf !<


    END SUBROUTINE gust_data_output_2d


!------------------------------------------------------------------------------!
!
! Description:
! ------------
!> Subroutine defining 3D output variables
!------------------------------------------------------------------------------!
    SUBROUTINE gust_data_output_3d( av, variable, found, local_pf, nzb_do, nzt_do )


       IMPLICIT NONE

       CHARACTER (LEN=*) ::  variable !<

       INTEGER(iwp) ::  av    !<
       INTEGER(iwp) ::  nzb_do !< lower limit of the data output (usually 0)
       INTEGER(iwp) ::  nzt_do !< vertical upper limit of the data output (usually nz_do3d)

       LOGICAL      ::  found !<

       REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf !<


    END SUBROUTINE gust_data_output_3d


!------------------------------------------------------------------------------!
! Description:
! ------------
!> This routine computes profile and timeseries data for the gust module.
!------------------------------------------------------------------------------!
    SUBROUTINE gust_statistics( mode, sr, tn, dots_max )


       IMPLICIT NONE

       CHARACTER (LEN=*) ::  mode   !<

       INTEGER(iwp) ::  sr   !<
       INTEGER(iwp) ::  tn   !<
       INTEGER(iwp) ::  dots_max   !<


    END SUBROUTINE gust_statistics


!------------------------------------------------------------------------------!
! Description:
! ------------
!> This routine reads the respective restart data for the gust module.
!------------------------------------------------------------------------------!
    SUBROUTINE gust_read_restart_data


       IMPLICIT NONE


    END SUBROUTINE gust_read_restart_data


!------------------------------------------------------------------------------!
! Description:
! ------------
!> This routine writes the respective restart data for the gust module.
!------------------------------------------------------------------------------!
    SUBROUTINE gust_write_restart_data


       IMPLICIT NONE


    END SUBROUTINE gust_write_restart_data



 END MODULE gust_mod