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

NAMELIST

PARAM_OCEAN_PHY_FILE

name	type	default value	comment
OCEAN_PHY_FILE_BASENAME	character(len=H_LONG)	''
OCEAN_PHY_FILE_ENABLE_PERIODIC_YEAR	logical	.false.
OCEAN_PHY_FILE_ENABLE_PERIODIC_MONTH	logical	.false.
OCEAN_PHY_FILE_ENABLE_PERIODIC_DAY	logical	.false.
OCEAN_PHY_FILE_STEP_FIXED	integer	0
OCEAN_PHY_FILE_OFFSET	real(RP)	0.0_RP
OCEAN_PHY_FILE_DEFVAL	real(RP)		= UNDEF
OCEAN_PHY_FILE_CHECK_COORDINATES	logical	.true.
OCEAN_PHY_FILE_STEP_LIMIT	integer	0

History Output: No history output

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_const::const_undef, scale_external_input::extin_regist(), scale_grid_index::ia, scale_grid_index::ie, scale_stdio::io_fid_conf, scale_stdio::io_fid_log, scale_stdio::io_fid_nml, scale_stdio::io_l, scale_stdio::io_nml, 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
     use scale_external_input, only: &
        extin_regist
     use scale_const, only: &
        undef => const_undef
     implicit none
 
     character(len=*), intent(in) :: OCEAN_TYPE
 
     character(len=H_LONG) :: OCEAN_PHY_FILE_basename              = ''
     logical               :: OCEAN_PHY_FILE_enable_periodic_year  = .false.
     logical               :: OCEAN_PHY_FILE_enable_periodic_month = .false.
     logical               :: OCEAN_PHY_FILE_enable_periodic_day   = .false.
     integer               :: OCEAN_PHY_FILE_step_fixed            = 0
     real(RP)              :: OCEAN_PHY_FILE_offset                = 0.0_rp
     real(RP)              :: OCEAN_PHY_FILE_defval              ! = UNDEF
     logical               :: OCEAN_PHY_FILE_check_coordinates     = .true.
     integer               :: OCEAN_PHY_FILE_step_limit            = 0
 
     namelist / param_ocean_phy_file / &
        ocean_phy_file_basename,              &
        ocean_phy_file_enable_periodic_year,  &
        ocean_phy_file_enable_periodic_month, &
        ocean_phy_file_enable_periodic_day,   &
        ocean_phy_file_step_fixed,            &
        ocean_phy_file_offset,                &
        ocean_phy_file_defval,                &
        ocean_phy_file_check_coordinates,     &
        ocean_phy_file_step_limit
 
     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
 
     ocean_phy_file_defval = undef
 
     !--- read namelist
     rewind(io_fid_conf)
     read(io_fid_conf,nml=param_ocean_phy_file,iostat=ierr)
     if( ierr < 0 ) then !--- missing
        if( io_l ) write(io_fid_log,*) '*** Not found namelist. Default used.'
     elseif( ierr > 0 ) then !--- fatal error
        write(*,*) 'xxx Not appropriate names in namelist PARAM_OCEAN_PHY_FILE. Check!'
        call prc_mpistop
     endif
     if( io_nml ) write(io_fid_nml,nml=param_ocean_phy_file)
 
     if ( ocean_phy_file_basename == '' ) then
        write(*,*) 'xxx OCEAN_PHY_FILE_basename is necessary'
        call prc_mpistop
     end if
 
     call extin_regist( ocean_phy_file_basename,              & ! [IN]
                        'OCEAN_TEMP',                         & ! [IN]
                        'XY',                                 & ! [IN]
                        ocean_phy_file_enable_periodic_year,  & ! [IN]
                        ocean_phy_file_enable_periodic_month, & ! [IN]
                        ocean_phy_file_enable_periodic_day,   & ! [IN]
                        ocean_phy_file_step_fixed,            & ! [IN]
                        ocean_phy_file_offset,                & ! [IN]
                        ocean_phy_file_defval,                & ! [IN]
                        ocean_phy_file_check_coordinates,     & ! [IN]
                        ocean_phy_file_step_limit             ) ! [IN]
 
     ! 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 150 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.