Index: /issm/trunk/src/m/classes/public/loadmultipleresultsfromcluster.m =================================================================== --- /issm/trunk/src/m/classes/public/loadmultipleresultsfromcluster.m (revision 2889) +++ /issm/trunk/src/m/classes/public/loadmultipleresultsfromcluster.m (revision 2889) @@ -0,0 +1,39 @@ +function md_list=loadmultipleresultsfromcluster(md_list) +%LOADMULTIPLERESULTSFROMCLUSTER - load multiple results of solution sequences from cluster +% +% Usage: +% md_list=loadresultsfromcluster(md_list); + +nummodels=length(md_list); + +%Get cluster settings +cluster=md_list{1}.cluster; +name=md_list{1}.name; +cluster_rc_location=which('cluster.rc'); +[codepath,executionpath]=ClusterParameters(cluster,cluster_rc_location); + +%Remote tar: +disp('tarring results'); +issmssh(cluster,['"cd ' executionpath ' && rm -rf file_list.txt ModelResults.tar.gz && find -iname ''*_*vs*.outbin'' > file_list.txt && tar zcvf ModelResults.tar.gz --files-from file_list.txt && rm -rf file_list.txt "']); + +%copy results from cluster to present directory +scpin(cluster, executionpath, {'ModelResults.tar.gz'}); + +%untar: +!tar -zxvf ModelResults.tar.gz + +%ok, go through list and load results from disk: +for i=1:nummodels, + %load results for this model + %md_list{i}.analysis_type=DiagnosticAnalysisEnum(); + %md_list{i}.sub_analysis_type=NoneAnalysisEnum(); + %md_list{i}=loadresultsfromdisk(md_list{i},[md_list{i}.name '_' num2str(i) 'vs' num2str(nummodels) '.outbin']); + md_list{i}=loadresultsfromdisk(md_list{i},[md_list{i}.name '.outbin']); + + %post solve phase + md_list{i}=postsolveparallel(md_list{i}); + delete([name '.outbin']); +end + +%erase files +delete('ModelResults.tar.gz'); Index: /issm/trunk/src/m/classes/public/modelsextract.m =================================================================== --- /issm/trunk/src/m/classes/public/modelsextract.m (revision 2889) +++ /issm/trunk/src/m/classes/public/modelsextract.m (revision 2889) @@ -0,0 +1,84 @@ +function md_list=modelsextract(md,flags,minel,varargin) +%modelsextract - extract several self contained models according to a list of element flags. +% +% The difference between this routine and the modelextract.m routine (without an 's') is that +% as many models are extracted as there are closed contours defined in area. +% This routine is needed for example when doing data assimilation of ice shelves in Antarctica. +% Many independent ice shelves are present, and we don't want data assimilation on one ice shelf +% to be hindered by another totally independent ice shelf. +% +% Usage: +% md_list=modelsextract(md,elementfalgs,minel); +% +% Examples: +% md_list=modelsextract(md,md.elementoniceshelf,1000); +% +% See also: EXTRUDE, COLLAPSE, MODELEXTRACT + +disp('selecting pools of elements'); +%go through flags and build as many independent element flags as there are groups of connected 1s +%in flags. + +%2D or 3D? +if strcmpi(md.type,'3d'), + numberofelements=md.numberofelements2d; %this will be forgotten when we get out. + flags=project2d(md,flags,1); +else + numberofelements=md.numberofelements; +end + +%recover extra arguments: +distance=0; +if nargin==4, + distance=varargin{1}; +end + +flag_list=cell(0,1); + +for i=1:size(flags,1), + + if (flags(i)), + + %ok, we are sure element i is part of a new pool. + pool=zeros(numberofelements,1); + pool=PropagateFlagsFromConnectivity(md.elementconnectivity,pool,i,flags); + flag_list{end+1,1}=pool; + + %speed up rest of computation by taking pool out of flags: + pos=find(pool);flags(pos)=0; + + end +end + +%go through flag_list and discard any pool of less than minel elements: +ex_pos=[]; +for i=1:length(flag_list), + if length(find(flag_list{i}))