- Timestamp:
- 07/16/21 17:19:58 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/parameterization/setflowequation.py
r25609 r26352 6 6 7 7 def setflowequation(md, *args): 8 """ 9 SETFLOWEQUATION - associate a solution type to each element 10 11 This routine works like plotmodel: it works with an even number of inputs 12 'SIA', 'SSA', 'HO', 'L1L2', 'MLHO', 'FS' and 'fill' are the possible options 13 that must be followed by the corresponding exp file or flags list 14 It can either be a domain file (argus type, .exp extension), or an array of element flags. 15 If user wants every element outside the domain to be 16 setflowequationd, add '~' to the name of the domain file (ex: '~HO.exp') 17 an empty string '' will be considered as an empty domain 18 a string 'all' will be considered as the entire domain 19 You can specify the type of coupling, 'penalties' or 'tiling', to use with the input 'coupling' 20 21 Usage: 22 md = setflowequation(md, varargin) 23 24 Example: 25 md = setflowequation(md, 'HO', 'HO.exp', fill', 'SIA', 'coupling', 'tiling') 8 """SETFLOWEQUATION - associate a solution type to each element 9 10 This routine works like plotmodel: it works with an even number of inputs 11 'SIA', 'SSA', 'HO', 'L1L2', 'MLHO', 'FS' and 'fill' are the possible 12 options that must be followed by the corresponding exp file or flags list. 13 It can either be a domain file (argus type, .exp extension), or an array of 14 element flags. 15 If user wants every element outside the domain to be setflowequationd, add 16 '~' to the name of the domain file (ex: '~HO.exp') an empty string '' will 17 be considered as an empty domain a string 'all' will be considered as the 18 entire domain. 19 You can specify the type of coupling, 'penalties' or 'tiling', to use with 20 the input 'coupling'. 21 22 Usage: 23 md = setflowequation(md, varargin) 24 25 Example: 26 md = setflowequation(md, 'HO', 'HO.exp', fill', 'SIA', 'coupling', 'tiling') 26 27 """ 27 28 … … 61 62 #check that each element has only one flag 62 63 if any(SIAflag + SSAflag + L1L2flag + MLHOflag + HOflag + FSflag > 1): 63 print( "setflowequation warning message: some elements have several types, higher order type is used for them")64 print('Warning: setflowequation.py: some elements have several types, higher order type is used for them') 64 65 SIAflag[np.where(np.logical_and(SIAflag, SSAflag))] = False 65 66 SIAflag[np.where(np.logical_and(SIAflag, HOflag))] = False
Note:
See TracChangeset
for help on using the changeset viewer.