scale_ocean_phy_file Module Reference

module OCEAN / Physics File More...

Functions/Subroutines
subroutine, public	ocean_phy_file_setup (OCEAN_TYPE)
	Setup. More...
subroutine, public	ocean_phy_file (OCEAN_TEMP_t, OCEAN_TEMP, OCEAN_SFLX_WH, OCEAN_SFLX_prec, OCEAN_SFLX_evap, dt)
	Slab ocean model. More...

Detailed Description

module OCEAN / Physics File

Description

ocean physics module, external file input

Author

Team SCALE

Function/Subroutine Documentation

ocean_phy_file_setup()

subroutine, public scale_ocean_phy_file::ocean_phy_file_setup

(

character(len=*), intent(in)

OCEAN_TYPE

)

Setup.

Definition at line 50 of file scale_ocean_phy_file.F90.

References scale_grid_index::ia, scale_grid_index::ie, scale_stdio::io_fid_log, scale_stdio::io_l, scale_grid_index::is, scale_grid_index::ja, scale_grid_index::je, scale_grid_index::js, scale_landuse::landuse_fact_ocean, and scale_process::prc_mpistop().

Referenced by scale_ocean_phy::ocean_phy_setup().

    use scale_process, only: &
       prc_mpistop
    use scale_landuse, only: &
       landuse_fact_ocean
    implicit none

    character(len=*), intent(in) :: ocean_type

    integer :: i, j
    integer :: ierr
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '++++++ Module[FILE] / Categ[OCEAN PHY] / Origin[SCALElib]'

    if( ocean_type /= 'FILE' ) then
       write(*,*) 'xxx wrong OCEAN_TYPE. Check!'
       call prc_mpistop
    end if

    ! judge to run slab ocean model
    allocate( is_ocn(ia,ja) )

    do j = js, je
    do i = is, ie
       if( landuse_fact_ocean(i,j) > 0.0_rp ) then
          is_ocn(i,j) = .true.
       else
          is_ocn(i,j) = .false.
       end if
    end do
    end do

    return

Here is the call graph for this function:

Here is the caller graph for this function:

ocean_phy_file()

subroutine, public scale_ocean_phy_file::ocean_phy_file	(	real(rp), dimension (ia,ja), intent(out)	OCEAN_TEMP_t,
		real(rp), dimension (ia,ja), intent(in)	OCEAN_TEMP,
		real(rp), dimension (ia,ja), intent(in)	OCEAN_SFLX_WH,
		real(rp), dimension(ia,ja), intent(in)	OCEAN_SFLX_prec,
		real(rp), dimension(ia,ja), intent(in)	OCEAN_SFLX_evap,
		real(dp), intent(in)	dt
	)

Slab ocean model.

Definition at line 95 of file scale_ocean_phy_file.F90.

References scale_grid_index::ie, scale_stdio::io_fid_log, scale_stdio::io_l, scale_grid_index::is, scale_grid_index::je, scale_grid_index::js, scale_process::prc_mpistop(), and scale_time::time_nowdaysec.

Referenced by scale_ocean_phy::ocean_phy_setup().

    use scale_grid_index
    use scale_time, only: &
         nowdaysec => time_nowdaysec
    use scale_process, only: &
         prc_mpistop
    use scale_external_input, only: &
         extin_update
    implicit none

    real(RP), intent(out) :: ocean_temp_t   (ia,ja)
    real(RP), intent(in)  :: ocean_temp     (ia,ja)
    real(RP), intent(in)  :: ocean_sflx_wh  (ia,ja)
    real(RP), intent(in)  :: ocean_sflx_prec(ia,ja)
    real(RP), intent(in)  :: ocean_sflx_evap(ia,ja)
    real(DP), intent(in)  :: dt

    real(RP) :: ocean_temp_new(ia,ja)

    logical :: error

    integer :: i, j
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*) '*** Ocean step: File'

    call extin_update( &
         ocean_temp_new, & ! (out)
         'OCEAN_TEMP',   & ! (in)
         'XY',           & ! (in)
         nowdaysec,      & ! (in)
         error           ) ! (out)
    if ( error ) then
       write(*,*) 'xxx Failed to read "OCEAN_TEMP" data from file!'
       write(*,*) 'xxx set namelist EXTITEM, for example'
       write(*,*) ' &EXTITEM'
       write(*,*) '  basename   = "filename",'
       write(*,*) '  varname    = "OCEAN_TEMP",'
       write(*,*) ' /'
       call prc_mpistop
    end if

    do j = js, je
    do i = is, ie
       if( is_ocn(i,j) ) then
          ocean_temp_t(i,j) = ( ocean_temp_new(i,j) - ocean_temp(i,j) ) / dt
       else
          ocean_temp_t(i,j) = 0.0_rp
       endif
    enddo
    enddo

    return

Here is the call graph for this function:

Here is the caller graph for this function: