source:
issm/oecreview/Archive/23390-24306/ISSM-24178-24179.diff@
24307
Last change on this file since 24307 was 24307, checked in by , 5 years ago | |
---|---|
File size: 1.8 KB |
-
../trunk-jpl/src/m/partition/AreaAverageOntoPartition.py
31 31 32 32 #finally, project vector: 33 33 vector = project2d(md3d,vector,layer) 34 md.qmu. partition = project2d(md3d,md3d.qmu.partition,layer)34 md.qmu.vpartition = project2d(md3d,md3d.qmu.vpartition,layer) 35 35 36 36 37 37 #ok, first check that part is Matlab indexed 38 part = (md.qmu. partition).copy()38 part = (md.qmu.vpartition).copy() 39 39 part = part.flatten() + 1 40 40 41 41 #some check: 42 42 if md.qmu.numberofpartitions != max(part): 43 raise RuntimeError('AreaAverageOntoPartition error message: ''npart'' should be equal to max(md.qmu. partition)')43 raise RuntimeError('AreaAverageOntoPartition error message: ''npart'' should be equal to max(md.qmu.vpartition)') 44 44 45 45 46 46 #initialize output -
../trunk-jpl/src/m/partition/AreaAverageOntoPartition.m
28 28 29 29 %finally, project vector: 30 30 vector=project2d(md3d,vector,layer); 31 md.qmu. partition=project2d(md3d,md3d.qmu.partition,layer);31 md.qmu.vpartition=project2d(md3d,md3d.qmu.vpartition,layer); 32 32 end 33 33 34 34 %ok, first check that part is Matlab indexed 35 part=md.qmu. partition+1;35 part=md.qmu.vpartition+1; 36 36 37 37 %some check: 38 38 if md.qmu.numberofpartitions~=max(part), 39 error('AreaAverageOntoPartition error message: ''npart'' should be equal to max(md.qmu. partition)');39 error('AreaAverageOntoPartition error message: ''npart'' should be equal to max(md.qmu.vpartition)'); 40 40 end 41 41 42 42 %initialize output
Note:
See TracBrowser
for help on using the repository browser.