write_forcing_file.c

超强的大尺度水文模拟工具

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

static char vcid[] = "$Id: write_forcing_file.c,v 1.2.2.1 2004/05/06 00:38:00 tbohn Exp $";

#if OUTPUT_FORCE
void write_forcing_file(atmos_data_struct *atmos,
			int                nrecs,
			outfiles_struct   *outfiles) 
/**********************************************************************
  write_forcing_file          Keith Cherkauer           July 19, 2000

  This routine writes the complete forcing data files for use in 
  future simulations.
**********************************************************************/
{
  extern option_struct options;

  int                 rec, j;
  short int          *tmp_siptr;
  unsigned short int *tmp_usiptr;

  for ( rec = 0; rec < nrecs; rec++ ) {
    for ( j = 0; j < NF; j++ ) {
      if(options.BINARY_OUTPUT) {
	/* Write a binary forcing file */
	tmp_siptr = (short int *)calloc(1,sizeof(short int));
	tmp_usiptr = (unsigned short int *)calloc(1,sizeof(unsigned short int));
	/* precipitation * 40 */
	tmp_usiptr[0] = (unsigned short int)(atmos[rec].prec[j]*40.);
	fwrite(tmp_usiptr,1,sizeof(unsigned short int),outfiles->fluxes);
	/* air temperature * 100 */
	tmp_siptr[0] = (short int)(atmos[rec].air_temp[j]*100.);
	fwrite(tmp_siptr,1,sizeof(short int),outfiles->fluxes);
	/* shortwave * 50 */
	tmp_usiptr[0] = (unsigned short int)(atmos[rec].shortwave[j]*50.);
	fwrite(tmp_usiptr,1,sizeof(unsigned short int),outfiles->fluxes);
	/* longwave * 80 */
	tmp_usiptr[0] = (unsigned short int)(atmos[rec].longwave[j]*80.);
	fwrite(tmp_usiptr,1,sizeof(unsigned short int),outfiles->fluxes);
	/* density * 100 */
	tmp_usiptr[0] = (unsigned short int)(atmos[rec].density[j]*100.);
	fwrite(tmp_usiptr,1,sizeof(unsigned short int),outfiles->fluxes);
	/* pressure * 100 */
	tmp_usiptr[0] = (unsigned short int)(atmos[rec].pressure[j]*100.);
	fwrite(tmp_usiptr,1,sizeof(unsigned short int),outfiles->fluxes);
	/* vp * 100 */
	tmp_siptr[0] = (short int)(atmos[rec].vp[j]*100.);
	fwrite(tmp_siptr,1,sizeof(short int),outfiles->fluxes);
	/* wind * 100 */
	tmp_usiptr[0] = (unsigned short int)(atmos[rec].wind[j]*100.);
	fwrite(tmp_usiptr,1,sizeof(unsigned short int),outfiles->fluxes);   

	free((char *)tmp_siptr);
	free((char *)tmp_usiptr);
      }
      else {
	/* Write an ASCII forcing file */
	fprintf(outfiles->fluxes,"%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\n", 
		atmos[rec].prec[j], atmos[rec].air_temp[j], 
		atmos[rec].shortwave[j], atmos[rec].longwave[j], 
		atmos[rec].density[j], atmos[rec].pressure[j], 
		atmos[rec].vp[j], atmos[rec].wind[j]);
      }
    }
  }
}
#endif /* OUTPUT_FORCE */