plate()

int plate

(

struct geometry *

geom

)

Plate geometry generation.

Parameters

[in,out]

geom

geometry structure

                                {
  DOUBLE l1,dl1;
  DOUBLE l2,dl2;
  int i,j;
  int inodb,ielmb,nelmb,ivbase;
  struct plate_geom* plategeom ;

  int i1,i2;

  logger(LOG_DEBUG,"building plate\n");
  plategeom = (struct plate_geom*) geom->ptrgeom;

  // steps
  dl1   = plategeom->lx/(plategeom->nx);
  dl2   = plategeom->ly/(plategeom->ny);

  inodb=0;
  l1 = -plategeom->lx/2;
  /* building geometry */
  for(j=0;j<(plategeom->nx+1);j++){
    l2 = plategeom->ly/2;
    for(i=0;i<(plategeom->ny+1);i++){
      geometry->nodes[inodb + _g_inodb].x = l1;
      geometry->nodes[inodb + _g_inodb].y = l2;
      geometry->nodes[inodb + _g_inodb].z = 0.0;
    
      l2 -= dl2;
      inodb++;
    }
    l1 += dl1;
  }

  /* building topology */
  nelmb = (plategeom->nx)*(plategeom->ny);
  ivbase=0;
  for(ielmb=0;ielmb<nelmb;ielmb++){
    i1 = (ielmb)/(plategeom->ny);
    i2 = i1+1;

    geometry->jnodb[ivbase+0 + _g_ielmb] = ielmb + i1 + 1;
    geometry->jnodb[ivbase+1 + _g_ielmb] = ielmb + i1 + 1 + plategeom->ny + 1;
    geometry->jnodb[ivbase+2 + _g_ielmb] = ielmb + i1     + plategeom->ny + 1;
    geometry->jnodb[ivbase+3 + _g_ielmb] = ielmb + i1;

    ivbase +=4;
  }

   geom_fill(geom); // filling geometry

   // adding to global vars
  _g_icntr += geom->ncntr;
  _g_inodb += geom->nnodb;
  _g_ielmb += geom->nelmb;

  logger(LOG_DEBUG,"plate built\n");

 return 1;

}