Changeset 25157
- Timestamp:
- 06/25/20 09:27:18 (5 years ago)
- Location:
- issm/trunk-jpl/src/m
- Files:
-
- 9 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/classes/giaivins.py
r25137 r25157 44 44 45 45 def checkconsistency(self, md, solution, analyses): #{{{ 46 if not 'GiaAnalysis'in analyses:46 if 'GiaAnalysis' not in analyses: 47 47 return md 48 48 -
issm/trunk-jpl/src/m/classes/giamme.py
r25130 r25157 42 42 43 43 def checkconsistency(self, md, solution, analyses): # {{{ 44 if not 'SealevelriseAnalysis' in analyses or solution == 'TransientSolution' and md.transient.istidallove == 0:44 if ('SealevelriseAnalysis' not in analyses) or (solution == 'TransientSolution' and md.transient.istidallove == 0): 45 45 return md 46 46 -
issm/trunk-jpl/src/m/classes/lovenumbers.py
r25156 r25157 1 1 import numpy as np 2 2 3 from checkfield import checkfield 3 4 from fielddisplay import fielddisplay … … 9 10 LOVENUMBERS numbers class definition 10 11 11 Usage:12 lovenumbers = lovenumbers() #will setup love numbers deg 1001 by default13 lovenumbers = lovenumbers('maxdeg', 10001); #supply numbers of degrees required (here, 10001)12 Usage: 13 lovenumbers = lovenumbers() #will setup love numbers deg 1001 by default 14 lovenumbers = lovenumbers('maxdeg', 10001); #supply numbers of degrees required (here, 10001) 14 15 ''' 15 16 16 17 def __init__(self, *args): #{{{ 17 18 18 #regular love numbers: 19 19 self.h = [] #provided by PREM model … … 26 26 self.tl = [] 27 27 self.tk2secular = 0 #deg 2 secular number. 28 # 28 #}}} 29 29 30 def setdefaultparameters(self, maxdeg, referenceframe): # {{{ 31 30 def setdefaultparameters(self, maxdeg, referenceframe): #{{{ 32 31 #initialize love numbers: 33 32 self.h = getlovenumbers('type', 'loadingverticaldisplacement', 'referenceframe', referenceframe, 'maxdeg', maxdeg) … … 40 39 #secular fluid love number: 41 40 self.tk2secular = 0.942 42 # 41 #}}} 43 42 44 def checkconsistency(self, md, solution, analyses): # {{{ 45 46 if 'SealevelriseAnalysis' not in analyses or (solution == 'TransientSolution' and md.transient.isslr == 0): 43 def checkconsistency(self, md, solution, analyses): #{{{ 44 if ('SealevelriseAnalysis' not in analyses) or (solution == 'TransientSolution' and md.transient.isslr == 0): 47 45 return 48 46 … … 58 56 if (size(self.h, 0) != size(self.k, 0) | size(self.h, 0) != size(self.l, 0)): 59 57 raise ValueError('lovenumbers error message: love numbers should be provided at the same level of accuracy') 60 # 58 #}}} 61 59 62 def defaultoutputs(self, md): #{{{60 def defaultoutputs(self, md): #{{{ 63 61 return[] 64 # 62 #}}} 65 63 66 def __repr__(self): #{{{64 def __repr__(self): #{{{ 67 65 s = ' lovenumbers parameters:' 68 66 … … 73 71 s += '{}\n'.format(fielddisplay(self, 'tk', 'tidal load Love number (deg 2)')) 74 72 s += '{}\n'.format(fielddisplay(self, 'tk2secular', 'secular fluid Love number')) 75 # 73 #}}} 76 74 77 def marshall(self, prefix, md, fid): # {{{ 78 75 def marshall(self, prefix, md, fid): #{{{ 79 76 WriteData(fid, prefix, 'object', self, 'fieldname', 'h', 'name', 'md.solidearth.lovenumbers.h', 'format', 'DoubleMat', 'mattype', 1) 80 77 WriteData(fid, prefix, 'object', self, 'fieldname', 'k', 'name', 'md.solidearth.lovenumbers.k', 'format', 'DoubleMat', 'mattype', 1) … … 85 82 WriteData(fid, prefix, 'object', self, 'fieldname', 'tl', 'name', 'md.solidearth.lovenumbers.tl', 'format', 'DoubleMat', 'mattype', 1) 86 83 WriteData(fid, prefix, 'object', self, 'data', self.tk2secular, 'fieldname', 'lovenumbers.tk2secular', 'format', 'Double') 87 # 84 #}}} 88 85 89 def extrude(self, md): #{{{86 def extrude(self, md): #{{{ 90 87 return 91 # 88 #}}} -
issm/trunk-jpl/src/m/classes/mesh3dsurface.py
r25129 r25157 114 114 md = checkfield(md, 'fieldname', 'mesh.average_vertex_connectivity', '>=', 9, 'message', '"mesh.average_vertex_connectivity" should be at least 9 in 2d') 115 115 116 if (solution == 'ThermalSolution'):116 if solution == 'ThermalSolution': 117 117 md = md.checkmessage('thermal not supported for 2d mesh') 118 118 -
issm/trunk-jpl/src/m/classes/rotational.py
r25125 r25157 48 48 49 49 def checkconsistency(self, md, solution, analyses): # {{{ 50 if not 'SealevelriseAnalysis' in analyses or solution == 'TransientSolution' and md.transient.isrotational == 0:50 if ('SealevelriseAnalysis' not in analyses) or (solution == 'TransientSolution' and md.transient.isrotational == 0): 51 51 return md 52 52 -
issm/trunk-jpl/src/m/classes/solidearth.py
r25153 r25157 57 57 58 58 def checkconsistency(self, md, solution, analyses): # {{{ 59 if 'SealevelriseAnalysis' not in analysesor (solution == 'TransientSolution' and md.transient.isslr == 0):59 if ('SealevelriseAnalysis' not in analyses) or (solution == 'TransientSolution' and md.transient.isslr == 0): 60 60 return md 61 61 -
issm/trunk-jpl/src/m/classes/solidearthsettings.py
r25125 r25157 66 66 67 67 def checkconsistency(self, md, solution, analyses): # {{{ 68 if not 'SealevelriseAnalysis' in analyses or solution == 'TransientSolution' and md.transient.issolidearthsettings == 0:68 if ('SealevelriseAnalysis' not in analyses) or (solution == 'TransientSolution' and md.transient.issolidearthsettings == 0): 69 69 return md 70 70 -
issm/trunk-jpl/src/m/classes/surfaceload.py
r25125 r25157 41 41 42 42 def checkconsistency(self, md, solution, analyses): # {{{ 43 if not 'SealevelriseAnalysis' in analyses or solution == 'TransientSolution' and md.transient.issurfaceload == 0:43 if ('SealevelriseAnalysis' not in analyses) or (solution == 'TransientSolution' and md.transient.issurfaceload == 0): 44 44 return md 45 45 -
issm/trunk-jpl/src/m/miscellaneous/fielddisplay.py
r25136 r25157 1 1 import numpy as np 2 2 3 import MatlabFuncs as m 3 4 4 5 5 6 def fielddisplay(md, name, comment): 6 """7 ''' 7 8 FIELDDISPLAY - display model field 8 9 9 10 Usage: 10 11 fielddisplay(md, name, comment) 11 """12 ''' 12 13 13 14 #get field … … 18 19 19 20 20 def parsedisplay(offset, name, field, comment): # {{{ 21 21 def parsedisplay(offset, name, field, comment): #{{{ 22 22 #string 23 23 if isinstance(field, str): … … 59 59 60 60 return string 61 # 61 #}}} 62 62 63 63 64 def dict_display(offset, name, field, comment): # {{{ 65 64 def dict_display(offset, name, field, comment): #{{{ 66 65 if field: 67 66 string = displayunit(offset, name, '{dictionary}', comment) + '\n' … … 78 77 79 78 return string 80 # 79 #}}} 81 80 82 81 83 def list_display(offset, name, field, comment): # {{{ 84 82 def list_display(offset, name, field, comment): #{{{ 85 83 #initialization 86 84 if isinstance(field, list): … … 110 108 #call displayunit 111 109 return displayunit(offset, name, string, comment) 112 # 110 #}}} 113 111 114 112 115 def displayunit(offset, name, characterization, comment): # {{{ 116 113 def displayunit(offset, name, characterization, comment): #{{{ 117 114 #take care of name 118 115 if len(name) > 23: … … 140 137 141 138 return string 142 # 139 #}}}
Note:
See TracChangeset
for help on using the changeset viewer.