#include  <stdio.h>
#include  <string.h>
#include <math.h>
#include <stdlib.h>
#include <time.h>
#include  <ctype.h>
#include <errno.h>

//---------------------------------------------------------
int main (int argc, char *argv[]) {

int nx=11, ny=11, nz=301;  
float x,z;

float v[500]={0.0}, curx[500]={0.0}; 
float dz=0.1;  

float alfa=0.1, u0=0.3, gama=0.05;  

FILE *stream_out;
char*outfile_out="a01_u03_gama005.rsf"; 
char*outfilex="a01_u03_gama005_x.rsf";


//--------------------  c(z)=a*z, z>0  ------------------------
for( int j=0; j<nz; j++) {
   z=j*dz;     
   v[j]=z*alfa;

   curx[j]=u0;   
   curx[j]=u0*(1.0+gama*z);  
}


if ( (stream_out  = fopen( outfile_out, "wb" )) == NULL  )  fprintf(stderr,"%s\n","Can not open output file!");
else {
	for( int i=0; i<ny; i++) 
	for( int j=0; j<nx; j++) 
         nz=fwrite(v,sizeof(float),nz,stream_out);
}
fclose( stream_out );  


if ( (stream_out  = fopen( outfilex, "wb" )) == NULL  )  fprintf(stderr,"%s\n","Can not open output file!");
else {
	for( int i=0; i<ny; i++) 
	for( int j=0; j<nx; j++) 
    	nz=fwrite(curx,sizeof(float),nz,stream_out);  
}
fclose( stream_out );

exit(0); 
}