mergemod.c 2.73 KB
Newer Older
tilman.metz's avatar
tilman.metz committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97
/*------------------------------------------------------------------------
 * Copyright (C) 2011 For the list of authors, see file AUTHORS.
 *
 * This file is part of SOFI3D.
 * 
 * SOFI3D is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, version 2.0 of the License only.
 * 
 * SOFI3D is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 * 
 * You should have received a copy of the GNU General Public License
 * along with SOFI3D. See file COPYING and/or 
  * <http://www.gnu.org/licenses/gpl-2.0.html>.
--------------------------------------------------------------------------*/
/*------------------------------------------------------------------------
 *   merge model files written by the different processes to 
 *   a single file                                 
 *
 *  ----------------------------------------------------------------------*/

#include "fd.h"


void mergemod(char modfile[STRING_SIZE], int format){



	extern int NXG, NYG, MYID, NPROCX, NPROCY, NPROCZ;
	extern int NX, NY, NZ, NPROC, IDX, IDY, IDZ;
	extern FILE *FP;


	char file[STRING_SIZE];
	FILE *fp[NPROCY_MAX][NPROCX_MAX][NPROCZ_MAX], *fpout;
	int i, j, k, ip, jp, kp;
	float a;


	if ((NPROCX>NPROCX_MAX)||(NPROCY>NPROCY_MAX)||(NPROCZ>NPROCZ_MAX))
		err(" merge.c: constant expression NPROC?_MAX < NPROC? ");


	printf("\n **Message from mergemod (printed by PE %d):\n",MYID);
	printf(" PE %d starts merge of %d model files \n",MYID,NPROC);	

	fprintf(FP,"\n writing merged model file to  %s \n",modfile);
	fpout=fopen(modfile,"w");


	for (kp=0;kp<=NPROCZ-1; kp++)
		for (ip=0;ip<=NPROCX-1; ip++)
			for (jp=0;jp<=NPROCY-1; jp++){
				sprintf(file,"%s.%i%i%i",modfile,ip,jp,kp);
				fp[jp][ip][kp]=fopen(file,"r");
				if (fp[jp][ip][kp]==NULL) err("merge: can't read modfile !");
			}

	fprintf(FP," ... finished. \n");



	fprintf(FP," Copying...");

	for (kp=0;kp<=NPROCZ-1; kp++)
		for (k=1;k<=NZ;k+=IDZ)
			for (ip=0;ip<=NPROCX-1; ip++)
				for (i=1;i<=NX;i+=IDX)
					for (jp=0;jp<=NPROCY-1; jp++)
						for (j=1;j<=NY;j+=IDY){
							a=readdsk(fp[jp][ip][kp],format);
							writedsk(fpout,a,format);
						}

	fprintf(FP," ... finished. \n");

	for (kp=0;kp<=NPROCZ-1; kp++)
		for (ip=0;ip<=NPROCX-1; ip++)
			for (jp=0;jp<=NPROCY-1; jp++){
				fclose(fp[jp][ip][kp]);
			}


	fclose(fpout);

	fprintf(FP," Use \n");
	fprintf(FP," xmovie n1=%d n2=%d < %s loop=1 label1=\"Y (vertical)\" label2=\"X (horizontal)\" title=%%g \n",
			((NYG-1)/IDY)+1,((NXG-1)/IDX)+1,modfile);
	fprintf(FP," to visualize model. \n");


}