segy.c

阐述了建立物质本构关系时应满足的基本公理, 较详细地推导了各种类型粘 弹性物质的本构关系, 论述了变率理论与泛函理论之间的关系, 推导了对于粘弹性问题 求解具有实用意义的对应定理, 并以简单的例子

#include "segy.h"

int bhdrupd_(bhdr,sr,spt,ntr)
     struct reel_binhdr_struct  *bhdr;
     int *sr;
     int *spt;
     int *ntr;
{
  int status;

  bhdr->job_id = 1;
  bhdr->line_number = 1;
  bhdr->reel_number = 1;
  bhdr->data_traces_per_rec = *ntr;
  bhdr->aux_traces_per_rec = 0;
  bhdr->sample_interval = *sr;
  bhdr->sample_interval_orig = *sr;
  bhdr->data_sample_format = 1;
  bhdr->CDP_fold = 1;
  bhdr->samples_per_trace = *spt;
  bhdr->samples_per_trace_orig = *spt;

  status = 0;
  return status;
}
int bhdrget_(bhdr,sr,spt,ntr)
     struct reel_binhdr_struct  *bhdr;
     int *sr;
     int *spt;
     int *ntr;
{
  int status;

  *ntr = bhdr->data_traces_per_rec;
  *sr =  bhdr->sample_interval;
  *spt =  bhdr->samples_per_trace;

  status = 0;
  return status;
}

int thdrupd_(thdr,orec,sr,spt)
     struct trace_header_struct  *thdr;
     int *orec;
     int *sr;
     int *spt;
{
  int status;

  thdr->reel_trace_number = *orec;
  thdr->line_trace_number = *orec;
  thdr->orig_field_rec_num = *orec;
  thdr->data_use = 1;
  thdr->num_samples_in_trace = *spt;
  thdr->sample_interval = *sr;

  status = 0;
  return status;
}

int thdrput_(thdr,orec,otrc,offset,sr,spt)
     struct trace_header_struct  *thdr;
     int *orec;
     int *otrc;
     int *offset;
     int *sr;
     int *spt;
{
  int status;

  thdr->reel_trace_number = *orec;
  thdr->line_trace_number = *orec;
  thdr->source_to_rec_grp_dist = *offset;
  thdr->orig_field_rec_num = *orec;
  thdr->orig_trace_num_in_rec = *otrc;
  thdr->CDP_ensemble_num = *orec;     
  thdr->trace_num_in_CDP= *otrc;
  thdr->trace_id_code = 1;
  thdr->src_coord_x = 0;
  thdr->src_coord_y = 0;
  thdr->grp_coord_x = *offset;
  thdr->grp_coord_y = 0;
  thdr->data_use = 1;
  thdr->num_samples_in_trace = *spt;
  thdr->sample_interval = *sr;

  status = 0;
  return status;
}
int thdrget_(thdr,orec,sr,spt)
     struct trace_header_struct  *thdr;
     int *orec;
     int *sr;
     int *spt;
{
  int status;

  *orec = thdr->reel_trace_number;
  *spt = thdr->num_samples_in_trace;
  *sr = thdr->sample_interval;

  status = 0;
  return status;
}
int txyget_(thdr,sx,rx,sz,rz)
     struct trace_header_struct  *thdr;
     float *sx;
     float *rx;
     float *sz;
     float *rz;

{
  int status;

  *sx = thdr->src_coord_x;
  *rx = thdr->grp_coord_x;
  *sz = thdr->src_coord_y;
  *rz = thdr->grp_coord_y;

  status = 0;
  return status;
}
int txyput_(thdr,sx,rx,sz,rz)
     struct trace_header_struct  *thdr;
     float *sx;
     float *rx;
     float *sz;
     float *rz;

{
  int status;

thdr->src_coord_x = *sx;
thdr->grp_coord_x = *rx;
thdr->src_coord_y = *sz;
thdr->grp_coord_y = *rz;

  status = 0;
  return status;
}