Changeset 9421 for issm/trunk

Timestamp:

08/25/11 14:41:40 (14 years ago)

Author:

jschierm

Message:

Flaim: fixed problem with directory specification in md.fm_tracks and improved error handling.

Location:

issm/trunk/src/m/solutions

Files:

• flaim.m (modified) (3 diffs)
• issm.m (modified) (1 diff)

issm/trunk/src/m/solutions/flaim.m

@@ -32 +32 @@
 %  process file names
 
+[pathstr,name,ext,versn] = fileparts(md.fm_tracks);
+fm_tracks=fullfile('',[name ext versn]);
 [pathstr,name,ext,versn] = fileparts(options.fmfile);
 if isempty(ext)
@@ -81 +83 @@
 display('Calling flaim.');
 FLAIM_DIR=flaimdir();
-%system([FLAIM_DIR '/clasp/swathModule/models/issm/issm.exe -kml ' md.fm_tracks ' -kml ' filekml ' -gridsatequator 200000 -opt 1 -hst 2017-276T02:57:00 -hdur 4d -sparsepoints -usevalueordering -split-antimeridian >& ' filelog]);
-flaim_cmd=[FLAIM_DIR '/clasp/swathModule/models/issm/issm.exe -kml ' md.fm_tracks ' -kml ' filekml ' -gridsatequator ' int2str(md.fm_gridsatequator) ' -opt 1 -hst 2017-276T02:57:00 -hdur 4d -sparsepoints'];
+%system([FLAIM_DIR '/clasp/swathModule/models/issm/issm.exe -kml ' fm_tracks ' -kml ' filekml ' -gridsatequator 200000 -opt 1 -hst 2017-276T02:57:00 -hdur 4d -sparsepoints -usevalueordering -split-antimeridian >& ' filelog]);
+flaim_cmd=[FLAIM_DIR '/clasp/swathModule/models/issm/issm.exe -kml ' fm_tracks ' -kml ' filekml ' -gridsatequator ' int2str(md.fm_gridsatequator) ' -opt 1 -hst 2017-276T02:57:00 -hdur 4d -sparsepoints'];
 if (md.fm_usevalueordering)
         flaim_cmd=[flaim_cmd ' -usevalueordering'];
@@ -100 +102 @@
 if exist(filelog,'file')
         [status,fline]=system(['grep -i quality ' filelog]);
-        strings=textscan(fline,'%s','delimiter',' =,','MultipleDelimsAsOne',1);
-        md.results.FlaimSolution.quality =sscanf(strings{1}{2},'%f');
+        if ~isempty(fline)
+                strings=textscan(fline,'%s','delimiter',' =,','MultipleDelimsAsOne',1);
+                md.results.FlaimSolution.quality =sscanf(strings{1}{2},'%f');
+        else
+                warning('Could not find ''quality'' in the ''%s/%s'' flaim log file.',options.fmdir,filelog);
+                md.results.FlaimSolution.quality =NaN;
+        end
+end
+if ~exist(md.results.FlaimSolution.solution,'file')
+        md.results.FlaimSolution.solution='';
 end
 display(md.results.FlaimSolution);

issm/trunk/src/m/solutions/issm.m

@@ -35 +35 @@
         else
                 %launch dakota driver for diagnostic core solution
+                display('calling Qmu');
                 Qmu(femmodel,femmodel.parameters);
+                display('finished with Qmu');
         end
+
+        dtest=Response(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,'MinVx',true);
+        disp(sprintf('MinVx=%g\n',dtest));
+        dtest=Response(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,'MaxVx',true);
+        disp(sprintf('MaxVx=%g\n',dtest));
+        dtest=Response(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,'MaxAbsVx',true);
+        disp(sprintf('MaxAbsVx=%g\n',dtest));
+        dtest=Response(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,'MinVy',true);
+        disp(sprintf('MinVy=%g\n',dtest));
+        dtest=Response(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,'MaxVy',true);
+        disp(sprintf('MaxVy=%g\n',dtest));
+        dtest=Response(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,'MaxAbsVy',true);
+        disp(sprintf('MaxAbsVy=%g\n',dtest));
+%       dtest=Response(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,'MinVel',true);
+%       disp(sprintf('MinVel=%g\n',dtest));
+%       dtest=Response(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,'MaxVel',true);
+%       disp(sprintf('MaxVel=%g\n',dtest));
 
         %stop timing
         t2=clock;
         disp(['solution converged in ' num2str(etime(t2,t1)) ' seconds']);
+
+        display('writing femmodel to femmodel.mat');
+        save femmodel femmodel
+