UGA Stratigraphy Lab

The data is in the strata

Patching Strata for use with Biostrat

Anil Deshpande and Peter Flemings of Penn State University have graciously provided a routine that allows Strata to produce water depth histories that can be used by Biostrat.

Instructions

Copy source code file mdump.c (see below) to the simbasin directory of Strata

Make the following two edits to the basmod.c file, located in the simbasin directory:

1) line 37, insert:

extern void mdump(void);

2) line 410, insert, following the line moviesave(sock): } and preceding the first of three lines beginning if(strcasecmp(.:

mdump();

Make one edit to the Makefile.in file, located in the simbasin directory:

1)line 48, insert, following the line isotope.o \ and preceding the line moviesave.o:

mdump.o \

Re-make simbasin.

mdump.c code:

/* ---------------------------------------------------------------
Function writes the elev[ ] array to the biostrat.wd file at every timeloop.

timeline[i][j][k] where:
i is 0 for depth to horizon, and 2 for bathymetry and is positive
  downwards.
j is timeline number where the deepest horizon has a j of 0
  and shallowest horizon has a j of timeloop[0]
  k is node number (lateral position)

kNodeSpacing determines spacing of outcrops; e.g., 10 means a water depth history for each 10th node across basin will be written to file
--------------------------------------------------------------- */

#include <stdio.h>
#include <stdlib.h>
#include "header.h"
#define kNodeSpacing 10

void mdump()
{
  int j,k;
  FILE *fpl;
  fpl=fopen("biostrat.wd","w");
  for (k=0;k<nnodes;k+=kNodeSpacing)
  {
    for (j=0; j<timeloop[0]; ++j)
    {
      fprintf(fpl, "%d %9.0lf %9.2lf %9.2lf\n",k,(j+1)*tlstep*
      timedelt, timeline[0][j][k], timeline[2][j][k]);
    }
  }
  fclose(fpl);
}

Correcting compile warnings for strata

In the setbasin directory, edit the following files:

1) In start.c, line 110, change void main(int argc, char **argv) to int main(int argc, char **argv)

2) In start.c, after line 597 XtAppMainLoop(app);, insert return 0;

3) In goodbye.c, after line 7 #include <stdio.h>, insert #include <stdlib.h>

4) In misc.c, after line 7 #include <stdio.h>, insert #include <stdlib.h>

In the setbasin directory, edit this file:

1) In openfile.c, after line 8 #include <stdlib.h>, insert #include <strings.h>

In the plotbasin directory, edit the following files:

1) In gif_here.c, after line 57 #include <stdio.h>, insert #include <stdlib.h>

2) In xwd_here.c, after line 78 #include <stdio.h>, insert #include <stdlib.h>

3) In xwd_dsimple.c, after line 6 #include <stdio.h>, insert #include <stdlib.h>

4) In start.c, line 123, change from void main(int argc, char **argv) to int main(int argc, char **argv)

5) In start.c, after line 606 XtAppMainLoop(app);, insert return 0;

In the filmbasin directory, make the following edits:

1) In film.c, line 60, change from void main(int argc, char **argv) to int main(int argc, char **argv)

2) In film.c, after lines 169-170 myexec(temp); }, insert return 0;