Changeset 14874

Timestamp:

05/03/13 15:55:08 (12 years ago)

Author:

Mathieu Morlighem

Message:

CHG: reverted back (but removed check on np>1 if no petsc

File:

• issm/trunk-jpl/src/m/solve/solve.py (modified) (1 diff)

issm/trunk-jpl/src/m/solve/solve.py

@@ +1 @@
+import datetime
+import os
+import shutil
+from pairoptions import *
+from process_solve_options import *
+from EnumDefinitions import *
+from EnumToString import EnumToString
+from ismodelselfconsistent import *
+from marshall import *
+from ispetsc import *
+from waitonlock import *
+from loadresultsfromcluster import *
+
+def solve(md,solutionenum,*args):
+        """
+        SOLVE - apply solution sequence for this model
+ 
+           Usage:
+              md=solve(md,solutionenum,varargin)
+              where varargin is a list of paired arguments of string OR enums
+ 
+           solution types available comprise:
+              - DiagnosticSolutionEnum
+              - PrognosticSolutionEnum
+              - ThermalSolutionEnum
+              - SteadystateSolutionEnum
+              - TransientSolutionEnum...
+              - BalancethicknessSolutionEnum
+              - BedSlopeSolutionEnum
+              - SurfaceSlopeSolutionEnum
+              - HydrologySolutionEnum
+              - FlaimSolutionEnum
+ 
+           extra options:
+              - loadonly : does not solve. only load results
+ 
+           Examples:
+              md=solve(md,DiagnosticSolutionEnum);
+        """
+
+        #recover and process solve options
+        options=pairoptions('solution_type',solutionenum,*args)
+        options=process_solve_options(options)
+
+        #recover some fields
+        md.private.solution=options['solution_type']
+        cluster=md.cluster
+
+        #check model consistency
+        print "checking model consistency"
+        if solutionenum == FlaimSolutionEnum():
+                md.private.isconsistent=True
+                md.mesh.checkconsistency(md,solutionenum)
+                md.flaim.checkconsistency(md,solutionenum)
+                if not md.private.isconsistent:
+                        raise RuntimeError("Model not consistent, see messages above.")
+        else:
+                ismodelselfconsistent(md)
+
+        #First, build a runtime name that is unique
+        c=datetime.datetime.now()
+        md.private.runtimename="%s-%02i-%02i-%04i-%02i-%02i-%02i-%i" % (md.miscellaneous.name,c.month,c.day,c.year,c.hour,c.minute,c.second,os.getpid())
+
+        #if running qmu analysis, some preprocessing of dakota files using models
+        #fields needs to be carried out. 
+        if md.qmu.isdakota:
+                md=preqmu(md,options)
+
+        #flaim analysis
         if options['solution_type'] == FlaimSolutionEnum():
                 md=flaim_sol(md,options)