scale_external_input_rm.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
module scale_external_input_rm
  !-----------------------------------------------------------------------------
  !
  !++ Used modules
  !
  use scale_precision
  use scale_stdio
  use scale_prof
  use scale_grid_index
  use scale_land_grid_index
  use scale_urban_grid_index
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedures
  !
  public :: extin_rm_get_dims_1d
  public :: extin_rm_get_dims_2d
  public :: extin_rm_get_dims_3d

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private procedures
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine extin_rm_get_dims_1d( &
       dim1_max, &
       dim1_S,   &
       dim1_E,   &
       varname,  &
       axistype  )
    use scale_process, only: &
       prc_mpistop
    implicit none
    integer,          intent(out) :: dim1_max
    integer,          intent(out) :: dim1_s
    integer,          intent(out) :: dim1_e
    character(len=*), intent(in)  :: varname
    character(len=*), intent(in)  :: axistype     ! axis type (Z/X/Y)

    select case ( axistype )
    case ( 'Z' )
       dim1_max = kmax
       dim1_s   = ks
       dim1_e   = ke
    case ( 'X' )
       dim1_max = imaxb
       dim1_s   = isb
       dim1_e   = ieb
    case ( 'Y' )
       dim1_max = jmaxb
       dim1_s   = jsb
       dim1_e   = jeb
    case default
       write(*,*) 'xxx [EXTIN_RM_get_dims_1D] unsupported axis type. Check! axistype:', trim(axistype), ', item:',trim(varname)
       call prc_mpistop
    end select

    return
  end subroutine extin_rm_get_dims_1d

  subroutine extin_rm_get_dims_2d( &
       dim1_max,  &
       dim1_S,    &
       dim1_E,    &
       dim2_max,  &
       dim2_S,    &
       dim2_E,    &
       transpose, &
       varname,   &
       axistype   )
    use scale_process, only: &
       prc_mpistop
    implicit none
    integer,          intent(out) :: dim1_max
    integer,          intent(out) :: dim1_s
    integer,          intent(out) :: dim1_e
    integer,          intent(out) :: dim2_max
    integer,          intent(out) :: dim2_s
    integer,          intent(out) :: dim2_e
    logical,          intent(out) :: transpose
    character(len=*), intent(in)  :: varname
    character(len=*), intent(in)  :: axistype     ! axis type (XY/XZ/ZX)

    select case ( axistype )
    case ( 'XY' )
       dim1_max = imaxb
       dim2_max = jmaxb
       dim1_s   = isb
       dim1_e   = ieb
       dim2_s   = jsb
       dim2_e   = jeb
       transpose = .false.
    case ( 'ZX' )
       dim1_max = kmax
       dim2_max = imaxb
       dim1_s   = ks
       dim1_e   = ke
       dim2_s   = isb
       dim2_e   = ieb
       transpose = .false.
    case ( 'XZ' )
       dim1_max = imaxb
       dim2_max = kmax
       dim1_s   = isb
       dim1_e   = ieb
       dim2_s   = ks
       dim2_e   = ke
       transpose = .true.
    case default
       write(*,*) 'xxx [EXTIN_RM_get_dims_2D] unsupported axis type. Check! axistype:', trim(axistype), ', item:',trim(varname)
       call prc_mpistop
    end select

    return
  end subroutine extin_rm_get_dims_2d

  subroutine extin_rm_get_dims_3d( &
       dim1_max,  &
       dim1_S,    &
       dim1_E,    &
       dim2_max,  &
       dim2_S,    &
       dim2_E,    &
       dim3_max,  &
       dim3_S,    &
       dim3_E,    &
       transpose, &
       varname,   &
       axistype   )
    use scale_process, only: &
       prc_mpistop
    implicit none
    integer,          intent(out) :: dim1_max
    integer,          intent(out) :: dim1_s
    integer,          intent(out) :: dim1_e
    integer,          intent(out) :: dim2_max
    integer,          intent(out) :: dim2_s
    integer,          intent(out) :: dim2_e
    integer,          intent(out) :: dim3_max
    integer,          intent(out) :: dim3_s
    integer,          intent(out) :: dim3_e
    logical,          intent(out) :: transpose
    character(len=*), intent(in)  :: varname
    character(len=*), intent(in)  :: axistype     ! axis type (ZXY/XYZ/LXY/XYL/UXY/XYU)

    select case ( axistype )
    case ( 'ZXY' )
       dim1_max = kmax
       dim2_max = imaxb
       dim3_max = jmaxb
       dim1_s   = ks
       dim1_e   = ke
       dim2_s   = isb
       dim2_e   = ieb
       dim3_s   = jsb
       dim3_e   = jeb
       transpose = .false.
    case ( 'XYZ' )
       dim1_max = imaxb
       dim2_max = jmaxb
       dim3_max = kmax
       dim1_s   = isb
       dim1_e   = ieb
       dim2_s   = jsb
       dim2_e   = jeb
       dim3_s   = ks
       dim3_e   = ke
       transpose = .true.
    case ( 'LXY' ) ! LAND
       dim1_max = lkmax
       dim2_max = imaxb
       dim3_max = jmaxb
       dim1_s   = lks
       dim1_e   = lke
       dim2_s   = isb
       dim2_e   = ieb
       dim3_s   = jsb
       dim3_e   = jeb
       transpose = .false.
    case ( 'XYL' ) ! LAND
       dim1_max = imaxb
       dim2_max = jmaxb
       dim3_max = lkmax
       dim1_s   = isb
       dim1_e   = ieb
       dim2_s   = jsb
       dim2_e   = jeb
       dim3_s   = lks
       dim3_e   = lke
       transpose = .true.
    case ( 'UXY' )
       dim1_max = ukmax
       dim2_max = imaxb
       dim3_max = jmaxb
       dim1_s   = uks
       dim1_e   = uke
       dim2_s   = isb
       dim2_e   = ieb
       dim3_s   = jsb
       dim3_e   = jeb
       transpose = .false.
    case ( 'XYU' )
       dim1_max = imaxb
       dim2_max = jmaxb
       dim3_max = ukmax
       dim1_s   = isb
       dim1_e   = ieb
       dim2_s   = jsb
       dim2_e   = jeb
       dim3_s   = uks
       dim3_e   = uke
       transpose = .true.
    case default
       write(*,*) 'xxx [EXTIN_RM_get_dims_3D] unsupported axis type. Check! axistype:', trim(axistype), ', item:',trim(varname)
       call prc_mpistop
    end select

    return
  end subroutine extin_rm_get_dims_3d

end module scale_external_input_rm