rain_inp.c File Reference

Establish linkage between two grids, map rainfall input data. More...

#include "rain.h"

Include dependency graph for rain_inp.c:

Go to the source code of this file.

Functions
int	rain_data_wmm (float *rainSME)

---

Detailed Description

Establish linkage between two grids, map rainfall input data.

This has a single function.

Note: documented with Doxygen, which expects specific syntax within special comments.

The Everglades Landscape Model (ELM).
last updated: Dec 2004

Definition in file rain_inp.c.

---

Function Documentation

int rain_data_wmm

(

float *

)

rainfall data array

Definition at line 26 of file rain_inp.c.

References conv_inTOtenths_mm, applicationStruct::dataELM, applicationStruct::day, debug, elm_OG_map, elm_wmm_map, initDataStruct(), mapGrids(), ModelPath, msgStr, printGridMap(), processData(), ProjName, rain_binfilename, rain_struct, applicationStruct::recRead, returnData(), s0, s1, SimTime, T, simTime::TIME, usrErr(), and usrErr0().

Referenced by cell_dyn1().

{

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

  if(SimTime.TIME==0)  {      

     // cleanUp(*elm_wmm_map, rainSME);
      
                /* 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 rain data...");
          usrErr0(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,"rain_data_wmm==> Finished mapping grids");
          usrErr (msgStr);
      }

      sprintf(rain_binfilename, "%s%s/Data/rain.BIN", ModelPath, ProjName );

/* initializing data structures, move pointer to initial date (gridmap.c) */
      stat = initDataStruct(rain_binfilename,&rain_struct);

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

  } /* end of SimTime.TIME=0 */
  
  
  if(rain_struct.day >= rain_struct.recRead) {      /* process the data in batch */
      sprintf(msgStr,"\nProcessing batch of rain data...");
      usrErr (msgStr); 
      stat = processData(rain_binfilename,&rain_struct);

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


  if (rain_struct.day < rain_struct.recRead) {      /* pass the data day by day */
    returnData(rainSME,&rain_struct);   
    /* change the unit here */
    for(i = 0; i < s0; i++) {
      for(j = 0; j < s1; j++) {
        k = i*s1+j;
        rainSME[T((i+1),(j+1))] = rain_struct.dataELM[k] * conv_inTOtenths_mm;  /* convert data from inches to tenths of mm */
      }
    } 

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

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

Here is the call graph for this function: