-
Notifications
You must be signed in to change notification settings - Fork 19
/
EurostagEchExport.java
445 lines (383 loc) · 23 KB
/
EurostagEchExport.java
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
/**
* Copyright (c) 2016, All partners of the iTesla project (http://www.itesla-project.eu/consortium)
* This Source Code Form is subject to the terms of the Mozilla Public
* License, v. 2.0. If a copy of the MPL was not distributed with this
* file, You can obtain one at http://mozilla.org/MPL/2.0/.
*/
package eu.itesla_project.iidm.eurostag.export;
import eu.itesla_project.commons.ITeslaException;
import eu.itesla_project.eurostag.network.*;
import eu.itesla_project.eurostag.network.io.EsgWriter;
import eu.itesla_project.iidm.network.*;
import eu.itesla_project.iidm.network.util.Identifiables;
import org.slf4j.Logger;
import org.slf4j.LoggerFactory;
import java.io.IOException;
import java.io.Writer;
import java.nio.charset.StandardCharsets;
import java.nio.file.Files;
import java.nio.file.Path;
import java.util.ArrayList;
import java.util.List;
import java.util.Objects;
import java.util.Properties;
/**
* @author Geoffroy Jamgotchian <geoffroy.jamgotchian at rte-france.com>
*/
public class EurostagEchExport {
private static final Logger LOGGER = LoggerFactory.getLogger(EurostagEchExport.class);
/**
* epsilon value for conductance
*/
public static final float G_EPSILON = 0.00001f;
/**
* epsilon value for susceptance
*/
public static final float B_EPSILON = 0.00001f;
private static final String XNODE_V_PROPERTY = "xnode_v";
private static final String XNODE_ANGLE_PROPERTY = "xnode_angle";
private final Network network;
private final EurostagEchExportConfig config;
private final BranchParallelIndexes parallelIndexes;
private final EurostagDictionary dictionary;
public EurostagEchExport(Network network, EurostagEchExportConfig config, BranchParallelIndexes parallelIndexes, EurostagDictionary dictionary) {
this.network = Objects.requireNonNull(network);
this.config = Objects.requireNonNull(config);
this.parallelIndexes = Objects.requireNonNull(parallelIndexes);
this.dictionary = Objects.requireNonNull(dictionary);
}
public EurostagEchExport(Network network, EurostagEchExportConfig config) {
this.network = Objects.requireNonNull(network);
this.config = config;
this.parallelIndexes = BranchParallelIndexes.build(network, config);
this.dictionary = EurostagDictionary.create(network, parallelIndexes, config);
}
public EurostagEchExport(Network network) {
this(network, new EurostagEchExportConfig());
}
private void createAreas(EsgNetwork esgNetwork) {
esgNetwork.addArea(new EsgArea(new Esg2charName(EchUtil.FAKE_AREA), EsgArea.Type.AC));
for (Country c : network.getCountries() ) {
esgNetwork.addArea(new EsgArea(new Esg2charName(c.toString()), EsgArea.Type.AC));
}
}
private EsgNode createNode(String busId, String countryIsoCode, float nominalV, float v, float angle, boolean slackBus) {
return new EsgNode(new Esg2charName(countryIsoCode),
new Esg8charName(dictionary.getEsgId(busId)),
nominalV,
Float.isNaN(v) ? 1f : v / nominalV,
Float.isNaN(angle) ? 0f : angle,
slackBus);
}
private EsgNode createNode(String busId, VoltageLevel vl, float v, float angle, boolean slackBus) {
return createNode(busId, vl.getSubstation().getCountry().name(), vl.getNominalV(), v, angle, slackBus);
}
private void createNodes(EsgNetwork esgNetwork) {
esgNetwork.addNode(createNode(EchUtil.FAKE_NODE_NAME1, EchUtil.FAKE_AREA, 380f, 380f, 0f, false));
esgNetwork.addNode(createNode(EchUtil.FAKE_NODE_NAME2, EchUtil.FAKE_AREA, 380f, 380f, 0f, false));
Bus sb = EchUtil.selectSlackbus(network, config);
if (sb == null) {
throw new RuntimeException("Stack bus not found");
}
LOGGER.debug("Slack bus: {} ({})", sb, sb.getVoltageLevel().getId());
for (Bus b : Identifiables.sort(EchUtil.getBuses(network, config))) {
esgNetwork.addNode(createNode(b.getId(), b.getVoltageLevel(), b.getV(), b.getAngle(), sb == b));
}
for (DanglingLine dl : Identifiables.sort(network.getDanglingLines())) {
Properties properties = dl.getProperties();
String strV = properties.getProperty(XNODE_V_PROPERTY);
String strAngle = properties.getProperty(XNODE_ANGLE_PROPERTY);
float v = strV != null ? Float.parseFloat(strV) : Float.NaN;
float angle = strAngle != null ? Float.parseFloat(strAngle) : Float.NaN;
esgNetwork.addNode(createNode(EchUtil.getBusId(dl), dl.getTerminal().getVoltageLevel(), v, angle, false));
}
}
private static EsgBranchConnectionStatus getStatus(ConnectionBus bus1, ConnectionBus bus2) {
if (!bus1.isConnected() && !bus2.isConnected()) {
return EsgBranchConnectionStatus.OPEN_AT_BOTH_SIDES;
} else if (bus1.isConnected() && bus2.isConnected()) {
return EsgBranchConnectionStatus.CLOSED_AT_BOTH_SIDE;
} else {
return bus1.isConnected() ? EsgBranchConnectionStatus.OPEN_AT_RECEIVING_SIDE
: EsgBranchConnectionStatus.OPEN_AT_SENDING_SIDE;
}
}
private void createCouplingDevices(EsgNetwork esgNetwork) {
for (VoltageLevel vl : Identifiables.sort(network.getVoltageLevels())) {
for (Switch sw : Identifiables.sort(EchUtil.getSwitches(vl, config))) {
Bus bus1 = EchUtil.getBus1(vl, sw.getId(), config);
Bus bus2 = EchUtil.getBus2(vl, sw.getId(), config);
esgNetwork.addCouplingDevice(new EsgCouplingDevice(new EsgBranchName(new Esg8charName(dictionary.getEsgId(bus1.getId())),
new Esg8charName(dictionary.getEsgId(bus2.getId())),
parallelIndexes.getParallelIndex(sw.getId()))
, sw.isOpen() ? EsgCouplingDevice.ConnectionStatus.OPEN : EsgCouplingDevice.ConnectionStatus.CLOSED));
}
}
}
private EsgLine createLine(String id, ConnectionBus bus1, ConnectionBus bus2, float nominalV, float r, float x, float g,
float b, EsgGeneralParameters parameters) {
EsgBranchConnectionStatus status = getStatus(bus1, bus2);
float rate = parameters.getSnref();
float vnom2 = (float) Math.pow(nominalV, 2);
float rb = r * parameters.getSnref() / vnom2;
float rxb = x * parameters.getSnref() / vnom2;
float gs = g / parameters.getSnref() * vnom2;
float bs = b / parameters.getSnref() * vnom2;
return new EsgLine(new EsgBranchName(new Esg8charName(dictionary.getEsgId(bus1.getId())),
new Esg8charName(dictionary.getEsgId(bus2.getId())),
parallelIndexes.getParallelIndex(id)),
status, rb, rxb, gs, bs, rate);
}
private EsgDissymmetricalBranch createDissymmetricalBranch(String id, ConnectionBus bus1, ConnectionBus bus2,
float nominalV, float r, float x, float g1, float b1, float g2, float b2,
EsgGeneralParameters parameters) {
EsgBranchConnectionStatus status = getStatus(bus1, bus2);
float rate = parameters.getSnref();
float vnom2 = (float) Math.pow(nominalV, 2);
float rb = (r * parameters.getSnref()) / vnom2;
float rxb = (x * parameters.getSnref()) / vnom2;
float gs1 = (g1 / parameters.getSnref()) * vnom2;
float bs1 = (b1 / parameters.getSnref()) * vnom2;
float gs2 = (g2 / parameters.getSnref()) * vnom2;
float bs2 = (b2 / parameters.getSnref()) * vnom2;
return new EsgDissymmetricalBranch(new EsgBranchName(new Esg8charName(bus1.getId()),
new Esg8charName(bus2.getId()),
parallelIndexes.getParallelIndex(id)),
status, rb / 2, rxb / 2, gs1, bs1, rate, rb / 2, rxb / 2, gs2, bs2);
}
private void createLines(EsgNetwork esgNetwork, EsgGeneralParameters parameters) {
for (Line l : Identifiables.sort(network.getLines())) {
ConnectionBus bus1 = ConnectionBus.fromTerminal(l.getTerminal1(), config, EchUtil.FAKE_NODE_NAME1);
ConnectionBus bus2 = ConnectionBus.fromTerminal(l.getTerminal2(), config, EchUtil.FAKE_NODE_NAME2);
// if the admittance are the same in the both side of PI line model
if (Math.abs(l.getG1() - l.getG2()) < G_EPSILON && Math.abs(l.getB1() - l.getB2()) < B_EPSILON) {
//...create a simple line
esgNetwork.addLine(createLine(l.getId(), bus1, bus2, l.getTerminal1().getVoltageLevel().getNominalV(),
l.getR(), l.getX(), l.getG1(), l.getB1(), parameters));
}
else {
// create a dissymmetrical branch
esgNetwork.addDissymmetricalBranch(createDissymmetricalBranch(l.getId(), bus1, bus2, l.getTerminal1().getVoltageLevel().getNominalV(),
l.getR(), l.getX(), l.getG1(), l.getB1(), l.getG2(), l.getB2(), parameters));
}
}
for (DanglingLine dl : Identifiables.sort(network.getDanglingLines())) {
ConnectionBus bus1 = ConnectionBus.fromTerminal(dl.getTerminal(), config, EchUtil.FAKE_NODE_NAME1);
ConnectionBus bus2 = new ConnectionBus(true, EchUtil.getBusId(dl));
esgNetwork.addLine(createLine(dl.getId(), bus1, bus2, dl.getTerminal().getVoltageLevel().getNominalV(),
dl.getR(), dl.getX(), dl.getG() / 2, dl.getB() / 2, parameters));
}
}
private EsgDetailedTwoWindingTransformer.Tap createTap(TwoWindingsTransformer twt, int iplo, float rho, float dr, float dx,
float dephas, float rate, EsgGeneralParameters parameters) {
float nomiU2 = twt.getTerminal2().getVoltageLevel().getNominalV();
float rho_i = twt.getRatedU2() / twt.getRatedU1() * rho;
float uno1 = nomiU2 / rho_i;
float uno2 = nomiU2;
float r = twt.getR() * (1 + dr / 100.0f);
float x = twt.getX() * (1 + dx / 100.0f);
//...mTrans.getR() = Get the nominal series resistance specified in Ω at the secondary voltage side.
float zb2 = (float) (Math.pow(nomiU2, 2) / parameters.getSnref());
float rpu2 = r / zb2; //...total line resistance [p.u.](Base snref)
float xpu2 = x / zb2; //...total line reactance [p.u.](Base snref)
//...leakage impedance [%] (base rate)
float ucc;
if (xpu2 < 0) {
ucc = xpu2 * 100f * rate / parameters.getSnref();
} else {
float zpu2 = (float) Math.hypot(rpu2, xpu2);
ucc = zpu2 * 100f * rate / parameters.getSnref();
}
return new EsgDetailedTwoWindingTransformer.Tap(iplo, dephas, uno1, uno2, ucc);
}
private void createTransformers(EsgNetwork esgNetwork, EsgGeneralParameters parameters) {
for (TwoWindingsTransformer twt : Identifiables.sort(network.getTwoWindingsTransformers())) {
ConnectionBus bus1 = ConnectionBus.fromTerminal(twt.getTerminal1(), config, EchUtil.FAKE_NODE_NAME1);
ConnectionBus bus2 = ConnectionBus.fromTerminal(twt.getTerminal2(), config, EchUtil.FAKE_NODE_NAME2);
EsgBranchConnectionStatus status = getStatus(bus1, bus2);
//...IIDM gives no rate value. we take rate = 100 MVA But we have to find the corresponding pcu, pfer ...
float rate = 100.f;
//**************************
//*** LOSSES COMPUTATION *** (Record 1)
//**************************
float nomiU2 = twt.getTerminal2().getVoltageLevel().getNominalV();
//...mTrans.getR() = Get the nominal series resistance specified in Ω at the secondary voltage side.
float Rpu2 = (twt.getR() * parameters.getSnref()) / nomiU2 / nomiU2; //...total line resistance [p.u.](Base snref)
float Gpu2 = (twt.getG() / parameters.getSnref()) * nomiU2 * nomiU2; //...semi shunt conductance [p.u.](Base snref)
float Bpu2 = (twt.getB() / parameters.getSnref()) * nomiU2 * nomiU2; //...semi shunt susceptance [p.u.](Base snref)
//...changing base snref -> base rate to compute losses
float pcu = Rpu2 * rate * 100f / parameters.getSnref(); //...base rate (100F -> %)
float pfer = 10000f * ( Gpu2 / rate) * (parameters.getSnref() / 100f) ; //...base rate
float modgb = (float) Math.sqrt(Math.pow(Gpu2,2.f) + Math.pow(Bpu2, 2.f) );
float cmagn = 10000 * ( modgb / rate) * (parameters.getSnref() / 100f); //...magnetizing current [% base rate]
float esat = 1.f;
//***************************
// *** TAP TRANSFORMATION *** (Record 2)
//***************************
EsgDetailedTwoWindingTransformer.RegulatingMode regulatingMode = EsgDetailedTwoWindingTransformer.RegulatingMode.NOT_REGULATING;
Esg8charName zbusr = null; //...regulated node name (if empty, no tap change)
float voltr = Float.NaN;
int ktpnom; //...nominal tap number is not available in IIDM. Take th median plot by default
int ktap8; //...initial tap position (tap number) (Ex: 10)
List<EsgDetailedTwoWindingTransformer.Tap> taps = new ArrayList<>();
RatioTapChanger rtc = twt.getRatioTapChanger();
PhaseTapChanger ptc = twt.getPhaseTapChanger();
if (rtc != null && ptc == null) {
if (rtc.isRegulating()) {
ConnectionBus regulatingBus = ConnectionBus.fromTerminal(rtc.getRegulationTerminal(), config, null);
if (regulatingBus.getId() != null) {
regulatingMode = EsgDetailedTwoWindingTransformer.RegulatingMode.VOLTAGE;
zbusr = new Esg8charName(dictionary.getEsgId(regulatingBus.getId()));
}
}
voltr = rtc.getTargetV();
ktap8 = rtc.getTapPosition() - rtc.getLowTapPosition() + 1;
ktpnom = rtc.getStepCount() / 2 + 1;
for (int p = rtc.getLowTapPosition(); p <= rtc.getHighTapPosition(); p++) {
int iplo = p - rtc.getLowTapPosition() + 1;
taps.add(createTap(twt, iplo, rtc.getStep(p).getRho(), rtc.getStep(p).getR(), rtc.getStep(p).getX(), 0f, rate, parameters));
}
} else if (ptc != null && rtc == null) {
if (ptc.getRegulationMode() == PhaseTapChanger.RegulationMode.CURRENT_LIMITER && ptc.isRegulating()) {
String regulbus = EchUtil.getBus(ptc.getRegulationTerminal(), config).getId();
if (regulbus.equals(bus1.getId())) {
regulatingMode = EsgDetailedTwoWindingTransformer.RegulatingMode.ACTIVE_FLUX_SIDE_1;
}
if (regulbus.equals(bus2.getId())) {
regulatingMode = EsgDetailedTwoWindingTransformer.RegulatingMode.ACTIVE_FLUX_SIDE_2;
}
if (regulatingMode == EsgDetailedTwoWindingTransformer.RegulatingMode.NOT_REGULATING) {
throw new ITeslaException("Phase transformer " + twt.getId() + " has an unknown regulated node");
}
}
ktap8 = ptc.getTapPosition() - ptc.getLowTapPosition() + 1;
ktpnom = ptc.getStepCount() / 2 + 1;
for (int p = ptc.getLowTapPosition(); p <= ptc.getHighTapPosition(); p++) {
int iplo = p - ptc.getLowTapPosition() + 1;
taps.add(createTap(twt, iplo, ptc.getStep(p).getRho(), ptc.getStep(p).getR(), ptc.getStep(p).getX(), ptc.getStep(p).getAlpha(), rate, parameters));
}
} else if (rtc == null && ptc == null) {
ktap8 = 1;
ktpnom = 1;
taps.add(createTap(twt, 1, 1f, 0f, 0f, 0f, rate, parameters));
} else {
throw new RuntimeException("Transformer " + twt.getId() + " with voltage and phase tap changer not yet supported");
}
float pregmin = Float.NaN; //...?
float pregmax = Float.NaN; //...?
EsgDetailedTwoWindingTransformer esgTransfo = new EsgDetailedTwoWindingTransformer(
new EsgBranchName(new Esg8charName(dictionary.getEsgId(bus1.getId())),
new Esg8charName(dictionary.getEsgId(bus2.getId())),
parallelIndexes.getParallelIndex(twt.getId())),
status,
cmagn,
rate,
pcu,
pfer,
esat,
ktpnom,
ktap8,
zbusr,
voltr,
pregmin,
pregmax,
regulatingMode);
//***************************
// *** TAP TRANSFORMATION *** (Record 3)
//***************************
esgTransfo.getTaps().addAll(taps);
esgNetwork.addDetailedTwoWindingTransformer(esgTransfo);
}
for (ThreeWindingsTransformer twt : network.getThreeWindingsTransformers()) {
throw new AssertionError("TODO");
}
}
private EsgLoad createLoad(ConnectionBus bus, String loadId, float p0, float q0) {
EsgLoad.ConnectionStatus status = bus.isConnected() ? EsgLoad.ConnectionStatus.CONNECTED : EsgLoad.ConnectionStatus.NOT_CONNECTED;
return new EsgLoad(status, new Esg8charName(dictionary.getEsgId(loadId)),
new Esg8charName(dictionary.getEsgId(bus.getId())),
0f, 0f, p0, 0f, 0f, q0);
}
private void createLoads(EsgNetwork esgNetwork) {
for (Load l : Identifiables.sort(network.getLoads())) {
ConnectionBus bus = ConnectionBus.fromTerminal(l.getTerminal(), config, EchUtil.FAKE_NODE_NAME1);
esgNetwork.addLoad(createLoad(bus, l.getId(), l.getP0(), l.getQ0()));
}
for (DanglingLine dl : Identifiables.sort(network.getDanglingLines())) {
ConnectionBus bus = new ConnectionBus(true, EchUtil.getBusId(dl));
esgNetwork.addLoad(createLoad(bus, EchUtil.getLoadId(dl), dl.getP0(), dl.getQ0()));
}
}
private void createGenerators(EsgNetwork esgNetwork) {
for (Generator g : Identifiables.sort(network.getGenerators())) {
ConnectionBus bus = ConnectionBus.fromTerminal(g.getTerminal(), config, EchUtil.FAKE_NODE_NAME1);
EsgGenerator.ConnectionStatus status = bus.isConnected() ? EsgGenerator.ConnectionStatus.CONNECTED : EsgGenerator.ConnectionStatus.NOT_CONNECTED;
float pgen = g.getTargetP();
float qgen = g.getTargetQ();
float pgmin = g.getMinP();
float pgmax = g.getMaxP();
float qgmin = config.isNoGeneratorMinMaxQ() ? -9999 : g.getReactiveLimits().getMinQ(pgen);
float qgmax = config.isNoGeneratorMinMaxQ() ? 9999 : g.getReactiveLimits().getMaxQ(pgen);
EsgGenerator.RegulatingMode mode = g.isVoltageRegulatorOn() && g.getTargetV() >= 0.1
? EsgGenerator.RegulatingMode.REGULATING : EsgGenerator.RegulatingMode.NOT_REGULATING;
float vregge = g.isVoltageRegulatorOn() ? g.getTargetV() : Float.NaN;
float qgensh = 1.f;
Bus regulatingBus = g.getRegulatingTerminal().getBusBreakerView().getConnectableBus();
esgNetwork.addGenerator(new EsgGenerator(new Esg8charName(dictionary.getEsgId(g.getId())),
new Esg8charName(dictionary.getEsgId(bus.getId())),
pgmin, pgen, pgmax, qgmin, qgen, qgmax, mode, vregge,
new Esg8charName(dictionary.getEsgId(regulatingBus.getId())),
qgensh, status));
}
}
private void createBanks(EsgNetwork esgNetwork) {
for (ShuntCompensator sc : Identifiables.sort(network.getShunts())) {
ConnectionBus bus = ConnectionBus.fromTerminal(sc.getTerminal(), config, EchUtil.FAKE_NODE_NAME1);
//...number of steps in service
int ieleba = bus.isConnected() ? sc.getCurrentSectionCount() : 0; // not really correct, because it can be connected with zero section, EUROSTAG should be modified...
float plosba = 0.f; // no active lost in the iidm shunt compensator
float vnom = sc.getTerminal().getVoltageLevel().getNominalV();
float rcapba = vnom * vnom * sc.getbPerSection();
int imaxba = sc.getMaximumSectionCount();
EsgCapacitorOrReactorBank.RegulatingMode xregba = EsgCapacitorOrReactorBank.RegulatingMode.NOT_REGULATING;
esgNetwork.addCapacitorsOrReactorBanks(new EsgCapacitorOrReactorBank(new Esg8charName(dictionary.getEsgId(sc.getId())),
new Esg8charName(dictionary.getEsgId(bus.getId())),
ieleba, plosba, rcapba, imaxba, xregba));
}
}
public EsgNetwork createNetwork(EsgGeneralParameters parameters) {
EsgNetwork esgNetwork = new EsgNetwork();
// areas
createAreas(esgNetwork);
// nodes
createNodes(esgNetwork);
// coupling devices
createCouplingDevices(esgNetwork);
// lines
createLines(esgNetwork, parameters);
// transformers
createTransformers(esgNetwork, parameters);
// loads
createLoads(esgNetwork);
// generators
createGenerators(esgNetwork);
// shunts
createBanks(esgNetwork);
return esgNetwork;
}
public void write(Writer writer, EsgGeneralParameters parameters) throws IOException {
EsgNetwork esgNetwork = createNetwork(parameters);
new EsgWriter(esgNetwork, parameters).write(writer, network.getId() + "/" + network.getStateManager().getWorkingStateId());
}
public void write(Writer writer) throws IOException {
write(writer, new EsgGeneralParameters());
}
public void write(Path file, EsgGeneralParameters parameters) throws IOException {
try (Writer writer = Files.newBufferedWriter(file, StandardCharsets.UTF_8)) {
write(writer, parameters);
}
}
public void write(Path file) throws IOException {
write(file, new EsgGeneralParameters());
}
}