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; }