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 72 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_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(:,:) = 1.0_rp
     landuse_fact_land(:,:) = 0.0_rp
     landuse_fact_urban(:,:) = 0.0_rp
 
 
     if    ( 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 149 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 248 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