evap_inp.c

Go to the documentation of this file.

/* General notes on revisions to this source file. 
       Nov/Dec 2004 v2.3.2: documentation upgrade 
                - Doxygen tags added 
        
*/
#include "evap.h"


/* \brief Get spatial input time series of Potential ET 

        Using the grid_map functions, this obtains potential 
        evapotranspiration (ET), at ELM (i.e., SME) grid scale. 
        \param evapSME potential ET data array (tenths mm) at the ELM/SME grid scale
        \return success         */




int evap_data_wmm(float* evapSME)           /* array of float of size s0*s1 */
{

  int i,j,k;
  int success = 1, fail = -1;
  int stat = success;
  char gridmapfilename[335];

  if(SimTime.TIME==0)  {      
/* elm_OG_map is data structure containing the mapping attributes at two scales */
    elm_wmm_map = elm_OG_map;

    if(elm_wmm_map == NULL) {
      sprintf(msgStr, "Mapping grids and setting up pET data...");
      usrErr (msgStr);

      sprintf(gridmapfilename, "%s%s/Data/gridmapping.txt", ModelPath, ProjName );   
      stat = mapGrids(gridmapfilename);
      elm_wmm_map = elm_OG_map;
    }

    if(debug > 4) {
      printGridMap();
      sprintf(msgStr,"evap_data_wmm==> Finished mapping grids");
      usrErr (msgStr);
    }

    sprintf(evap_binfilename, "%s%s/Data/ETp.BIN", ModelPath, ProjName );
    /* initializing data structures, move pointer to initial date (gridmap.c) */
    stat = initDataStruct(evap_binfilename,&evap_struct);  

    if(debug > 4) {
      /*printELM2Grid_io(); */
      /*drawELM2Grid_io(); */
      sprintf(msgStr,"evap_data_wmm==> Finished initializing");
      usrErr (msgStr); 
    }

  } /* end of SimTime.TIME=0 */
  
  
  if(evap_struct.day >= evap_struct.recRead) {      /* process the data in batch */
      sprintf(msgStr,"Processing batch of pET data...");
      usrErr (msgStr); 
      stat = processData(evap_binfilename,&evap_struct);

      if(debug > 4 ) {
        /*printBatchData(evapWMM,gridio_batch_len,widCnt);*/ /* TODO: remove this printBatchData function when sure is no longer needed */
        sprintf(msgStr,"evap_data_wmm==> Finished processing data");
        usrErr (msgStr); 
      }
  } /* end of if */


  if(evap_struct.day < evap_struct.recRead) {      /* pass the data day by day */
    returnData(evapSME,&evap_struct); 

    for(i = 0; i < s0; i++) {
      for(j = 0; j < s1; j++) {
        k = i*s1+j;
        evapSME[T((i+1),(j+1))] = evap_struct.dataELM[k] * conv_inTOtenths_mm; /* convert data from inches to tenths of mm */  
      }
    } 

    if(debug > 4) {
      sprintf(msgStr,"evap_data_wmm==> Finished returning data");
      usrErr (msgStr); 
    }

  } /* end of if */
    
  return success;
}