# $Id$ from .solver import Solver class Marc(Solver): #-------------------------- def __init__(self): #-------------------------- self.solver = 'Marc' self.releases = { \ '2012': ['linux64',''], '2011': ['linux64',''], '2010.2':['linux64',''], '2010': ['linux64',''], '2008r1':[''], '2007r1':[''], '2005r3':[''], } #-------------------------- def version(self,rootRelation = ''): #-------------------------- import os,damask.environment MSCpath = damask.environment.Environment(rootRelation).pathInfo['MSC_ROOT'] for release,subdirs in sorted(self.releases.items(),reverse=True): for subdir in subdirs: path = '%s/mentat%s/shlib/%s'%(MSCpath,release,subdir) if os.path.exists(path): return release else: continue return '' #-------------------------- def libraryPath(self,rootRelation = '',releases = []): #-------------------------- import os,damask.environment MSCpath = damask.environment.Environment(rootRelation).pathInfo['MSC_ROOT'] if len(releases) == 0: releases = self.releases.keys() if type(releases) is not list: releases = [releases] for release in sorted(releases,reverse=True): if release not in self.releases: continue for subdir in self.releases[release]: libPath = '%s/mentat%s/shlib/%s'%(MSCpath,release,subdir) if os.path.exists(libPath): return libPath else: continue return '' #-------------------------- def toolsPath(self,rootRelation = '',release = ''): #-------------------------- import os,damask.environment MSCpath = damask.environment.Environment(rootRelation).pathInfo['MSC_ROOT'] if len(release) == 0: release = self.version(rootRelation) path = '%s/marc%s/tools'%(MSCpath,release) if os.path.exists(path): return path else: return '' #-------------------------- def submit_job(self, #-------------------------- rootRelation = '', release = '', model = 'model', job = 'job1', logfile = None, compile = False, optimization ='', openMP = False ): import os,damask.environment import subprocess,shlex,shutil if len(release) == 0: release = self.version(rootRelation) damaskEnv = damask.environment.Environment(rootRelation) user = os.path.join(damaskEnv.relPath('code/'),'DAMASK_marc') # Define options [see Marc Installation and Operation Guide, pp 23] if int(release) >= 2013: script = 'run_damask%s'%({False:'',True:'_'}[optimization!='' or openMP]) script = script+'%s%s'%({False:'',True:optimization}[optimization!=''],{False:'',True:'mp'}[openMP]) else: script = 'run_marc' cmd = os.path.join(self.toolsPath(rootRelation,release),script) + \ ' -jid ' + model + '_' + job + \ ' -nprocd 1 -autorst 0 -ci n -cr n -dcoup 0 -b no -v no' if compile: cmd += ' -u ' + user+'.f90' + ' -save y' else: cmd += ' -prog ' + user print('job submission with%s compilation: %s'%({False:'out',True:''}[compile],user)) if logfile: log = open(logfile, 'w') print(cmd) self.p = subprocess.Popen(shlex.split(cmd),stdout = log,stderr = subprocess.STDOUT) log.close() self.p.wait() #-------------------------- def exit_number_from_outFile(self,outFile=None): #-------------------------- import string exitnumber = -1 fid_out = open(outFile,'r') for ln in fid_out: if (string.find(ln,'tress iteration') is not -1): print(ln) elif (string.find(ln,'Exit number') is not -1): substr = ln[string.find(ln,'Exit number'):len(ln)] exitnumber = int(substr[12:16]) fid_out.close() return exitnumber