source: issm/trunk/src/m/classes/public/ismodelselfconsistent.m@ 5360

Last change on this file since 5360 was 5360, checked in by Eric.Larour, 15 years ago

reverted back

File size: 17.2 KB
RevLine 
[2326]1function ismodelselfconsistent(md),
[1]2%ISMODELSELFCONSISTENT - check that model forms a closed form solvable problem.
3%
4% Usage:
[2326]5% ismodelselfconsistent(md),
[1]6
[2366]7%tolerance we use in litmus checks for the consistency of the model
[1]8tolerance=10^-12;
[4756]9%check usage {{{1
[2326]10if nargin~=1,
[27]11 help ismodelselfconsistent
12 error('ismodelselfconsistent error message: wrong usage');
[1]13end
[4756]14%}}}
[1]15
[4756]16%recursive call for TRANSIENT {{{1
[4303]17if (md.analysis_type==Transient2DSolutionEnum | md.analysis_type==Transient3DSolutionEnum),
[2326]18 if md.dt<=0,
19 error('model not consistent: field dt must be positive for a transient run')
20 end
21
22 %recursive call to ismodelselfconsistent
[3994]23 if (md.dim==2),
[4303]24 analysis=[DiagnosticSolutionEnum PrognosticSolutionEnum];
[2528]25 else
[4303]26 analysis=[DiagnosticSolutionEnum PrognosticSolutionEnum ThermalSolutionEnum];
[2528]27 end
28
[2326]29 for i=1:length(analysis),
[2520]30 md.analysis_type=analysis(i);
[2326]31 ismodelselfconsistent(md);
32 end
33end
[4756]34%}}}
[2326]35
[4756]36% Common checks
37%COUNTER {{{1
[27]38if md.counter<3,
[2326]39 error(['model ' md.name ' is not correctly configured. You forgot one step in the following sequence (mesh, geography, parameterize,setelementstype)!']);
[1]40end
[4756]41%}}}
42%NAME{{{1
[1]43if isempty(md.name),
[2326]44 error(['model is not correctly configured: missing name!']);
[1]45end
[4756]46%}}}
47%ELEMENTS{{{1
48fields={'elements'};
49if (md.dim==2),
50 checksize(md,fields,[md.numberofelements 3]);
51else
52 checksize(md,fields,[md.numberofelements 6]);
53end
54%}}}
55%ELEMENTSTYPE{{{1
[4891]56%Check the size of elements_type
57fields={'elements_type'};
58checksize(md,fields,[md.numberofelements 2]);
59%Check the values of elements_type(1)
[5134]60checkvalues(md,{'elements_type(:,1)'},[MacAyealApproximationEnum() HutterApproximationEnum() PattynApproximationEnum() MacAyealPattynApproximationEnum()]);
[4891]61%Check the values of elements_type(2)
[5126]62checkvalues(md,{'elements_type(:,2)'},[StokesApproximationEnum() NoneApproximationEnum()]);
[3994]63if (md.dim==2),
[5126]64 checkvalues(md,{'elements_type(:,1)'},[MacAyealApproximationEnum() HutterApproximationEnum()]);
[1]65end
[358]66if (md.ismacayealpattyn==0 && md.ishutter==0 && md.isstokes==0),
[2326]67 error(['model not consistent: no elements type set for this model. at least one of ismacayealpattyn, ishutter and isstokes need to be =1']);
[358]68end
[4756]69%}}}
[5071]70%VERTICESTYPE{{{1
71%Check the size of verticess_type
72fields={'vertices_type'};
73checksize(md,fields,[md.numberofgrids 2]);
74%Check the values of elements_type(1)
[5126]75checkvalues(md,{'elements_type(:,1)'},[MacAyealApproximationEnum() HutterApproximationEnum() PattynApproximationEnum() MacAyealPattynApproximationEnum()]);
[5071]76%Check the values of elements_type(2)
[5126]77checkvalues(md,{'elements_type(:,2)'},[StokesApproximationEnum() NoneApproximationEnum()]);
[5071]78if (md.dim==2),
[5126]79 checkvalues(md,{'elements_type(:,1)'},[MacAyealApproximationEnum() HutterApproximationEnum()]);
[5071]80end
81if (md.ismacayealpattyn==0 && md.ishutter==0 && md.isstokes==0),
82 error(['model not consistent: no elements type set for this model. at least one of ismacayealpattyn, ishutter and isstokes need to be =1']);
83end
84%}}}
[4756]85%DG {{{1
[4764]86if md.prognostic_DG==1;
[4756]87 %CHECK THE COMPATIBILITY OF THE EDGES
88 fields={'edges'};
89 checksize(md,fields,[NaN 4]);
90 checknan(md,fields);
91end
92%}}}
93%PRESSURELOAD{{{1
[3994]94if (md.dim==2),
[2326]95 fields={'pressureload'};
[3099]96 checksize(md,fields,[NaN 4]);
[3994]97elseif md.dim==3,
[2326]98 fields={'pressureload'};
[3099]99 checksize(md,fields,[NaN 6]);
[3356]100else
[4796]101 error('dim should be either 2 3');
[1796]102end
[4870]103checkvalues(md,{'pressureload(:,end)'},[WaterEnum() AirEnum()]);
[4756]104%}}}
105%NO NAN {{{1
[3758]106fields={'numberofelements','numberofgrids','x','y','z','drag_coefficient','drag_type','drag_p','drag_q',...
107 'rho_ice','rho_water','rheology_B','elementoniceshelf','surface','thickness','bed','g','lowmem','sparsity','nsteps','maxiter',...
108 'tolx','np','eps_res','max_nonlinear_iterations','exclusive','rheology_n','gridonbed','gridonsurface','elementonbed','elementonsurface','deltaH','DeltaH','timeacc','timedec'};
[2366]109checknan(md,fields);
[4756]110%}}}}
111%FIELDS >= 0 {{{1
[3758]112fields={'numberofelements','numberofgrids','elements','drag_coefficient','drag_type','drag_p','drag_q',...
113 'rho_ice','rho_water','rheology_B','elementoniceshelf','thickness','g','eps_res','max_nonlinear_iterations','eps_rel','eps_abs','nsteps','maxiter','tolx','exclusive',...
114 'sparsity','lowmem','rheology_n','gridonbed','gridonsurface','elementonbed','elementonsurface','deltaH','DeltaH','timeacc','timedec'};
[2366]115checkgreater(md,fields,0);
[4756]116%}}}
117%FIELDS > 0 {{{1
[3758]118fields={'numberofelements','numberofgrids','elements','drag_type','drag_p',...
119 'rho_ice','rho_water','rheology_B','thickness','g','max_nonlinear_iterations','eps_res','eps_rel','eps_abs','maxiter','tolx',...
[27]120 'sparsity','deltaH','DeltaH','timeacc','timedec'};
[2366]121checkgreaterstrict(md,fields,0);
[4756]122%}}}
123%SIZE NUMBEROFELEMENTS {{{1
[3759]124fields={'drag_p','drag_q','elementoniceshelf','rheology_n','elementonbed'};
[2991]125checksize(md,fields,[md.numberofelements 1]);
[4756]126%}}}
127%SIZE NUMBEROFGRIDS {{{1
[3758]128fields={'x','y','z','rheology_B','drag_coefficient','melting_rate','accumulation_rate','surface','thickness','bed','gridonbed','gridonsurface'};
[2991]129checksize(md,fields,[md.numberofgrids 1]);
[4756]130%}}}
131%OTHER SIZES {{{1
[2405]132fields={'spcvelocity'};
133checksize(md,fields,[md.numberofgrids 6]);
[4756]134%}}}
135%THICKNESS = SURFACE - BED {{{1
[681]136if any((md.thickness-md.surface+md.bed)>tolerance),
[2326]137 error(['model not consistent: model ' md.name ' violates the equality thickness=surface-bed!']);
[27]138end
[4756]139%}}}
140%RIFTS{{{1
[27]141if md.numrifts,
[3994]142 if ~(md.dim==2),
[2326]143 error(['model not consistent: models with rifts are only supported in 2d for now!']);
[1]144 end
[2660]145 if ~isstruct(md.rifts),
146 error(['model not consistent: md.rifts should be a structure!']);
147 end
148 if ~isempty(find(md.segmentmarkers>=2)),
149 %We have segments with rift markers, but no rift structure!
150 error(['model not consistent: model ' md.name ' should be processed for rifts (run meshprocessrifts)!']);
151 end
152 %Check that rifts are filled with proper material
153 checkvalues(md,{'rifts.fill'},[WaterEnum() AirEnum() IceEnum() MelangeEnum()]);
154else
[2748]155 if ~isnans(md.rifts),
[2660]156 error(['model not consistent: md.rifts shoud be NaN since md.numrifts is 0!']);
[27]157 end
158end
[4756]159%}}}
160%FLAGS (0 or 1){{{1
[2326]161if ~ismember(md.artificial_diffusivity,[0 1]),
162 error('model not consistent: artificial_diffusivity should be a scalar (1 or 0)');
[1]163end
[4756]164if ~ismember(md.prognostic_DG,[0 1]),
165 error('model not consistent: prognostic_DG should be a scalar (1 or 0)');
166end
167if ~ismember(md.lowmem,[0 1]),
168 error(['model not consistent: model ' md.name ' lowmem field should be 0 or 1']);
169end
170%}}}
171%PARAMETEROUTPUT {{{1
[2307]172if md.numoutput~=length(md.parameteroutput),
[2326]173 error('model not consistent: numoutput should be the same size as parameteroutput');
[2307]174end
[4756]175%}}}
176%CONNECTIVITY {{{1
[3994]177if (md.dim==2),
[2326]178 if md.connectivity<9,
179 error('model not consistent: connectivity should be at least 9 for 2d models');
180 end
181end
[3994]182if md.dim==3,
[2326]183 if md.connectivity<24,
184 error('model not consistent: connectivity should be at least 24 for 3d models');
185 end
186end
[4756]187%}}}
188%PARALLEL{{{1
[2326]189if ~strcmpi(md.cluster,'none'),
190
191 %NAN VALUES
192 fields={'time','np'};
[2366]193 checknan(md,fields);
[2326]194
195 %FIELD > 0
196 fields={'time','np'};
[2366]197 checkgreaterstrict(md,fields,0);
[2326]198
199end
[4756]200%}}}
[2326]201
[4756]202% Solution checks
203%DIAGNOSTIC{{{1
[4303]204if md.analysis_type==DiagnosticSolutionEnum,
[27]205
206 %HUTTER ON ICESHELF WARNING
[5126]207 if any(md.elements_type(:,1)==HutterApproximationEnum & md.elementoniceshelf),
[27]208 disp(sprintf('\n !!! Warning: Hutter''s model is not consistent on ice shelves !!!\n'));
209 end
210
211 %SINGULAR
[2282]212 if ~any(sum(md.spcvelocity(:,1:2),2)==2),
[2326]213 error(['model not consistent: model ' md.name ' is not well posed (singular). You need at least one grid with fixed velocity!'])
[1]214 end
215
[27]216 %DIRICHLET IF THICKNESS <= 0
[681]217 if any(md.thickness<=0),
[27]218 pos=find(md.thickness<=0);
[1755]219 if any(find(md.spcthickness(pos,1)==0)),
[2326]220 error(['model not consistent: model ' md.name ' has some grids with 0 thickness']);
[27]221 end
222 end
[2357]223
[2367]224 %INITIAL VELOCITY
225 if length(md.vx)==md.numberofgrids & length(md.vy)==md.numberofgrids,
226 fields={'vx','vy'};
227 checknan(md,fields);
228 end
[1]229end
[4756]230%}}}
231%PROGNOSTIC{{{1
[4303]232if md.analysis_type==PrognosticSolutionEnum,
[1]233
[2357]234 %INITIAL VELOCITIES
235 fields={'vx','vy'};
[2366]236 checksize(md,fields,[md.numberofgrids 1]);
237 checknan(md,fields);
[2357]238
[2519]239 %CHECK THAT WE ARE NOT FULLY CONSTRAINED
[3994]240 if (md.dim==2),
[2519]241 if isempty(find(~md.spcthickness(:,1))),
242 error(['model not consistent: model ' md.name ' is totally constrained for prognostic, no need to solve!']);
243 end
244 end
245
[1]246end
[4756]247%}}}
248%STEADYSTATE{{{1
[4303]249if md.analysis_type==SteadystateSolutionEnum,
[2326]250
[1901]251 %NDT
252 if md.dt~=0,
[2326]253 error(['model not consistent: for a steadystate computation, dt must be zero.']);
[1901]254 end
[2326]255
[1787]256 %PRESSURE
257 if isnans(md.pressure),
[2326]258 error(['model not consistent: for a steadystate computation, the model must have an initial pressure, even lithostatic will do.']);
[1787]259 end
260
261 %eps:
262 if isnan(md.eps_rel),
[2326]263 error(['model not consistent: for a steadystate computation, eps_rel (relative convergence criterion) must be defined!']);
[1787]264 end
265
266 %dim:
[3994]267 if (md.dim==2),
[2326]268 error(['model not consistent: for a steadystate computation, model needs to be 3d']);
[1787]269 end
270end
[4756]271%}}}
272%THERMAL {{{1
[27]273%THERMAL STEADY AND THERMAL TRANSIENT
[4303]274if md.analysis_type==ThermalSolutionEnum,
[27]275
276 %EXTRUSION
[3994]277 if (md.dim==2),
[2326]278 error(['model not consistent: for a ' md.analysis_type ' computation, the model must be 3d, extrude it first!'])
[1]279 end
[27]280
[2519]281 %CHECK THAT WE ARE NOT FULLY CONSTRAINED
[3994]282 if isempty(find(~md.spctemperature(:,1))),
283 error(['model not consistent: model ' md.name ' is totally constrained for temperature, no need to solve!']);
[2519]284 end
285
[27]286 %VELOCITIES AND PRESSURE
[2357]287 fields={'vx','vy','vz','pressure'};
[2366]288 checksize(md,fields,[md.numberofgrids 1]);
289 checknan(md,fields);
[2357]290
[1]291end
292
[1311]293%THERMAL TRANSIENT
[4303]294if md.analysis_type==ThermalSolutionEnum & md.dt~=0,
[27]295
296 %DT and NDT
297 fields={'dt','ndt'};
[2366]298 checkgreaterstrict(md,fields,0);
[27]299
[1311]300 %INITIAL TEMPERATURE, MELTING AND ACCUMULATION
[3758]301 fields={'temperature','accumulation_rate','melting_rate'};
[2366]302 checksize(md,fields,[md.numberofgrids 1]);
[2440]303 checknan(md,fields);
[1]304
[2714]305 %INITIAL TEMPERATURE
[2669]306 fields={'temperature','spctemperature(:,2)','observed_temperature'};
307 checkgreater(md,fields,0)
308
[1]309end
[4756]310%}}}
311%BALANCEDTHICKNESS{{{1
[4303]312if md.analysis_type==BalancedthicknessSolutionEnum
[2714]313
314 %VELOCITIES MELTING AND ACCUMULATION
[4756]315 fields={'vx','vy','accumulation_rate','melting_rate','dhdt'};
[3589]316 checksize(md,fields,[md.numberofgrids 1]);
317 checknan(md,fields);
318
319 %SPC
[4965]320 if ~md.prognostic_DG,
321 if any(md.spcthickness(find(md.gridonboundary))~=1),
322 error(['model not consistent: model ' md.name ' should have all the nodes on boundary constrained in field spcthickness']);
323 end
324 end
[3589]325end
[4756]326%}}}
327%BALANCEDVELOCITIES{{{1
[4303]328if md.analysis_type==BalancedvelocitiesSolutionEnum
[2723]329
330 %VELOCITIES MELTING AND ACCUMULATION
[3758]331 fields={'vx','vy','accumulation_rate','melting_rate'};
[2723]332 checksize(md,fields,[md.numberofgrids 1]);
333 checknan(md,fields);
334
335 %SPC
336 if any(md.spcvelocity(find(md.gridonboundary),[1:2])~=1),
337 error(['model not consistent: model ' md.name ' should have all the nodes on boundary constrained in field spcvelocity']);
338 end
339end
[4756]340%}}}
341%CONTROL{{{1
[1907]342if md.control_analysis,
[1]343
[27]344 %CONTROL TYPE
[5359]345 if ~isnumeric(md.control_type),
346 error('model not consistent: control_type should be an enum');
[1]347 end
[5359]348 checkvalues(md,{'control_type'},[DhDtEnum DragCoefficientEnum RheologyBbarEnum]);
[1]349
[27]350 %LENGTH CONTROL FIELDS
[5360]351 fields={'maxiter','optscal','cm_responses','cm_jump'};
[2991]352 checksize(md,fields,[md.nsteps 1]);
[27]353
[5282]354 %RESPONSES
[5289]355 checkvalues(md,{'cm_responses'},[SurfaceAbsVelMisfitEnum SurfaceRelVelMisfitEnum SurfaceLogVelMisfitEnum SurfaceLogVxVyMisfitEnum SurfaceAverageVelMisfitEnum ThicknessAbsMisfitEnum]);
[1]356
[3203]357 %WEIGHTS
358 fields={'weights'};
359 checksize(md,fields,[md.numberofgrids 1]);
360 checkgreater(md,fields,0);
361
[27]362 %OBSERVED VELOCITIES
[5324]363 if md.analysis_type==BalancedthicknessSolutionEnum
364 fields={'thickness_obs'};
365 checksize(md,fields,[md.numberofgrids 1]);
366 checknan(md,fields);
367 else
368 fields={'vx_obs','vy_obs'};
369 checksize(md,fields,[md.numberofgrids 1]);
370 checknan(md,fields);
371 end
[1361]372
373 %DIRICHLET IF THICKNESS <= 0
374 if any(md.thickness<=0),
375 pos=find(md.thickness<=0);
[1755]376 if any(find(md.spcthickness(pos,1)==0)),
[2326]377 error(['model not consistent: model ' md.name ' has some grids with 0 thickness']);
[1361]378 end
379 end
[2686]380
381 %parameters
382 fields={'cm_noisedmp'};
383 checknan(md,fields);
[1]384end
[4756]385%}}}
386%QMU {{{1
387if md.qmu_analysis,
388 if md.qmu_params.evaluation_concurrency~=1,
389 error(['model not consistent: concurrency should be set to 1 when running dakota in library mode']);
390 end
391 if ~isempty(md.part),
392 if numel(md.part)~=md.numberofgrids,
393 error(['model not consistent: user supplied partition for qmu analysis should have size md.numberofgrids x 1 ']);
394 end
395 if find(md.part)>=md.numberofgrids,
396 error(['model not consistent: user supplied partition should be indexed from 0 (c-convention)']);
397 end
398 if min(md.part)~=0,
399 error(['model not consistent: partition vector not indexed from 0 on']);
400 end
401 if max(md.part)>=md.numberofgrids,
402 error(['model not consistent: partition vector cannot have maximum index larger than number of grids']);
403 end
404 if ~isempty(find(md.part<0)),
405 error(['model not consistent: partition vector cannot have values less than 0']);
406 end
407 if ~isempty(find(md.part>=md.npart)),
408 error(['model not consistent: partition vector cannot have values more than md.npart-1']);
409 end
410 if max(md.part)>=md.npart,
411 error(['model not consistent: for qmu analysis, partitioning vector cannot go over npart, number of partition areas']);
412 end
413 end
414 if md.eps_rel>1.1*10^-5,
415 error(['model not consistent: for qmu analysis, eps_rel should be least than 10^-5, 10^-15 being a better value']);
416 end
417end
[1]418
[465]419if strcmpi(md.analysis_type,'qmu'),
[27]420 if ~strcmpi(md.cluster,'none'),
421 if md.waitonlock==0,
[2326]422 error(['model is not correctly configured: waitonlock should be activated when running qmu in parallel mode!']);
[1]423 end
424 end
425end
[4756]426%}}}
[1]427
[4756]428end
[1]429
[4756]430%checks additional functions
431%checklength {{{1
[2366]432function checklength(md,fields,fieldlength)
433 %CHECKSIZE - check length of a field
[2326]434 for i=1:length(fields),
435 if length(eval(['md.' fields{i}]))~=fieldlength,
436 error(['model not consistent: field ' fields{i} ' length should be ' num2str(fieldlength)]);
[2164]437 end
438 end
439end
[4756]440%}}}
441%checksize {{{1
[2366]442function checksize(md,fields,fieldsize)
443 %CHECKSIZE - check size of a field
[2326]444 for i=1:length(fields),
445 if isnan(fieldsize(1)),
446 if (size(eval(['md.' fields{i}]),2)~=fieldsize(2)),
[3203]447 %LOG warnings:
[3104]448 if strcmpi(fields{i},'pressureload'),
449 disp(' ');
450 disp(' If pressureload does not have the right size (one column missing), this is due to a change in md.pressure.');
451 disp(' To update your model, use the following line');
452 disp(' md.pressureload=[md.pressureload WaterEnum*md.elementoniceshelf(md.pressureload(:,end))+AirEnum*md.elementonicesheet(md.pressureload(:,end))];');
453 disp(' ');
454 end
[2326]455 error(['model not consistent: field ' fields{i} ' should have ' num2str(fieldsize(2)) ' columns']);
456 end
457 elseif isnan(fieldsize(2)),
458 if (size(eval(['md.' fields{i}]),1)~=fieldsize(1)),
459 error(['model not consistent: field ' fields{i} ' should have ' num2str(fieldsize(1)) ' rows']);
460 end
461 else
462 if ((size(eval(['md.' fields{i}]),1)~=fieldsize(1)) | (size(eval(['md.' fields{i}]),2)~=fieldsize(2)))
[3219]463 %LOG warnings:
464 if strcmpi(fields{i},'weights'),
465 disp(' ');
466 disp(' ''weights'' is a new field that should be of length numberofgrids.');
467 disp(' To update your model, use the following line');
468 disp(' md.weights=ones(md.numberofgrids,1)');
469 disp(' ');
470 end
[2326]471 error(['model not consistent: field ' fields{i} ' size should be ' num2str(fieldsize(1)) ' x ' num2str(fieldsize(2))]);
472 end
[2164]473 end
474 end
475end
[4756]476%}}}
477%checknan {{{1
[2366]478function checknan(md,fields)
479 %CHECKNAN - check nan values of a field
[2326]480 for i=1:length(fields),
481 if any(isnan(eval(['md.' fields{i}]))),
482 error(['model not consistent: NaN values in field ' fields{i}]);
[2164]483 end
484 end
485end
[4756]486%}}}
487%checkreal{{{1
[2366]488function checkreal(md,fields)
489 %CHECKREAL - check real values of a field
[2326]490 for i=1:length(fields),
491 if any(eval(['~isreal(md.' fields{i} ')'])),
492 error(['model not consistent: complex values in field ' fields{i}]);
493 end
494 end
[2164]495end
[4756]496%}}}
497%checkgreaterstrict{{{1
[2366]498function checkgreaterstrict(md,fields,lowerbound)
499 %CHECKGREATERSTRICT - check values of a field
[2326]500 for i=1:length(fields),
501 if any(eval(['md.' fields{i} '<=' num2str(lowerbound) ])),
502 error(['model not consistent: field ' fields{i} ' should have values stricly above ' num2str(lowerbound)]);
503 end
[2164]504 end
505end
[4756]506%}}}
507%checkgreater{{{1
[2366]508function checkgreater(md,fields,lowerbound)
509 %CHECKGREATER - check values of a field
[27]510 for i=1:length(fields),
[2326]511 if any(eval(['md.' fields{i} '<' num2str(lowerbound) ])),
512 error(['model not consistent: field ' fields{i} ' should have values above ' num2str(lowerbound)]);
[27]513 end
514 end
[2326]515end
[4756]516%}}}
517%checklessstrict{{{1
[2366]518function checklessstrict(md,fields,upperbound)
519 %CHECKLESSSTRICT - check values of a field
[27]520 for i=1:length(fields),
[2326]521 if any(eval(['md.' fields{i} '>=' num2str(upperbound) ])),
522 error(['model not consistent: field ' fields{i} ' should have values stricly below ' num2str(upperbound)]);
[27]523 end
524 end
[2326]525end
[4756]526%}}}
527%checkless{{{1
[2366]528function checkless(md,fields,upperbound)
529 %CHECKLESS - check values of a field
[27]530 for i=1:length(fields),
[2326]531 if any(eval(['md.' fields{i} '>' num2str(upperbound) ])),
532 error(['model not consistent: field ' fields{i} ' should have values below ' num2str(upperbound)]);
[27]533 end
[1]534 end
535end
[4756]536%}}}
537%checkvalues {{{1
[2660]538function checkvalues(md,fields,values)
[4870]539 %CHECKVALUE - check that a field has specified values
[2660]540 for i=1:length(fields),
[4870]541 if eval(['any(~ismember( md.' fields{i} ',values))']),
[2660]542 error(['model not consistent: field ' fields{i} ' should have values in ' num2str(values)]);
543 end
544 end
545end
[4756]546%}}}
Note: See TracBrowser for help on using the repository browser.