![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)
Changeset 12662
- Timestamp:
- 07/19/12 15:32:03 (13 years ago)
- Location:
- issm/trunk-jpl/src/m/utils
- Files:
-
- 2 deleted
- 2 edited
-
DataProcessing/pkriging.m (modified) (2 diffs)
-
consistency/checkfield.m (modified) (15 diffs)
-
consistency/checkmessage.m (deleted)
-
consistency/modelconsistency.m (deleted)
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/utils/DataProcessing/pkriging.m
r12384 r12662 8 8 cluster=getfieldvalue(options,'cluster',generic('np',10)); 9 9 options=removefield(options,'cluster',0); 10 name = 'krig';10 name = ['krig' num2str(feature('GetPid'))]; 11 11 12 12 % ========================================= MARSHALL.m ================================================= … … 38 38 LaunchQueueJob(cluster,name,name,{[name '.bin'] [name '.queue']}); 39 39 toc 40 choice=input('Is the job successfully completed? (y/n)','s'); 40 41 Download(cluster,name,{[name '.outbin']}); 42 structure=parseresultsfromdisk([name '.outbin'],0); 43 delete([name '.outlog']); 44 delete([name '.errlog']); 45 delete([name '.outbin']); 46 delete([name '.bin']); 47 if ~ispc, 48 delete([name '.tar.gz']); 49 end 41 50 42 51 %Process results 43 structure=parseresultsfromdisk([name '.outbin'],0);44 52 B=structure.AutodiffForward; 45 B=reshape(B,size(x_interp, 1),size(x_interp,2))';53 B=reshape(B,size(x_interp,2),size(x_interp,1))'; 46 54 E=structure.AutodiffIsautodiff; 47 E=reshape(E,size(x_interp, 1),size(x_interp,2))';55 E=reshape(E,size(x_interp,2),size(x_interp,1))'; -
issm/trunk-jpl/src/m/utils/consistency/checkfield.m
r11264 r12662 1 function checkfield(md,fieldname,varargin)1 function md = checkfield(md,fieldname,varargin) 2 2 %CHECKFIELD - check field consistency 3 3 % … … 19 19 % 20 20 % Usage: 21 % checkfield(md,fieldname,options);21 % md = checkfield(md,fieldname,options); 22 22 % 23 23 % Example: 24 % checkfield(md,'mesh.elementonbed','size',[md.mesh.numberofelements 1],'values',[0 1]);25 % checkfield(md,'diagnostic.icefront','size',[NaN 4],'NaN',1);26 % checkfield(md,'diagnostic.icefront(:,end)','values',[0 1 2]);24 % md = checkfield(md,'mesh.elementonbed','size',[md.mesh.numberofelements 1],'values',[0 1]); 25 % md = checkfield(md,'diagnostic.icefront','size',[NaN 4],'NaN',1); 26 % md = checkfield(md,'diagnostic.icefront(:,end)','values',[0 1 2]); 27 27 28 28 %get options … … 35 35 if exist(options,'empty') 36 36 if isempty(field), 37 checkmessage(getfieldvalue(options,'message',...37 md = checkmessage(md,getfieldvalue(options,'message',... 38 38 ['field ''' fieldname ''' is empty'])); 39 39 end … … 45 45 if isnan(fieldsize(1)), 46 46 if (size(field,2)~=fieldsize(2)), 47 checkmessage(getfieldvalue(options,'message',...47 md = checkmessage(md,getfieldvalue(options,'message',... 48 48 ['field ''' fieldname ''' should have ' num2str(fieldsize(2)) ' columns'])); 49 49 end 50 50 elseif isnan(fieldsize(2)), 51 51 if (size(field,1)~=fieldsize(1)), 52 checkmessage(getfieldvalue(options,'message',...52 md = checkmessage(md,getfieldvalue(options,'message',... 53 53 ['field ''' fieldname ''' should have ' num2str(fieldsize(1)) ' lines'])); 54 54 end 55 55 else 56 56 if ((size(field)~=fieldsize(1)) | (size(field,2)~=fieldsize(2))) 57 checkmessage(getfieldvalue(options,'message',...57 md = checkmessage(md,getfieldvalue(options,'message',... 58 58 ['field ''' fieldname ''' size should be ' num2str(fieldsize(1)) ' x ' num2str(fieldsize(2))])); 59 59 end … … 66 66 if ~ismember(numel(field),fieldnumel), 67 67 if length(fieldnumel)==1 68 checkmessage(getfieldvalue(options,'message',...68 md = checkmessage(md,getfieldvalue(options,'message',... 69 69 ['field ''' fieldname ''' size should be ' sprintf('%g ',fieldnumel) ])); 70 70 elseif length(fieldnumel)==2 71 checkmessage(getfieldvalue(options,'message',...71 md = checkmessage(md,getfieldvalue(options,'message',... 72 72 ['field ''' fieldname ''' size should be ' num2str(fieldnumel(1)) ' or ' num2str(fieldnumel(2)) ])); 73 73 else 74 checkmessage(getfieldvalue(options,'message',...74 md = checkmessage(md,getfieldvalue(options,'message',... 75 75 ['field ''' fieldname ''' size should be ' sprintf('%g, ',fieldnumel(1:end-1)) ' or ' num2str(fieldnumel(end)) ])); 76 76 end … … 82 82 field2=reshape(field,prod(size(field)),1); 83 83 if any(isnan(field2)), 84 checkmessage(getfieldvalue(options,'message',...84 md = checkmessage(md,getfieldvalue(options,'message',... 85 85 ['NaN values found in field ''' fieldname ''''])); 86 86 end … … 90 90 if getfieldvalue(options,'cell',0); 91 91 if ~iscell(field), 92 checkmessage(getfieldvalue(options,'message',...92 md = checkmessage(md,getfieldvalue(options,'message',... 93 93 ['field ''' fieldname ''' should be a cell'])); 94 94 end … … 102 102 if any(~ismember(field,fieldvalues)), 103 103 if length(fieldvalues)==1 104 checkmessage(getfieldvalue(options,'message',...104 md = checkmessage(md,getfieldvalue(options,'message',... 105 105 ['field ''' fieldname ''' value should be ' fieldvalues{1} ])); 106 106 elseif length(fieldvalues)==2 107 checkmessage(getfieldvalue(options,'message',...107 md = checkmessage(md,getfieldvalue(options,'message',... 108 108 ['field ''' fieldname ''' values should be ' fieldvalues{1} ' or ' fieldvalues{2} ])); 109 109 else 110 checkmessage(getfieldvalue(options,'message',...110 md = checkmessage(md,getfieldvalue(options,'message',... 111 111 ['field ''' fieldname ''' should have values in ' sprintf('''%s'', ',fieldvalues{1:end-1}) 'or ''' fieldvalues{end} ''''])); 112 112 end 113 113 end 114 114 else 115 checkmessage(getfieldvalue(options,'message',...115 md = checkmessage(md,getfieldvalue(options,'message',... 116 116 ['field ''' fieldname ''' should be one of the following strings: ' sprintf('''%s'', ',fieldvalues{1:end-1}) 'or ''' fieldvalues{end} ''''])); 117 117 end … … 120 120 if isnumeric(field), 121 121 if any(~ismember(field2,fieldvalues)), 122 checkmessage(getfieldvalue(options,'message',...122 md = checkmessage(md,getfieldvalue(options,'message',... 123 123 ['field ''' fieldname ''' should have values in [' num2str(fieldvalues) ']'])); 124 124 end 125 125 else 126 checkmessage(getfieldvalue(options,'message',...126 md = checkmessage(md,getfieldvalue(options,'message',... 127 127 ['field ''' fieldname ''' should be a number in [' num2str(fieldvalues) ']'])); 128 128 end … … 135 135 field2=reshape(field,prod(size(field)),1); 136 136 if any(field2<lowerbound), 137 checkmessage(getfieldvalue(options,'message',...137 md = checkmessage(md,getfieldvalue(options,'message',... 138 138 ['field ''' fieldname ''' should have values above ' num2str(lowerbound)])); 139 139 end … … 143 143 field2=reshape(field,prod(size(field)),1); 144 144 if any(field2<=lowerbound), 145 checkmessage(getfieldvalue(options,'message',...145 md = checkmessage(md,getfieldvalue(options,'message',... 146 146 ['field ''' fieldname ''' should have values above ' num2str(lowerbound)])); 147 147 end … … 153 153 field2=reshape(field,prod(size(field)),1); 154 154 if any(field2>upperbound), 155 checkmessage(getfieldvalue(options,'message',...155 md = checkmessage(md,getfieldvalue(options,'message',... 156 156 ['field ''' fieldname ''' should have values below ' num2str(upperbound)])); 157 157 end … … 161 161 field2=reshape(field,prod(size(field)),1); 162 162 if any(field2>=upperbound), 163 checkmessage(getfieldvalue(options,'message',...163 md = checkmessage(md,getfieldvalue(options,'message',... 164 164 ['field ''' fieldname ''' should have values below ' num2str(upperbound(1))])); 165 165 end … … 169 169 if getfieldvalue(options,'file',0), 170 170 if ~exist(field,'file') 171 checkmessage(['file profided in ''' fieldname ''': ''' field ''' does not exist']);171 md = checkmessage(md,['file profided in ''' fieldname ''': ''' field ''' does not exist']); 172 172 end 173 173 end … … 177 177 if size(field,1)==md.mesh.numberofvertices, 178 178 if ~size(field,2)==1, 179 checkmessage(getfieldvalue(options,'message',...179 md = checkmessage(md,getfieldvalue(options,'message',... 180 180 ['field ''' fieldname ''' should have only one column as there are md.mesh.numberofvertices lines'])); 181 181 end 182 182 elseif size(field,1)==md.mesh.numberofvertices+1 183 183 if any(field(end,:)~=sort(field(end,:))), 184 checkmessage(getfieldvalue(options,'message',...184 md = checkmessage(md,getfieldvalue(options,'message',... 185 185 ['field ''' fieldname ''' columns should be chronological'])); 186 186 end 187 187 if any(field(end,1:end-1)==field(end,2:end)), 188 checkmessage(getfieldvalue(options,'message',...188 md = checkmessage(md,getfieldvalue(options,'message',... 189 189 ['field ''' fieldname ''' columns must not contain duplicate timesteps'])); 190 190 end 191 191 else 192 checkmessage(getfieldvalue(options,'message',...192 md = checkmessage(md,getfieldvalue(options,'message',... 193 193 ['field ''' fieldname ''' should have md.mesh.numberofvertices or md.mesh.numberofvertices+1 lines'])); 194 194 end
Note:
See TracChangeset
for help on using the changeset viewer.
