/*!\file: transient.cpp * \brief: transient solution */ #ifdef HAVE_CONFIG_H #include "config.h" #else #error "Cannot compile with HAVE_CONFIG_H symbol! run configure first!" #endif #include "../objects/objects.h" #include "../shared/shared.h" #include "../DataSet/DataSet.h" #include "../EnumDefinitions/EnumDefinitions.h" #include "../include/globals.h" #include "./parallel.h" int main(int argc,char* *argv){ /*I/O: */ FILE* fid=NULL; char* inputfilename=NULL; char* outputfilename=NULL; char* lockname=NULL; int numberofnodes; bool qmu_analysis=false; bool waitonlock=false; /*Model: */ Model* model=NULL; int dim=-1; /*Results: */ DataSet* results=NULL; DataSet* processed_results=NULL; Result* result=NULL; Param* param=NULL; /*time*/ double start, finish; double start_core, finish_core; double start_init, finish_init; MODULEBOOT(); #if !defined(_PARALLEL_) || (defined(_PARALLEL_) && !defined(_HAVE_PETSC_)) ISSMERROR(" parallel executable was compiled without support of parallel libraries!"); #endif /*Initialize Petsc and get start time*/ PetscInitialize(&argc,&argv,(char *)0,""); MPI_Barrier(MPI_COMM_WORLD); start=MPI_Wtime(); /*Size and rank: */ MPI_Comm_rank(MPI_COMM_WORLD,&my_rank); MPI_Comm_size(MPI_COMM_WORLD,&num_procs); _printf_("recover , input file name and output file name:\n"); inputfilename=argv[2]; outputfilename=argv[3]; lockname=argv[4]; /*Open handle to data on disk: */ fid=pfopen(inputfilename,"rb"); /*Initialize model structure: */ MPI_Barrier(MPI_COMM_WORLD); start_init=MPI_Wtime(); model=new Model(); _printf_("read and create finite element model:\n"); _printf_("\n reading diagnostic horiz model data:\n"); model->AddFormulation(fid,DiagnosticAnalysisEnum,HorizAnalysisEnum); _printf_("\n reading diagnostic vert model data:\n"); model->AddFormulation(fid,DiagnosticAnalysisEnum,VertAnalysisEnum); _printf_("\n reading diagnostic stokes model data:\n"); model->AddFormulation(fid,DiagnosticAnalysisEnum,StokesAnalysisEnum); _printf_("\n reading diagnostic hutter model data:\n"); model->AddFormulation(fid,DiagnosticAnalysisEnum,HutterAnalysisEnum); _printf_("\n reading surface and bed slope computation model data:\n"); model->AddFormulation(fid,SlopecomputeAnalysisEnum); _printf_("\n reading prognositc model data:\n"); model->AddFormulation(fid,PrognosticAnalysisEnum); /*Do we run in 3d?, in which case we need thermal and melting also:*/ model->FindParam(&dim,DimEnum); if(dim==3){ _printf_("read and create thermal finite element model:\n"); model->AddFormulation(fid,ThermalAnalysisEnum,TransientAnalysisEnum); _printf_("read and create melting finite element model:\n"); model->AddFormulation(fid,MeltingAnalysisEnum,TransientAnalysisEnum); } /*recover parameters: */ model->FindParam(&waitonlock,WaitOnLockEnum); model->FindParam(&qmu_analysis,QmuAnalysisEnum); _printf_("initialize results:\n"); results=new DataSet(ResultsEnum); MPI_Barrier(MPI_COMM_WORLD); finish_init=MPI_Wtime(); /*are we running the solution sequence, or a qmu wrapper around it? : */ if(!qmu_analysis){ /*run diagnostic analysis: */ _printf_("call computational core:\n"); MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( ); transient_core(results,model); MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( ); /*Add analysis_type to results: */ result=new Result(results->Size()+1,0,1,"analysis_type","transient"); results->AddObject(result); _printf_("process results:\n"); ProcessResults(&processed_results,results,model,TransientAnalysisEnum); _printf_("write results to disk:\n"); OutputResults(processed_results,outputfilename); } else{ /*run qmu analysis: */ _printf_("calling qmu analysis on transient core:\n"); #ifdef _HAVE_DAKOTA_ MPI_Barrier(MPI_COMM_WORLD); start_core=MPI_Wtime( ); Qmux(model,TransientAnalysisEnum,NoneAnalysisEnum); MPI_Barrier(MPI_COMM_WORLD); finish_core=MPI_Wtime( ); #else ISSMERROR(" Dakota not present, cannot do qmu!"); #endif } if (waitonlock>0){ _printf_("write lock file:\n"); WriteLockFile(lockname); } /*Free ressources:*/ delete results; delete processed_results; delete model; /*Get finish time and close*/ MPI_Barrier(MPI_COMM_WORLD); finish = MPI_Wtime( ); _printf_("\n %-34s %f seconds \n","Model initialization elapsed time:",finish_init-start_init); _printf_(" %-34s %f seconds \n","Core solution elapsed time:",finish_core-start_core); _printf_("\n %s %i hrs %i min %i sec\n\n","Total elapsed time:",int((finish-start)/3600),int(int(finish-start)%3600/60),int(finish-start)%60); _printf_("closing MPI and Petsc\n"); PetscFinalize(); /*end module: */ MODULEEND(); return 0; //unix success return; }