Initialize Soil Variables

Subroutine Insoil (cupinp.f)

This subroutine reads in soil data.

READS IN:

  CLODSZ     Soil clod size index

  NDSOIL     Number of depths in soil
             (number of soil layers)

  REFDPT     Reference depth (m)
  
  IRESDU     Is there residue on the surface? (1=T, 0=F)
  
  Z0SOIL     Roughness length in meters
  	       0.003 is typical for clear ground
  	       0.01  is typical for ground with residue
  
  BD         Bulk density

  SANDFC     Sand fraction

  SILTFC     Silt fraction

  CLAYFC     Clay fraction

  QRTZFC     Quartz fraction in fractions of total weight

  DPE        Air entry potential (J/kg) (optional). If 0 is entered 
             it is calculated by G. Cambbell's equations.

  DBX        Exponent in the soil moisture release curve (optional).

  DAKS       Saturating water conductivity (optional).
  
  ZSOIL      Depth of nodes in soil (m) (ZOSOIL in output)
  
  IDOROC     Does rock layer exist? (1=T, 0=F)
  
  IROCLY     The NDSOIL layer where rock starts
  
  AKROCK     Rock thermal conductivity
  
  CPROCK     Rock heat capacity
  
  DELD       Convergence criteria of diff. equation    
             solution
             

             In code: "DELD is the tolerance on sum D(JZ)
             in newt-raphson soil eqs--1.e-5-1.e-7
             Should be at least 1.e-5 to get storage
             change to match fluxes."

  DRHSFC     Value to match soil water flow equations  
             with canopy vapor flow equations
             

             In code: "Tolerance on matching sfc
             rel hum from cpy vapor eqs with rel hum
             from soil eqs.  Determines how close ecpy and
             wcpys are.  0.001 is small and 0.01 is usually
             acceptable."

  RHSLOP     Value to match soil water flow equations
             with canopy vapor flow equations
             

  	     In code: "RHSLOP is the change of rhsfc
  	     below value from cpy vapor eqs to get slope
  	     of rhsfc vs soil evap for newt-raph iteration
  	     in soil eqs."
  	     

COMPUTES:

  DGEOM      Geometric diameter
  
  SGEOM	     Geometric standard deviation
  
  A-ESOIL    Coefficients for soil thermal conductivity calc.
  
  PE	     Air entry potential (j/kg)
  
  AKS        Saturated hydraulic conductivity (kg s m-3)
  
  FC1-5
  
  AKSMMH     Saturated hydraulic conductivity (mm/hr)



Return to cupidg2