module REAL input More...

Functions/Subroutines
subroutine, public	realinput_atmos

subroutine, public	realinput_surface

subroutine	land_interporation (tg, strg, lst, albg, ust, albu, tg_org, strg_org, smds_org, lst_org, albg_org, ust_org, sst_org, lmask_org, lsmask_nest, topo_org, lz_org, llon_org, llat_org, LCZ, LON, LAT, ldims, odims, maskval_tg, maskval_strg, init_landwater_ratio, use_file_landwater, use_waterratio, soilwater_ds2vc_flag, elevation_collection, intrp_iter_max, ol_interp, URBAN_do)

subroutine	make_mask (gmask, data, nx, ny, landdata)

subroutine	replace_misval_const (data, maskval, frac_land)

subroutine	replace_misval_map (data, maskval, nx, ny, elem)

Variables
integer, parameter, public	iscale = 1

integer, parameter, public	iwrfarw = 2

integer, parameter, public	igrads = 4

Detailed Description

module REAL input

Description: read data from file for real atmospheric simulations

Author: Team SCALE

NAMELIST

PARAM_MKINIT_REAL_ATMOS

name	type	default value	comment
NUMBER_OF_FILES	integer	1
NUMBER_OF_TSTEPS	integer	1	num of time steps in one file
NUMBER_OF_SKIP_TSTEPS	integer	0	num of skipped first several data
SERIAL_PROC_READ	logical	.true.	read by one MPI process and broadcast
FILETYPE_ORG	character(len=H_LONG)	''
BASENAME_ORG	character(len=H_LONG)	''
BASENAME_ADD_NUM	logical	.false.
BASENAME_BOUNDARY	character(len=H_LONG)	''
BOUNDARY_POSTFIX_TIMELABEL	logical	.false.
BOUNDARY_TITLE	character(len=H_LONG)	'SCALE-RM BOUNDARY CONDITION for REAL CASE'
BOUNDARY_DTYPE	character(len=H_SHORT)	'DEFAULT'
BOUNDARY_UPDATE_DT	real(DP)	0.0_DP	inteval time of boudary data update [s]
FILTER_ORDER	integer	8	order of the hyper-diffusion (must be even)
FILTER_NITER	integer	5	times for hyper-diffusion iteration
USE_FILE_DENSITY	logical	.false.	use density data from files
SAME_MP_TYPE	logical	.false.	microphysics type of the parent model is same as it in this model

PARAM_MKINIT_REAL_LAND

name	type	default value	comment
NUMBER_OF_FILES	integer	1
NUMBER_OF_TSTEPS	integer	1	num of time steps in one file
NUMBER_OF_SKIP_TSTEPS	integer	0	num of skipped first several data
FILETYPE_ORG	character(len=H_LONG)	''
BASENAME_ORG	character(len=H_LONG)	''
BASENAME_ADD_NUM	logical	.false.
BASENAME_BOUNDARY	character(len=H_LONG)	''
BOUNDARY_POSTFIX_TIMELABEL	logical	.false.
BOUNDARY_TITLE	character(len=H_LONG)	'SCALE-RM BOUNDARY CONDITION for REAL CASE'
BOUNDARY_UPDATE_DT	real(DP)	0.0_DP	inteval time of boudary data update [s]
USE_FILE_LANDWATER	logical
INIT_LANDWATER_RATIO	real(RP)
INTRP_LAND_TEMP	character(len=*)
INTRP_LAND_WATER	character(len=*)
INTRP_LAND_SFC_TEMP	character(len=*)
INTRP_ITER_MAX	integer
FILTER_ORDER	integer	8	order of the hyper-diffusion (must be even)
FILTER_NITER	integer	5	times for hyper-diffusion iteration
SOILWATER_DS2VC	character(len=H_SHORT)	'limit'
ELEVATION_COLLECTION	logical
SERIAL_PROC_READ	logical	.true.	read by one MPI process and broadcast

PARAM_MKINIT_REAL_OCEAN

name	type	default value	comment
NUMBER_OF_FILES	integer	1
NUMBER_OF_TSTEPS	integer	1	num of time steps in one file
NUMBER_OF_SKIP_TSTEPS	integer	0	num of skipped first several data
FILETYPE_ORG	character(len=H_LONG)	''
BASENAME_ORG	character(len=H_LONG)	''
BASENAME_ADD_NUM	logical	.false.
BASENAME_BOUNDARY	character(len=H_LONG)	''
BOUNDARY_POSTFIX_TIMELABEL	logical	.false.
BOUNDARY_TITLE	character(len=H_LONG)	'SCALE-RM BOUNDARY CONDITION for REAL CASE'
BOUNDARY_UPDATE_DT	real(DP)	0.0_DP	inteval time of boudary data update [s]
INIT_OCEAN_ALB_LW	real(RP)
INIT_OCEAN_ALB_SW	real(RP)
INIT_OCEAN_Z0W	real(RP)
INTRP_OCEAN_TEMP	character(len=*)
INTRP_OCEAN_SFC_TEMP	character(len=*)
INTRP_ITER_MAX	integer
FILTER_ORDER	integer	8	order of the hyper-diffusion (must be even)
FILTER_NITER	integer	5	times for hyper-diffusion iteration
SERIAL_PROC_READ	logical	.true.	read by one MPI process and broadcast

History Output: No history output

Function/Subroutine Documentation

◆ realinput_atmos()

subroutine, public mod_realinput::realinput_atmos ( )

Definition at line 167 of file mod_realinput.F90.

References mod_atmos_admin::atmos_phy_mp_type, mod_atmos_vars::dens, scale_atmos_grid_cartesc_index::ia, scale_io::io_fid_conf, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::ka, mod_atmos_vars::momx, mod_atmos_vars::momy, mod_atmos_vars::momz, scale_prc::prc_abort(), scale_tracer::qa, mod_atmos_vars::qtrc, mod_atmos_vars::rhot, and scale_time::time_gettimelabel().

Referenced by mod_mkinit::read_sounding().

     use scale_time, only: &
        time_gettimelabel
     use mod_atmos_vars, only: &
        dens, &
        momz, &
        momx, &
        momy, &
        rhot, &
        qtrc
     use mod_atmos_admin, only: &
        atmos_phy_mp_type
     implicit none
 
     namelist / param_mkinit_real_atmos / &
        number_of_files,            &
        number_of_tsteps,           &
        number_of_skip_tsteps,      &
        serial_proc_read,           &
        filetype_org,               &
        basename_org,               &
        basename_add_num,           &
        basename_boundary,          &
        boundary_postfix_timelabel, &
        boundary_title,             &
        boundary_dtype,             &
        boundary_update_dt,         &
        filter_order,               &
        filter_niter,               &
        use_file_density,           &
        same_mp_type
 
     character(len=H_LONG) :: basename_mod
     character(len=H_LONG) :: basename_out_mod
     character(len=19)     :: timelabel
 
     integer  :: dims(6) ! dims 1-3: normal, 4-6: staggerd
     integer  :: timelen
 
     integer  :: fid_atmos
     integer  :: vid_atmos(5+qa)
 
     real(RP) :: dens_in(ka,ia,ja)
     real(RP) :: momz_in(ka,ia,ja) ! staggered point
     real(RP) :: momx_in(ka,ia,ja) ! staggered point
     real(RP) :: momy_in(ka,ia,ja) ! staggered point
     real(RP) :: rhot_in(ka,ia,ja)
     real(RP) :: qtrc_in(ka,ia,ja,qa)
 
     real(RP) :: velz_in(ka,ia,ja) ! staggered point
     real(RP) :: velx_in(ka,ia,ja) ! staggered point
     real(RP) :: vely_in(ka,ia,ja) ! staggered point
     real(RP) :: pott_in(ka,ia,ja)
 
     integer  :: ifile, istep, t, tall
     integer  :: k, i, j, iq
     integer  :: ierr
     !---------------------------------------------------------------------------
 
     log_newline
     log_info('REALINPUT_atmos',*) 'Setup'
 
     !--- read namelist
     rewind(io_fid_conf)
     read(io_fid_conf,nml=param_mkinit_real_atmos,iostat=ierr)
     if( ierr < 0 ) then !--- missing
        log_info("REALINPUT_atmos",*) 'Not found namelist. Default used.'
     elseif( ierr > 0 ) then !--- fatal error
        log_error("REALINPUT_atmos",*) 'Not appropriate names in namelist PARAM_MKINIT_REAL_ATMOS. Check!'
        call prc_abort
     endif
     log_nml(param_mkinit_real_atmos)
 
     if ( boundary_update_dt <= 0.0_dp ) then
        log_error("REALINPUT_atmos",*) 'BOUNDARY_UPDATE_DT is necessary in real case preprocess'
        call prc_abort
     endif
 
     if ( filetype_org == 'GrADS' ) then
        basename_mod = trim(basename_org) ! namelist file name
     else
        if ( number_of_files > 1 .OR. basename_add_num ) then
           basename_mod = trim(basename_org)//'_00000'
        else
           basename_mod = trim(basename_org)
        endif
     endif
 
     call parentatmossetup( filetype_org,     & ![IN]
                            basename_mod,     & ![IN]
                            serial_proc_read, & ![IN]
                            use_file_density, & ![IN]
                            dims(:),          & ![OUT]
                            timelen           ) ![OUT]
 
     if ( timelen > 0 ) then
        number_of_tsteps = timelen ! read from file
     endif
 
     log_newline
     log_info("REALINPUT_atmos",*) 'Number of temporal data in each file : ', number_of_tsteps
 
     do ifile = 1, number_of_files
 
        if ( filetype_org == 'GrADS' ) then
           if ( number_of_files > 1 .OR. basename_add_num ) then
              write(basename_mod,'(A,I5.5)') '_', ifile-1 ! only the number postfix
           else
              basename_mod = ''
           endif
        else
           if ( number_of_files > 1 .OR. basename_add_num ) then
              write(basename_mod,'(A,A,I5.5)') trim(basename_org), '_', ifile-1
           else
              basename_mod = trim(basename_org)
           endif
        endif
 
        log_newline
        log_info("REALINPUT_atmos",*) 'read external data from : ', trim(basename_mod)
 
        call parentatmosopen( filetype_org, & ![IN]
                              basename_mod, & ![IN]
                              dims(:)       ) ![IN]
 
        do istep = 1, number_of_tsteps
 
           tall = number_of_tsteps * (ifile-1) + istep ! consecutive time step (input)
           t    = tall - number_of_skip_tsteps         ! time step (output)
 
           if ( t <= 0 ) then
              log_progress('(1x,A,I4,A,I5,A,I6,A)') &
                           '[file,step,cons.] = [', ifile, ',', istep, ',', tall, '] ...skip.'
              cycle
           endif
 
           if ( t == 1 .OR. basename_boundary /= '' ) then
 
              log_progress('(1x,A,I4,A,I5,A,I6,A)') &
                           '[file,step,cons.] = [', ifile, ',', istep, ',', tall, ']'
 
              ! read prepared data
              call parentatmosinput( filetype_org,     & ! [IN]
                                     basename_mod,     & ! [IN]
                                     dims(:),          & ! [IN]
                                     istep,            & ! [IN]
                                     same_mp_type,     & ! [IN]
                                     dens_in(:,:,:),   & ! [OUT]
                                     momz_in(:,:,:),   & ! [OUT]
                                     momx_in(:,:,:),   & ! [OUT]
                                     momy_in(:,:,:),   & ! [OUT]
                                     rhot_in(:,:,:),   & ! [OUT]
                                     qtrc_in(:,:,:,:), & ! [OUT]
                                     velz_in(:,:,:),   & ! [OUT]
                                     velx_in(:,:,:),   & ! [OUT]
                                     vely_in(:,:,:),   & ! [OUT]
                                     pott_in(:,:,:)    ) ! [OUT]
           else
              log_progress('(1x,A,I4,A,I5,A,I6,A)') &
                           '[file,step,cons.] = [', ifile, ',', istep, ',', tall, '] ...skip.'
           endif
 
           !--- store prognostic variables as initial
           if ( t == 1 ) then
              log_newline
              log_info("REALINPUT_atmos",*) 'store initial state.'
 
              do j = 1, ja
              do i = 1, ia
              do k = 1, ka
                 dens(k,i,j) = dens_in(k,i,j)
                 momz(k,i,j) = momz_in(k,i,j)
                 momx(k,i,j) = momx_in(k,i,j)
                 momy(k,i,j) = momy_in(k,i,j)
                 rhot(k,i,j) = rhot_in(k,i,j)
              enddo
              enddo
              enddo
 
              do iq = 1, qa
              do j  = 1, ja
              do i  = 1, ia
              do k  = 1, ka
                 qtrc(k,i,j,iq) = qtrc_in(k,i,j,iq)
              enddo
              enddo
              enddo
              enddo
 
           endif
 
           !--- output boundary data
           if ( basename_boundary /= '' ) then
 
              if ( t == 1 ) then
                 if ( boundary_postfix_timelabel ) then
                    call time_gettimelabel( timelabel )
                    basename_out_mod = trim(basename_boundary)//'_'//trim(timelabel)
                 else
                    basename_out_mod = trim(basename_boundary)
                 endif
 
                 call boundaryatmossetup( basename_out_mod,   & ! [IN]
                                          boundary_title,     & ! [IN]
                                          boundary_dtype,     & ! [IN]
                                          boundary_update_dt, & ! [IN]
                                          fid_atmos,          & ! [OUT]
                                          vid_atmos(:)        ) ! [OUT]
              endif
 
              call boundaryatmosoutput( dens_in(:,:,:),     & ! [IN]
                                        velz_in(:,:,:),     & ! [IN]
                                        velx_in(:,:,:),     & ! [IN]
                                        vely_in(:,:,:),     & ! [IN]
                                        pott_in(:,:,:),     & ! [IN]
                                        qtrc_in(:,:,:,:),   & ! [IN]
                                        fid_atmos,          & ! [IN]
                                        vid_atmos(:),       & ! [IN]
                                        boundary_update_dt, & ! [IN]
                                        t                   ) ! [IN]
           endif
 
        enddo ! istep loop
     enddo ! ifile loop
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ realinput_surface()

subroutine, public mod_realinput::realinput_surface ( )

Definition at line 396 of file mod_realinput.F90.

References scale_atmos_grid_cartesc_metric::atmos_grid_cartesc_metric_rotc, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_cz, scale_atmos_hydrometeor::atmos_hydrometeor_dry, mod_atmos_admin::atmos_phy_ch_type, mod_atmos_phy_mp_driver::atmos_phy_mp_driver_qhyd2qtrc(), mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_albedo, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_temp, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_z0e, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_z0h, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_z0m, scale_comm_cartesc_nest::comm_cartesc_nest_interp_level, scale_const::const_eps, scale_const::const_i_lw, scale_const::const_i_sw, scale_const::const_laps, scale_const::const_undef, mod_atmos_vars::dens, scale_file_cartesc::file_cartesc_create(), scale_file_cartesc::file_cartesc_def_var(), scale_file_cartesc::file_cartesc_enddef(), scale_atmos_hydrometeor::i_qv, scale_cpl_sfc_index::i_r_diffuse, scale_cpl_sfc_index::i_r_direct, scale_cpl_sfc_index::i_r_ir, scale_cpl_sfc_index::i_r_nir, scale_cpl_sfc_index::i_r_vis, scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ie, igrads, scale_interp::interp_domain_compatibility(), scale_interp::interp_factor2d(), scale_interp::interp_factor3d(), scale_interp::interp_interp2d(), scale_interp::interp_interp3d(), scale_io::io_fid_conf, scale_atmos_grid_cartesc_index::is, iscale, iwrfarw, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::js, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, mod_land_admin::land_do, scale_land_grid_cartesc::land_grid_cartesc_cz, land_interporation(), mod_land_vars::land_sfc_albedo, mod_land_vars::land_sfc_temp, mod_land_vars::land_temp, mod_land_vars::land_water, scale_landuse::landuse_fact_land, scale_landuse::landuse_fact_ocean, scale_landuse::landuse_fact_urban, scale_landuse::landuse_frac_land, scale_land_grid_cartesc_index::lkmax, make_mask(), mod_atmos_vars::momx, mod_atmos_vars::momy, mod_atmos_vars::momz, scale_atmos_hydrometeor::n_hyd, scale_cpl_sfc_index::n_rad_dir, scale_cpl_sfc_index::n_rad_rgn, mod_ocean_admin::ocean_do, mod_ocean_vars::ocean_ice_mass, mod_ocean_vars::ocean_ice_temp, mod_ocean_vars::ocean_ocn_z0m, scale_ocean_phy_ice_simple::ocean_phy_ice_freezetemp, mod_ocean_vars::ocean_salt, mod_ocean_vars::ocean_sfc_albedo, mod_ocean_vars::ocean_sfc_temp, mod_ocean_vars::ocean_sfc_z0e, mod_ocean_vars::ocean_sfc_z0h, mod_ocean_vars::ocean_sfc_z0m, mod_ocean_vars::ocean_temp, mod_ocean_vars::ocean_uvel, mod_ocean_vars::ocean_vvel, scale_ocean_grid_cartesc_index::okmax, scale_ocean_grid_cartesc_index::oks, mod_realinput_grads::parentatmosinputgrads(), mod_realinput_scale::parentatmosinputscale(), mod_realinput_wrfarw::parentatmosinputwrfarw(), mod_realinput_grads::parentatmosopengrads(), mod_realinput_scale::parentatmosopenscale(), mod_realinput_wrfarw::parentatmosopenwrfarw(), mod_realinput_grads::parentatmossetupgrads(), mod_realinput_scale::parentatmossetupscale(), mod_realinput_wrfarw::parentatmossetupwrfarw(), mod_realinput_grads::parentlandinputgrads(), mod_realinput_scale::parentlandinputscale(), mod_realinput_wrfarw::parentlandinputwrfarw(), mod_realinput_grads::parentlandsetupgrads(), mod_realinput_scale::parentlandsetupscale(), mod_realinput_wrfarw::parentlandsetupwrfarw(), mod_realinput_grads::parentoceaninputgrads(), mod_realinput_scale::parentoceaninputscale(), mod_realinput_wrfarw::parentoceaninputwrfarw(), mod_realinput_grads::parentoceanopengrads(), mod_realinput_scale::parentoceanopenscale(), mod_realinput_wrfarw::parentoceanopenwrfarw(), mod_realinput_grads::parentoceansetupgrads(), mod_realinput_scale::parentoceansetupscale(), mod_realinput_wrfarw::parentoceansetupwrfarw(), scale_prc::prc_abort(), scale_prc::prc_ismaster, scale_prof::prof_rapend(), scale_prof::prof_rapstart(), scale_tracer::qa, mod_atmos_phy_ch_vars::qe_ch, mod_atmos_phy_mp_vars::qe_mp, scale_atmos_hydrometeor::qle, scale_atmos_hydrometeor::qls, mod_atmos_phy_ch_vars::qs_ch, mod_atmos_phy_mp_vars::qs_mp, mod_atmos_vars::qtrc, replace_misval_map(), mod_atmos_vars::rhot, scale_time::time_gettimelabel(), scale_time::time_nowdate, scale_topography::topo_zsfc, scale_tracer::tracer_cp, scale_tracer::tracer_cv, scale_tracer::tracer_mass, scale_tracer::tracer_name, scale_tracer::tracer_r, scale_urban_grid_cartesc_index::ukmax, mod_urban_admin::urban_do, mod_urban_vars::urban_qc, mod_urban_vars::urban_rainb, mod_urban_vars::urban_raing, mod_urban_vars::urban_rainr, mod_urban_vars::urban_roff, mod_urban_vars::urban_sfc_albedo, mod_urban_vars::urban_sfc_temp, mod_urban_vars::urban_tb, mod_urban_vars::urban_tbl, mod_urban_vars::urban_tc, mod_urban_vars::urban_tg, mod_urban_vars::urban_tgl, mod_urban_vars::urban_tr, mod_urban_vars::urban_trl, and mod_urban_vars::urban_uc.

Referenced by mod_mkinit::read_sounding().

     use scale_time, only: &
        time_gettimelabel
     use scale_landuse, only: &
        fact_ocean => landuse_fact_ocean, &
        fact_land  => landuse_fact_land,  &
        fact_urban => landuse_fact_urban
     use mod_atmos_phy_sf_vars, only: &
        atmos_phy_sf_sfc_temp,   &
        atmos_phy_sf_sfc_albedo, &
        atmos_phy_sf_sfc_z0m,    &
        atmos_phy_sf_sfc_z0h,    &
        atmos_phy_sf_sfc_z0e
     use mod_ocean_admin, only: &
        ocean_do
     use scale_ocean_phy_ice_simple, only: &
        ocean_phy_ice_freezetemp
     use mod_ocean_vars, only: &
        ocean_temp,       &
        ocean_salt,       &
        ocean_uvel,       &
        ocean_vvel,       &
        ocean_ocn_z0m,    &
        ocean_ice_temp,   &
        ocean_ice_mass,   &
        ocean_sfc_temp,   &
        ocean_sfc_albedo, &
        ocean_sfc_z0m,    &
        ocean_sfc_z0h,    &
        ocean_sfc_z0e
     use mod_land_admin, only: &
        land_do
     use mod_land_vars, only: &
        land_temp,       &
        land_water,      &
        land_sfc_temp,   &
        land_sfc_albedo
     use mod_urban_admin, only: &
        urban_do
     use mod_urban_vars, only: &
        urban_tc,        &
        urban_qc,        &
        urban_uc,        &
        urban_tr,        &
        urban_tb,        &
        urban_tg,        &
        urban_trl,       &
        urban_tbl,       &
        urban_tgl,       &
        urban_rainr,     &
        urban_rainb,     &
        urban_raing,     &
        urban_roff,      &
        urban_sfc_temp,  &
        urban_sfc_albedo
     implicit none
 
     logical                  :: use_file_landwater   = .true.    ! use land water data from files
     real(RP)                 :: init_landwater_ratio = 0.5_rp    ! Ratio of land water to storage is constant, if USE_FILE_LANDWATER is ".false."
     real(RP)                 :: init_ocean_alb_lw    = 0.04_rp   ! initial LW albedo on the ocean
     real(RP)                 :: init_ocean_alb_sw    = 0.10_rp   ! initial SW albedo on the ocean
     real(RP)                 :: init_ocean_z0w       = 1.0e-3_rp ! initial surface roughness on the ocean
     character(len=H_SHORT)   :: intrp_land_temp      = 'off'
     character(len=H_SHORT)   :: intrp_land_water     = 'off'
     character(len=H_SHORT)   :: intrp_land_sfc_temp  = 'off'
     character(len=H_SHORT)   :: intrp_ocean_temp     = 'off'
     character(len=H_SHORT)   :: intrp_ocean_sfc_temp = 'off'
     integer                  :: intrp_iter_max       = 100
     character(len=H_SHORT)   :: soilwater_ds2vc      = 'limit'
     logical                  :: soilwater_ds2vc_flag           ! true: 'critical', false: 'limit'
     logical                  :: elevation_collection = .true.
     logical                  :: elevation_collection_land
     logical                  :: elevation_collection_ocean
 
     namelist / param_mkinit_real_land / &
        number_of_files,            &
        number_of_tsteps,           &
        number_of_skip_tsteps,      &
        filetype_org,               &
        basename_org,               &
        basename_add_num,           &
        basename_boundary,          &
        boundary_postfix_timelabel, &
        boundary_title,             &
        boundary_update_dt,         &
        use_file_landwater,         &
        init_landwater_ratio,       &
        intrp_land_temp,            &
        intrp_land_water,           &
        intrp_land_sfc_temp,        &
        intrp_iter_max,             &
        filter_order,               &
        filter_niter,               &
        soilwater_ds2vc,            &
        elevation_collection,       &
        serial_proc_read
 
     namelist / param_mkinit_real_ocean / &
        number_of_files,            &
        number_of_tsteps,           &
        number_of_skip_tsteps,      &
        filetype_org,               &
        basename_org,               &
        basename_add_num,           &
        basename_boundary,          &
        boundary_postfix_timelabel, &
        boundary_title,             &
        boundary_update_dt,         &
        init_ocean_alb_lw,          &
        init_ocean_alb_sw,          &
        init_ocean_z0w,             &
        intrp_ocean_temp,           &
        intrp_ocean_sfc_temp,       &
        intrp_iter_max,             &
        filter_order,               &
        filter_niter,               &
        serial_proc_read
 
     character(len=H_LONG) :: filetype_land
     character(len=H_LONG) :: filetype_ocean
     character(len=H_LONG) :: basename_land
     character(len=H_LONG) :: basename_ocean
     character(len=5)      :: num               = ''
 
     ! land
     real(RP), allocatable :: land_temp_org      (:,:,:,:)
     real(RP), allocatable :: land_water_org     (:,:,:,:)
     real(RP), allocatable :: land_sfc_temp_org  (:,:,:)
     real(RP), allocatable :: land_sfc_albedo_org(:,:,:,:,:)
 
     ! urban
     real(RP) :: urban_tc_org(ia,ja)
     real(RP) :: urban_qc_org(ia,ja)
     real(RP) :: urban_uc_org(ia,ja)
     real(RP) :: urban_sfc_temp_org(ia,ja)
     real(RP) :: urban_sfc_albedo_org(ia,ja,n_rad_dir,n_rad_rgn)
 
     ! ocean
     real(RP), allocatable :: ocean_temp_org      (:,:,:,:)
     real(RP), allocatable :: ocean_sfc_temp_org  (:,:,:)
     real(RP), allocatable :: ocean_sfc_albedo_org(:,:,:,:,:)
     real(RP), allocatable :: ocean_sfc_z0_org    (:,:,:)
 
     integer :: number_of_files_land        = 1
     integer :: number_of_files_ocean       = 1
     integer :: number_of_tsteps_land       = 1       ! num of time steps in one file
     integer :: number_of_tsteps_ocean      = 1       ! num of time steps in one file
     integer :: number_of_skip_tsteps_land  = 0       ! num of skipped first several data
     integer :: number_of_skip_tsteps_ocean = 0       ! num of skipped first several data
 
     character(len=H_LONG) :: basename_boundary_land           = ''
     character(len=H_LONG) :: basename_boundary_ocean          = ''
     logical               :: boundary_postfix_timelabel_land  = .false.
     logical               :: boundary_postfix_timelabel_ocean = .false.
     character(len=H_LONG) :: boundary_title_land              = 'SCALE-RM BOUNDARY CONDITION for REAL CASE'
     character(len=H_LONG) :: boundary_title_ocean             = 'SCALE-RM BOUNDARY CONDITION for REAL CASE'
     real(DP)              :: boundary_update_dt_land          = 0.0_dp  ! inteval time of boudary data update [s]
     real(DP)              :: boundary_update_dt_ocean         = 0.0_dp  ! inteval time of boudary data update [s]
 
     integer :: mdlid_land, mdlid_ocean
     integer :: ldims(3), odims(2)
 
     integer :: totaltimesteps = 1
     integer :: timelen
     integer :: skip_steps
     integer :: ierr
 
     character(len=H_LONG) :: basename_out_mod
     character(len=19)     :: timelabel
 
     logical :: boundary_flag = .false.
     logical :: land_flag
 
     integer :: k, i, j, n, ns, ne, idir, irgn
     !---------------------------------------------------------------------------
 
     if ( land_do .or. urban_do ) then
        land_flag = .true.
     else
        land_flag = .false.
     end if
 
     if ( .not. land_flag .or. .not. ocean_do ) then
        log_error("REALINPUT_surface",*) 'OCEAN_ and LAND_DYN_TYPE must be set'
     end if
 
 
     log_newline
     log_info('REALINPUT_surface',*) 'Setup LAND'
 
     ! LAND/URBAN
 
     !--- read namelist
     rewind(io_fid_conf)
     read(io_fid_conf,nml=param_mkinit_real_land,iostat=ierr)
     if( ierr < 0 ) then !--- missing
        log_info("REALINPUT_surface",*) 'Not found namelist. Default used.'
     elseif( ierr > 0 ) then !--- fatal error
        log_error("REALINPUT_surface",*) 'Not appropriate names in namelist PARAM_MKINIT_REAL_LAND. Check!'
        call prc_abort
     endif
     log_nml(param_mkinit_real_land)
 
     number_of_files_land            = number_of_files
     number_of_tsteps_land           = number_of_tsteps
     number_of_skip_tsteps_land      = number_of_skip_tsteps
     filetype_land                   = filetype_org
     basename_boundary_land          = basename_boundary
     boundary_postfix_timelabel_land = boundary_postfix_timelabel
     boundary_title_land             = boundary_title
     boundary_update_dt_land         = boundary_update_dt
     elevation_collection_land       = elevation_collection
 
     if ( filetype_land .ne. "GrADS" .and. ( number_of_files > 1 .OR. basename_add_num ) ) then
        basename_land = trim(basename_org)//"_00000"
     else
        basename_land = trim(basename_org)
     endif
 
     select case( soilwater_ds2vc )
     case( 'critical' )
        soilwater_ds2vc_flag = .true.
     case('limit' )
        soilwater_ds2vc_flag = .false.
     case default
        log_error("REALINPUT_surface",*) 'Unsupported SOILWATER_DS2CV TYPE:', trim(soilwater_ds2vc)
        call prc_abort
     end select
 
     serial_land = serial_proc_read
 
     log_newline
     log_info('REALINPUT_surface',*) 'Setup OCEAN'
 
     !--- read namelist
     rewind(io_fid_conf)
     read(io_fid_conf,nml=param_mkinit_real_ocean,iostat=ierr)
     if( ierr < 0 ) then !--- missing
        log_info("REALINPUT_surface",*) 'Not found namelist. Default used.'
     elseif( ierr > 0 ) then !--- fatal error
        log_error("REALINPUT_surface",*) 'Not appropriate names in namelist PARAM_MKINIT_REAL_OCEAN. Check!'
        call prc_abort
     endif
     log_nml(param_mkinit_real_ocean)
 
     number_of_files_ocean            = number_of_files
     number_of_tsteps_ocean           = number_of_tsteps
     number_of_skip_tsteps_ocean      = number_of_skip_tsteps
     filetype_ocean                   = filetype_org
     basename_boundary_ocean          = basename_boundary
     boundary_postfix_timelabel_ocean = boundary_postfix_timelabel
     boundary_title_ocean             = boundary_title
     boundary_update_dt_ocean         = boundary_update_dt
     elevation_collection_ocean       = elevation_collection
 
     if ( filetype_ocean .ne. "GrADS" .and. ( number_of_files > 1 .OR. basename_add_num ) ) then
        basename_ocean = trim(basename_org)//"_00000"
     else
        basename_ocean = trim(basename_org)
     endif
 
     serial_ocean = serial_proc_read
 
     ! check land/ocean parameters
     if( number_of_files_land            .NE.   number_of_files_ocean            .OR. &
         number_of_tsteps_land           .NE.   number_of_tsteps_ocean           .OR. &
         number_of_skip_tsteps_land      .NE.   number_of_skip_tsteps_ocean      .OR. &
         basename_boundary_land          .NE.   basename_boundary_ocean          .OR. &
         boundary_postfix_timelabel_land .NEQV. boundary_postfix_timelabel_ocean .OR. &
         boundary_title_land             .NE.   boundary_title_ocean             .OR. &
         boundary_update_dt_land         .NE.   boundary_update_dt_ocean              ) then
        log_error("REALINPUT_surface",*) 'The following LAND/OCEAN parameters must be consistent due to technical problem:'
        log_error_cont(*) '           NUMBER_OF_FILES, NUMBER_OF_TSTEPS, NUMBER_OF_SKIP_TSTEPS,'
        log_error_cont(*) '           BASENAME_BOUNDARY, BOUNDARY_POSTFIX_TIMELABEL, BOUNDARY_TITLE, BOUNDARY_UPDATE_DT.'
        call prc_abort
     end if
 
     call parentsurfacesetup( ldims, odims,           & ![OUT]
                              mdlid_land,             & ![OUT]
                              mdlid_ocean,            & ![OUT]
                              timelen,                & ![OUT]
                              basename_land,          & ![IN]
                              basename_ocean,         & ![IN]
                              filetype_land,          & ![IN]
                              filetype_ocean,         & ![IN]
                              use_file_landwater,     & ![IN]
                              intrp_land_temp,        & ![IN]
                              intrp_land_water,       & ![IN]
                              intrp_land_sfc_temp,    & ![IN]
                              intrp_ocean_temp,       & ![IN]
                              intrp_ocean_sfc_temp    ) ![IN]
 
     if ( timelen > 0 ) then
        number_of_tsteps = timelen ! read from file
     endif
 
     totaltimesteps = number_of_files * number_of_tsteps
 
     allocate( land_temp_org(lkmax,ia,ja,                    1+number_of_skip_tsteps:totaltimesteps) )
     allocate( land_water_org(lkmax,ia,ja,                    1+number_of_skip_tsteps:totaltimesteps) )
     allocate( land_sfc_temp_org(      ia,ja,                    1+number_of_skip_tsteps:totaltimesteps) )
     allocate( land_sfc_albedo_org(      ia,ja,n_rad_dir,n_rad_rgn,1+number_of_skip_tsteps:totaltimesteps) )
 
     allocate( ocean_temp_org(okmax,ia,ja,                    1+number_of_skip_tsteps:totaltimesteps) )
     allocate( ocean_sfc_temp_org(      ia,ja,                    1+number_of_skip_tsteps:totaltimesteps) )
     allocate( ocean_sfc_albedo_org(      ia,ja,n_rad_dir,n_rad_rgn,1+number_of_skip_tsteps:totaltimesteps) )
     allocate( ocean_sfc_z0_org(      ia,ja,                    1+number_of_skip_tsteps:totaltimesteps) )
 
     if ( mdlid_ocean == igrads ) then
        basename_org = ""
     endif
 
     if ( basename_boundary /= '' ) then
        boundary_flag = .true.
     endif
 
     !--- read external file
     do n = 1, number_of_files
 
        if ( number_of_files > 1 .OR. basename_add_num ) then
           write(num,'(I5.5)') n-1
           basename_land  = trim(basename_org)//"_"//num
           basename_ocean = trim(basename_org)//"_"//num
        else
           basename_land  = trim(basename_org)
           basename_ocean = trim(basename_org)
        endif
 
        log_newline
        log_info("REALINPUT_surface",*) 'Target File Name (Land) : ', trim(basename_land)
        log_info("REALINPUT_surface",*) 'Target File Name (Ocean): ', trim(basename_ocean)
        log_info("REALINPUT_surface",*) 'Time Steps in One File  : ', number_of_tsteps
 
        ns = number_of_tsteps * (n - 1) + 1
        ne = ns + (number_of_tsteps - 1)
 
        if ( ne <= number_of_skip_tsteps ) then
           log_info("REALINPUT_surface",*) '    SKIP'
           cycle
        endif
 
        skip_steps = max(number_of_skip_tsteps - ns + 1, 0)
        ns = max(ns, number_of_skip_tsteps+1)
 
        ! read all prepared data
        call parentsurfaceinput( land_temp_org(:,:,:,  ns:ne), &
                                 land_water_org(:,:,:,  ns:ne), &
                                 land_sfc_temp_org(:,:,    ns:ne), &
                                 land_sfc_albedo_org(:,:,:,:,ns:ne), &
                                 urban_tc_org,         &
                                 urban_qc_org,         &
                                 urban_uc_org,         &
                                 urban_sfc_temp_org,   &
                                 urban_sfc_albedo_org, &
                                 ocean_temp_org(oks,:,:,    ns:ne), &
                                 ocean_sfc_temp_org(    :,:,    ns:ne), &
                                 ocean_sfc_albedo_org(    :,:,:,:,ns:ne), &
                                 ocean_sfc_z0_org(    :,:,    ns:ne), &
                                 basename_land,                &
                                 basename_ocean,               &
                                 mdlid_land, mdlid_ocean,      &
                                 ldims, odims,                 &
                                 use_file_landwater,           &
                                 init_landwater_ratio,         &
                                 init_ocean_alb_lw,            &
                                 init_ocean_alb_sw,            &
                                 init_ocean_z0w,               &
                                 intrp_iter_max,               &
                                 soilwater_ds2vc_flag,         &
                                 elevation_collection_land,    &
                                 elevation_collection_ocean,   &
                                 boundary_flag,                &
                                 number_of_tsteps, skip_steps, &
                                 urban_do                      )
 
        ! required one-step data only
        if( basename_boundary == '' ) exit
 
     enddo
 
 
     !--- input initial data
     ns = number_of_skip_tsteps + 1  ! skip first several data
 
     do j = 1, ja
     do i = 1, ia
        ocean_sfc_temp(i,j) = ocean_sfc_temp_org(i,j,ns)
        ocean_sfc_z0m(i,j) = ocean_sfc_z0_org(i,j,ns)
        ocean_sfc_z0h(i,j) = ocean_sfc_z0_org(i,j,ns)
        ocean_sfc_z0e(i,j) = ocean_sfc_z0_org(i,j,ns)
        do irgn = i_r_ir, i_r_vis
        do idir = i_r_direct, i_r_diffuse
           ocean_sfc_albedo(i,j,idir,irgn) = ocean_sfc_albedo_org(i,j,idir,irgn,ns)
        enddo
        enddo
        do k = 1, okmax
           ocean_temp(k,i,j) = ocean_temp_org(oks,i,j,ns)
           ocean_salt(k,i,j) = 0.0_rp
           ocean_uvel(k,i,j) = 0.0_rp
           ocean_vvel(k,i,j) = 0.0_rp
        enddo
        ocean_ocn_z0m(i,j) = ocean_sfc_z0_org(i,j,ns)
        ocean_ice_temp(i,j) = min( ocean_sfc_temp_org(i,j,ns), ocean_phy_ice_freezetemp )
        ocean_ice_mass(i,j) = 0.0_rp
 
        land_sfc_temp(i,j)      = land_sfc_temp_org(i,j,     ns)
        do irgn = i_r_ir, i_r_vis
        do idir = i_r_direct, i_r_diffuse
           land_sfc_albedo(i,j,idir,irgn) = land_sfc_albedo_org(i,j,idir,irgn,ns)
        enddo
        enddo
        do k = 1, lkmax
           land_temp(k,i,j) = land_temp_org(k,i,j,ns)
           land_water(k,i,j) = land_water_org(k,i,j,ns)
        enddo
 
        if ( urban_do ) then
           urban_sfc_temp(i,j)      = urban_sfc_temp_org(i,j)
           do irgn = i_r_ir, i_r_vis
           do idir = i_r_direct, i_r_diffuse
              urban_sfc_albedo(i,j,idir,irgn) = urban_sfc_albedo_org(i,j,idir,irgn)
           enddo
           enddo
           do k = 1, ukmax
              urban_trl(k,i,j) = urban_sfc_temp_org(i,j)
              urban_tbl(k,i,j) = urban_sfc_temp_org(i,j)
              urban_tgl(k,i,j) = urban_sfc_temp_org(i,j)
           enddo
           urban_tc(i,j) = urban_tc_org(i,j)
           urban_qc(i,j) = urban_qc_org(i,j)
           urban_uc(i,j) = urban_uc_org(i,j)
           urban_tr(i,j) = urban_sfc_temp_org(i,j)
           urban_tb(i,j) = urban_sfc_temp_org(i,j)
           urban_tg(i,j) = urban_sfc_temp_org(i,j)
           urban_rainr(i,j) = 0.0_rp
           urban_rainb(i,j) = 0.0_rp
           urban_raing(i,j) = 0.0_rp
           urban_roff(i,j) = 0.0_rp
        end if
 
        atmos_phy_sf_sfc_z0m(i,j) = ocean_sfc_z0m(i,j)
        atmos_phy_sf_sfc_z0h(i,j) = ocean_sfc_z0h(i,j)
        atmos_phy_sf_sfc_z0e(i,j) = ocean_sfc_z0e(i,j)
 
        if ( urban_do ) then
           atmos_phy_sf_sfc_temp(i,j) = fact_ocean(i,j) * ocean_sfc_temp(i,j) &
                                      + fact_land(i,j) * land_sfc_temp(i,j) &
                                      + fact_urban(i,j) * urban_sfc_temp(i,j)
           do irgn = i_r_ir, i_r_vis
           do idir = i_r_direct, i_r_diffuse
              atmos_phy_sf_sfc_albedo(i,j,idir,irgn) = fact_ocean(i,j) * ocean_sfc_albedo(i,j,idir,irgn) &
                                                     + fact_land(i,j) * land_sfc_albedo(i,j,idir,irgn) &
                                                     + fact_urban(i,j) * urban_sfc_albedo(i,j,idir,irgn)
           enddo
           enddo
        else
           atmos_phy_sf_sfc_temp(i,j) = fact_ocean(i,j) * ocean_sfc_temp(i,j) &
                                      + fact_land(i,j) * land_sfc_temp(i,j)
           do irgn = i_r_ir, i_r_vis
           do idir = i_r_direct, i_r_diffuse
              atmos_phy_sf_sfc_albedo(i,j,idir,irgn) = fact_ocean(i,j) * ocean_sfc_albedo(i,j,idir,irgn) &
                                                     + fact_land(i,j) * land_sfc_albedo(i,j,idir,irgn)
           enddo
           enddo
        endif
     enddo
     enddo
 
 
     !--- output boundary data
     if( basename_boundary /= '' ) then
        totaltimesteps = totaltimesteps - number_of_skip_tsteps ! skip first several data
        if ( totaltimesteps > 1 ) then
           if ( boundary_update_dt <= 0.0_dp ) then
              log_error("REALINPUT_surface",*) 'BOUNDARY_UPDATE_DT is necessary in real case preprocess'
              call prc_abort
           endif
 
           if ( boundary_postfix_timelabel ) then
              call time_gettimelabel( timelabel )
              basename_out_mod = trim(basename_boundary)//'_'//trim(timelabel)
           else
              basename_out_mod = trim(basename_boundary)
           endif
 
           call parentsurfaceboundary( land_temp_org(:,:,:,ns:ne), &
                                       land_water_org(:,:,:,ns:ne), &
                                       land_sfc_temp_org(  :,:,ns:ne), &
                                       ocean_temp_org(:,:,:,ns:ne), &
                                       ocean_sfc_temp_org(  :,:,ns:ne), &
                                       ocean_sfc_z0_org(  :,:,ns:ne), &
                                       totaltimesteps,                  &
                                       boundary_update_dt,              &
                                       basename_out_mod,                &
                                       boundary_title                   )
 
        endif
     endif
 
     deallocate( land_temp_org        )
     deallocate( land_water_org       )
     deallocate( land_sfc_temp_org    )
     deallocate( land_sfc_albedo_org  )
     deallocate( ocean_temp_org       )
     deallocate( ocean_sfc_temp_org   )
     deallocate( ocean_sfc_albedo_org )
     deallocate( ocean_sfc_z0_org     )
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ land_interporation()

subroutine mod_realinput::land_interporation	(	real(rp), dimension(lkmax,ia,ja), intent(out)	tg,
		real(rp), dimension(lkmax,ia,ja), intent(out)	strg,
		real(rp), dimension(ia,ja), intent(out)	lst,
		real(rp), dimension(ia,ja,n_rad_dir,n_rad_rgn), intent(out)	albg,
		real(rp), dimension(ia,ja), intent(out)	ust,
		real(rp), dimension(ia,ja,n_rad_dir,n_rad_rgn), intent(out)	albu,
		real(rp), dimension(:,:,:), intent(inout)	tg_org,
		real(rp), dimension(:,:,:), intent(inout)	strg_org,
		real(rp), dimension(:,:,:), intent(inout)	smds_org,
		real(rp), dimension(:,:), intent(inout)	lst_org,
		real(rp), dimension(:,:,:,:), intent(inout)	albg_org,
		real(rp), dimension(:,:), intent(inout)	ust_org,
		real(rp), dimension(:,:), intent(inout)	sst_org,
		real(rp), dimension(:,:), intent(in)	lmask_org,
		real(rp), dimension(:,:), intent(in)	lsmask_nest,
		real(rp), dimension(:,:), intent(in)	topo_org,
		real(rp), dimension(:), intent(in)	lz_org,
		real(rp), dimension(:,:), intent(in)	llon_org,
		real(rp), dimension(:,:), intent(in)	llat_org,
		real(rp), dimension(lkmax), intent(in)	LCZ,
		real(rp), dimension(ia,ja), intent(in)	LON,
		real(rp), dimension(ia,ja), intent(in)	LAT,
		integer, dimension(3), intent(in)	ldims,
		integer, dimension(2), intent(in)	odims,
		real(rp), intent(in)	maskval_tg,
		real(rp), intent(in)	maskval_strg,
		real(rp), intent(in)	init_landwater_ratio,
		logical, intent(in)	use_file_landwater,
		logical, intent(in)	use_waterratio,
		logical, intent(in)	soilwater_ds2vc_flag,
		logical, intent(in)	elevation_collection,
		integer, intent(in)	intrp_iter_max,
		logical, intent(in)	ol_interp,
		logical, intent(in)	URBAN_do
	)

Definition at line 2909 of file mod_realinput.F90.

References scale_const::const_eps, scale_const::const_i_lw, scale_const::const_i_sw, scale_const::const_laps, scale_const::const_undef, mod_land_vars::convert_ws2vwc(), scale_cpl_sfc_index::i_r_diffuse, scale_cpl_sfc_index::i_r_direct, scale_cpl_sfc_index::i_r_ir, scale_cpl_sfc_index::i_r_nir, scale_cpl_sfc_index::i_r_vis, scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ie, scale_interp::interp_factor2d(), scale_interp::interp_factor3d(), scale_interp::interp_interp2d(), scale_interp::interp_interp3d(), scale_atmos_grid_cartesc_index::is, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::js, scale_land_grid_cartesc_index::lke, scale_land_grid_cartesc_index::lkmax, scale_land_grid_cartesc_index::lks, make_mask(), scale_prc::prc_abort(), replace_misval_const(), replace_misval_map(), and scale_topography::topo_zsfc.

Referenced by realinput_surface().

     use scale_prc, only: &
          prc_abort
     use scale_const, only: &
          undef => const_undef, &
          eps   => const_eps, &
          i_sw => const_i_sw, &
          i_lw => const_i_lw, &
          laps => const_laps
     use scale_interp, only: &
          interp_factor2d, &
          interp_factor3d, &
          interp_interp2d, &
          interp_interp3d
     use scale_comm_cartesc, only: &
          comm_vars8, &
          comm_wait
     use scale_filter, only: &
          filter_hyperdiff
     use scale_topography, only: &
          topo_zsfc
     use mod_land_vars, only: &
          convert_ws2vwc
     implicit none
     real(RP), intent(out)   :: tg(lkmax,ia,ja)
     real(RP), intent(out)   :: strg(lkmax,ia,ja)
     real(RP), intent(out)   :: lst(ia,ja)
     real(RP), intent(out)   :: albg(ia,ja,n_rad_dir,n_rad_rgn)
     real(RP), intent(out)   :: ust(ia,ja)
     real(RP), intent(out)   :: albu(ia,ja,n_rad_dir,n_rad_rgn)
     real(RP), intent(inout) :: tg_org(:,:,:)
     real(RP), intent(inout) :: strg_org(:,:,:)
     real(RP), intent(inout) :: smds_org(:,:,:)
     real(RP), intent(inout) :: lst_org(:,:)
     real(RP), intent(inout) :: albg_org(:,:,:,:)
     real(RP), intent(inout) :: ust_org(:,:)
     real(RP), intent(inout) :: sst_org(:,:)
     real(RP), intent(in)    :: lmask_org(:,:)
     real(RP), intent(in)    :: lsmask_nest(:,:)
     real(RP), intent(in)    :: topo_org(:,:)
     real(RP), intent(in)    :: lz_org(:)
     real(RP), intent(in)    :: llon_org(:,:)
     real(RP), intent(in)    :: llat_org(:,:)
     real(RP), intent(in)    :: lcz(lkmax)
     real(RP), intent(in)    :: lon(ia,ja)
     real(RP), intent(in)    :: lat(ia,ja)
     integer,  intent(in)    :: ldims(3)
     integer,  intent(in)    :: odims(2)
     real(RP), intent(in)    :: maskval_tg
     real(RP), intent(in)    :: maskval_strg
     real(RP), intent(in)    :: init_landwater_ratio
     logical,  intent(in)    :: use_file_landwater
     logical,  intent(in)    :: use_waterratio
     logical,  intent(in)    :: soilwater_ds2vc_flag
     logical,  intent(in)    :: elevation_collection
     integer,  intent(in)    :: intrp_iter_max
     logical,  intent(in)    :: ol_interp
     logical,  intent(in)    :: urban_do
 
     real(RP) :: lmask(ldims(2), ldims(3))
     real(RP) :: smds(lkmax,ia,ja)
 
     ! data for interporation
     real(RP) :: hfact_l(ldims(2), ldims(3), itp_nh)
     integer  :: igrd_l (ldims(2), ldims(3), itp_nh)
     integer  :: jgrd_l (ldims(2), ldims(3), itp_nh)
     real(RP) :: vfactl(lkmax,ia,ja,itp_nh,itp_nv)
     integer  :: kgrdl (lkmax,ia,ja,itp_nh,itp_nv)
 
     real(RP) :: sst_land(ldims(2), ldims(3))
     real(RP) :: work(ldims(2), ldims(3))
 
     real(RP) :: lz3d_org(ldims(1),ldims(2),ldims(3))
     real(RP) :: lcz_3d(lkmax,ia,ja)
 
     ! elevation collection
     real(RP) :: topo(ia,ja)
     real(RP) :: tdiff
 
     real(RP) :: one2d(ia,ja)
     real(RP) :: one3d(ka,ia,ja)
 
     integer :: k, i, j
 
 
     ! Surface skin temp: interpolate over the ocean
     if ( i_intrp_land_sfc_temp .ne. i_intrp_off ) then
        select case( i_intrp_land_sfc_temp )
        case( i_intrp_mask )
           lmask = lmask_org
        case( i_intrp_fill )
           call make_mask( lmask, lst_org, ldims(2), ldims(3), landdata=.true.)
        case default
           log_error("land_interporation",*) 'INTRP_LAND_SFC_TEMP is invalid.'
           call prc_abort
        end select
        call interp_oceanland_data(lst_org, lmask, ldims(2), ldims(3), .true., intrp_iter_max)
     end if
 
     ! Urban surface temp: interpolate over the ocean
     ! if ( i_INTRP_URB_SFC_TEMP .ne. i_intrp_off ) then
     !   select case( i_INTRP_URB_SFC_TEMP )
     !   case( i_intrp_mask )
     !      lmask = lmask_org
     !   case( i_intrp_fill )
     !      call make_mask( lmask, ust_org, ldims(2), ldims(3), landdata=.true.)
     !   case default
     !      LOG_ERROR("land_interporation",*) 'INTRP_URB_SFC_TEMP is invalid.'
     !      call PRC_abort
     !   end select
     !   call interp_OceanLand_data(ust_org, lmask, ldims(2), ldims(3), .true., intrp_iter_max)
     !end if
 
     if ( ol_interp ) then
        ! interpolation facter between outer land grid and ocean grid
        call interp_factor2d( itp_nh,             & ! [IN]
                              odims(1), odims(2), & ! [IN]
                              olon_org(:,:),      & ! [IN]
                              olat_org(:,:),      & ! [IN]
                              ldims(2), ldims(3), & ! [IN]
                              llon_org(:,:),      & ! [IN]
                              llat_org(:,:),      & ! [IN]
                              igrd_l(:,:,:),    & ! [OUT]
                              jgrd_l(:,:,:),    & ! [OUT]
                              hfact_l(:,:,:)     ) ! [OUT]
 
        ! sst on land grid
        call interp_interp2d( itp_nh,             & ! [IN]
                              odims(1), odims(2), & ! [IN]
                              ldims(2), ldims(3), & ! [IN]
                              igrd_l(:,:,:),   & ! [IN]
                              jgrd_l(:,:,:),   & ! [IN]
                              hfact_l(:,:,:),   & ! [IN]
                              sst_org(:,:),     & ! [IN]
                              sst_land(:,:)      ) ! [OUT]
     else
        sst_land(:,:) = sst_org(:,:)
     end if
 
     do j = 1, ldims(3)
     do i = 1, ldims(2)
        if ( topo_org(i,j) > undef + eps ) then ! ignore UNDEF value
           sst_land(i,j) = sst_land(i,j) - topo_org(i,j) * laps
        end if
     end do
     end do
 
     call replace_misval_map( lst_org, sst_land, ldims(2), ldims(3), "SKINT" )
 
     ! replace missing value
     do j = 1, ldims(3)
     do i = 1, ldims(2)
 !       if ( skinw_org(i,j) == UNDEF ) skinw_org(i,j) = 0.0_RP
 !       if ( snowq_org(i,j) == UNDEF ) snowq_org(i,j) = 0.0_RP
 !       if ( snowt_org(i,j) == UNDEF ) snowt_org(i,j) = TEM00
        if( albg_org(i,j,i_r_direct ,i_r_ir ) == undef ) albg_org(i,j,i_r_direct ,i_r_ir ) = 0.03_rp
        if( albg_org(i,j,i_r_diffuse,i_r_ir ) == undef ) albg_org(i,j,i_r_diffuse,i_r_ir ) = 0.03_rp  ! emissivity of general ground surface : 0.95-0.98
        if( albg_org(i,j,i_r_direct ,i_r_nir) == undef ) albg_org(i,j,i_r_direct ,i_r_nir) = 0.22_rp
        if( albg_org(i,j,i_r_diffuse,i_r_nir) == undef ) albg_org(i,j,i_r_diffuse,i_r_nir) = 0.22_rp
        if( albg_org(i,j,i_r_direct ,i_r_vis) == undef ) albg_org(i,j,i_r_direct ,i_r_vis) = 0.22_rp
        if( albg_org(i,j,i_r_diffuse,i_r_vis) == undef ) albg_org(i,j,i_r_diffuse,i_r_vis) = 0.22_rp
     end do
     end do
     if ( urban_do ) then
        do j = 1, ldims(3)
        do i = 1, ldims(2)
           if ( ust_org(i,j) == undef ) ust_org(i,j) = lst_org(i,j)
        end do
        end do
     end if
 
     ! Land temp: interpolate over the ocean
     if ( i_intrp_land_temp .ne. i_intrp_off ) then
        do k = 1, ldims(1)
           work(:,:) = tg_org(k,:,:)
           select case( i_intrp_land_temp )
           case( i_intrp_mask )
              lmask = lmask_org
           case( i_intrp_fill )
              call make_mask( lmask, work, ldims(2), ldims(3), landdata=.true.)
           end select
           call interp_oceanland_data( work, lmask, ldims(2), ldims(3), .true., intrp_iter_max )
           !replace land temp using skin temp
           call replace_misval_map( work, lst_org, ldims(2),  ldims(3),  "STEMP")
           tg_org(k,:,:) = work(:,:)
        end do
     end if
 
 
     ! fill grid data
     do j = 1, ldims(3)
     do i = 1, ldims(2)
        lz3d_org(:,i,j) = lz_org(:)
     end do
     end do
 
     do j = 1, ja
     do i = 1, ia
        lcz_3d(:,i,j) = lcz(:)
     enddo
     enddo
 
     call interp_factor3d( itp_nh,                & ! [IN]
                           ldims(1), 1, ldims(1), & ! [IN]
                           ldims(2), ldims(3),    & ! [IN]
                           llon_org(:,:),         & ! [IN]
                           llat_org(:,:),         & ! [IN]
                           lz3d_org(:,:,:),       & ! [IN]
                           lkmax, lks, lke,       & ! [IN]
                           ia, ja,                & ! [IN]
                           lon(:,:),         & ! [IN]
                           lat(:,:),         & ! [IN]
                           lcz_3d(:,:,:),       & ! [IN]
                           igrd(    :,:,:),   & ! [OUT]
                           jgrd(    :,:,:),   & ! [OUT]
                           hfact(    :,:,:),   & ! [OUT]
                           kgrdl(:,:,:,:,:),   & ! [OUT]
                           vfactl(:,:,:,:,:)    ) ! [OUT]
 
     call interp_interp2d( itp_nh,             & ! [IN]
                           ldims(2), ldims(3), & ! [IN]
                           ia, ja,             & ! [IN]
                           igrd(:,:,:),    & ! [IN]
                           jgrd(:,:,:),    & ! [IN]
                           hfact(:,:,:),    & ! [IN]
                           lst_org(:,:),      & ! [IN]
                           lst(:,:)       ) ! [OUT]
     if ( filter_niter > 0 ) then
        call filter_hyperdiff( ia, is, ie, ja, js, je, &
                               lst(:,:), filter_order, filter_niter )
        call comm_vars8( lst(:,:), 1 )
        call comm_wait ( lst(:,:), 1 )
     end if
 
     if ( urban_do ) then
        call interp_interp2d( itp_nh,             & ! [IN]
                              ldims(2), ldims(3), & ! [IN]
                              ia, ja,             & ! [IN]
                              igrd(:,:,:),    & ! [IN]
                              jgrd(:,:,:),    & ! [IN]
                              hfact(:,:,:),    & ! [IN]
                              ust_org(:,:),      & ! [IN]
                              ust(:,:)       ) ! [OUT]
        if ( filter_niter > 0 ) then
           call filter_hyperdiff( ia, is, ie, ja, js, je, &
                                  ust(:,:), filter_order, filter_niter )
        call comm_vars8( ust(:,:), 1 )
        call comm_wait ( ust(:,:), 1 )
        end if
     end if
 
     call interp_interp2d( itp_nh,                            & ! [IN]
                           ldims(2), ldims(3),                & ! [IN]
                           ia, ja,                            & ! [IN]
                           igrd(:,:,:),                   & ! [IN]
                           jgrd(:,:,:),                   & ! [IN]
                           hfact(:,:,:),                   & ! [IN]
                           albg_org(:,:,i_r_direct ,i_r_ir ), & ! [IN]
                           albg(:,:,i_r_direct ,i_r_ir )  ) ! [OUT]
     if ( filter_niter > 0 ) then
        one2d(:,:) = 1.0_rp
        call filter_hyperdiff( ia, is, ie, ja, js, je, &
                               albg(:,:,i_r_direct,i_r_ir), filter_order, filter_niter, &
                               limiter_sign = one2d(:,:) )
        call comm_vars8( albg(:,:,i_r_direct,i_r_ir), 1 )
        call comm_wait ( albg(:,:,i_r_direct,i_r_ir), 1 )
     end if
 
     call interp_interp2d( itp_nh,                            & ! [IN]
                           ldims(2), ldims(3),                & ! [IN]
                           ia, ja,                            & ! [IN]
                           igrd(:,:,:),                   & ! [IN]
                           jgrd(:,:,:),                   & ! [IN]
                           hfact(:,:,:),                   & ! [IN]
                           albg_org(:,:,i_r_diffuse,i_r_ir ), & ! [IN]
                           albg(:,:,i_r_diffuse,i_r_ir )  ) ! [OUT]
     if ( filter_niter > 0 ) then
        call filter_hyperdiff( ia, is, ie, ja, js, je, &
                               albg(:,:,i_r_diffuse,i_r_ir), filter_order, filter_niter, &
                               limiter_sign = one2d(:,:) )
        call comm_vars8( albg(:,:,i_r_diffuse,i_r_ir), 1 )
        call comm_wait ( albg(:,:,i_r_diffuse,i_r_ir), 1 )
     end if
 
     call interp_interp2d( itp_nh,                            & ! [IN]
                           ldims(2), ldims(3),                & ! [IN]
                           ia, ja,                            & ! [IN]
                           igrd(:,:,:),                   & ! [IN]
                           jgrd(:,:,:),                   & ! [IN]
                           hfact(:,:,:),                   & ! [IN]
                           albg_org(:,:,i_r_direct ,i_r_nir), & ! [IN]
                           albg(:,:,i_r_direct ,i_r_nir)  ) ! [OUT]
     if ( filter_niter > 0 ) then
        call filter_hyperdiff( ia, is, ie, ja, js, je, &
                               albg(:,:,i_r_direct,i_r_nir), filter_order, filter_niter, &
                               limiter_sign = one2d(:,:) )
        call comm_vars8( albg(:,:,i_r_direct,i_r_nir), 1 )
        call comm_wait ( albg(:,:,i_r_direct,i_r_nir), 1 )
     end if
 
     call interp_interp2d( itp_nh,                            & ! [IN]
                           ldims(2), ldims(3),                & ! [IN]
                           ia, ja,                            & ! [IN]
                           igrd(:,:,:),                   & ! [IN]
                           jgrd(:,:,:),                   & ! [IN]
                           hfact(:,:,:),                   & ! [IN]
                           albg_org(:,:,i_r_diffuse,i_r_nir), & ! [IN]
                           albg(:,:,i_r_diffuse,i_r_nir)  ) ! [OUT]
     if ( filter_niter > 0 ) then
        call filter_hyperdiff( ia, is, ie, ja, js, je, &
                               albg(:,:,i_r_diffuse,i_r_nir), filter_order, filter_niter, &
                               limiter_sign = one2d(:,:) )
        call comm_vars8( albg(:,:,i_r_diffuse,i_r_nir), 1 )
        call comm_wait ( albg(:,:,i_r_diffuse,i_r_nir), 1 )
     end if
 
     call interp_interp2d( itp_nh,                            & ! [IN]
                           ldims(2), ldims(3),                & ! [IN]
                           ia, ja,                            & ! [IN]
                           igrd(:,:,:),                   & ! [IN]
                           jgrd(:,:,:),                   & ! [IN]
                           hfact(:,:,:),                   & ! [IN]
                           albg_org(:,:,i_r_direct ,i_r_vis), & ! [IN]
                           albg(:,:,i_r_direct ,i_r_vis)  ) ! [OUT]
     if ( filter_niter > 0 ) then
        call filter_hyperdiff( ia, is, ie, ja, js, je, &
                               albg(:,:,i_r_direct,i_r_vis), filter_order, filter_niter, &
                               limiter_sign = one2d(:,:) )
        call comm_vars8( albg(:,:,i_r_direct,i_r_vis), 1 )
        call comm_wait ( albg(:,:,i_r_direct,i_r_vis), 1 )
     end if
 
     call interp_interp2d( itp_nh,                            & ! [IN]
                           ldims(2), ldims(3),                & ! [IN]
                           ia, ja,                            & ! [IN]
                           igrd(:,:,:),                   & ! [IN]
                           jgrd(:,:,:),                   & ! [IN]
                           hfact(:,:,:),                   & ! [IN]
                           albg_org(:,:,i_r_diffuse,i_r_vis), & ! [IN]
                           albg(:,:,i_r_diffuse,i_r_vis)  ) ! [OUT]
     if ( filter_niter > 0 ) then
        call filter_hyperdiff( ia, is, ie, ja, js, je, &
                               albg(:,:,i_r_diffuse,i_r_vis), filter_order, filter_niter, &
                               limiter_sign = one2d(:,:) )
        call comm_vars8( albg(:,:,i_r_diffuse,i_r_vis), 1 )
        call comm_wait ( albg(:,:,i_r_diffuse,i_r_vis), 1 )
     end if
 
     call interp_interp3d( itp_nh,                       & ! [IN]
                           ldims(1), ldims(2), ldims(3), & ! [IN]
                           lkmax, lks, lke,              & ! [IN]
                           ia, ja,                       & ! [IN]
                           igrd(    :,:,:),            & ! [IN]
                           jgrd(    :,:,:),            & ! [IN]
                           hfact(    :,:,:),            & ! [IN]
                           kgrdl(:,:,:,:,:),            & ! [IN]
                           vfactl(:,:,:,:,:),            & ! [IN]
                           tg_org(:,:,:),                & ! [IN]
                           tg(:,:,:)                 ) ! [OUT]
 
     do j = 1, ja
     do i = 1, ia
        tg(lkmax,i,j) = tg(lkmax-1,i,j)
     enddo
     enddo
 
     ! replace values over the ocean
     do k = 1, lkmax
        call replace_misval_const( tg(k,:,:), maskval_tg, lsmask_nest )
     enddo
     if ( filter_niter > 0 ) then
        call filter_hyperdiff( lkmax, 1, lkmax, ia, is, ie, ja, js, je, &
                               tg(:,:,:), filter_order, filter_niter )
        call comm_vars8( tg(:,:,:), 1 )
        call comm_wait ( tg(:,:,:), 1 )
     end if
 
 
     ! elevation collection
     if ( elevation_collection ) then
        call interp_interp2d( itp_nh,             & ! [IN]
                              ldims(2), ldims(3), & ! [IN]
                              ia, ja,             & ! [IN]
                              igrd(:,:,:),    & ! [IN]
                              jgrd(:,:,:),    & ! [IN]
                              hfact(:,:,:),    & ! [IN]
                              topo_org(:,:),      & ! [IN]
                              topo(:,:)       ) ! [OUT]
        if ( filter_niter > 0 ) then
           call filter_hyperdiff( ia, is, ie, ja, js, je, &
                                  topo(:,:), filter_order, filter_niter )
           call comm_vars8( topo(:,:), 1 )
           call comm_wait ( topo(:,:), 1 )
        end if
 
        do j = 1, ja
        do i = 1, ia
           if ( topo(i,j) > undef + eps ) then ! ignore UNDEF value
              tdiff = ( topo_zsfc(i,j) - topo(i,j) ) * laps
              lst(i,j) = lst(i,j) - tdiff
              do k = 1, lkmax
                 tg(k,i,j) = tg(k,i,j) - tdiff
              end do
           end if
        end do
        end do
 
        if ( urban_do ) then
           do j = 1, ja
           do i = 1, ia
              if ( topo(i,j) > 0.0_rp ) then ! ignore UNDEF value
                 tdiff = ( topo_zsfc(i,j) - topo(i,j) ) * laps
                 ust(i,j) = ust(i,j) - tdiff
              end if
           end do
           end do
        end if
 
     end if
 
 
 
     ! Land water: interpolate over the ocean
     if( use_file_landwater )then
 
        if ( use_waterratio ) then
 
           if ( i_intrp_land_water .ne. i_intrp_off ) then
              do k = 1, ldims(1)
                 work(:,:) = smds_org(k,:,:)
                 select case( i_intrp_land_water )
                 case( i_intrp_mask )
                    lmask = lmask_org
                 case( i_intrp_fill )
                    call make_mask( lmask, work, ldims(2), ldims(3), landdata=.true.)
                 end select
                 call interp_oceanland_data(work, lmask, ldims(2), ldims(3), .true., intrp_iter_max)
                 lmask(:,:) = init_landwater_ratio
                 !replace missing value to init_landwater_ratio
                 call replace_misval_map( work, lmask, ldims(2), ldims(3),  "SMOISDS")
                 smds_org(k,:,:) = work(:,:)
              enddo
           end if
 
           call interp_interp3d( itp_nh,                       & ! [IN]
                                 ldims(1), ldims(2), ldims(3), & ! [IN]
                                 lkmax, lks, lke,              & ! [IN]
                                 ia, ja,                       & ! [IN]
                                 igrd(    :,:,:),          & ! [IN]
                                 jgrd(    :,:,:),          & ! [IN]
                                 hfact(    :,:,:),          & ! [IN]
                                 kgrdl(:,:,:,:,:),          & ! [IN]
                                 vfactl(:,:,:,:,:),          & ! [IN]
                                 smds_org(:,:,:),              & ! [IN]
                                 smds(:,:,:)               ) ! [OUT]
 
           do k = 1, lkmax-1
              strg(k,:,:) = convert_ws2vwc( smds(k,:,:), critical=soilwater_ds2vc_flag )
           end do
 
        else
 
           if ( i_intrp_land_water .ne. i_intrp_off ) then
              do k = 1, ldims(1)
                 work(:,:) = strg_org(k,:,:)
                 select case( i_intrp_land_water )
                 case( i_intrp_mask )
                    lmask = lmask_org
                 case( i_intrp_fill )
                    call make_mask( lmask, work, ldims(2), ldims(3), landdata=.true.)
                 end select
                 call interp_oceanland_data(work, lmask, ldims(2), ldims(3), .true., intrp_iter_max)
                 lmask(:,:) = maskval_strg
                 !replace missing value to init_landwater_ratio
                 call replace_misval_map( work, lmask, ldims(2), ldims(3),  "SMOIS")
                 strg_org(k,:,:) = work(:,:)
              enddo
           end if
 
           call interp_interp3d( itp_nh,                       & ! [IN]
                                 ldims(1), ldims(2), ldims(3), & ! [IN]
                                 lkmax, lks, lke,              & ! [IN]
                                 ia, ja,                       & ! [IN]
                                 igrd(    :,:,:),          & ! [IN]
                                 jgrd(    :,:,:),          & ! [IN]
                                 hfact(    :,:,:),          & ! [IN]
                                 kgrdl(:,:,:,:,:),          & ! [IN]
                                 vfactl(:,:,:,:,:),          & ! [IN]
                                 strg_org(:,:,:),              & ! [IN]
                                 strg(:,:,:)               ) ! [OUT]
        end if
 
        ! replace values over the ocean
        do k = 1, lkmax-1
           call replace_misval_const( strg(k,:,:), maskval_strg, lsmask_nest )
        enddo
 
        if ( filter_niter > 0 ) then
           one3d(:,:,:) = 1.0_rp
           call filter_hyperdiff( lkmax, 1, lkmax-1, ia, is, ie, ja, js, je, &
                                  strg(:,:,:), filter_order, filter_niter, &
                                  limiter_sign = one3d(:,:,:) )
           call comm_vars8( strg(:,:,:), 1 )
           call comm_wait ( strg(:,:,:), 1 )
        end if
 
        do j = 1, ja
        do i = 1, ia
           strg(lkmax,i,j) = strg(lkmax-1,i,j)
        enddo
        enddo
 
     else  ! not read from boundary file
 
        smds(:,:,:) = init_landwater_ratio
        ! conversion from water saturation [fraction] to volumetric water content [m3/m3]
        do k = 1, lkmax
           strg(k,:,:) = convert_ws2vwc( smds(k,:,:), critical=.true. )
        end do
 
     endif ! use_file_waterratio
 
 
     if ( urban_do ) then
        ! copy albedo of land to urban
        do j = 1, ja
        do i = 1, ia
           albu(i,j,:,:) = albg(i,j,:,:)
        enddo
        enddo
     end if
 
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ make_mask()

subroutine mod_realinput::make_mask	(	real(rp), dimension(:,:), intent(out)	gmask,
		real(rp), dimension(:,:), intent(in)	data,
		integer, intent(in)	nx,
		integer, intent(in)	ny,
		logical, intent(in)	landdata
	)

Definition at line 3451 of file mod_realinput.F90.

References scale_const::const_eps, and scale_const::const_undef.

Referenced by land_interporation(), and realinput_surface().

     use scale_const, only: &
        eps => const_eps, &
        undef => const_undef
     implicit none
     real(RP), intent(out)  :: gmask(:,:)
     real(RP), intent(in)   :: data(:,:)
     integer,  intent(in)   :: nx
     integer,  intent(in)   :: ny
     logical,  intent(in)   :: landdata   ! .true. => land data , .false. => ocean data
 
     real(RP)               :: dd
     integer                :: i,j
 
     if( landdata )then
        gmask(:,:) = 1.0_rp  ! gmask=1 will be skip in "interp_OceanLand_data"
        dd         = 0.0_rp
     else
        gmask(:,:) = 0.0_rp  ! gmask=0 will be skip in "interp_OceanLand_data"
        dd         = 1.0_rp
     endif
 
     do j = 1, ny
     do i = 1, nx
        if( abs(data(i,j) - undef) < sqrt(eps) )then
           gmask(i,j) = dd
        endif
     enddo
     enddo
 
     return

Here is the caller graph for this function:

◆ replace_misval_const()

subroutine mod_realinput::replace_misval_const	(	real(rp), dimension(:,:), intent(inout)	data,
		real(rp), intent(in)	maskval,
		real(rp), dimension(:,:), intent(in)	frac_land
	)

Definition at line 3596 of file mod_realinput.F90.

References scale_const::const_eps, scale_atmos_grid_cartesc_index::ia, and scale_atmos_grid_cartesc_index::ja.

Referenced by land_interporation().

     use scale_const, only: &
        eps => const_eps
     implicit none
     real(RP), intent(inout) :: data(:,:)
     real(RP), intent(in)    :: maskval
     real(RP), intent(in)    :: frac_land(:,:)
     integer                 :: i, j
 
     do j = 1, ja
     do i = 1, ia
        if( abs(frac_land(i,j)-0.0_rp) < eps )then ! ocean grid
           data(i,j) = maskval
        endif
     enddo
     enddo
 

Here is the caller graph for this function:

◆ replace_misval_map()

subroutine mod_realinput::replace_misval_map	(	real(rp), dimension(:,:), intent(inout)	data,
		real(rp), dimension(:,:), intent(in)	maskval,
		integer, intent(in)	nx,
		integer, intent(in)	ny,
		character(len=*), intent(in)	elem
	)

Definition at line 3616 of file mod_realinput.F90.

References scale_const::const_eps, scale_const::const_undef, and scale_prc::prc_abort().

Referenced by land_interporation(), and realinput_surface().

     use scale_const, only: &
        eps => const_eps, &
        undef => const_undef
     implicit none
 
     real(RP),         intent(inout) :: data(:,:)
     real(RP),         intent(in)    :: maskval(:,:)
     integer,          intent(in)    :: nx, ny
     character(len=*), intent(in)    :: elem
 
     integer :: i, j
 
     do j = 1, ny
     do i = 1, nx
        if( abs(data(i,j) - undef) < sqrt(eps) )then
           if( abs(maskval(i,j) - undef) < sqrt(eps) )then
              log_error("replace_misval_map",*) "data for mask of "//trim(elem)//"(",i,",",j,") includes missing value."
              log_error_cont(*) "Please check input data of SKINTEMP or SST. "
              call prc_abort
           else
              data(i,j) = maskval(i,j)
           endif
        endif
     enddo
     enddo
 