void
ABec4::compFlux (D_DECL(MultiFab &xflux, MultiFab &yflux, MultiFab &zflux),
MultiFab& in, const BC_Mode& bc_mode,
int src_comp, int dst_comp, int num_comp, int bnd_comp)
{
compFlux(D_DECL(xflux, yflux, zflux), in, true, bc_mode, src_comp, dst_comp, num_comp, bnd_comp);
}
int main(int argc, char *argv[])
{
FILE *inPtr, *alphaPtr;
int i, k, m, nin, ei, ti;
char *dir, sysCmd[512], filename[64], *NS, alphaFile[128];
float L, *arr, *arrarr[128]; // <- array of pointers to precip arrays
float J[100][100];
double Jext;
// Get directory name from command line
if(argc != 3) {
printf("\n\aWrong number of input arguments\n");
exit(0);
}
dir = argv[1];
L_TARG = atof(argv[2]);
NUM_TIMES = (int)ceil( RES_FINT / RES_DT );
sprintf( alphaFile, "%s/alpha_%g_%s", dir, L_TARG, "N" );
arr = getArr(NUM_E, NUM_TIMES);
if( (alphaPtr=fopen(alphaFile, "r"))!=NULL ) {
nin = fread(arr,sizeof(float),(NUM_E*NUM_TIMES),alphaPtr);
fclose(alphaPtr);
arrarr[0] = arr;
sprintf( alphaFile, "%s/alpha_%g_%s", dir, L_TARG, "S" );
if((alphaPtr=fopen(alphaFile, "r"))==NULL)
printf("\n\aprob opening %s\n", alphaFile);
arr = getArr(NUM_E, NUM_TIMES);
nin = fread(arr,sizeof(float),(NUM_E*NUM_TIMES),alphaPtr);
fclose(alphaPtr);
arrarr[1] = arr;
} else {
printf("\n%s cannot be opened\n", alphaFile);
} // if alphaFile exists
for(k=0; k<2; k++) {
compFlux(arrarr[k], L_TARG, k, dir);
} // N/S - hemisphere
return 0;
}
