mod_realinput Module Reference

module REAL input More...

Functions/Subroutines
subroutine, public	realinput_atmos (flg_intrp)
subroutine, public	realinput_surface
subroutine	parentatomsetup (dims, timelen, mdlid, basename_org, filetype, use_file_density_in, serial_in)
	Atmos Setup. More...
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)
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)

Detailed Description

module REAL input

Description

read data from file for real atmospheric simulations

Author

Team SCALE

History

• 2014-04-28 (R.Yoshida) [new]

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
  FILETYPE_ORG	character(len=H_LONG)	''
  BASENAME_ORG	character(len=*)
  BASENAME_ADD_NUM	logical	.false.
  BASENAME_BOUNDARY	character(len=H_LONG)	'boundary_atmos'
  BOUNDARY_TITLE	character(len=H_LONG)	'SCALE-RM BOUNDARY CONDITION for REAL CASE'
  BOUNDARY_UPDATE_DT	real(RP)	0.0_RP	inteval time of boudary data update [s]
  PARENT_MP_TYPE	integer	6	microphysics type of the parent model (number of classes)
  SERIAL_PROC_READ	logical	.true.	read by one MPI process and broadcast
  USE_FILE_DENSITY	logical	.false.	use density data from files

  
  
• 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=*)
  BASENAME_ADD_NUM	logical	.false.
  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
  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=*)
  BASENAME_ADD_NUM	logical	.false.
  BASENAME_BOUNDARY	character(len=H_LONG)	'boundary_atmos'
  BOUNDARY_TITLE	character(len=H_LONG)	'SCALE-RM BOUNDARY CONDITION for REAL CASE'
  BOUNDARY_UPDATE_DT	real(RP)	0.0_RP	inteval time of boudary data update [s]
  INTRP_OCEAN_TEMP	character(len=*)
  INTRP_OCEAN_SFC_TEMP	character(len=*)
  INTRP_ITER_MAX	integer
  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

(

logical, intent(in)

flg_intrp

)

Definition at line 164 of file mod_realinput.f90.

References mod_atmos_vars::dens, scale_grid_index::ia, scale_stdio::io_fid_conf, scale_stdio::io_fid_log, scale_stdio::io_l, scale_stdio::io_lnml, scale_grid_index::ja, scale_grid_index::ka, scale_grid_index::ke, scale_grid_index::ks, mod_atmos_vars::momx, mod_atmos_vars::momy, mod_atmos_vars::momz, parentatomsetup(), scale_process::prc_mpistop(), scale_tracer::qa, mod_atmos_vars::qtrc, mod_atmos_vars::rhot, and scale_tracer::tracer_type.

Referenced by mod_mkinit::interporation_fact().

    use mod_atmos_vars, only: &
         dens, &
         momz, &
         momx, &
         momy, &
         rhot, &
         qtrc
    implicit none

    logical, intent(in)  :: flg_intrp ! flag for interpolation of SBM(S10) from outer bulk-MP model


    namelist / param_mkinit_real_atmos / &
         number_of_files,        &
         number_of_tsteps,       &
         number_of_skip_tsteps,  &
         filetype_org,           &
         basename_org,           &
         basename_add_num,       &
         basename_boundary,      &
         boundary_title,         &
         boundary_update_dt,     &
         parent_mp_type,         &
         serial_proc_read,       &
         use_file_density

    character(len=H_LONG) :: basename_atmos    = ''
    character(len=H_LONG) :: basename_withnum  = ''
    character(len=5)      :: num               = ''

    ! atmos
    real(RP), allocatable :: dens_org(:,:,:,:)
    real(RP), allocatable :: momz_org(:,:,:,:)
    real(RP), allocatable :: momx_org(:,:,:,:)
    real(RP), allocatable :: momy_org(:,:,:,:)
    real(RP), allocatable :: rhot_org(:,:,:,:)
    real(RP), allocatable :: qtrc_org(:,:,:,:,:)

    integer :: mdlid
    integer :: dims(6) ! dims 1-3: normal, 4-6: staggerd

    integer :: totaltimesteps = 1
    integer :: timelen
    integer :: skip_steps
    integer :: ierr
    logical :: flg_bin          ! flag for SBM(S10) is used or not 0-> not used, 1-> used

    integer :: k, i, j, iq, n, ns, ne
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '+++ Module[RealCaseAtmos]/Categ[MKINIT]'

    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_mkinit_real_atmos,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_MKINIT_REAL_ATMOS. Check!'
       call prc_mpistop
    endif
    if( io_lnml ) write(io_fid_log,nml=param_mkinit_real_atmos)

    if( tracer_type == 'SUZUKI10' ) then
       flg_bin = .true.
    else
       flg_bin = .false.
    endif

    if ( filetype_org == "GrADS" ) then
       basename_withnum = basename_org ! namelist file name
       basename_atmos = ""
    else
       if ( number_of_files > 1 .or. basename_add_num ) then
          basename_withnum = trim(basename_org)//"_00000"
       else
          basename_withnum = trim(basename_org)
       end if
       basename_atmos = basename_org
    end if

    call parentatomsetup( dims(:), timelen, mdlid, & ![OUT]
                          basename_withnum,        & ![IN]
                          filetype_org,            & ![IN]
                          use_file_density,        & ![IN]
                          serial_proc_read         ) ![IN]

    if ( boundary_update_dt <= 0.0_rp ) then
       write(*,*) 'xxx BOUNDARY_UPDATE_DT is necessary in real case preprocess'
       call prc_mpistop
    endif

    if ( timelen > 0 ) then
       number_of_tsteps = timelen ! read from file
    endif

    totaltimesteps = number_of_files * number_of_tsteps

    allocate( dens_org(ka,ia,ja,1+number_of_skip_tsteps:totaltimesteps   ) )
    allocate( momz_org(ka,ia,ja,1+number_of_skip_tsteps:totaltimesteps   ) )
    allocate( momx_org(ka,ia,ja,1+number_of_skip_tsteps:totaltimesteps   ) )
    allocate( momy_org(ka,ia,ja,1+number_of_skip_tsteps:totaltimesteps   ) )
    allocate( rhot_org(ka,ia,ja,1+number_of_skip_tsteps:totaltimesteps   ) )
    allocate( qtrc_org(ka,ia,ja,qa,1+number_of_skip_tsteps:totaltimesteps) )

    !--- 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_withnum = trim(basename_atmos)//"_"//num
       else
          basename_withnum = trim(basename_atmos)
       end if

       if( io_l ) write(io_fid_log,*) ' '
       if( io_l ) write(io_fid_log,*) '+++ Target File Name: ',trim(basename_withnum)
       if( io_l ) write(io_fid_log,*) '    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
          if( io_l ) write(io_fid_log,*) '    SKIP'
          cycle
       end if

       skip_steps = max(number_of_skip_tsteps - ns + 1, 0)
       ns = max(ns, number_of_skip_tsteps+1)

       ! read all prepared data
       call parentatominput( dens_org(:,:,:,ns:ne),   &
                             momz_org(:,:,:,ns:ne),   &
                             momx_org(:,:,:,ns:ne),   &
                             momy_org(:,:,:,ns:ne),   &
                             rhot_org(:,:,:,ns:ne),   &
                             qtrc_org(:,:,:,:,ns:ne), &
                             basename_withnum,        &
                             dims(:),                 &
                             mdlid,                   &
                             flg_bin, flg_intrp,      &
                             parent_mp_type,          &
                             number_of_tsteps,        &
                             skip_steps               )
    enddo

    !--- input initial data
    ns = number_of_skip_tsteps + 1  ! skip first several data
    do j = 1, ja
    do i = 1, ia
    do k = ks, ke
       dens(k,i,j) = dens_org(k,i,j,ns)
       momz(k,i,j) = momz_org(k,i,j,ns)
       momx(k,i,j) = momx_org(k,i,j,ns)
       momy(k,i,j) = momy_org(k,i,j,ns)
       rhot(k,i,j) = rhot_org(k,i,j,ns)

       do iq = 1, qa
          qtrc(k,i,j,iq) = qtrc_org(k,i,j,iq,ns)
       enddo
    enddo
    enddo
    enddo

    !--- output boundary data
    totaltimesteps = totaltimesteps - number_of_skip_tsteps ! skip first several data
    call parentatomboundary( dens_org(:,:,:,ns:ne),   &
                             momz_org(:,:,:,ns:ne),   &
                             momx_org(:,:,:,ns:ne),   &
                             momy_org(:,:,:,ns:ne),   &
                             rhot_org(:,:,:,ns:ne),   &
                             qtrc_org(:,:,:,:,ns:ne), &
                             totaltimesteps,          &
                             boundary_update_dt,      &
                             basename_boundary,       &
                             boundary_title           )

    deallocate( dens_org )
    deallocate( momz_org )
    deallocate( momx_org )
    deallocate( momy_org )
    deallocate( rhot_org )
    deallocate( qtrc_org )

    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 354 of file mod_realinput.f90.

References 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_const::const_i_lw, scale_const::const_i_sw, mod_atmos_vars::dens, scale_grid_index::ia, scale_external_io::igrads, scale_stdio::io_fid_conf, scale_stdio::io_fid_log, scale_stdio::io_l, scale_stdio::io_lnml, scale_grid_index::ja, 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_land_grid_index::lkmax, mod_atmos_vars::momx, mod_atmos_vars::momy, mod_atmos_vars::momz, 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, scale_process::prc_mpistop(), mod_atmos_vars::qtrc, mod_atmos_vars::rhot, scale_urban_grid_index::uke, scale_urban_grid_index::uks, 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::interporation_fact().

    use scale_const, only: &
         i_sw => const_i_sw, &
         i_lw => const_i_lw
    use mod_atmos_vars, only: &
         dens, &
         momz, &
         momx, &
         momy, &
         rhot, &
         qtrc
    use mod_land_vars, only: &
         land_temp, &
         land_water, &
         land_sfc_temp, &
         land_sfc_albedo
    use mod_urban_vars, only: &
         urban_sfc_temp, &
         urban_sfc_albedo, &
         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
    use mod_ocean_vars, only: &
         ocean_temp, &
         ocean_sfc_temp, &
         ocean_sfc_albedo, &
         ocean_sfc_z0m, &
         ocean_sfc_z0h, &
         ocean_sfc_z0e
    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 scale_landuse, only: &
         fact_ocean  => landuse_fact_ocean, &
         fact_land   => landuse_fact_land, &
         fact_urban  => landuse_fact_urban
    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."
    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       = 20
    character(len=H_SHORT)   :: soilwater_ds2vc      = 'limit'
    logical                  :: soilwater_ds2vc_flag           ! true: 'critical', false: 'limit'
    logical                  :: elevation_collection = .true.

    namelist / param_mkinit_real_land / &
         number_of_files,        &
         number_of_tsteps,       &
         number_of_skip_tsteps,  &
         filetype_org,           &
         basename_org,           &
         basename_add_num,       &
         use_file_landwater,     &
         init_landwater_ratio,   &
         intrp_land_temp,        &
         intrp_land_water,       &
         intrp_land_sfc_temp,    &
         intrp_iter_max,         &
         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_title,         &
         boundary_update_dt,     &
         intrp_ocean_temp,       &
         intrp_ocean_sfc_temp,   &
         intrp_iter_max,         &
         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=H_LONG) :: basename_withnum  = ''
    character(len=5)      :: num               = ''
    logical               :: serial_proc_read_land
    logical               :: serial_proc_read_ocean

    ! land
    real(RP) :: land_temp_org(lkmax,ia,ja)
    real(RP) :: land_water_org(lkmax,ia,ja)
    real(RP) :: land_sfc_temp_org(ia,ja)
    real(RP) :: land_sfc_albedo_org(ia,ja,2)

    ! 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,2)

    ! 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 :: mdlid_land, mdlid_ocean
    integer :: ldims(3), odims(2)

    integer :: totaltimesteps = 1
    integer :: timelen
    integer :: skip_steps
    integer :: lit
    integer :: lfn
    integer :: ierr

    integer :: k, i, j, n, ns, ne
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '+++ Module[RealCaseSurface]/Categ[MKINIT]'


    ! LAND/URBAN

    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_mkinit_real_land,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_MKINIT_REAL_LAND. Check!'
       call prc_mpistop
    endif
    if( io_lnml ) write(io_fid_log,nml=param_mkinit_real_land)

    filetype_land = filetype_org

    lfn = number_of_skip_tsteps / number_of_tsteps
    if ( filetype_land .ne. "GrADS" .and. ( number_of_files > 1 .or. basename_add_num ) ) then
       write(num,'(I5.5)') lfn
       basename_land = trim(basename_org)//"_"//num
    else
       basename_land = trim(basename_org)
    end if

    serial_land = serial_proc_read

    lit = mod(number_of_skip_tsteps,number_of_tsteps)+1

    !--- read external file
    if( io_l ) write(io_fid_log,*) ' '
    if( io_l ) write(io_fid_log,*) '+++ Target File Name (Land): ',trim(basename_land)
    if( io_l ) write(io_fid_log,*) '    Time Steps in One File: ', number_of_tsteps
    if( io_l ) write(io_fid_log,*) '    Time Step to read: ', lit



    ! OCEAN

    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_mkinit_real_ocean,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_MKINIT_REAL_OCEAN. Check!'
       call prc_mpistop
    endif
    if( io_lnml ) write(io_fid_log,nml=param_mkinit_real_ocean)

    filetype_ocean = filetype_org

    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)
    end if

    select case ( soilwater_ds2vc )
    case ( 'critical' )
       soilwater_ds2vc_flag = .true.
    case ('limit' )
       soilwater_ds2vc_flag = .false.
    case default
      write(*,*) ' xxx Unsupported SOILWATER_DS2CV TYPE:', trim(soilwater_ds2vc)
      call prc_mpistop
    end select

    serial_ocean = serial_proc_read


    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( ocean_temp_org(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,2,1+number_of_skip_tsteps:totaltimesteps   ) )
    allocate( ocean_sfc_z0_org(ia,ja,1+number_of_skip_tsteps:totaltimesteps   ) )

    if ( mdlid_land == igrads .and. ( number_of_files > 1 .or. basename_add_num ) ) then
       write(num,'(I5.5)') lfn
       basename_land = "_"//num
    end if

    if ( mdlid_ocean == igrads ) then
       basename_org = ""
    end if

    !--- 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_ocean = trim(basename_org)//"_"//num
       else
          basename_ocean = trim(basename_org)
       end if

       if( io_l ) write(io_fid_log,*) ' '
       if( io_l ) write(io_fid_log,*) '+++ Target File Name (Ocean): ', trim(basename_ocean)
       if( io_l ) write(io_fid_log,*) '    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
          if( io_l ) write(io_fid_log,*) '    SKIP'
          cycle
       end if

       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,        &
                                land_water_org,       &
                                land_sfc_temp_org,    &
                                land_sfc_albedo_org,  &
                                urban_tc_org,         &
                                urban_qc_org,         &
                                urban_uc_org,         &
                                urban_sfc_temp_org,   &
                                urban_sfc_albedo_org, &
                                ocean_temp_org(:,:,ns:ne),         &
                                ocean_sfc_temp_org(:,:,ns:ne),     &
                                ocean_sfc_albedo_org(:,:,:,ns:ne), &
                                ocean_sfc_z0_org(:,:,ns:ne),       &
                                dens, &
                                momz, &
                                momx, &
                                momy, &
                                rhot, &
                                qtrc, &
                                basename_land,           &
                                basename_ocean,          &
                                mdlid_land, mdlid_ocean, &
                                ldims, odims,            &
                                use_file_landwater,      &
                                init_landwater_ratio,    &
                                intrp_iter_max,          &
                                soilwater_ds2vc_flag,    &
                                elevation_collection,    &
                                number_of_tsteps,        &
                                skip_steps, lit          )

    enddo


    !--- input initial data
    do j = 1, ja
    do i = 1, ia
       land_sfc_temp(i,j) = land_sfc_temp_org(i,j)
       land_sfc_albedo(i,j,i_lw) = land_sfc_albedo_org(i,j,i_lw)
       land_sfc_albedo(i,j,i_sw) = land_sfc_albedo_org(i,j,i_sw)
       do k = 1, lkmax
          land_temp(k,i,j) = land_temp_org(k,i,j)
          land_water(k,i,j) = land_water_org(k,i,j)
       end do

       urban_sfc_temp(i,j) = urban_sfc_temp_org(i,j)
       urban_sfc_albedo(i,j,i_lw) = urban_sfc_albedo_org(i,j,i_lw)
       urban_sfc_albedo(i,j,i_sw) = urban_sfc_albedo_org(i,j,i_sw)
       do k = uks, uke
          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)
       end do
       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
    enddo
    enddo


    ns = number_of_skip_tsteps + 1  ! skip first several data
    do j = 1, ja
    do i = 1, ia
       ocean_temp(i,j)     = ocean_temp_org(i,j,ns)
       ocean_sfc_temp(i,j) = ocean_sfc_temp_org(i,j,ns)
       ocean_sfc_albedo(i,j,i_lw) = ocean_sfc_albedo_org(i,j,i_lw,ns)
       ocean_sfc_albedo(i,j,i_sw) = ocean_sfc_albedo_org(i,j,i_sw,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)
    enddo
    enddo

    do j = 1, ja
    do i = 1, ia
       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)
       atmos_phy_sf_sfc_albedo(i,j,i_lw) = fact_ocean(i,j) * ocean_sfc_albedo(i,j,i_lw) &
                                         + fact_land(i,j) * land_sfc_albedo(i,j,i_lw) &
                                         + fact_urban(i,j) * urban_sfc_albedo(i,j,i_lw)
       atmos_phy_sf_sfc_albedo(i,j,i_sw) = fact_ocean(i,j) * ocean_sfc_albedo(i,j,i_sw) &
                                         + fact_land(i,j) * land_sfc_albedo(i,j,i_sw) &
                                         + fact_urban(i,j) * urban_sfc_albedo(i,j,i_sw)
       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)
    end do
    end do


    !--- output boundary data
    totaltimesteps = totaltimesteps - number_of_skip_tsteps ! skip first several data
    if ( totaltimesteps > 1 ) then
       if ( boundary_update_dt <= 0.0_rp ) then
          write(*,*) 'xxx BOUNDARY_UPDATE_DT is necessary in real case preprocess'
          call prc_mpistop
       endif

       call parentoceanboundary( ocean_temp_org(:,:,ns:ne),         &
                                 ocean_sfc_temp_org(:,:,ns:ne),     &
                                 ocean_sfc_albedo_org(:,:,:,ns:ne), &
                                 ocean_sfc_z0_org(:,:,ns:ne),       &
                                 totaltimesteps,     &
                                 boundary_update_dt, &
                                 basename_boundary,  &
                                 boundary_title      )

    end if

    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:

parentatomsetup()

subroutine mod_realinput::parentatomsetup	(	integer, dimension(6), intent(out)	dims,
		integer, intent(out)	timelen,
		integer, intent(out)	mdlid,
		character(len=*), intent(in)	basename_org,
		character(len=*), intent(in)	filetype,
		logical, intent(in)	use_file_density_in,
		logical, intent(in)	serial_in
	)

Atmos Setup.

Definition at line 759 of file mod_realinput.f90.

References scale_tracer::aq_name, scale_tracer::aq_unit, scale_const::const_eps, scale_const::const_i_lw, scale_const::const_i_sw, scale_const::const_undef, scale_land_grid::grid_lcz, scale_gridtrans::gtrans_rotc, scale_tracer::i_qc, scale_tracer::i_qv, scale_grid_index::ia, scale_external_io::igrads, scale_external_io::inicam, scale_interpolation_nest::intrpnest_domain_compatibility(), scale_interpolation_nest::intrpnest_interp_2d, scale_interpolation_nest::intrpnest_interp_3d, scale_interpolation_nest::intrpnest_interp_fact_latlon(), scale_interpolation_nest::intrpnest_interp_fact_llz(), scale_stdio::io_fid_log, scale_stdio::io_l, scale_external_io::iscale, scale_external_io::iwrfarw, scale_grid_index::ja, scale_grid_index::ka, scale_grid_index::ke, scale_grid_index::ks, land_interporation(), scale_landuse::landuse_frac_land, scale_land_grid_index::lkmax, dc_log::log(), make_mask(), scale_grid_nest::nest_interp_level, mod_realinput_grads::parentatominputgrads(), mod_realinput_nicam::parentatominputnicam(), mod_realinput_scale::parentatominputscale(), mod_realinput_wrfarw::parentatominputwrfarw(), mod_realinput_grads::parentatomopengrads(), mod_realinput_nicam::parentatomopennicam(), mod_realinput_scale::parentatomopenscale(), mod_realinput_wrfarw::parentatomopenwrfarw(), mod_realinput_grads::parentatomsetupgrads(), mod_realinput_nicam::parentatomsetupnicam(), mod_realinput_scale::parentatomsetupscale(), mod_realinput_wrfarw::parentatomsetupwrfarw(), mod_realinput_grads::parentlandinputgrads(), mod_realinput_nicam::parentlandinputnicam(), mod_realinput_scale::parentlandinputscale(), mod_realinput_wrfarw::parentlandinputwrfarw(), mod_realinput_grads::parentlandsetupgrads(), mod_realinput_nicam::parentlandsetupnicam(), mod_realinput_scale::parentlandsetupscale(), mod_realinput_wrfarw::parentlandsetupwrfarw(), mod_realinput_grads::parentoceaninputgrads(), mod_realinput_nicam::parentoceaninputnicam(), mod_realinput_scale::parentoceaninputscale(), mod_realinput_wrfarw::parentoceaninputwrfarw(), mod_realinput_grads::parentoceanopengrads(), mod_realinput_nicam::parentoceanopennicam(), mod_realinput_scale::parentoceanopenscale(), mod_realinput_wrfarw::parentoceanopenwrfarw(), mod_realinput_grads::parentoceansetupgrads(), mod_realinput_nicam::parentoceansetupnicam(), mod_realinput_scale::parentoceansetupscale(), mod_realinput_wrfarw::parentoceansetupwrfarw(), scale_process::prc_ismaster, scale_process::prc_mpistop(), scale_tracer::qa, scale_tracer::qwe, scale_tracer::qws, replace_misval_map(), scale_time::time_nowdate, and scale_tracer::tracer_type.

Referenced by realinput_atmos().

    use scale_external_io, only: &
         iscale, &
         iwrfarw, &
         inicam, &
         igrads
    use mod_realinput_scale, only: &
         parentatomsetupscale
    use mod_realinput_wrfarw, only: &
         parentatomsetupwrfarw
    use mod_realinput_nicam, only: &
         parentatomsetupnicam
    use mod_realinput_grads, only: &
         parentatomsetupgrads
    implicit none

    integer,          intent(out) :: dims(6)
    integer,          intent(out) :: timelen
    integer,          intent(out) :: mdlid
    character(len=*), intent(in)  :: basename_org
    character(len=*), intent(in)  :: filetype
    logical,          intent(in)  :: serial_in             ! read by a serial process
    logical,          intent(in)  :: use_file_density_in   ! use density data from files

    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '+++ ScaleLib/IO[RealinputAtmos]/Categ[Setup]'

    serial = serial_in
    if( serial ) then
       if( prc_ismaster ) then
          do_read_atom = .true.
       else
          do_read_atom = .false.
       endif
    else
       do_read_atom = .true.
    endif

    select case(trim(filetype))
    case('SCALE-RM')

       mdlid = iscale
       do_read_atom = .true.
       serial = .false.
       call parentatomsetupscale( dims ) ! (out)
       update_coord = .false.
       use_file_density = use_file_density_in
       use_temp = .false.
       rotate = .false.
       timelen = -1

    case('WRF-ARW')

       mdlid = iwrfarw
       if ( do_read_atom ) call parentatomsetupwrfarw( dims, timelen, & ! (out)
                                                       basename_org   ) ! (in)
       update_coord = .true.
       use_file_density = .false.
       use_temp = .true.
       rotate = .true.

    case('NICAM-NETCDF')

       mdlid = inicam
       if ( do_read_atom ) call parentatomsetupnicam( dims, timelen, & ! (out)
                                                      basename_org   ) ! (in)
       update_coord = .false.
       use_file_density = .false.
       use_temp = .true.
       rotate = .true.

    case('GrADS')

       mdlid = igrads
       if ( do_read_atom ) call parentatomsetupgrads( dims,        & ! (out)
                                                      basename_org ) ! (in)
       update_coord = .true.
       use_file_density = .false.
       use_temp = .true.
       rotate = .true.
       timelen = -1

    case default

       write(*,*) ' xxx Unsupported FILE TYPE:', trim(filetype)
       call prc_mpistop

    endselect

    if( serial ) then
       call comm_bcast( dims(:), 6 )
       call comm_bcast( timelen )
    endif

    if( io_l ) write(io_fid_log,*) '+++ Horizontal Interpolation Level:', &
                                    nest_interp_level
    itp_nh = int( nest_interp_level )
    itp_nv = 2

    allocate( hfact(        ia, ja, itp_nh         ) )
    allocate( vfact( ka,    ia, ja, itp_nh, itp_nv ) )
    allocate( igrd(        ia, ja, itp_nh         ) )
    allocate( jgrd(        ia, ja, itp_nh         ) )
    allocate( kgrd( ka,    ia, ja, itp_nh, itp_nv ) )
    allocate( ncopy(        ia, ja, itp_nh         ) )

    allocate( lon_org(            dims(2), dims(3) ) )
    allocate( lat_org(            dims(2), dims(3) ) )
    allocate( cz_org( dims(1)+2, dims(2), dims(3) ) )

    allocate( velz_org( dims(1)+2, dims(2), dims(3) ) )
    allocate( velx_org( dims(1)+2, dims(2), dims(3) ) )
    allocate( vely_org( dims(1)+2, dims(2), dims(3) ) )
    allocate( pott_org( dims(1)+2, dims(2), dims(3) ) )
    allocate( temp_org( dims(1)+2, dims(2), dims(3) ) )
    allocate( pres_org( dims(1)+2, dims(2), dims(3) ) )
    allocate( qtrc_org( dims(1)+2, dims(2), dims(3), qa ) )
    allocate( dens_org( dims(1)+2, dims(2), dims(3) ) )

    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,2), intent(out)	albg,
		real(rp), dimension(ia,ja), intent(out)	ust,
		real(rp), dimension(ia,ja,2), 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
	)

Definition at line 2301 of file mod_realinput.f90.

References scale_const::const_i_lw, scale_const::const_i_sw, scale_const::const_laps, scale_const::const_undef, mod_land_vars::convert_ws2vwc(), scale_grid_index::ia, scale_interpolation_nest::intrpnest_interp_2d, scale_interpolation_nest::intrpnest_interp_3d, scale_interpolation_nest::intrpnest_interp_fact_latlon(), scale_interpolation_nest::intrpnest_interp_fact_llz(), scale_grid_index::ja, scale_land_grid_index::lkmax, make_mask(), scale_process::prc_mpistop(), replace_misval_const(), replace_misval_map(), and scale_topography::topo_zsfc.

Referenced by parentatomsetup().

    use scale_process, only: &
         prc_mpistop
    use scale_const, only: &
         undef => const_undef, &
         i_sw => const_i_sw, &
         i_lw => const_i_lw, &
         laps => const_laps
    use scale_interpolation_nest, only: &
         intrpnest_interp_fact_llz, &
         intrpnest_interp_fact_latlon, &
         intrpnest_interp_3d, &
         intrpnest_interp_2d
    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,2)
    real(RP), intent(out)   :: ust(ia,ja)
    real(RP), intent(out)   :: albu(ia,ja,2)
    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

    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

    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
          write(*,*) 'xxx INTRP_LAND_SFC_TEMP is invalid.'
          call prc_mpistop
       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
    !      write(*,*) 'xxx INTRP_URB_SFC_TEMP is invalid.'
    !      call PRC_MPIstop
    !   end select
    !   call interp_OceanLand_data(ust_org, lmask, ldims(2), ldims(3), .true., intrp_iter_max)
    !end if

    ! interpolation facter between outer land grid and ocean grid
    call intrpnest_interp_fact_latlon( hfact_l(:,:,:),               & ! [OUT]
                                       igrd_l(:,:,:), jgrd_l(:,:,:), & ! [OUT]
                                       llat_org(:,:), llon_org(:,:), & ! [IN]
                                       ldims(2), ldims(3),           & ! [IN]
                                       olat_org(:,:), olon_org(:,:), & ! [IN]
                                       odims(1), odims(2)            ) ! [IN]


    ! sst on land grid
    call intrpnest_interp_2d( sst_land(:,:), sst_org(:,:), hfact_l(:,:,:),      &
                              igrd_l(:,:,:), jgrd_l(:,:,:), ldims(2), ldims(3) )
    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 ( ust_org(i,j) == undef ) ust_org(i,j) = lst_org(i,j)
!       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_lw) == undef ) albg_org(i,j,i_lw) = 0.03_rp  ! emissivity of general ground surface : 0.95-0.98
       if ( albg_org(i,j,i_sw) == undef ) albg_org(i,j,i_sw) = 0.22_rp
    end do
    end do

    ! 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 intrpnest_interp_fact_llz( hfact(:,:,:),               & ! [OUT]
                                    vfactl(:,:,:,:,:),           & ! [OUT]
                                    kgrdl(:,:,:,:,:),           & ! [OUT]
                                    igrd(:,:,:), jgrd(:,:,:),  & ! [OUT]
                                    ncopy(:,:,:),               & ! [OUT]
                                    lcz_3d(:,:,:),               & ! [IN]
                                    lat(:,:), lon(:,:),   & ! [IN]
                                    1, lkmax, ia, ja,             & ! [IN]
                                    lz3d_org(:,:,:),              & ! [IN]
                                    llat_org(:,:), llon_org(:,:), & ! [IN]
                                    ldims(1), ldims(2), ldims(3), & ! [IN]
                                    landgrid=.true.               ) ! [IN]

    call intrpnest_interp_2d( lst(:,:),   lst_org(:,:),   hfact(:,:,:), &
                              igrd(:,:,:), jgrd(:,:,:), ia, ja )
    call intrpnest_interp_2d( ust(:,:),   ust_org(:,:),   hfact(:,:,:), &
                              igrd(:,:,:), jgrd(:,:,:), ia, ja )
!    call INTRPNEST_interp_2d( skinw(:,:), skinw_org(:,:), hfact(:,:,:), &
!                              igrd(:,:,:), jgrd(:,:,:), IA, JA )
!    call INTRPNEST_interp_2d( snowq(:,:), snowq_org(:,:), hfact(:,:,:), &
!                              igrd(:,:,:), jgrd(:,:,:), IA, JA )
!    call INTRPNEST_interp_2d( snowt(:,:), snowt_org(:,:), hfact(:,:,:), &
!                              igrd(:,:,:), jgrd(:,:,:), IA, JA )
    call intrpnest_interp_2d( albg(:,:,i_lw), albg_org(:,:,i_lw), hfact(:,:,:), &
                              igrd(:,:,:), jgrd(:,:,:), ia, ja )
    call intrpnest_interp_2d( albg(:,:,i_sw), albg_org(:,:,i_sw), hfact(:,:,:), &
                              igrd(:,:,:), jgrd(:,:,:), ia, ja )

    call intrpnest_interp_3d( tg(:,:,:),     &
                              tg_org(:,:,:),     &
                              hfact(:,:,:),     &
                              vfactl(:,:,:,:,:), &
                              kgrdl(:,:,:,:,:), &
                              igrd(:,:,:),     &
                              jgrd(:,:,:),     &
                              ia, ja, 1, lkmax-1 )

    do j = 1, ja
    do i = 1, ia
       tg(lkmax,i,j) = tg(lkmax-1,i,j)
    enddo ! i
    enddo ! j

    ! replace values over the ocean
    do k = 1, lkmax
       call replace_misval_const( tg(k,:,:), maskval_tg, lsmask_nest )
    enddo


    ! elevation collection
    if ( elevation_collection ) then
       call intrpnest_interp_2d( topo(:,:),   topo_org(:,:),   hfact(:,:,:), &
                                 igrd(:,:,:), jgrd(:,:,:), ia, ja )

       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
             lst(i,j) = lst(i,j) - tdiff
             ust(i,j) = ust(i,j) - tdiff
             do k = 1, lkmax
                tg(k,i,j) = tg(k,i,j) - tdiff
             end do
          end if
       end do
       end do
    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 intrpnest_interp_3d( smds(:,:,:),     &
                                    smds_org(:,:,:),     &
                                    hfact(:,:,:),     &
                                    vfactl(:,:,:,:,:), &
                                    kgrdl(:,:,:,:,:), &
                                    igrd(:,:,:),     &
                                    jgrd(:,:,:),     &
                                    ia, ja, 1, lkmax-1   )
          do k = 1, lkmax
             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 intrpnest_interp_3d( strg(:,:,:),     &
                                    strg_org(:,:,:),     &
                                    hfact(:,:,:),     &
                                    vfactl(:,:,:,:,:), &
                                    kgrdl(:,:,:,:,:), &
                                    igrd(:,:,:),     &
                                    jgrd(:,:,:),     &
                                    ia, ja, 1, lkmax-1   )
          ! interpolation
          do j = 1, ja
          do i = 1, ia
             strg(lkmax,i,j) = strg(lkmax-1,i,j)
          enddo
          enddo

       end if

       ! replace values over the ocean
       do k = 1, lkmax
          call replace_misval_const( strg(k,:,:), maskval_strg, lsmask_nest )
       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


    ! copy albedo of land to urban
    do j = 1, ja
    do i = 1, ia
       albu(i,j,:) = albg(i,j,:)
    enddo
    enddo


    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 2633 of file mod_realinput.f90.

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

Referenced by land_interporation(), and parentatomsetup().

    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 2774 of file mod_realinput.f90.

References scale_const::const_eps, scale_grid_index::ia, and scale_grid_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(*), intent(in)	elem
	)

Definition at line 2794 of file mod_realinput.f90.

References scale_const::const_eps, scale_const::const_pi, scale_const::const_undef, scale_tracer::i_nc, scale_tracer::i_ng, scale_tracer::i_ni, scale_tracer::i_nr, scale_tracer::i_ns, scale_tracer::i_qc, scale_tracer::i_qg, scale_tracer::i_qi, scale_tracer::i_qr, scale_tracer::i_qs, and scale_process::prc_mpistop().

Referenced by land_interporation(), and parentatomsetup().

    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(*), 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
             write(*,*) "Data for mask has missing value. ",trim(elem),i,j
             call prc_mpistop
          else
             data(i,j) = maskval(i,j)
          endif
       endif
    enddo
    enddo

Here is the call graph for this function:

Here is the caller graph for this function: