



//======================================================================
//Parameters

const char a_record = '0';

//const double beta = .02*.1*64/4; // Basal
//const double gama = .1; // Side
//const double zeta = 1.0/8;


const double dx = .001;

const int T_end = 100000;
const int T_record = 1000;
const int T_display = 1000;

// Flags to define the conditions of the simulation
const int Sheet_condition = 2; // 1 for bent sheet, 2 for flat sheet

const int fixed_y = 1; // flag for fixing y (1 fixed, 0 not fixed)
const int fixed_first_last_cell_z = 0;// flag for fixing first 2 and last 2 cells along z
const int fixed_first_last_vert_xyz = 0; // flag for fixing the first/last verts along xyz 
const int fixed_top = 1; // flag for introducingthe apical membrane

const int yolk_present = 0; //  flag for the yolk
const int plane_constraint_present = 1; // flag that the plane is contraining

const int neighbor_apical_active = 0; // flag for having the neighboring cells also constrict
const int neighbor_basal_strong = 0; // flag for having the neighboring cells have strong base

const int neighbor_cell_N = 2; // Number of neighbor cells

const int gaussian_KA_flag = 1; // flag for having a gaussian distribution of 

const int midline_wall_flag = 1; // flag for having the midline act like a wall

const int Percent_stop_flag = 0;
const double Percent_stop = .0001;

//======================================================================




