#include <stdio.h>
#include <stdlib.h>
#include <vicNl.h>

static char vcid[] = "$Id: write_layer.c,v 5.2 2003/07/31 17:42:13 cherkaue Exp $";

void write_layer(layer_data_struct *layer,
                 int                veg,
                 int                Nlayer,
#if SPATIAL_FROST
                 double            *frost_fract,
#endif
		 double            *depth)
/**********************************************************************
	write_soilvar		Keith Cherkauer		July 17, 1997

  This routine writes soil variables to stdout.  It creates a table 
  of soil moisture values which shows how much liquid water and ice 
  are contained in the thawed, frozen and unfrozen sublayers of each 
  soil layer.  It also gives the total soil moisture for each layer.

  xx-xx-01 Modified to handle spatial soil frost.                 KAC

**********************************************************************/
{
  extern option_struct options;

  int index;
  double layer_moist;
  double sum_moist;
#if SPATIAL_FROST
  int    frost_area;
  double avg_ice;
#endif

  printf("Layer Data for Vegetation Type #%i\n",veg);
  printf("Layer:\t");
  for(index=0;index<Nlayer;index++) printf("\t\t%i",index+1);
  printf("\nEvaporation:\t");
  for(index=0;index<Nlayer;index++) printf("\t%f",layer[index].evap);
  printf("\n      Kappa:\t");
  for(index=0;index<Nlayer;index++) printf("\t%f",layer[index].kappa);
  printf("\n         Cs:\t");
  for(index=0;index<Nlayer;index++) printf("\t%f",layer[index].Cs);
  printf("\n\nMoisture Table\n---------------------------------------------------------------------------\n Moist:\t");
  for(index=0;index<Nlayer;index++) printf("\t%f",layer[index].moist);
  printf("\n        Ice:\t");
#if SPATIAL_FROST
  for(index=0;index<Nlayer;index++) {
    avg_ice = 0;
    for ( frost_area = 0; frost_area < FROST_SUBAREAS; frost_area++ )
      avg_ice += layer[index].ice[frost_area] * frost_fract[frost_area];
    printf("\t%f",avg_ice);
  }
#else
  for(index=0;index<Nlayer;index++) printf("\t%f",layer[index].ice);
#endif
  printf("\n---------------------------------------------------------------------------\nLayer Moist:\t");
  sum_moist = 0.;
  for(index=0;index<Nlayer;index++) {
    layer_moist = layer[index].moist;
    sum_moist += layer_moist;
    printf("\t%f",layer_moist);
  }
  printf("\n\n-----> Total Moisture = %f\n\n",sum_moist);
}