Modules

@file @brief calcdecomp This module computes the IO decomposition when generated by PO internally.

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

! END OF INIT !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

TODO: this should be a shallow object, with a simple init() and get() function, and the contents should go back to being private. TODO: same comment as above.

@file pio_kinds.F90 @brief basic data types

@file pio_support.F90 @brief internal code for compiler workarounds, aborts and debug functions

! NOTE: If all atmospheres are identical in number of processes, ! number of threads, and grid layout, we should check that the ! user-provided number of processes and threads are consistent ! (or else, only accept one entry for these quantities when reading ! the namelist). ARE OTHER PROTECTIONS/CHECKS NEEDED???

PUBLIC MEMBER FUNCTIONS

! END OF INIT !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

$ subroutine init_soil_resistance() $ ! $ !DESCRIPTIONS $ ! initialize method for soil resis calculation $ ! $ ! !USES: $ use abortutils , only : endrun $ use fileutils , only : getavu, relavu $ use spmdMod , only : mpicom, masterproc $ use shr_mpi_mod , only : shr_mpi_bcast $ use clm_varctl , only : iulog, use_bedrock $ use controlMod , only : NLFilename $ use clm_nlUtilsMod , only : find_nlgroup_name $ $ ! !ARGUMENTS: $ !------------------------------------------------------------------------------ $ implicit none $ integer :: nu_nml ! unit for namelist file $ integer :: nml_error ! namelist i/o error flag $ character(), parameter :: subName = "('init_soil_resistance')" $ $ !----------------------------------------------------------------------- $ $! MUST agree with name in namelist and read statement $ namelist /soil_resis_inparm/ soil_resis_method $ $ ! Default values for namelist $ $ soil_resis_method = sl_14 $ $ ! Read soil_resis namelist $ if (masterproc) then $ nu_nml = getavu() $ open( nu_nml, file=trim(NLFilename), status='old', iostat=nml_error ) $ call find_nlgroup_name(nu_nml, 'soil_resis_inparm', status=nml_error) $ if (nml_error == 0) then $ read(nu_nml, nml=soil_resis_inparm,iostat=nml_error) $ if (nml_error /= 0) then $ call endrun(subname // ':: ERROR reading soil_resis namelist') $ end if $ end if $ close(nu_nml) $ call relavu( nu_nml ) $ $ endif $ $ call shr_mpi_bcast(soil_resis_method, mpicom) $ $ if (masterproc) then $ write(iulog,) ' ' $ write(iulog,) 'soil_resis settings:' $ write(iulog,) ' soil_resis_method = ',soil_resis_method $ endif $!scs $! soil_resis_method = leepielke_1992 $! soil_resis_method = sl_14 $!scs $ $ end subroutine init_soil_resistance