Changeset 24313 for issm/trunk/src/m/inversions/marshallcostfunctions.py
- Timestamp:
- 11/01/19 12:01:57 (5 years ago)
- Location:
- issm/trunk
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk
-
issm/trunk/src
- Property svn:mergeinfo changed
-
issm/trunk/src/m/inversions/marshallcostfunctions.py
r21341 r24313 1 import copy2 1 3 2 def marshallcostfunctions(cost_functions): 4 3 5 cfDict={101:'SurfaceAbsVelMisfit',6 102:'SurfaceRelVelMisfit',7 103:'SurfaceLogVelMisfit',8 104:'SurfaceLogVxVyMisfit',9 105:'SurfaceAverageVelMisfit',10 201:'ThicknessAbsMisfit',11 501:'DragCoefficientAbsGradient',12 502:'RheologyBbarAbsGradient',13 503:'ThicknessAbsGradient',14 504:'ThicknessAlongGradient',15 505:'ThicknessAcrossGradient',}4 cfDict = {101: 'SurfaceAbsVelMisfit', 5 102: 'SurfaceRelVelMisfit', 6 103: 'SurfaceLogVelMisfit', 7 104: 'SurfaceLogVxVyMisfit', 8 105: 'SurfaceAverageVelMisfit', 9 201: 'ThicknessAbsMisfit', 10 501: 'DragCoefficientAbsGradient', 11 502: 'RheologyBbarAbsGradient', 12 503: 'ThicknessAbsGradient', 13 504: 'ThicknessAlongGradient', 14 505: 'ThicknessAcrossGradient'} 16 15 17 data=[cfDict[cf] for cf in cost_functions]18 ##copy list first19 # data=copy.deepcopy(cost_functions)16 data = [cfDict[cf] for cf in cost_functions] 17 # #copy list first 18 # data = copy.deepcopy(cost_functions) 20 19 21 # #convert to strings 22 # pos=[i for i,x in enumerate(cost_functions) if x==101]; 23 # for i in pos: data[i]='SurfaceAbsVelMisfit' 24 # pos=[i for i,x in enumerate(cost_functions) if x==102]; 25 # for i in pos: data[i]='SurfaceRelVelMisfit' 26 # pos=[i for i,x in enumerate(cost_functions) if x==103]; 27 # for i in pos: data[i]='SurfaceLogVelMisfit' 28 # pos=[i for i,x in enumerate(cost_functions) if x==104]; 29 # for i in pos: data[i]='SurfaceLogVxVyMisfit' 30 # pos=[i for i,x in enumerate(cost_functions) if x==105]; 31 # for i in pos: data[i]='SurfaceAverageVelMisfit' 32 # pos=[i for i,x in enumerate(cost_functions) if x==201]; 33 # for i in pos: data[i]='ThicknessAbsMisfit' 34 # pos=[i for i,x in enumerate(cost_functions) if x==501]; 35 # for i in pos: data[i]='DragCoefficientAbsGradient' 36 # pos=[i for i,x in enumerate(cost_functions) if x==502]; 37 # for i in pos: data[i]='RheologyBbarAbsGradient' 38 # pos=[i for i,x in enumerate(cost_functions) if x==503]; 39 # for i in pos: data[i]='ThicknessAbsGradient' 40 # pos=[i for i,x in enumerate(cost_functions) if x==504]; 41 # for i in pos: data[i]='ThicknessAlongGradient' 42 # pos=[i for i,x in enumerate(cost_functions) if x==505]; 43 # for i in pos: data[i]='ThicknessAcrossGradient' 20 # #convert to strings 21 # pos = [i for i, x in enumerate(cost_functions) if x == 101] 22 # for i in pos: data[i] = 'SurfaceAbsVelMisfit' 23 # pos = [i for i, x in enumerate(cost_functions) if x == 102] 24 # for i in pos: data[i] = 'SurfaceRelVelMisfit' 25 # pos = [i for i, x in enumerate(cost_functions) if x == 103] 26 # for i in pos: data[i] = 'SurfaceLogVelMisfit' 27 # pos = [i for i, x in enumerate(cost_functions) if x == 104] 28 # for i in pos: data[i] = 'SurfaceLogVxVyMisfit' 29 # pos = [i for i, x in enumerate(cost_functions) if x == 105] 30 # for i in pos: data[i] = 'SurfaceAverageVelMisfit' 31 # pos = [i for i, x in enumerate(cost_functions) if x == 201] 32 # for i in pos: data[i] = 'ThicknessAbsMisfit' 33 # pos = [i for i, x in enumerate(cost_functions) if x == 501] 34 # for i in pos: data[i] = 'DragCoefficientAbsGradient' 35 # pos = [i for i, x in enumerate(cost_functions) if x == 502] 36 # for i in pos: data[i] = 'RheologyBbarAbsGradient' 37 # pos = [i for i, x in enumerate(cost_functions) if x == 503] 38 # for i in pos: data[i] = 'ThicknessAbsGradient' 39 # pos = [i for i, x in enumerate(cost_functions) if x == 504] 40 # for i in pos: data[i] = 'ThicknessAlongGradient' 41 # pos = [i for i, x in enumerate(cost_functions) if x == 505] 42 # for i in pos: data[i] = 'ThicknessAcrossGradient' 44 43 45 44 return data
Note:
See TracChangeset
for help on using the changeset viewer.