module Convert LandUseIndex More...

Functions/Subroutines
subroutine, public	cnvlanduse_setup
	Setup. More...

subroutine, public	cnvlanduse
	Driver. More...

Variables
logical, public	cnvlanduse_donothing

logical, public	cnvlanduse_useglccv2 = .false.

logical, public	cnvlanduse_uselu100m = .false.

logical, public	cnvlanduse_usejibis = .false.

Detailed Description

module Convert LandUseIndex

Description: subroutines for preparing landuse index data (convert from external file)

Author: Team SCALE

NAMELIST

PARAM_CNVLANDUSE

name	type	default value	comment
CNVLANDUSE_NAME	character(len=H_SHORT)	'NONE'	keep backward compatibility
CNVLANDUSE_USEGLCCV2	logical	.false.
CNVLANDUSE_USELU100M	logical	.false.
CNVLANDUSE_USEJIBIS	logical	.false.
CNVLANDUSE_UNITTILE_DDEG	real(RP)	0.0_RP	dx for unit tile [deg]
CNVLANDUSE_OVERSAMPLING_FACTOR	real(RP)	2.0_RP	factor of min. dx against the unit tile

PARAM_CNVLANDUSE_GLCCv2

name	type	default value	comment
GLCCV2_IN_CATALOGUE	character(len=H_LONG)	''	metadata files for GLCCv2
GLCCV2_IN_DIR	character(len=H_LONG)	''	directory contains GLCCv2 files (GrADS format)
LIMIT_URBAN_FRACTION	real(RP)	1.0_RP	fraction limiter for urban area

PARAM_CNVLANDUSE_LU100M

name	type	default value	comment
LU100M_IN_CATALOGUE	character(len=H_LONG)	''	metadata files for LU100M
LU100M_IN_DIR	character(len=H_LONG)	''	directory contains LU100M files (GrADS format)
LIMIT_URBAN_FRACTION	real(RP)	1.0_RP	fraction limiter for urban area

History Output: No history output

Function/Subroutine Documentation

◆ cnvlanduse_setup()

subroutine, public mod_cnvlanduse::cnvlanduse_setup ( )

Setup.

Definition at line 59 of file mod_cnvlanduse.f90.

References cnvlanduse_donothing, cnvlanduse_useglccv2, cnvlanduse_usejibis, cnvlanduse_uselu100m, scale_const::const_d2r, scale_stdio::io_fid_conf, scale_stdio::io_fid_log, scale_stdio::io_fid_nml, scale_stdio::io_l, scale_stdio::io_nml, scale_process::prc_mpistop(), scale_grid_real::real_dlat, and scale_grid_real::real_dlon.

Referenced by mod_convert::convert_setup().

     use scale_process, only: &
        prc_mpistop
     use scale_const, only: &
        d2r => const_d2r
     use scale_comm, only: &
        comm_horizontal_min
     use scale_grid_real, only: &
        real_dlat, &
        real_dlon
     implicit none
 
     character(len=H_SHORT) :: CNVLANDUSE_name = 'NONE' ! keep backward compatibility
 
     namelist / param_cnvlanduse / &
        cnvlanduse_name,                &
        cnvlanduse_useglccv2,           &
        cnvlanduse_uselu100m,           &
        cnvlanduse_usejibis,            &
        cnvlanduse_unittile_ddeg,       &
        cnvlanduse_oversampling_factor
 
     real(RP) :: drad(IA,JA)
     real(RP) :: drad_min
 
     integer  :: ierr
     !---------------------------------------------------------------------------
 
     if( io_l ) write(io_fid_log,*)
     if( io_l ) write(io_fid_log,*) '++++++ Module[convert landuseindex] / Categ[preprocess] / Origin[SCALE-RM]'
     !--- read namelist
     rewind(io_fid_conf)
     read(io_fid_conf,nml=param_cnvlanduse,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_CNVLANDUSE. Check!'
        call prc_mpistop
     endif
     if( io_nml ) write(io_fid_nml,nml=param_cnvlanduse)
 
     select case(cnvlanduse_name)
     case('NONE')
        ! do nothing
     case('GLCCv2')
        cnvlanduse_useglccv2 = .true.
        cnvlanduse_uselu100m = .false.
        cnvlanduse_usejibis  = .false.
     case('LU100M')
        cnvlanduse_useglccv2 = .false.
        cnvlanduse_uselu100m = .true.
        cnvlanduse_usejibis  = .false.
     case('COMBINE')
        cnvlanduse_useglccv2 = .true.
        cnvlanduse_uselu100m = .true.
        cnvlanduse_usejibis  = .true.
     case default
        write(*,*) 'xxx Unsupported TYPE:', trim(cnvlanduse_name)
        call prc_mpistop
     endselect
 
     cnvlanduse_donothing = .true.
 
     if ( cnvlanduse_useglccv2 ) then
        cnvlanduse_donothing = .false.
        if( io_l ) write(io_fid_log,*) '*** Use GLCC ver.2, global 30 arcsec. data'
        if ( cnvlanduse_uselu100m ) then
           if( io_l ) write(io_fid_log,*) '*** Use KSJ landuse 100m data for Japan region'
           if( io_l ) write(io_fid_log,*) '*** Overwrite Japan region'
           if ( cnvlanduse_usejibis ) then
              if( io_l ) write(io_fid_log,*) '*** Use J-IBIS map 100m data for Japan region'
              if( io_l ) write(io_fid_log,*) '*** Overwrite Japan region (PFT only)'
           endif
        endif
     elseif ( cnvlanduse_uselu100m ) then
        cnvlanduse_donothing = .false.
        if( io_l ) write(io_fid_log,*) '*** Use KSJ landuse 100m data, Japan region only'
        if ( cnvlanduse_usejibis ) then
           if( io_l ) write(io_fid_log,*) '*** Use J-IBIS map 100m data for Japan region'
           if( io_l ) write(io_fid_log,*) '*** Overwrite Japan region (PFT only)'
        endif
     endif
 
     if ( cnvlanduse_donothing ) then
        if( io_l ) write(io_fid_log,*) '*** Do nothing for landuse index'
     else
        drad(:,:) = min( real_dlat(:,:), real_dlon(:,:) )
        call comm_horizontal_min( drad_min, drad(:,:) )
 
        if ( cnvlanduse_unittile_ddeg > 0.0_rp ) then
           cnvlanduse_oversampling_factor = ( drad_min / d2r ) / cnvlanduse_unittile_ddeg
        endif
        cnvlanduse_oversampling_factor = max( 1.0_rp, cnvlanduse_oversampling_factor )
        cnvlanduse_unittile_ddeg       = ( drad_min / d2r ) / cnvlanduse_oversampling_factor
 
        if( io_l ) write(io_fid_log,*) '*** The size of tile [deg] = ', cnvlanduse_unittile_ddeg
        if( io_l ) write(io_fid_log,*) '*** oversampling factor    = ', cnvlanduse_oversampling_factor
     endif
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ cnvlanduse()

subroutine, public mod_cnvlanduse::cnvlanduse ( )

Driver.

Definition at line 163 of file mod_cnvlanduse.f90.

Referenced by mod_convert::convert().

     use scale_process, only: &
        prc_mpistop
     use scale_landuse, only: &
        landuse_calc_fact,  &
        landuse_write
     implicit none
     !---------------------------------------------------------------------------
 
     if ( cnvlanduse_donothing ) then
        if( io_l ) write(io_fid_log,*)
        if( io_l ) write(io_fid_log,*) '++++++ SKIP  CONVERT LANDUSE DATA ++++++'
     else
        if( io_l ) write(io_fid_log,*)
        if( io_l ) write(io_fid_log,*) '++++++ START CONVERT LANDUSE DATA ++++++'
 
        if ( cnvlanduse_useglccv2 ) then
           call cnvlanduse_glccv2
        endif
 
        if ( cnvlanduse_uselu100m ) then
           call cnvlanduse_lu100m
        endif
 
        if ( cnvlanduse_usejibis ) then
           call cnvlanduse_jibis
        endif
 
        ! calculate landuse factors
        call landuse_calc_fact
 
        if( io_l ) write(io_fid_log,*) '++++++ END   CONVERT LANDUSE DATA ++++++'
 
        ! output landuse file
        call landuse_write
     endif
 
     return