mod_land_vars.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
#include "scalelib.h"
module mod_land_vars
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_io
  use scale_prof
  use scale_debug
  use scale_land_grid_cartesc_index
  use scale_tracer
  use scale_cpl_sfc_index
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: land_vars_setup
  public :: land_vars_restart_read
  public :: land_vars_restart_write
  public :: land_vars_history
  public :: land_vars_total

  public :: land_vars_restart_create
  public :: land_vars_restart_open
  public :: land_vars_restart_def_var
  public :: land_vars_restart_enddef
  public :: land_vars_restart_close

  public :: convert_ws2vwc

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  logical,               public :: land_restart_output                 = .false.         

  character(len=H_LONG),  public :: land_restart_in_basename           = ''
  logical,                public :: land_restart_in_aggregate
  logical,                public :: land_restart_in_postfix_timelabel  = .false.         
  character(len=H_LONG),  public :: land_restart_out_basename          = ''
  logical,                public :: land_restart_out_aggregate
  logical,                public :: land_restart_out_postfix_timelabel = .true.          
  character(len=H_MID),   public :: land_restart_out_title             = 'LAND restart'
  character(len=H_SHORT), public :: land_restart_out_dtype             = 'DEFAULT'

  ! prognostic variables
  real(RP), public, allocatable :: land_temp      (:,:,:)
  real(RP), public, allocatable :: land_water     (:,:,:)
  real(RP), public, allocatable :: land_sfc_temp  (:,:)
  real(RP), public, allocatable :: land_sfc_albedo(:,:,:,:)

  ! for snow model
  real(RP), public, allocatable :: snow_sfc_temp  (:,:)
  real(RP), public, allocatable :: snow_swe       (:,:)
  real(RP), public, allocatable :: snow_depth     (:,:)
  real(RP), public, allocatable :: snow_dzero     (:,:)
  real(RP), public, allocatable :: snow_nosnowsec (:,:)

  ! tendency variables
  real(RP), public, allocatable :: land_temp_t    (:,:,:)
  real(RP), public, allocatable :: land_water_t   (:,:,:)

  ! surface variables for restart
  real(RP), public, allocatable :: land_sflx_mw   (:,:)
  real(RP), public, allocatable :: land_sflx_mu   (:,:)
  real(RP), public, allocatable :: land_sflx_mv   (:,:)
  real(RP), public, allocatable :: land_sflx_sh   (:,:)
  real(RP), public, allocatable :: land_sflx_lh   (:,:)
  real(RP), public, allocatable :: land_sflx_qtrc (:,:,:)
  real(RP), public, allocatable :: land_sflx_gh   (:,:)
  real(RP), public, allocatable :: land_sflx_water(:,:)
  real(RP), public, allocatable :: land_sflx_ice  (:,:)

  ! diagnostic variables
  real(RP), public, allocatable :: land_u10       (:,:)
  real(RP), public, allocatable :: land_v10       (:,:)
  real(RP), public, allocatable :: land_t2        (:,:)
  real(RP), public, allocatable :: land_q2        (:,:)

  ! recieved atmospheric variables
  real(RP), public, allocatable :: atmos_temp       (:,:)
  real(RP), public, allocatable :: atmos_pres       (:,:)
  real(RP), public, allocatable :: atmos_w          (:,:)
  real(RP), public, allocatable :: atmos_u          (:,:)
  real(RP), public, allocatable :: atmos_v          (:,:)
  real(RP), public, allocatable :: atmos_dens       (:,:)
  real(RP), public, allocatable :: atmos_qv         (:,:)
  real(RP), public, allocatable :: atmos_pbl        (:,:)
  real(RP), public, allocatable :: atmos_sfc_dens   (:,:)
  real(RP), public, allocatable :: atmos_sfc_pres   (:,:)
  real(RP), public, allocatable :: atmos_sflx_rad_dn(:,:,:,:)
  real(RP), public, allocatable :: atmos_cossza     (:,:)
  real(RP), public, allocatable :: atmos_sflx_rain  (:,:)
  real(RP), public, allocatable :: atmos_sflx_snow  (:,:)

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

  character(len=H_LONG), public :: land_property_in_filename  = ''

  integer,  public, parameter   :: land_property_nmax = 11
  integer,  public, parameter   :: i_waterlimit       =  1 ! maximum  soil moisture           [m3/m3]
  integer,  public, parameter   :: i_watercritical    =  2 ! critical soil moisture           [m3/m3]
  integer,  public, parameter   :: i_stomataresist    =  3 ! stomata resistance               [1/s]
  integer,  public, parameter   :: i_thermalcond      =  4 ! thermal conductivity for soil    [W/K/m]
  integer,  public, parameter   :: i_heatcapacity     =  5 ! heat capacity for soil           [J/K/m3]
  integer,  public, parameter   :: i_waterdiff        =  6 ! moisture diffusivity in the soil [m2/s]
  integer,  public, parameter   :: i_alblw            =  7 ! surface albedo for long  wave    [1]
  integer,  public, parameter   :: i_albsw            =  8 ! surface albedo for short wave    [1]
  integer,  public, parameter   :: i_z0m              =  9 ! roughness length for momemtum    [m]
  integer,  public, parameter   :: i_z0h              = 10 ! roughness length for heat        [m]
  integer,  public, parameter   :: i_z0e              = 11 ! roughness length for vapor       [m]

  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  private :: land_param_read

  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  logical,                private :: land_vars_checkrange      = .false.

  integer,                private, parameter :: vmax              = 19
  integer,                private, parameter :: i_temp            =  1
  integer,                private, parameter :: i_water           =  2
  integer,                private, parameter :: i_waterds         =  3
  integer,                private, parameter :: i_sfc_temp        =  4
  integer,                private, parameter :: i_sfc_alb_ir_dir  =  5
  integer,                private, parameter :: i_sfc_alb_ir_dif  =  6
  integer,                private, parameter :: i_sfc_alb_nir_dir =  7
  integer,                private, parameter :: i_sfc_alb_nir_dif =  8
  integer,                private, parameter :: i_sfc_alb_vis_dir =  9
  integer,                private, parameter :: i_sfc_alb_vis_dif = 10
  integer,                private, parameter :: i_sflx_mw         = 11
  integer,                private, parameter :: i_sflx_mu         = 12
  integer,                private, parameter :: i_sflx_mv         = 13
  integer,                private, parameter :: i_sflx_sh         = 14
  integer,                private, parameter :: i_sflx_lh         = 15
  integer,                private, parameter :: i_sflx_evap       = 16
  integer,                private, parameter :: i_sflx_gh         = 17
  integer,                private, parameter :: i_sflx_water      = 18
  integer,                private, parameter :: i_sflx_ice        = 19

  character(len=H_SHORT), private            :: var_name(vmax)
  character(len=H_MID),   private            :: var_desc(vmax)
  character(len=H_MID),   private            :: var_stdn(vmax)
  character(len=H_SHORT), private            :: var_unit(vmax)
  integer,                private            :: var_id(vmax)
  integer,                private            :: restart_fid = -1  ! file ID

  data var_name / 'LAND_TEMP',            &
                  'LAND_WATER',           &
                  'LAND_DSAT',            &
                  'LAND_SFC_TEMP',        &
                  'LAND_SFC_ALB_IR_dir',  &
                  'LAND_SFC_ALB_IR_dif',  &
                  'LAND_SFC_ALB_NIR_dir', &
                  'LAND_SFC_ALB_NIR_dif', &
                  'LAND_SFC_ALB_VIS_dir', &
                  'LAND_SFC_ALB_VIS_dif', &
                  'LAND_SFLX_MW',         &
                  'LAND_SFLX_MU',         &
                  'LAND_SFLX_MV',         &
                  'LAND_SFLX_SH',         &
                  'LAND_SFLX_LH',         &
                  'LAND_SFLX_evap',       &
                  'LAND_SFLX_GH',         &
                  'LAND_SFLX_water',      &
                  'LAND_SFLX_ice'         /
  data var_desc / 'temperature at each soil layer',            &
                  'moisture at each soil layer',               &
                  'degree of saturation at each soil layer',   &
                  'land surface skin temperature',             &
                  'land surface albedo for IR (direct)',       &
                  'land surface albedo for IR (diffuse)',      &
                  'land surface albedo for NIR (direct)',      &
                  'land surface albedo for NIR (diffuse)',     &
                  'land surface albedo for VIS (direct)',      &
                  'land surface albedo for VIS (diffuse)',     &
                  'land surface w-momentum flux (upward)',     &
                  'land surface u-momentum flux (upward)',     &
                  'land surface v-momentum flux (upward)',     &
                  'land surface sensible heat flux (upward)',  &
                  'land surface latent heat flux (upward)',    &
                  'land surface water vapor flux (upward)',    &
                  'land subsurface heat flux (downward)',      &
                  'land surface liquid water flux (downward)', &
                  'land surface ice water flux (downward)'     /
  data var_stdn / 'soil_temperature', &
                  'volume_fraction_of_condensed_water_in_soil', &
                  'volume_fraction_of_condensed_water_in_soil_at_field_capacity', &
                  'surface_temperature_where_land', &
                  '', &
                  '', &
                  '', &
                  '', &
                  '', &
                  '', &
                  '', &
                  '', &
                  '', &
                  '', &
                  '', &
                  '', &
                  '', &
                  '', &
                  '' /
  data var_unit / 'K',       &
                  'm3/m3',   &
                  '1',       &
                  'K',       &
                  '1',       &
                  '1',       &
                  '1',       &
                  '1',       &
                  '1',       &
                  '1',       &
                  'kg/m2/s', &
                  'kg/m2/s', &
                  'kg/m2/s', &
                  'J/m2/s',  &
                  'J/m2/s',  &
                  'kg/m2/s', &
                  'J/m2/s',  &
                  'kg/m2/s', &
                  'kg/m2/s'  /

  real(RP), private, allocatable :: land_property_table(:,:)

  logical,  private :: land_restart_in_check_coordinates = .true.

  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine land_vars_setup
    use scale_prc, only: &
       prc_abort
    use scale_const, only: &
       undef => const_undef
    use scale_comm_cartesc, only: &
       comm_vars8, &
       comm_wait
    use scale_landuse, only: &
       landuse_index_pft, &
       landuse_pft_nmin, &
       landuse_pft_nmax
    implicit none

    namelist / param_land_vars / &
       land_restart_in_basename,           &
       land_restart_in_aggregate,          &
       land_restart_in_postfix_timelabel,  &
       land_restart_in_check_coordinates,  &
       land_restart_output,                &
       land_restart_out_basename,          &
       land_restart_out_aggregate,         &
       land_restart_out_postfix_timelabel, &
       land_restart_out_title,             &
       land_restart_out_dtype,             &
       land_vars_checkrange

    integer :: ierr
    integer :: i, j, iv, p
    !---------------------------------------------------------------------------

    log_newline
    log_info("LAND_vars_setup",*) 'Setup'

    allocate( land_temp(lkmax,lia,lja)               )
    allocate( land_water(lkmax,lia,lja)               )
    allocate( land_sfc_temp(lia,lja)                     )
    allocate( land_sfc_albedo(lia,lja,n_rad_dir,n_rad_rgn) )
    land_temp(:,:,:)   = undef
    land_water(:,:,:)   = undef
    land_sfc_temp(:,:)     = undef
    land_sfc_albedo(:,:,:,:) = undef

    allocate( snow_sfc_temp(lia,lja) )
    allocate( snow_swe(lia,lja) )
    allocate( snow_depth(lia,lja) )
    allocate( snow_dzero(lia,lja) )
    allocate( snow_nosnowsec(lia,lja) )
    snow_sfc_temp(:,:)     = undef
    snow_swe(:,:)     = undef
    snow_depth(:,:)     = undef
    snow_dzero(:,:)     = undef
    snow_nosnowsec(:,:)     = undef

    allocate( land_temp_t(lkmax,lia,lja) )
    allocate( land_water_t(lkmax,lia,lja) )
    land_temp_t(:,:,:)   = undef
    land_water_t(:,:,:)   = undef

    allocate( land_sflx_mw(lia,lja) )
    allocate( land_sflx_mu(lia,lja) )
    allocate( land_sflx_mv(lia,lja) )
    allocate( land_sflx_sh(lia,lja) )
    allocate( land_sflx_lh(lia,lja) )
    allocate( land_sflx_gh(lia,lja) )
    allocate( land_sflx_qtrc(lia,lja,qa) )
    allocate( land_sflx_water(lia,lja) )
    allocate( land_sflx_ice(lia,lja) )
    land_sflx_mw(:,:)     = undef
    land_sflx_mu(:,:)     = undef
    land_sflx_mv(:,:)     = undef
    land_sflx_sh(:,:)     = undef
    land_sflx_lh(:,:)     = undef
    land_sflx_qtrc(:,:,:)   = undef
    land_sflx_gh(:,:)     = undef
    land_sflx_water(:,:)     = undef
    land_sflx_ice(:,:)     = undef

    allocate( land_u10(lia,lja) )
    allocate( land_v10(lia,lja) )
    allocate( land_t2(lia,lja) )
    allocate( land_q2(lia,lja) )
    land_u10(:,:)     = undef
    land_v10(:,:)     = undef
    land_t2(:,:)     = undef
    land_q2(:,:)     = undef

    allocate( atmos_temp(lia,lja)                     )
    allocate( atmos_pres(lia,lja)                     )
    allocate( atmos_w(lia,lja)                     )
    allocate( atmos_u(lia,lja)                     )
    allocate( atmos_v(lia,lja)                     )
    allocate( atmos_dens(lia,lja)                     )
    allocate( atmos_qv(lia,lja)                     )
    allocate( atmos_pbl(lia,lja)                     )
    allocate( atmos_sfc_dens(lia,lja)                     )
    allocate( atmos_sfc_pres(lia,lja)                     )
    allocate( atmos_sflx_rad_dn(lia,lja,n_rad_dir,n_rad_rgn) )
    allocate( atmos_cossza(lia,lja)                     )
    allocate( atmos_sflx_rain(lia,lja)                     )
    allocate( atmos_sflx_snow(lia,lja)                     )
    atmos_temp(:,:)     = undef
    atmos_pres(:,:)     = undef
    atmos_w(:,:)     = undef
    atmos_u(:,:)     = undef
    atmos_v(:,:)     = undef
    atmos_dens(:,:)     = undef
    atmos_qv(:,:)     = undef
    atmos_pbl(:,:)     = undef
    atmos_sfc_dens(:,:)     = undef
    atmos_sfc_pres(:,:)     = undef
    atmos_sflx_rad_dn(:,:,:,:) = undef
    atmos_cossza(:,:)     = undef
    atmos_sflx_rain(:,:)     = undef
    atmos_sflx_snow(:,:)     = undef

    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_land_vars,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("LAND_vars_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("LAND_vars_setup",*) 'Not appropriate names in namelist PARAM_LAND_VARS. Check!'
       call prc_abort
    endif
    log_nml(param_land_vars)

    log_newline
    log_info("LAND_vars_setup",*) 'List of prognostic variables (LAND) '
    log_info_cont('(1x,A,A24,A,A48,A,A12,A)') &
               '      |', 'VARNAME                 ','|', &
               'DESCRIPTION                                     ', '[', 'UNIT        ', ']'
    do iv = 1, vmax
       log_info_cont('(1x,A,I3,A,A24,A,A48,A,A12,A)') &
                  'NO.',iv,'|',var_name(iv),'|',var_desc(iv),'[',var_unit(iv),']'
    enddo

    log_newline
    if ( land_restart_in_basename /= '' ) then
       log_info("LAND_vars_setup",*) 'Restart input?  : YES, file = ', trim(land_restart_in_basename)
       log_info("LAND_vars_setup",*) 'Add timelabel?  : ', land_restart_in_postfix_timelabel
    else
       log_info("LAND_vars_setup",*) 'Restart input?  : NO'
    endif
    if (       land_restart_output             &
         .AND. land_restart_out_basename /= '' ) then
       log_info("LAND_vars_setup",*) 'Restart output? : YES, file = ', trim(land_restart_out_basename)
       log_info("LAND_vars_setup",*) 'Add timelabel?  : ', land_restart_out_postfix_timelabel
    else
       log_info("LAND_vars_setup",*) 'Restart output? : NO'
       land_restart_output = .false.
    endif

    ! Read land property table
    allocate( land_property_table(landuse_pft_nmin:landuse_pft_nmax,land_property_nmax) )
    land_property_table(:,:) = undef

    call land_param_read

    ! Apply land property to 2D map
    allocate( land_property(lia,lja,land_property_nmax) )

    ! tentative, mosaic is off
    do p = 1, land_property_nmax
    do j = ljs, lje
    do i = lis, lie
       land_property(i,j,p) = land_property_table( landuse_index_pft(i,j,1), p )
    enddo
    enddo
    enddo

    do p = 1, land_property_nmax
       call comm_vars8( land_property(:,:,p), p )
    enddo
    do p = 1, land_property_nmax
       call comm_wait ( land_property(:,:,p), p )
    enddo

    return
  end subroutine land_vars_setup

  !-----------------------------------------------------------------------------
  subroutine land_vars_restart_open
    use scale_time, only: &
       time_gettimelabel
    use scale_file_cartesc, only: &
       file_cartesc_open, &
       file_cartesc_check_coordinates
    use mod_land_admin, only: &
       land_do
    implicit none

    character(len=19)     :: timelabel
    character(len=H_LONG) :: basename
    !---------------------------------------------------------------------------

    log_newline
    log_info("LAND_vars_restart_open",*) 'Open restart file (LAND) '

    if ( land_do .and. land_restart_in_basename /= '' ) then

       if ( land_restart_in_postfix_timelabel ) then
          call time_gettimelabel( timelabel )
          basename = trim(land_restart_in_basename)//'_'//trim(timelabel)
       else
          basename = trim(land_restart_in_basename)
       endif

       log_info("LAND_vars_restart_open",*) 'basename: ', trim(basename)

       call file_cartesc_open( basename, restart_fid, aggregate=land_restart_in_aggregate )

       if ( land_restart_in_check_coordinates ) then
          call file_cartesc_check_coordinates( restart_fid, land=.true. )
       end if

    else
       log_info("LAND_vars_restart_open",*) 'restart file for land is not specified.'
    endif

    return
  end subroutine land_vars_restart_open

  !-----------------------------------------------------------------------------
  subroutine land_vars_restart_read
    use scale_prc, only: &
       prc_abort
    use scale_file, only: &
       file_get_aggregate
    use scale_file_cartesc, only: &
       file_cartesc_read, &
       file_cartesc_flush
    implicit none
    !---------------------------------------------------------------------------

    if ( restart_fid /= -1 ) then
       log_newline
       log_info("LAND_vars_restart_read",*) 'Read from restart file (LAND) '

       call file_cartesc_read( restart_fid, var_name(i_temp),            'LXY', & ! [IN]
                               land_temp(:,:,:)                           ) ! [OUT]
       call file_cartesc_read( restart_fid, var_name(i_water),           'LXY', & ! [IN]
                               land_water(:,:,:)                           ) ! [OUT]
       call file_cartesc_read( restart_fid, var_name(i_sfc_temp),        'XY',  & ! [IN]
                               land_sfc_temp(:,:)                             ) ! [OUT]
       call file_cartesc_read( restart_fid, var_name(i_sfc_alb_ir_dir),  'XY',  & ! [IN]
                               land_sfc_albedo(:,:,i_r_direct ,i_r_ir )         ) ! [OUT]
       call file_cartesc_read( restart_fid, var_name(i_sfc_alb_ir_dif),  'XY',  & ! [IN]
                               land_sfc_albedo(:,:,i_r_diffuse,i_r_ir )         ) ! [OUT]
       call file_cartesc_read( restart_fid, var_name(i_sfc_alb_nir_dir), 'XY',  & ! [IN]
                               land_sfc_albedo(:,:,i_r_direct ,i_r_nir)         ) ! [OUT]
       call file_cartesc_read( restart_fid, var_name(i_sfc_alb_nir_dif), 'XY',  & ! [IN]
                               land_sfc_albedo(:,:,i_r_diffuse,i_r_nir)         ) ! [OUT]
       call file_cartesc_read( restart_fid, var_name(i_sfc_alb_vis_dir), 'XY',  & ! [IN]
                               land_sfc_albedo(:,:,i_r_direct ,i_r_vis)         ) ! [OUT]
       call file_cartesc_read( restart_fid, var_name(i_sfc_alb_vis_dif), 'XY',  & ! [IN]
                               land_sfc_albedo(:,:,i_r_diffuse,i_r_vis)         ) ! [OUT]

       !call FILE_CARTESC_read( restart_fid, 'SNOW_SFC_TEMP',      'XY',  & ! [OUT]
       !                        SNOW_SFC_TEMP(:,:)                        ) ! [IN]
       !call FILE_CARTESC_read( restart_fid, 'SNOW_SWE',           'XY',  & ! [OUT]
       !                        SNOW_SWE(:,:)                             ) ! [IN]
       !call FILE_CARTESC_read( restart_fid, 'SNOW_Depth',         'XY',  & ! [OUT]
       !                        SNOW_Depth(:,:)                           ) ! [IN]
       !call FILE_CARTESC_read( restart_fid, 'SNOW_Dzero',         'XY',  & ! [OUT]
       !                        SNOW_Dzero(:,:)                           ) ! [IN]
       !call FILE_CARTESC_read( restart_fid, 'SNOW_nosnowsec',     'XY'   & ! [OUT]
       !                        SNOW_nosnowsec(:,:)                       ) ! [IN]

       !!!!! Tentative for snow model !!!!!
       snow_sfc_temp    = 273.15_rp
       snow_swe         = 0.0_rp
       snow_depth       = 0.0_rp
       snow_dzero       = 0.0_rp
       snow_nosnowsec   = 0.0_rp

       if( file_get_aggregate(restart_fid) ) call file_cartesc_flush( restart_fid ) ! commit all pending read requests

       call land_vars_total
    else
       log_error("LAND_vars_restart_read",*) 'invalid restart file ID for land.'
       call prc_abort
    endif

    return
  end subroutine land_vars_restart_read

  !-----------------------------------------------------------------------------
  subroutine land_vars_history
    use scale_file_history, only: &
       file_history_in
    use scale_atmos_hydrometeor, only: &
       i_qv
    implicit none

    real(RP) :: LAND_WATERDS(lkmax,lia,lja)
    integer :: k, i, j
    !---------------------------------------------------------------------------

    call prof_rapstart('LND_History', 1)


    if ( land_vars_checkrange ) then
       call valcheck( land_temp(lks:lke,lis:lie,ljs:lje),             0.0_rp, 1000.0_rp, &
                      var_name(i_temp),                                     __file__, __line__ )
       call valcheck( land_water(lks:lke,lis:lie,ljs:lje),             0.0_rp, 1000.0_rp, &
                      var_name(i_water),                                    __file__, __line__ )
       call valcheck( land_sfc_temp(lis:lie,ljs:lje),                     0.0_rp, 1000.0_rp, &
                      var_name(i_sfc_temp),                                 __file__, __line__ )
       call valcheck( land_sfc_albedo(lis:lie,ljs:lje,i_r_direct ,i_r_ir ), 0.0_rp,    2.0_rp, &
                      var_name(i_sfc_alb_ir_dir ),                          __file__, __line__ )
       call valcheck( land_sfc_albedo(lis:lie,ljs:lje,i_r_diffuse,i_r_ir ), 0.0_rp,    2.0_rp, &
                      var_name(i_sfc_alb_ir_dif ),                          __file__, __line__ )
       call valcheck( land_sfc_albedo(lis:lie,ljs:lje,i_r_direct ,i_r_nir), 0.0_rp,    2.0_rp, &
                      var_name(i_sfc_alb_nir_dir),                          __file__, __line__ )
       call valcheck( land_sfc_albedo(lis:lie,ljs:lje,i_r_diffuse,i_r_nir), 0.0_rp,    2.0_rp, &
                      var_name(i_sfc_alb_nir_dif),                          __file__, __line__ )
       call valcheck( land_sfc_albedo(lis:lie,ljs:lje,i_r_direct ,i_r_vis), 0.0_rp,    2.0_rp, &
                      var_name(i_sfc_alb_vis_dir),                          __file__, __line__ )
       call valcheck( land_sfc_albedo(lis:lie,ljs:lje,i_r_diffuse,i_r_vis), 0.0_rp,    2.0_rp, &
                      var_name(i_sfc_alb_vis_dif),                          __file__, __line__ )

       !call VALCHECK( SNOW_SFC_TEMP (IS:IE,JS:JE),      0.0_RP, 1000.0_RP, 'SNOW_SFC_TEMP',        &
       !              __FILE__, __LINE__ )
       !call VALCHECK( SNOW_SWE      (IS:IE,JS:JE),      0.0_RP, 1000.0_RP, 'SNOW_SWE',             &
       !              __FILE__, __LINE__ )
       !call VALCHECK( SNOW_Depth    (IS:IE,JS:JE),      0.0_RP, 1000.0_RP, 'SNOW_Depth',           &
       !              __FILE__, __LINE__ )
       !call VALCHECK( SNOW_Dzero    (IS:IE,JS:JE),      0.0_RP, 1000.0_RP, 'SNOW_Dzero',           &
       !              __FILE__, __LINE__ )
    endif

    call file_history_in( land_temp(:,:,:), var_name(i_temp),  var_desc(i_temp),  var_unit(i_temp),  dim_type='LXY', standard_name=var_stdn(i_temp) )
    call file_history_in( land_water(:,:,:), var_name(i_water), var_desc(i_water), var_unit(i_water), dim_type='LXY', standard_name=var_stdn(i_water) )
    do j = ljs, lje
    do i = lis, lie
    do k = 1, lkmax
       land_waterds(k,i,j) = land_water(k,i,j) / land_property(i,j,i_waterlimit)
    end do
    end do
    end do
    call file_history_in( land_waterds(:,:,:), var_name(i_waterds), var_desc(i_waterds), var_unit(i_waterds), dim_type='LXY', fill_halo=.true., standard_name=var_stdn(i_waterds) )


    call file_history_in( land_sfc_temp(:,:),                     var_name(i_sfc_temp),                                     &
                          var_desc(i_sfc_temp),        var_unit(i_sfc_temp),        standard_name=var_stdn(i_sfc_temp)        )
    call file_history_in( land_sfc_albedo(:,:,i_r_direct ,i_r_ir ), var_name(i_sfc_alb_ir_dir),                               &
                          var_desc(i_sfc_alb_ir_dir),  var_unit(i_sfc_alb_ir_dir),  standard_name=var_stdn(i_sfc_alb_ir_dir)  )
    call file_history_in( land_sfc_albedo(:,:,i_r_diffuse,i_r_ir ), var_name(i_sfc_alb_ir_dif),                               &
                          var_desc(i_sfc_alb_ir_dif),  var_unit(i_sfc_alb_ir_dif),  standard_name=var_stdn(i_sfc_alb_ir_dif)  )
    call file_history_in( land_sfc_albedo(:,:,i_r_direct ,i_r_nir), var_name(i_sfc_alb_nir_dir),                              &
                          var_desc(i_sfc_alb_nir_dir), var_unit(i_sfc_alb_nir_dir), standard_name=var_stdn(i_sfc_alb_nir_dir) )
    call file_history_in( land_sfc_albedo(:,:,i_r_diffuse,i_r_nir), var_name(i_sfc_alb_nir_dif),                              &
                          var_desc(i_sfc_alb_nir_dif), var_unit(i_sfc_alb_nir_dif), standard_name=var_stdn(i_sfc_alb_nir_dif) )
    call file_history_in( land_sfc_albedo(:,:,i_r_direct ,i_r_vis), var_name(i_sfc_alb_vis_dir),                              &
                          var_desc(i_sfc_alb_vis_dir), var_unit(i_sfc_alb_vis_dir), standard_name=var_stdn(i_sfc_alb_vis_dir) )
    call file_history_in( land_sfc_albedo(:,:,i_r_diffuse,i_r_vis), var_name(i_sfc_alb_vis_dif),                              &
                          var_desc(i_sfc_alb_vis_dif), var_unit(i_sfc_alb_vis_dif), standard_name=var_stdn(i_sfc_alb_vis_dif) )

    call file_history_in( land_sflx_mw(:,:),                     var_name(i_sflx_mw),                                      &
                          var_desc(i_sflx_mw),         var_unit(i_sflx_mw),         standard_name=var_stdn(i_sflx_mw)         )
    call file_history_in( land_sflx_mu(:,:),                     var_name(i_sflx_mu),                                      &
                          var_desc(i_sflx_mu),         var_unit(i_sflx_mu),         standard_name=var_stdn(i_sflx_mu)         )
    call file_history_in( land_sflx_mv(:,:),                     var_name(i_sflx_mv),                                      &
                          var_desc(i_sflx_mv),         var_unit(i_sflx_mv),         standard_name=var_stdn(i_sflx_mv)         )
    call file_history_in( land_sflx_sh(:,:),                     var_name(i_sflx_sh),                                      &
                          var_desc(i_sflx_sh),         var_unit(i_sflx_sh),         standard_name=var_stdn(i_sflx_sh)         )
    call file_history_in( land_sflx_lh(:,:),                     var_name(i_sflx_lh),                                      &
                          var_desc(i_sflx_lh),         var_unit(i_sflx_lh),         standard_name=var_stdn(i_sflx_lh)         )
    if ( i_qv > 0 ) then
    call file_history_in( land_sflx_qtrc(:,:,i_qv),                var_name(i_sflx_evap),                                    &
                          var_desc(i_sflx_evap),       var_unit(i_sflx_evap),       standard_name=var_stdn(i_sflx_evap)       )
    endif
    call file_history_in( land_sflx_gh(:,:),                     var_name(i_sflx_gh),                                      &
                          var_desc(i_sflx_gh),         var_unit(i_sflx_gh),         standard_name=var_stdn(i_sflx_gh)         )
    call file_history_in( land_sflx_water(:,:),                     var_name(i_sflx_water),                                   &
                          var_desc(i_sflx_water),      var_unit(i_sflx_water),      standard_name=var_stdn(i_sflx_water)      )
    call file_history_in( land_sflx_ice(:,:),                     var_name(i_sflx_ice),                                     &
                          var_desc(i_sflx_ice),        var_unit(i_sflx_ice),        standard_name=var_stdn(i_sflx_ice)        )

    ! snow model
    call file_history_in( snow_sfc_temp(:,:), 'SNOW_SFC_TEMP',  'Snow surface temperature',    'K'      )
    call file_history_in( snow_swe(:,:), 'SNOW_SWE',       'Snow water equivalent',       'kg/m2'  )
    call file_history_in( snow_depth(:,:), 'SNOW_Depth',     'Snow depth',                  'm'      )
    call file_history_in( snow_dzero(:,:), 'SNOW_Dzero',     'Snow depth at melting point', 'm'      )
    call file_history_in( snow_nosnowsec(:,:), 'SNOW_nosnowsec', 'Time duration without snow',  's'      )

    call prof_rapend  ('LND_History', 1)

    return
  end subroutine land_vars_history

  !-----------------------------------------------------------------------------
  subroutine land_vars_total
    use scale_statistics, only: &
       statistics_checktotal, &
       statistics_total
    use scale_land_grid_cartesc_real, only: &
       land_grid_cartesc_real_area,    &
       land_grid_cartesc_real_totarea, &
       land_grid_cartesc_real_vol,     &
       land_grid_cartesc_real_totvol
    implicit none

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

    if ( statistics_checktotal ) then

       ! 3D
       call statistics_total( lka, lks, lke, lia, lis, lie, lja, ljs, lje, &
                              land_temp(:,:,:), var_name(i_temp),  & ! (in)
                              land_grid_cartesc_real_vol(:,:,:),    & ! (in)
                              land_grid_cartesc_real_totvol         ) ! (in)
       call statistics_total( lka, lks, lke, lia, lis, lie, lja, ljs, lje, &
                              land_water(:,:,:), var_name(i_water), & ! (in)
                              land_grid_cartesc_real_vol(:,:,:),    & ! (in)
                              land_grid_cartesc_real_totvol         ) ! (in)

       ! 2D
       call statistics_total( lia, lis, lie, lja, ljs, lje,                                          & ! [IN]
                              land_sfc_temp(:,:),                     var_name(i_sfc_temp),        & ! [IN]
                              land_grid_cartesc_real_area(:,:),                                      & ! [IN]
                              land_grid_cartesc_real_totarea                                         ) ! [IN]
       call statistics_total( lia, lis, lie, lja, ljs, lje,                                          & ! [IN]
                              land_sfc_albedo(:,:,i_r_direct ,i_r_ir ), var_name(i_sfc_alb_ir_dir),  & ! [IN]
                              land_grid_cartesc_real_area(:,:),                                      & ! [IN]
                              land_grid_cartesc_real_totarea                                         ) ! [IN]
       call statistics_total( lia, lis, lie, lja, ljs, lje,                                          & ! [IN]
                              land_sfc_albedo(:,:,i_r_diffuse,i_r_ir ), var_name(i_sfc_alb_ir_dif),  & ! [IN]
                              land_grid_cartesc_real_area(:,:),                                      & ! [IN]
                              land_grid_cartesc_real_totarea                                         ) ! [IN]
       call statistics_total( lia, lis, lie, lja, ljs, lje,                                          & ! [IN]
                              land_sfc_albedo(:,:,i_r_direct ,i_r_nir), var_name(i_sfc_alb_nir_dir), & ! [IN]
                              land_grid_cartesc_real_area(:,:),                                      & ! [IN]
                              land_grid_cartesc_real_totarea                                         ) ! [IN]
       call statistics_total( lia, lis, lie, lja, ljs, lje,                                          & ! [IN]
                              land_sfc_albedo(:,:,i_r_diffuse,i_r_nir), var_name(i_sfc_alb_nir_dif), & ! [IN]
                              land_grid_cartesc_real_area(:,:),                                      & ! [IN]
                              land_grid_cartesc_real_totarea                                         ) ! [IN]
       call statistics_total( lia, lis, lie, lja, ljs, lje,                                          & ! [IN]
                              land_sfc_albedo(:,:,i_r_direct ,i_r_vis), var_name(i_sfc_alb_vis_dir), & ! [IN]
                              land_grid_cartesc_real_area(:,:),                                      & ! [IN]
                              land_grid_cartesc_real_totarea                                         ) ! [IN]
       call statistics_total( lia, lis, lie, lja, ljs, lje,                                          & ! [IN]
                              land_sfc_albedo(:,:,i_r_diffuse,i_r_vis), var_name(i_sfc_alb_vis_dif), & ! [IN]
                              land_grid_cartesc_real_area(:,:),                                      & ! [IN]
                              land_grid_cartesc_real_totarea                                         ) ! [IN]

       call statistics_total( lia, lis, lie, lja, ljs, lje, &
                              snow_sfc_temp(:,:),     'SNOW_SFC_TEMP', & ! (in)
                              land_grid_cartesc_real_area(:,:),          & ! (in)
                              land_grid_cartesc_real_totarea             ) ! (in)
       call statistics_total( lia, lis, lie, lja, ljs, lje, &
                              snow_swe(:,:),     'SNOW_SWE',      &
                              land_grid_cartesc_real_area(:,:),          & ! (in)
                              land_grid_cartesc_real_totarea             ) ! (in)
       call statistics_total( lia, lis, lie, lja, ljs, lje, &
                              snow_depth(:,:),     'SNOW_Depth',    & ! (in)
                              land_grid_cartesc_real_area(:,:),          & ! (in)
                              land_grid_cartesc_real_totarea             ) ! (in)
       call statistics_total( lia, lis, lie, lja, ljs, lje, &
                              snow_dzero(:,:),     'SNOW_Dzero',    & ! (in)
                              land_grid_cartesc_real_area(:,:),          & ! (in)
                              land_grid_cartesc_real_totarea             ) ! (in)
    endif

    return
  end subroutine land_vars_total

  !-----------------------------------------------------------------------------
  subroutine land_param_read
    use scale_prc, only: &
       prc_abort
    use scale_landuse, only: &
       landuse_pft_nmin, &
       landuse_pft_nmax
    implicit none

    integer              :: index
    character(len=H_MID) :: description
    real(RP)             :: STRGMAX
    real(RP)             :: STRGCRT
    real(RP)             :: RSTOMA
    real(RP)             :: TCS
    real(RP)             :: HCS
    real(RP)             :: DFW
    real(RP)             :: ALBLW
    real(RP)             :: ALBSW
    real(RP)             :: Z0M
    real(RP)             :: Z0H
    real(RP)             :: Z0E

    namelist / param_land_property / &
       land_property_in_filename

    namelist / param_land_data / &
       index,       &
       description, &
       strgmax,     &
       strgcrt,     &
       rstoma,      &
       tcs,         &
       hcs,         &
       dfw,         &
       alblw,       &
       albsw,       &
       z0m,         &
       z0h,         &
       z0e

    integer :: n
    integer :: ierr

    integer :: IO_FID_LAND_PROPERTY
    !---------------------------------------------------------------------------

    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_land_property,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("LAND_param_read",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("LAND_param_read",*) 'Not appropriate names in namelist PARAM_LAND_PROPERTY. Check!'
       call prc_abort
    endif
    log_nml(param_land_property)

    if( land_property_in_filename /= '' ) then
      !--- Open land parameter file
      io_fid_land_property = io_get_available_fid()
      open( io_fid_land_property,                     &
            file   = trim(land_property_in_filename), &
            form   = 'formatted',                     &
            status = 'old',                           &
            iostat = ierr                             )

      if ( ierr /= 0 ) then
         log_error("LAND_param_read",*) 'Failed to open land parameter file! :', &
                    trim(land_property_in_filename)
         call prc_abort
      else
        log_newline
        log_info("LAND_param_read",*) 'Properties for each plant functional type (PFT)'
        log_info_cont('(12(1x,A))') '                         PFT DESCRIPTION', &
                                    'Max Stg', &
                                    'CRT Stg', &
                                    'Stm.Res', &
                                    'T condu', &
                                    'H capac', &
                                    'DFC Wat', &
                                    'LW ALB',  &
                                    'SW ALB',  &
                                    ' Z0(m)',  &
                                    ' Z0(h)',  &
                                    ' Z0(e)'

        !--- read namelist
        rewind(io_fid_land_property)

        do n = landuse_pft_nmin, landuse_pft_nmax
           ! default value
           albsw   =  0.2_rp
           strgmax =  0.2_rp
           strgcrt =  0.1_rp
           rstoma  = 50.0_rp
           tcs     =  1.0_rp
           hcs     =  2.e+6_rp
           dfw     =  1.e-6_rp
           alblw   =  0.04_rp
           albsw   =  0.22_rp
           z0m     =  0.1_rp
           z0h     = -1.0_rp
           z0e     = -1.0_rp

           read(io_fid_land_property,nml=param_land_data,iostat=ierr)
           if ( ierr < 0 ) then !--- no more data
              exit
           elseif( ierr > 0 ) then !--- fatal error
              log_error("LAND_param_read",*) 'Not appropriate names in namelist PARAM_LAND_DATA. Check!'
              call prc_abort
           endif

           if( z0h < 0.0_rp ) then
             z0h = z0m / 7.4_rp ! defined by Garratt and Francey (1978)
           endif
           if( z0e < 0.0_rp ) then
             z0e = z0m / 7.4_rp ! defined by Garratt and Francey (1978)
           endif

           land_property_table(index,i_waterlimit   ) = strgmax
           land_property_table(index,i_watercritical) = strgcrt
           land_property_table(index,i_stomataresist) = rstoma
           land_property_table(index,i_thermalcond  ) = tcs
           land_property_table(index,i_heatcapacity ) = hcs
           land_property_table(index,i_waterdiff    ) = dfw
           land_property_table(index,i_alblw        ) = alblw
           land_property_table(index,i_albsw        ) = albsw
           land_property_table(index,i_z0m          ) = z0m
           land_property_table(index,i_z0h          ) = z0h
           land_property_table(index,i_z0e          ) = z0e

           log_info_cont('(1x,A4,I4.3,1x,A32,4(1x,F7.3),2(1x,ES7.1),5(1x,F6.3))') &
                                         'IDX=', index, &
                                         trim(description), &
                                         strgmax, &
                                         strgcrt, &
                                         rstoma,  &
                                         tcs,     &
                                         hcs,     &
                                         dfw,     &
                                         alblw,   &
                                         albsw,   &
                                         z0m,     &
                                         z0h,     &
                                         z0e
        enddo

      end if

      close( io_fid_land_property )

    endif

    return
  end subroutine land_param_read

  !-----------------------------------------------------------------------------
  function convert_ws2vwc( WS, critical ) result( VWC )
    implicit none

    real(RP), intent(in) :: WS(lia,lja) ! water saturation [fraction]
    logical,  intent(in) :: critical  ! is I_WaterCritical used?

    real(RP) :: VWC(lia,lja) ! volumetric water content [m3/m3]

    ! work
    integer :: i, j, num
    !---------------------------------------------------------------------------

    if( critical ) then
      num = i_watercritical
    else
      num = i_waterlimit
    end if

    do j = ljs, lje
    do i = lis, lie
      vwc(i,j) = max( min( ws(i,j)*land_property(i,j,num), land_property(i,j,num) ), 0.0_rp )
    end do
    end do

    return
  end function convert_ws2vwc

  !-----------------------------------------------------------------------------
  subroutine land_vars_restart_create
    use scale_time, only: &
       time_gettimelabel
    use scale_file_cartesc, only: &
       file_cartesc_create
    use mod_land_admin, only: &
       land_do
    implicit none

    character(len=19)     :: timelabel
    character(len=H_LONG) :: basename
    !---------------------------------------------------------------------------

    if ( land_do .and. land_restart_out_basename /= '' ) then

       log_newline
       log_info("LAND_vars_restart_create",*) 'Create restart file (LAND) '

       if ( land_restart_out_postfix_timelabel ) then
          call time_gettimelabel( timelabel )
          basename = trim(land_restart_out_basename)//'_'//trim(timelabel)
       else
          basename = trim(land_restart_out_basename)
       endif

       log_info("LAND_vars_restart_create",*) 'basename: ', trim(basename)

       call file_cartesc_create( &
            basename, land_restart_out_title, land_restart_out_dtype, & ! [IN]
            restart_fid,                                              & ! [OUT]
            aggregate=land_restart_out_aggregate                      ) ! [IN]

    endif

    return
  end subroutine land_vars_restart_create

  !-----------------------------------------------------------------------------
  subroutine land_vars_restart_enddef
    use scale_file_cartesc, only: &
       file_cartesc_enddef
    implicit none

    if ( restart_fid /= -1 ) then
       call file_cartesc_enddef( restart_fid ) ! [IN]
    endif

    return
  end subroutine land_vars_restart_enddef

  !-----------------------------------------------------------------------------
  subroutine land_vars_restart_close
    use scale_file_cartesc, only: &
       file_cartesc_close
    implicit none
    !---------------------------------------------------------------------------

    if ( restart_fid /= -1 ) then
       log_newline
       log_info("LAND_vars_restart_close",*) 'Close restart file (LAND) '

       call file_cartesc_close( restart_fid ) ! [IN]

       restart_fid = -1
    endif

    return
  end subroutine land_vars_restart_close

  !-----------------------------------------------------------------------------
  subroutine land_vars_restart_def_var
    use scale_file_cartesc, only: &
       file_cartesc_def_var
    implicit none
    integer :: i
    !---------------------------------------------------------------------------

    if ( restart_fid /= -1 ) then

       do i = i_temp, i_water
          if ( i == i_waterds ) cycle
          call file_cartesc_def_var( restart_fid,     & ! [IN]
               var_name(i), var_desc(i), var_unit(i), & ! [IN]
               'LXY', land_restart_out_dtype,         & ! [IN]
               var_id(i),                             & ! [OUT]
               standard_name=var_stdn(i)              ) ! [IN]
       end do
       do i = i_sfc_temp, i_sfc_alb_vis_dif
          call file_cartesc_def_var( restart_fid,     & ! [IN]
               var_name(i), var_desc(i), var_unit(i), & ! [IN]
               'XY', land_restart_out_dtype,          & ! [IN]
               var_id(i),                             & ! [OUT]
               standard_name=var_stdn(i)              ) ! [IN]
       end do

       !call FILE_CARTESC_def_var( restart_fid, 'SNOW_SFC_TEMP',  'Snow surface temperature',    &
       !                           'K',                   'XY',   LAND_RESTART_OUT_DTYPE, ?????  )
       !call FILE_CARTESC_def_var( restart_fid, 'SNOW_SWE',       'Snow water equivalent',       &
       !                           'kg/m2',               'XY',   LAND_RESTART_OUT_DTYPE, ?????  )
       !call FILE_CARTESC_def_var( restart_fid, 'SNOW_Depth',     'Snow depth',                  &
       !                           'm',                   'XY',   LAND_RESTART_OUT_DTYPE, ?????  )
       !call FILE_CARTESC_def_var( restart_fid, 'SNOW_Dzero',     'Snow depth at melting point', &
       !                           'm',                   'XY',   LAND_RESTART_OUT_DTYPE, ?????  )
       !call FILE_CARTESC_def_var( restart_fid, 'SNOW_nosnowsec', 'Time duration without snow', &
       !                           's',                 'XY',   LAND_RESTART_OUT_DTYPE,   ?????  )

    endif

    return
  end subroutine land_vars_restart_def_var

  !-----------------------------------------------------------------------------
  subroutine land_vars_restart_write
    use scale_file_cartesc, only: &
       file_cartesc_write_var
    implicit none
    !---------------------------------------------------------------------------

    if ( restart_fid /= -1 ) then

       call land_vars_total

       call file_cartesc_write_var( restart_fid, var_id(i_temp),                         & ! [IN]
                                    land_temp(:,:,:),                                    & ! [IN]
                                    var_name(i_temp),            'LXY', fill_halo=.true. ) ! [IN]
       call file_cartesc_write_var( restart_fid, var_id(i_water),                        & ! [IN]
                                    land_water(:,:,:),                                   & ! [IN]
                                    var_name(i_water),           'LXY', fill_halo=.true. ) ! [IN]
       call file_cartesc_write_var( restart_fid, var_id(i_sfc_temp),                     & ! [IN]
                                    land_sfc_temp(:,:),                                  & ! [IN]
                                    var_name(i_sfc_temp),        'XY',  fill_halo=.true. ) ! [IN]
       call file_cartesc_write_var( restart_fid, var_id(i_sfc_alb_ir_dir),               & ! [IN]
                                    land_sfc_albedo(:,:,i_r_direct ,i_r_ir ),            & ! [IN]
                                    var_name(i_sfc_alb_ir_dir),  'XY',  fill_halo=.true. ) ! [IN]
       call file_cartesc_write_var( restart_fid, var_id(i_sfc_alb_ir_dif),               & ! [IN]
                                    land_sfc_albedo(:,:,i_r_diffuse,i_r_ir ),            & ! [IN]
                                    var_name(i_sfc_alb_ir_dif),  'XY',  fill_halo=.true. ) ! [IN]
       call file_cartesc_write_var( restart_fid, var_id(i_sfc_alb_nir_dir),              & ! [IN]
                                    land_sfc_albedo(:,:,i_r_direct ,i_r_nir),            & ! [IN]
                                    var_name(i_sfc_alb_nir_dir), 'XY',  fill_halo=.true. ) ! [IN]
       call file_cartesc_write_var( restart_fid, var_id(i_sfc_alb_nir_dif),              & ! [IN]
                                    land_sfc_albedo(:,:,i_r_diffuse,i_r_nir),            & ! [IN]
                                    var_name(i_sfc_alb_nir_dif), 'XY',  fill_halo=.true. ) ! [IN]
       call file_cartesc_write_var( restart_fid, var_id(i_sfc_alb_vis_dir),              & ! [IN]
                                    land_sfc_albedo(:,:,i_r_direct ,i_r_vis),            & ! [IN]
                                    var_name(i_sfc_alb_vis_dir), 'XY',  fill_halo=.true. ) ! [IN]
       call file_cartesc_write_var( restart_fid, var_id(i_sfc_alb_vis_dif),              & ! [IN]
                                    land_sfc_albedo(:,:,i_r_diffuse,i_r_vis),            & ! [IN]
                                    var_name(i_sfc_alb_vis_dif), 'XY',  fill_halo=.true. ) ! [IN]

       !call FILE_CARTESC_write_var( restart_fid, VMAX+1, SNOW_SFC_TEMP(:,:),                    &
       !                           'SNOW_SFC_TEMP',       'XY',   fill_halo=.true.               )
       !call FILE_CARTESC_write_var( restart_fid, VMAX+2, SNOW_SWE(:,:),                         &
       !                           'SNOW_SWE',            'XY',   fill_halo=.true.               )
       !call FILE_CARTESC_write_var( restart_fid, VMAX+3, SNOW_Depth(:,:),                       &
       !                           'SNOW_Depth',          'XY',   fill_halo=.true.               )
       !call FILE_CARTESC_write_var( restart_fid, VMAX+4, SNOW_Dzero(:,:),                       &
       !                           'SNOW_Dzero',          'XY',   fill_halo=.true.               )
       !call FILE_CARTESC_write_var( restart_fid, VMAX+5, SNOW_nosnowsec(:,:),                   &
       !                           'SNOW_nosnowsec',      'XY',   fill_halo=.true.               )

    endif

    return
  end subroutine land_vars_restart_write

end module mod_land_vars