scale_ocean_grid_cartesC_real.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
#include "scalelib.h"
module scale_ocean_grid_cartesc_real
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_io
  use scale_prof
  use scale_atmos_grid_cartesc_index
  use scale_ocean_grid_cartesc_index
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: ocean_grid_cartesc_real_setup
  public :: ocean_grid_cartesc_real_finalize
  public :: ocean_grid_cartesc_real_set_areavol
 
  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  real(rp), public, allocatable :: ocean_grid_cartesc_real_area(:,:)
  real(rp), public              :: ocean_grid_cartesc_real_totarea
  real(rp), public, allocatable :: ocean_grid_cartesc_real_vol (:,:,:)
  real(rp), public              :: ocean_grid_cartesc_real_totvol
  !$acc declare create(OCEAN_GRID_CARTESC_REAL_TOTAREA,OCEAN_GRID_CARTESC_REAL_TOTVOL)
 
  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine ocean_grid_cartesc_real_setup
    implicit none
 
    ! at this moment, horizontal grid is identical to that of the atmosphere
    allocate( ocean_grid_cartesc_real_area(    oia,oja) )
    allocate( ocean_grid_cartesc_real_vol(oka,oia,oja) )
    !$acc enter data create(OCEAN_GRID_CARTESC_REAL_AREA,OCEAN_GRID_CARTESC_REAL_VOL)
 
    return
  end subroutine ocean_grid_cartesc_real_setup
 
  !-----------------------------------------------------------------------------
  subroutine ocean_grid_cartesc_real_finalize
    implicit none
 
    !$acc exit data delete(OCEAN_GRID_CARTESC_REAL_AREA,OCEAN_GRID_CARTESC_REAL_VOL)
    deallocate( ocean_grid_cartesc_real_area )
    deallocate( ocean_grid_cartesc_real_vol )
 
    return
  end subroutine ocean_grid_cartesc_real_finalize
 
  !-----------------------------------------------------------------------------
  subroutine ocean_grid_cartesc_real_set_areavol
    use scale_atmos_grid_cartesc_real, only: &
       atmos_grid_cartesc_real_area
    use scale_ocean_grid_cartesc, only: &
       ocean_grid_cartesc_cdz
    use scale_file_cartesc, only: &
       file_cartesc_set_coordinates_ocean
    use scale_landuse, only: &
       landuse_fact_ocean
    implicit none
 
    integer :: k, i, j
 
    do j = 1,   oja
    do i = 1,   oia
       ocean_grid_cartesc_real_area(i,j) = atmos_grid_cartesc_real_area(i,j) * landuse_fact_ocean(i,j)
    end do
    end do
 
    ocean_grid_cartesc_real_totarea = 0.0_rp
    do j = ojs, oje
    do i = ois, oie
       ocean_grid_cartesc_real_totarea = ocean_grid_cartesc_real_totarea + ocean_grid_cartesc_real_area(i,j)
    end do
    end do
 
    do j = 1,   oja
    do i = 1,   oia
    do k = oks, oke
       ocean_grid_cartesc_real_vol(k,i,j) = ocean_grid_cartesc_real_area(i,j) * ocean_grid_cartesc_cdz(k)
    enddo
    enddo
    enddo
 
    ocean_grid_cartesc_real_totvol = 0.0_rp
    do j = ojs, oje
    do i = ois, oie
    do k = oks, oke
       ocean_grid_cartesc_real_totvol = ocean_grid_cartesc_real_totvol + ocean_grid_cartesc_real_vol(k,i,j)
    end do
    end do
    end do
 
    call file_cartesc_set_coordinates_ocean( ocean_grid_cartesc_real_vol(:,:,:) ) ! [IN]
 
    !$acc update device(OCEAN_GRID_CARTESC_REAL_AREA,OCEAN_GRID_CARTESC_REAL_VOL)
    !$acc update device(OCEAN_GRID_CARTESC_REAL_TOTAREA,OCEAN_GRID_CARTESC_REAL_TOTVOL)
    return
  end subroutine ocean_grid_cartesc_real_set_areavol
 
end module scale_ocean_grid_cartesc_real