scale_landuse Module Reference

module LANDUSE More...

Functions/Subroutines
subroutine, public	landuse_setup
	Setup. More...
subroutine, public	landuse_calc_fact
subroutine, public	landuse_write
	Write landuse data. More...

Variables
real(rp), dimension(:,:), allocatable, public	landuse_fact_ocean
	ocean factor More...
real(rp), dimension(:,:), allocatable, public	landuse_fact_land
	land factor More...
real(rp), dimension(:,:), allocatable, public	landuse_fact_urban
	urban factor More...
real(rp), dimension(:,:), allocatable, public	landuse_frac_land
	land fraction More...
real(rp), dimension(:,:), allocatable, public	landuse_frac_lake
	lake fraction More...
real(rp), dimension(:,:), allocatable, public	landuse_frac_urban
	urban fraction More...
integer, public	landuse_pft_mosaic = 2
	number of PFT mosaic More...
integer, public	landuse_pft_nmax = 15
	number of plant functional type(PFT) More...
real(rp), dimension(:,:,:), allocatable, public	landuse_frac_pft
	fraction of PFT for each mosaic More...
integer, dimension(:,:,:), allocatable, public	landuse_index_pft
	index of PFT for each mosaic More...

Detailed Description

module LANDUSE

Description

Land use category module Manage land/lake/urban/PFT fraction and PFT index

Author

Team SCALE

Function/Subroutine Documentation

landuse_setup()

subroutine, public scale_landuse::landuse_setup

(

)

Setup.

Definition at line 73 of file scale_landuse.F90.

References 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, landuse_calc_fact(), landuse_fact_land, landuse_fact_ocean, landuse_fact_urban, landuse_frac_lake, landuse_frac_land, landuse_frac_pft, landuse_frac_urban, landuse_index_pft, landuse_pft_mosaic, landuse_pft_nmax, and scale_process::prc_mpistop().

Referenced by mod_rm_driver::scalerm(), and mod_rm_prep::scalerm_prep().

    use scale_process, only: &
       prc_mpistop
    implicit none

    namelist / param_landuse / &
       landuse_in_basename,  &
       landuse_out_basename, &
       landuse_out_dtype,    &
       landuse_pft_mosaic,   &
       landuse_pft_nmax,     &
       landuse_allocean,     &
       landuse_allland,      &
       landuse_allurban,     &
       landuse_mosaicworld

    integer :: ierr
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '++++++ Module[LANDUSE] / Categ[COUPLER] / Origin[SCALElib]'

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

    allocate( landuse_frac_land(ia,ja) )
    allocate( landuse_frac_lake(ia,ja) )
    allocate( landuse_frac_urban(ia,ja) )
    landuse_frac_land(:,:) = 0.0_rp
    landuse_frac_lake(:,:) = 0.0_rp
    landuse_frac_urban(:,:) = 0.0_rp

    allocate( landuse_index_pft(ia,ja,landuse_pft_mosaic) )
    allocate( landuse_frac_pft(ia,ja,landuse_pft_mosaic) )
    landuse_frac_pft(:,:,:) = 0.0_rp
    landuse_frac_pft(:,:,1) = 1.0_rp ! tentative, mosaic is off
    landuse_index_pft(:,:,:) = 1      ! default

    allocate( landuse_fact_ocean(ia,ja) )
    allocate( landuse_fact_land(ia,ja) )
    allocate( landuse_fact_urban(ia,ja) )
    landuse_fact_ocean(:,:) = 0.0_rp
    landuse_fact_land(:,:) = 0.0_rp
    landuse_fact_urban(:,:) = 0.0_rp


    if    ( landuse_allocean ) then
       if( io_l ) write(io_fid_log,*) '*** Assume all grids are ocean'
       call landuse_calc_fact
    elseif( landuse_allland ) then
       if( io_l ) write(io_fid_log,*) '*** Assume all grids are land'
       landuse_frac_land(:,:) = 1.0_rp
       call landuse_calc_fact
    elseif( landuse_allurban ) then
       if( io_l ) write(io_fid_log,*) '*** Assume all grids are land'
       landuse_frac_land(:,:) = 1.0_rp
       if( io_l ) write(io_fid_log,*) '*** Assume all lands are urban'
       landuse_frac_urban(:,:) = 1.0_rp
       call landuse_calc_fact
    elseif( landuse_mosaicworld ) then
       if( io_l ) write(io_fid_log,*) '*** Assume all grids have ocean, land, and urban'
       landuse_frac_land(:,:) = 0.5_rp
       landuse_frac_urban(:,:) = 0.5_rp
       call landuse_calc_fact
    else
       ! read from file
       call landuse_read
    endif

    return

Here is the call graph for this function:

Here is the caller graph for this function:

landuse_calc_fact()

subroutine, public scale_landuse::landuse_calc_fact

(

)

Definition at line 154 of file scale_landuse.F90.

References scale_stdio::io_fid_log, scale_stdio::io_l, landuse_fact_land, landuse_fact_ocean, landuse_fact_urban, landuse_frac_lake, landuse_frac_land, landuse_frac_pft, landuse_frac_urban, landuse_index_pft, and landuse_pft_mosaic.

Referenced by mod_cnvlanduse::cnvlanduse(), mod_mkinit::interporation_fact(), and landuse_setup().

    implicit none
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '+++ calculate landuse factor'

    ! tentative treatment for lake fraction
    landuse_frac_land(:,:) = landuse_frac_land(:,:) * ( 1.0_rp - landuse_frac_lake(:,:) )

    ! make factors
    landuse_fact_ocean(:,:) = ( 1.0_rp - landuse_frac_land(:,:) )
    landuse_fact_land(:,:) = (          landuse_frac_land(:,:) ) * ( 1.0_rp - landuse_frac_urban(:,:) )
    landuse_fact_urban(:,:) = (          landuse_frac_land(:,:) ) * (          landuse_frac_urban(:,:) )

    return

Here is the caller graph for this function:

landuse_write()

subroutine, public scale_landuse::landuse_write

(

)

Write landuse data.

Definition at line 253 of file scale_landuse.F90.

References scale_stdio::io_fid_log, scale_stdio::io_l, landuse_fact_land, landuse_fact_ocean, landuse_fact_urban, landuse_frac_lake, landuse_frac_land, landuse_frac_pft, landuse_frac_urban, and landuse_pft_mosaic.

Referenced by mod_cnvlanduse::cnvlanduse(), and mod_mkinit::mkinit().

    use scale_fileio, only: &
       fileio_write
    implicit none

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

    character(len=H_SHORT) :: varname
    integer :: p
    !---------------------------------------------------------------------------

    if ( landuse_out_basename /= '' ) then

       if( io_l ) write(io_fid_log,*)
       if( io_l ) write(io_fid_log,*) '*** Output landuse file ***'

       call fileio_write( landuse_frac_land(:,:), landuse_out_basename, landuse_out_title, & ! [IN]
                          'FRAC_LAND',  'LAND fraction',  '1', 'XY',   landuse_out_dtype, & ! [IN]
                          nozcoord=.true. )

       call fileio_write( landuse_frac_lake(:,:), landuse_out_basename, landuse_out_title, & ! [IN]
                          'FRAC_LAKE',  'LAKE fraction',  '1', 'XY',   landuse_out_dtype, & ! [IN]
                          nozcoord=.true. )
       call fileio_write( landuse_frac_urban(:,:), landuse_out_basename, landuse_out_title, & ! [IN]
                          'FRAC_URBAN', 'URBAN fraction', '1', 'XY',   landuse_out_dtype, & ! [IN]
                          nozcoord=.true. )

       do p = 1, landuse_pft_mosaic
          write(varname,'(A8,I1.1)') 'FRAC_PFT', p

          call fileio_write( landuse_frac_pft(:,:,p), landuse_out_basename, landuse_out_title, & ! [IN]
                             varname, 'PFT fraction', '1', 'XY',          landuse_out_dtype, & ! [IN]
                             nozcoord=.true. )

          write(varname,'(A9,I1.1)') 'INDEX_PFT', p
          temp(:,:) = real(LANDUSE_index_PFT(:,:,p),kind=rp)

          call fileio_write( temp(:,:), landuse_out_basename,   landuse_out_title, & ! [IN]
                             varname, 'PFT index', '1', 'XY', landuse_out_dtype, & ! [IN]
                             nozcoord=.true. )
       enddo

       call fileio_write( landuse_fact_ocean(:,:), landuse_out_basename, landuse_out_title, & ! [IN]
                          'FRAC_OCEAN_abs', 'absolute OCEAN fraction',  '0-1', 'XY',   landuse_out_dtype, & ! [IN]
                          nozcoord=.true. )

       call fileio_write( landuse_fact_land(:,:), landuse_out_basename, landuse_out_title, & ! [IN]
                          'FRAC_LAND_abs ', 'absolute LAND fraction',  '0-1', 'XY',   landuse_out_dtype, & ! [IN]
                          nozcoord=.true. )

       call fileio_write( landuse_fact_urban(:,:), landuse_out_basename, landuse_out_title, & ! [IN]
                          'FRAC_URBAN_abs', 'absolute URBAN fraction',  '0-1', 'XY',   landuse_out_dtype, & ! [IN]
                          nozcoord=.true. )
    endif

    return

Here is the caller graph for this function:

Variable Documentation

landuse_fact_ocean

real(rp), dimension(:,:), allocatable, public scale_landuse::landuse_fact_ocean

ocean factor

Definition at line 35 of file scale_landuse.F90.

Referenced by mod_atmos_phy_rd_driver::atmos_phy_rd_driver(), mod_cpl_admin::cpl_admin_setup(), mod_cpl_vars::cpl_getsfc_atm(), landuse_calc_fact(), landuse_setup(), landuse_write(), scale_ocean_phy_file::ocean_phy_file_setup(), scale_ocean_phy_slab::ocean_phy_slab_setup(), scale_ocean_sfc_slab::ocean_sfc_slab_setup(), and mod_realinput::realinput_surface().

  real(RP), public, allocatable :: landuse_fact_ocean(:,:)

landuse_fact_land

real(rp), dimension (:,:), allocatable, public scale_landuse::landuse_fact_land

land factor

Definition at line 36 of file scale_landuse.F90.

Referenced by mod_atmos_phy_rd_driver::atmos_phy_rd_driver(), mod_cpl_admin::cpl_admin_setup(), mod_cpl_vars::cpl_getsfc_atm(), scale_land_phy_slab::land_phy_slab_setup(), scale_land_sfc_slab::land_sfc_slab_setup(), landuse_calc_fact(), landuse_setup(), landuse_write(), and mod_realinput::realinput_surface().

  real(RP), public, allocatable :: landuse_fact_land (:,:)

landuse_fact_urban

real(rp), dimension(:,:), allocatable, public scale_landuse::landuse_fact_urban

urban factor

Definition at line 37 of file scale_landuse.F90.

Referenced by mod_atmos_phy_rd_driver::atmos_phy_rd_driver(), mod_cpl_admin::cpl_admin_setup(), mod_cpl_vars::cpl_getsfc_atm(), landuse_calc_fact(), landuse_setup(), landuse_write(), mod_realinput::realinput_surface(), and scale_urban_phy_slc::urban_phy_slc_setup().

  real(RP), public, allocatable :: landuse_fact_urban(:,:)

landuse_frac_land

real(rp), dimension (:,:), allocatable, public scale_landuse::landuse_frac_land

land fraction

Definition at line 39 of file scale_landuse.F90.

Referenced by mod_cnvlanduse::cnvlanduse(), scale_history::hist_switch(), mod_mkinit::interporation_fact(), landuse_calc_fact(), landuse_setup(), landuse_write(), and mod_realinput::parentatomsetup().

  real(RP), public, allocatable :: landuse_frac_land (:,:)

landuse_frac_lake

real(rp), dimension (:,:), allocatable, public scale_landuse::landuse_frac_lake

lake fraction

Definition at line 40 of file scale_landuse.F90.

Referenced by mod_cnvlanduse::cnvlanduse(), landuse_calc_fact(), landuse_setup(), and landuse_write().

  real(RP), public, allocatable :: landuse_frac_lake (:,:)

landuse_frac_urban

real(rp), dimension(:,:), allocatable, public scale_landuse::landuse_frac_urban

urban fraction

Definition at line 41 of file scale_landuse.F90.

Referenced by mod_cnvlanduse::cnvlanduse(), mod_mkinit::interporation_fact(), landuse_calc_fact(), landuse_setup(), and landuse_write().

  real(RP), public, allocatable :: landuse_frac_urban(:,:)

landuse_pft_mosaic

integer, public scale_landuse::landuse_pft_mosaic = 2

number of PFT mosaic

Definition at line 43 of file scale_landuse.F90.

Referenced by mod_cnvlanduse::cnvlanduse(), landuse_calc_fact(), landuse_setup(), and landuse_write().

  integer,  public              :: landuse_pft_mosaic = 2

landuse_pft_nmax

integer, public scale_landuse::landuse_pft_nmax = 15

number of plant functional type(PFT)

Definition at line 44 of file scale_landuse.F90.

Referenced by mod_cnvlanduse::cnvlanduse(), mod_land_vars::land_vars_external_in(), mod_land_vars::land_vars_setup(), and landuse_setup().

  integer,  public              :: landuse_pft_nmax   = 15

landuse_frac_pft

real(rp), dimension (:,:,:), allocatable, public scale_landuse::landuse_frac_pft

fraction of PFT for each mosaic

Definition at line 46 of file scale_landuse.F90.

Referenced by mod_cnvlanduse::cnvlanduse(), landuse_calc_fact(), landuse_setup(), and landuse_write().

  real(RP), public, allocatable :: landuse_frac_pft (:,:,:)

landuse_index_pft

integer, dimension(:,:,:), allocatable, public scale_landuse::landuse_index_pft

index of PFT for each mosaic

Definition at line 47 of file scale_landuse.F90.

Referenced by mod_cnvlanduse::cnvlanduse(), mod_land_vars::land_vars_setup(), landuse_calc_fact(), and landuse_setup().

  integer,  public, allocatable :: landuse_index_pft(:,:,:)