issm_ocean.cpp File Reference

: ISSM OCEAN main program. More...

#include "./issm.h"
#include <stdlib.h>

Go to the source code of this file.

Functions
int	main (int argc, char **argv)

Detailed Description

: ISSM OCEAN main program.

Definition in file issm_ocean.cpp.

Function Documentation

main()

int main	(	int	argc,
		char **	argv
	)

Definition at line 8 of file issm_ocean.cpp.

                              {
 
   /*diverse:*/
   int    icecommsize;
   int    my_rank,my_local_rank,my_size,my_local_size;
   ISSM_MPI_Comm worldcomm;
   ISSM_MPI_Comm modelcomm;
   ISSM_MPI_Comm frommitgcm;
   ISSM_MPI_Comm tomitgcmcomm;
   ISSM_MPI_Status status;
 
   /*Initialize exception trapping: */
   ExceptionTrapBegin();
 
   /*Initialize environment (MPI, PETSC, MUMPS, etc ...)*/
   worldcomm=EnvironmentInit(argc,argv);
 
   /*What is my rank?:*/
   ISSM_MPI_Comm_rank(worldcomm,&my_rank);
   ISSM_MPI_Comm_size(worldcomm,&my_size);
 
   /*First model is ice, second is ocean*/
   /*ice comm size: */
   icecommsize=(int) strtol(argv[2], (char **)NULL, 10);
 
   /*Split world into sub-communicators for each and every model:*/
   ISSM_MPI_Comm_split(worldcomm,0, my_rank, &modelcomm);
   ISSM_MPI_Comm_rank(modelcomm,&my_local_rank);
   ISSM_MPI_Comm_size(modelcomm,&my_local_size);
 
   ISSM_MPI_Intercomm_create( modelcomm, 0, worldcomm, my_local_size, 0, &tomitgcmcomm); 
 
   FemModel *femmodel = new FemModel(argc,argv,modelcomm);
 
   /*Now that the models are initialized, keep communicator information in the parameters datasets of each model: */
   femmodel->parameters->AddObject(new GenericParam<ISSM_MPI_Comm>(worldcomm,WorldCommEnum));
   femmodel->parameters->AddObject(new GenericParam<ISSM_MPI_Comm>(tomitgcmcomm,ToMITgcmCommEnum));
 
   /*Solve: */
   femmodel->Solve();
 
   /*Output results: */
   OutputResultsx(femmodel);
 
   /*Wrap up: */
   femmodel->CleanUp();
 
   /*Delete Model: */
   delete femmodel;
 
   /*Finalize environment:*/
   EnvironmentFinalize();
 
   /*Finalize exception trapping: */
   ExceptionTrapEnd();
 
   /*Free ressources:*/
 
   /*Return unix success: */
   return 0; 
}