summaryrefslogtreecommitdiffstats
path: root/benchmarks/sacoi1.dat
blob: 3bbd58c82ce151ce27753a8d477ddf0252b54c9f (plain)
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
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
1001
1002
1003
1004
1005
1006
1007
1008
1009
1010
1011
1012
1013
1014
1015
1016
1017
1018
1019
1020
1021
1022
1023
1024
1025
1026
1027
1028
1029
1030
1031
1032
1033
1034
1035
1036
1037
1038
1039
1040
1041
1042
1043
1044
1045
1046
1047
1048
1049
1050
1051
1052
1053
1054
1055
1056
1057
1058
1059
1060
1061
1062
1063
1064
1065
1066
1067
1068
1069
1070
1071
1072
1073
1074
1075
1076
1077
1078
1079
1080
1081
1082
1083
1084
1085
1086
1087
1088
1089
1090
1091
1092
1093
1094
1095
1096
1097
1098
1099
1100
1101
1102
1103
1104
1105
1106
1107
1108
1109
1110
1111
1112
1113
1114
1115
1116
1117
1118
1119
1120
1121
1122
1123
1124
1125
1126
1127
1128
1129
1130
1131
1132
1133
1134
1135
1136
1137
1138
1139
1140
1141
1142
1143
1144
1145
1146
1147
1148
1149
1150
1151
begin new data case
c 
c 
c ------------------------------------------------------------------------------
c miscellaneous data cards
c ------------------------------------------------------------------------------
c  dt  >< tmax >< xopt >< copt ><epsiln>
  2.5e-5     1.5
     500       1       1       1       1       0       0       1       0
c 	 
c ------------------------------------------------------------------------------
c blocco tacs
c ------------------------------------------------------------------------------
tacs hybrid
c ----- misura dell'angolo gamma, convertitore suvereto ---------------
90vysuva                                                                     10.
90vysuvb                                                                     10.
90vysuvc                                                                     10.
90vdsuva                                                                     10.
90vdsuvb                                                                     10.
90vdsuvc                                                                     10.
90elesuv                                                                     10.
90meio11                                                                     10.
90polosu                                                                     10.
c
98va1101  =suvera.lt.suverb
98va1102  =suvera.lt.suverc
98va1103  =suverb.lt.suverc
98va1104  =suverb.lt.suvera
98va1105  =suverc.lt.suvera
98va1106  =suverc.lt.suverb
98va1107  =suvera.lt.zero
98va1108  =suverc.gt.zero
98va1109  =suverb.lt.zero
98va1110  =suvera.gt.zero
98va1111  =suverc.lt.zero
98va1112  =suverb.gt.zero
c
98vd110153                +va1101                       2.5e-5
98vd110253                +va1102                       2.5e-5
98vd110353                +va1103                       2.5e-5
98vd110453                +va1104                       2.5e-5
98vd110553                +va1105                       2.5e-5
98vd110653                +va1106                       2.5e-5
98vd110753                +va1107                       2.5e-5
98vd110853                +va1108                       2.5e-5
98vd110953                +va1109                       2.5e-5
98vd111053                +va1110                       2.5e-5
98vd111153                +va1111                       2.5e-5
98vd111253                +va1112                       2.5e-5
c
98re1101  =(vysuva.lt.meio11).and.(timex.gt.0.02)
98re1102  =(elesuv.lt.vysuvc).and.(timex.gt.0.02)
98re1103  =(vysuvb.lt.meio11).and.(timex.gt.0.02)
98re1104  =(elesuv.lt.vysuva).and.(timex.gt.0.02)
98re1105  =(vysuvc.lt.meio11).and.(timex.gt.0.02)
98re1106  =(elesuv.lt.vysuvb).and.(timex.gt.0.02)
98re1107  =(vdsuva.lt.polosu).and.(timex.gt.0.02)
98re1108  =(meio11.lt.vdsuvc).and.(timex.gt.0.02)
98re1109  =(vdsuvb.lt.polosu).and.(timex.gt.0.02)
98re1110  =(meio11.lt.vdsuva).and.(timex.gt.0.02)
98re1111  =(vdsuvc.lt.polosu).and.(timex.gt.0.02)
98re1112  =(meio11.lt.vdsuvb).and.(timex.gt.0.02)
c
98tr1101  =(.not.va1101).and.vd1101
98tr1102  =(.not.va1102).and.vd1102
98tr1103  =(.not.va1103).and.vd1103
98tr1104  =(.not.va1104).and.vd1104
98tr1105  =(.not.va1105).and.vd1105
98tr1106  =(.not.va1106).and.vd1106
98tr1107  =(.not.va1107).and.vd1107
98tr1108  =(.not.va1108).and.vd1108
98tr1109  =(.not.va1109).and.vd1109
98tr1110  =(.not.va1110).and.vd1110
98tr1111  =(.not.va1111).and.vd1111
98tr1112  =(.not.va1112).and.vd1112
c
98ig110158        +re1101                          1.8e4          1.vd1101
98ig110258        +re1102                          1.8e4          1.vd1102
98ig110358        +re1103                          1.8e4          1.vd1103
98ig110458        +re1104                          1.8e4          1.vd1104
98ig110558        +re1105                          1.8e4          1.vd1105
98ig110658        +re1106                          1.8e4          1.vd1106
98ig110758        +re1107                          1.8e4          1.vd1107
98ig110858        +re1108                          1.8e4          1.vd1108
98ig110958        +re1109                          1.8e4          1.vd1109
98ig111058        +re1110                          1.8e4          1.vd1110
98ig111158        +re1111                          1.8e4          1.vd1111
98ig111258        +re1112                          1.8e4          1.vd1112
c
98gi110158        +vd1101                          1.8e4          1.ig1101
98gi110258        +vd1102                          1.8e4          1.ig1102
98gi110358        +vd1103                          1.8e4          1.ig1103
98gi110458        +vd1104                          1.8e4          1.ig1104
98gi110558        +vd1105                          1.8e4          1.ig1105
98gi110658        +vd1106                          1.8e4          1.ig1106
98gi110758        +vd1107                          1.8e4          1.ig1107
98gi110858        +vd1108                          1.8e4          1.ig1108
98gi110958        +vd1109                          1.8e4          1.ig1109
98gi111058        +vd1110                          1.8e4          1.ig1110
98gi111158        +vd1111                          1.8e4          1.ig1111
98gi111258        +vd1112                          1.8e4          1.ig1112
c
98gm110162        +gi1101                                           tr1101
98gm110262        +gi1102                                           tr1102
98gm110362        +gi1103                                           tr1103
98gm110462        +gi1104                                           tr1104
98gm110562        +gi1105                                           tr1105
98gm110662        +gi1106                                           tr1106
98gm110762        +gi1107                                           tr1107
98gm110862        +gi1108                                           tr1108
98gm110962        +gi1109                                           tr1109
98gm111062        +gi1110                                           tr1110
98gm111162        +gi1111                                           tr1111
98gm111262        +gi1112                                           tr1112
c
98tra11y  =(tr1101+tr1102+tr1103+tr1104+tr1105+tr1106)/deltat
98gm11_y61+gm1106 +gm1101 +gm1102 +gm1103 +gm1104                99.gm1105gam11y
98gam11y58                +tra11y                     1.          1.nt1106unity
98nt1106  =(.not.tr1106).eq.1
98tra11d  =(tr1107+tr1108+tr1109+tr1110+tr1111+tr1112)/deltat
98gm11_d61+gm1112 +gm1107 +gm1108 +gm1109 +gm1110                99.gm1111gam11d
98gam11d58                +tra11d                     1.          1.nt1112unity
98nt1112  =(.not.tr1112).eq.1
98alfar 63+gm11_y                         +gm11_d          -1.
c
90suvera
90suverb
90suverc
c
c -------  misura della tensione media (50hz), sbarre di commutazione ---------
98vacm11  =sqrt(suvera*suvera+suverb*suverb+suverc*suverc)/400.e3
33vacm11alfar
c
c -------  pll (per la gfu), lato suvereto  -----------------------------------
98omeg11  =2.*pi*   50.
98iga11   =(2*suvera-suverb-suverc)/(sqrt(6.)*400.e3)
98igb11   =(suverb-suverc)/(sqrt(2.)*400.e3)
88err11   =iga11*cos(p11pll)+igb11*sin(p11pll)
 1erin11  +err11                                     10.
        1.
                  1.
 0erpp11  +err11                                     50.
98res11   =phd11.le.(2*pi-2*w11*deltat)
98p11pll58                +w11                        1.          1.res11 zero
 0w11             +erpp11 +erin11 +omeg11             1. 301.6 452.4
 0phd11                   +p11pll                     1.
77p11pll        1.57
c 33p##pll
c 33w##
c ------------------  innesco dei tiristori ----------------------------------
90uctd11
11grd11   .017453292
98vad11   = sin(p11pll-(uctd11+30)*grd11)
98vbd11   = sin(p11pll-(uctd11+150)*grd11)
98vcd11   = sin(p11pll-(uctd11-90)*grd11)
c -------------------------------------------------
98vay11   = sin(p11pll-uctd11*grd11)
98vby11   = sin(p11pll-(uctd11+120)*grd11)
98vcy11   = sin(p11pll-(uctd11-120)*grd11)
98ds1101  = vay11.gt.vcy11.and.vay11.gt.vby11.and.timex.gt.   0.0
98ds1102  = vcy11.lt.vay11.and.vcy11.lt.vby11.and.timex.gt.   0.0
98ds1103  = vby11.gt.vcy11.and.vby11.gt.vay11.and.timex.gt.   0.0
98ds1104  = vay11.lt.vcy11.and.vay11.lt.vby11.and.timex.gt.   0.0
98ds1105  = vcy11.gt.vay11.and.vcy11.gt.vby11.and.timex.gt.   0.0
98ds1106  = vby11.lt.vcy11.and.vby11.lt.vay11.and.timex.gt.   0.0
c -------------------------------------------------
98ds1107  = vad11.gt.vcd11.and.vad11.gt.vbd11.and.timex.gt.   0.0
98ds1108  = vcd11.lt.vad11.and.vcd11.lt.vbd11.and.timex.gt.   0.0
98ds1109  = vbd11.gt.vcd11.and.vbd11.gt.vad11.and.timex.gt.   0.0
98ds1110  = vad11.lt.vcd11.and.vad11.lt.vbd11.and.timex.gt.   0.0
98ds1111  = vcd11.gt.vad11.and.vcd11.gt.vbd11.and.timex.gt.   0.0
98ds1112  = vbd11.lt.vcd11.and.vbd11.lt.vad11.and.timex.gt.   0.0
c
c
c ----- misura dell'angolo gamma, convertitore codrongianos --------------------
90vycoda                                                                     10.
90vycodb                                                                     10.
90vycodc                                                                     10.
90vdcoda                                                                     10.
90vdcodb                                                                     10.
90vdcodc                                                                     10.
90poloco                                                                     10.
90meio21                                                                     10.
90elecod                                                                     10.
c
98va2101  =codroa.lt.codrob
98va2102  =codroa.lt.codroc
98va2103  =codrob.lt.codroc
98va2104  =codrob.lt.codroa
98va2105  =codroc.lt.codroa
98va2106  =codroc.lt.codrob
98va2107  =codroa.lt.zero
98va2108  =codroc.gt.zero
98va2109  =codrob.lt.zero
98va2110  =codroa.gt.zero
98va2111  =codroc.lt.zero
98va2112  =codrob.gt.zero
c
98vd210153                +va2101                       2.5e-5
98vd210253                +va2102                       2.5e-5
98vd210353                +va2103                       2.5e-5
98vd210453                +va2104                       2.5e-5
98vd210553                +va2105                       2.5e-5
98vd210653                +va2106                       2.5e-5
98vd210753                +va2107                       2.5e-5
98vd210853                +va2108                       2.5e-5
98vd210953                +va2109                       2.5e-5
98vd211053                +va2110                       2.5e-5
98vd211153                +va2111                       2.5e-5
98vd211253                +va2112                       2.5e-5
c
98re2101  =(vycoda.lt.meio21).and.(timex.gt.0.02)
98re2102  =(poloco.lt.vycodc).and.(timex.gt.0.02)
98re2103  =(vycodb.lt.meio21).and.(timex.gt.0.02)
98re2104  =(poloco.lt.vycoda).and.(timex.gt.0.02)
98re2105  =(vycodc.lt.meio21).and.(timex.gt.0.02)
98re2106  =(poloco.lt.vycodb).and.(timex.gt.0.02)
98re2107  =(vdcoda.lt.elecod).and.(timex.gt.0.02)
98re2108  =(meio21.lt.vdcodc).and.(timex.gt.0.02)
98re2109  =(vdcodb.lt.elecod).and.(timex.gt.0.02)
98re2110  =(meio21.lt.vdcoda).and.(timex.gt.0.02)
98re2111  =(vdcodc.lt.elecod).and.(timex.gt.0.02)
98re2112  =(meio21.lt.vdcodb).and.(timex.gt.0.02)
c
98tr2101  =(.not.va2101).and.vd2101
98tr2102  =(.not.va2102).and.vd2102
98tr2103  =(.not.va2103).and.vd2103
98tr2104  =(.not.va2104).and.vd2104
98tr2105  =(.not.va2105).and.vd2105
98tr2106  =(.not.va2106).and.vd2106
98tr2107  =(.not.va2107).and.vd2107
98tr2108  =(.not.va2108).and.vd2108
98tr2109  =(.not.va2109).and.vd2109
98tr2110  =(.not.va2110).and.vd2110
98tr2111  =(.not.va2111).and.vd2111
98tr2112  =(.not.va2112).and.vd2112
c
98ig210158        +re2101                          1.8e4          1.vd2101
98ig210258        +re2102                          1.8e4          1.vd2102
98ig210358        +re2103                          1.8e4          1.vd2103
98ig210458        +re2104                          1.8e4          1.vd2104
98ig210558        +re2105                          1.8e4          1.vd2105
98ig210658        +re2106                          1.8e4          1.vd2106
98ig210758        +re2107                          1.8e4          1.vd2107
98ig210858        +re2108                          1.8e4          1.vd2108
98ig210958        +re2109                          1.8e4          1.vd2109
98ig211058        +re2110                          1.8e4          1.vd2110
98ig211158        +re2111                          1.8e4          1.vd2111
98ig211258        +re2112                          1.8e4          1.vd2112
c
98gi210158        +vd2101                          1.8e4          1.ig2101
98gi210258        +vd2102                          1.8e4          1.ig2102
98gi210358        +vd2103                          1.8e4          1.ig2103
98gi210458        +vd2104                          1.8e4          1.ig2104
98gi210558        +vd2105                          1.8e4          1.ig2105
98gi210658        +vd2106                          1.8e4          1.ig2106
98gi210758        +vd2107                          1.8e4          1.ig2107
98gi210858        +vd2108                          1.8e4          1.ig2108
98gi210958        +vd2109                          1.8e4          1.ig2109
98gi211058        +vd2110                          1.8e4          1.ig2110
98gi211158        +vd2111                          1.8e4          1.ig2111
98gi211258        +vd2112                          1.8e4          1.ig2112
c
98gm210162        +gi2101                                           tr2101
98gm210262        +gi2102                                           tr2102
98gm210362        +gi2103                                           tr2103
98gm210462        +gi2104                                           tr2104
98gm210562        +gi2105                                           tr2105
98gm210662        +gi2106                                           tr2106
98gm210762        +gi2107                                           tr2107
98gm210862        +gi2108                                           tr2108
98gm210962        +gi2109                                           tr2109
98gm211062        +gi2110                                           tr2110
98gm211162        +gi2111                                           tr2111
98gm211262        +gi2112                                           tr2112
c
98tra21y  =(tr2101+tr2102+tr2103+tr2104+tr2105+tr2106)/deltat
98gm21_y61+gm2106 +gm2101 +gm2102 +gm2103 +gm2104                99.gm2105gam21y
98gam21y58                +tra21y                     1.          1.nt2106unity
98nt2106  =(.not.tr2106).eq.1
98tra21d  =(tr2107+tr2108+tr2109+tr2110+tr2111+tr2112)/deltat
98gm21_d61+gm2112 +gm2107 +gm2108 +gm2109 +gm2110                99.gm2111gam21d
98gam21d58                +tra21d                     1.          1.nt2112unity
98nt2112  =(.not.tr2112).eq.1
98gammai63+gm21_y                         +gm21_d          -1.
c
90codroa
90codrob
90codroc
c
c -------  misura della tensione media (50hz), sbarre di commutazione ---------
98vacm21  =sqrt(codroa*codroa+codrob*codrob+codroc*codroc)/400.e3
33vacm21gammai
c
c -------  pll (per la gfu), lato codrongianos  -----------------------------------
98omeg21  =2.*pi*   50.
98iga21   =(2*codroa-codrob-codroc)/(sqrt(6.)*400.e3)
98igb21   =(codrob-codroc)/(sqrt(2.)*400.e3)
88err21   =iga21*cos(p21pll)+igb21*sin(p21pll)
 1erin21  +err21                                     10.
        1.
                  1.
 0erpp21  +err21                                     50.
98res21   =phd21.le.(2*pi-2*w21*deltat)
98p21pll58                +w21                        1.          1.res21 zero
 0w21             +erpp21 +erin21 +omeg21             1. 301.6 452.4
 0phd21                   +p21pll                     1.
77p21pll        1.57
c 33p##pll
c 33w##
c ------------------  innesco dei tiristori ----------------------------------
90uctd21
11grd21   .017453292
98vad21   = sin(p21pll-(uctd21+30)*grd21)
98vbd21   = sin(p21pll-(uctd21+150)*grd21)
98vcd21   = sin(p21pll-(uctd21-90)*grd21)
c -------------------------------------------------
98vay21   = sin(p21pll-uctd21*grd21)
98vby21   = sin(p21pll-(uctd21+120)*grd21)
98vcy21   = sin(p21pll-(uctd21-120)*grd21)
98ds2101  = vay21.gt.vcy21.and.vay21.gt.vby21.and.timex.gt.   0.0
98ds2102  = vcy21.lt.vay21.and.vcy21.lt.vby21.and.timex.gt.   0.0
98ds2103  = vby21.gt.vcy21.and.vby21.gt.vay21.and.timex.gt.   0.0
98ds2104  = vay21.lt.vcy21.and.vay21.lt.vby21.and.timex.gt.   0.0
98ds2105  = vcy21.gt.vay21.and.vcy21.gt.vby21.and.timex.gt.   0.0
98ds2106  = vby21.lt.vcy21.and.vby21.lt.vay21.and.timex.gt.   0.0
c -------------------------------------------------
98ds2107  = vad21.gt.vcd21.and.vad21.gt.vbd21.and.timex.gt.   0.0
98ds2108  = vcd21.lt.vad21.and.vcd21.lt.vbd21.and.timex.gt.   0.0
98ds2109  = vbd21.gt.vcd21.and.vbd21.gt.vad21.and.timex.gt.   0.0
98ds2110  = vad21.lt.vcd21.and.vad21.lt.vbd21.and.timex.gt.   0.0
98ds2111  = vcd21.gt.vad21.and.vcd21.gt.vbd21.and.timex.gt.   0.0
98ds2112  = vbd21.lt.vcd21.and.vbd21.lt.vad21.and.timex.gt.   0.0
c 
c ----- misura dell'angolo gamma, convertitore lucciana  -----------------------
90vyluca                                                                     10.
90vylucb                                                                     10.
90vylucc                                                                     10.
90vdluca                                                                     10.
90vdlucb                                                                     10.
90vdlucc                                                                     10.
90pololu                                                                     10.
90meio31                                                                     10.
90eleluc                                                                     10.
c
98va3101  =luccia.lt.luccib
98va3102  =luccia.lt.luccic
98va3103  =luccib.lt.luccic
98va3104  =luccib.lt.luccia
98va3105  =luccic.lt.luccia
98va3106  =luccic.lt.luccib
98va3107  =luccia.lt.zero
98va3108  =luccic.gt.zero
98va3109  =luccib.lt.zero
98va3110  =luccia.gt.zero
98va3111  =luccic.lt.zero
98va3112  =luccib.gt.zero
c
98vd310153                +va3101                       2.5e-5
98vd310253                +va3102                       2.5e-5
98vd310353                +va3103                       2.5e-5
98vd310453                +va3104                       2.5e-5
98vd310553                +va3105                       2.5e-5
98vd310653                +va3106                       2.5e-5
98vd310753                +va3107                       2.5e-5
98vd310853                +va3108                       2.5e-5
98vd310953                +va3109                       2.5e-5
98vd311053                +va3110                       2.5e-5
98vd311153                +va3111                       2.5e-5
98vd311253                +va3112                       2.5e-5
c
98re3101  =(vyluca.lt.meio31).and.(timex.gt.0.02)
98re3102  =(pololu.lt.vylucc).and.(timex.gt.0.02)
98re3103  =(vylucb.lt.meio31).and.(timex.gt.0.02)
98re3104  =(pololu.lt.vyluca).and.(timex.gt.0.02)
98re3105  =(vylucc.lt.meio31).and.(timex.gt.0.02)
98re3106  =(pololu.lt.vylucb).and.(timex.gt.0.02)
98re3107  =(vdluca.lt.eleluc).and.(timex.gt.0.02)
98re3108  =(meio31.lt.vdlucc).and.(timex.gt.0.02)
98re3109  =(vdlucb.lt.eleluc).and.(timex.gt.0.02)
98re3110  =(meio31.lt.vdluca).and.(timex.gt.0.02)
98re3111  =(vdlucc.lt.eleluc).and.(timex.gt.0.02)
98re3112  =(meio31.lt.vdlucb).and.(timex.gt.0.02)
c
98tr3101  =(.not.va3101).and.vd3101
98tr3102  =(.not.va3102).and.vd3102
98tr3103  =(.not.va3103).and.vd3103
98tr3104  =(.not.va3104).and.vd3104
98tr3105  =(.not.va3105).and.vd3105
98tr3106  =(.not.va3106).and.vd3106
98tr3107  =(.not.va3107).and.vd3107
98tr3108  =(.not.va3108).and.vd3108
98tr3109  =(.not.va3109).and.vd3109
98tr3110  =(.not.va3110).and.vd3110
98tr3111  =(.not.va3111).and.vd3111
98tr3112  =(.not.va3112).and.vd3112
c
98ig310158        +re3101                          1.8e4          1.vd3101
98ig310258        +re3102                          1.8e4          1.vd3102
98ig310358        +re3103                          1.8e4          1.vd3103
98ig310458        +re3104                          1.8e4          1.vd3104
98ig310558        +re3105                          1.8e4          1.vd3105
98ig310658        +re3106                          1.8e4          1.vd3106
98ig310758        +re3107                          1.8e4          1.vd3107
98ig310858        +re3108                          1.8e4          1.vd3108
98ig310958        +re3109                          1.8e4          1.vd3109
98ig311058        +re3110                          1.8e4          1.vd3110
98ig311158        +re3111                          1.8e4          1.vd3111
98ig311258        +re3112                          1.8e4          1.vd3112
c
98gi310158        +vd3101                          1.8e4          1.ig3101
98gi310258        +vd3102                          1.8e4          1.ig3102
98gi310358        +vd3103                          1.8e4          1.ig3103
98gi310458        +vd3104                          1.8e4          1.ig3104
98gi310558        +vd3105                          1.8e4          1.ig3105
98gi310658        +vd3106                          1.8e4          1.ig3106
98gi310758        +vd3107                          1.8e4          1.ig3107
98gi310858        +vd3108                          1.8e4          1.ig3108
98gi310958        +vd3109                          1.8e4          1.ig3109
98gi311058        +vd3110                          1.8e4          1.ig3110
98gi311158        +vd3111                          1.8e4          1.ig3111
98gi311258        +vd3112                          1.8e4          1.ig3112
c
98gm310162        +gi3101                                           tr3101
98gm310262        +gi3102                                           tr3102
98gm310362        +gi3103                                           tr3103
98gm310462        +gi3104                                           tr3104
98gm310562        +gi3105                                           tr3105
98gm310662        +gi3106                                           tr3106
98gm310762        +gi3107                                           tr3107
98gm310862        +gi3108                                           tr3108
98gm310962        +gi3109                                           tr3109
98gm311062        +gi3110                                           tr3110
98gm311162        +gi3111                                           tr3111
98gm311262        +gi3112                                           tr3112
c
98tra31y  =(tr3101+tr3102+tr3103+tr3104+tr3105+tr3106)/deltat
98gm31_y61+gm3106 +gm3101 +gm3102 +gm3103 +gm3104                99.gm3105gam31y
98gam31y58                +tra31y                     1.          1.nt3106unity
98nt3106  =(.not.tr3106).eq.1
98tra31d  =(tr3107+tr3108+tr3109+tr3110+tr3111+tr3112)/deltat
98gm31_d61+gm3112 +gm3107 +gm3108 +gm3109 +gm3110                99.gm3111gam31d
98gam31d58                +tra31d                     1.          1.nt3112unity
98nt3112  =(.not.tr3112).eq.1
98gammal63+gm31_y                         +gm31_d          -1.
c
90luccia
90luccib
90luccic
c
c -------  misura della tensione media (50hz), sbarre di commutazione ---------
98vacm31  =sqrt(luccia*luccia+luccib*luccib+luccic*luccic)/90.e3
33vacm31gammal
c
c -------  pll (per la gfu), lato lucciana  ------------------------------------
98omeg31  =2.*pi*   50.
98iga31   =(2*luccia-luccib-luccic)/(sqrt(6.)*90.e3)
98igb31   =(luccib-luccic)/(sqrt(2.)*90.e3)
88err31   =iga31*cos(p31pll)+igb31*sin(p31pll)
 1erin31  +err31                                     10.
        1.
                  1.
 0erpp31  +err31                                     50.
98res31   =phd31.le.(2*pi-2*w31*deltat)
98p31pll58                +w31                        1.          1.res31 zero
 0w31             +erpp31 +erin31 +omeg31             1. 301.6 452.4
 0phd31                   +p31pll                     1.
77p31pll        1.57
c ------------------  innesco dei tiristori ----------------------------------
90uctd31
11grd31   .017453292
98vad31   = sin(p31pll-(uctd31+30)*grd31)
98vbd31   = sin(p31pll-(uctd31+150)*grd31)
98vcd31   = sin(p31pll-(uctd31-90)*grd31)
c -------------------------------------------------
98vay31   = sin(p31pll-uctd31*grd31)
98vby31   = sin(p31pll-(uctd31+120)*grd31)
98vcy31   = sin(p31pll-(uctd31-120)*grd31)
98ds3101  = vay31.gt.vcy31.and.vay31.gt.vby31.and.timex.gt.   0.0
98ds3102  = vcy31.lt.vay31.and.vcy31.lt.vby31.and.timex.gt.   0.0
98ds3103  = vby31.gt.vcy31.and.vby31.gt.vay31.and.timex.gt.   0.0
98ds3104  = vay31.lt.vcy31.and.vay31.lt.vby31.and.timex.gt.   0.0
98ds3105  = vcy31.gt.vay31.and.vcy31.gt.vby31.and.timex.gt.   0.0
98ds3106  = vby31.lt.vcy31.and.vby31.lt.vay31.and.timex.gt.   0.0
c -------------------------------------------------
98ds3107  = vad31.gt.vcd31.and.vad31.gt.vbd31.and.timex.gt.   0.0
98ds3108  = vcd31.lt.vad31.and.vcd31.lt.vbd31.and.timex.gt.   0.0
98ds3109  = vbd31.gt.vcd31.and.vbd31.gt.vad31.and.timex.gt.   0.0
98ds3110  = vad31.lt.vcd31.and.vad31.lt.vbd31.and.timex.gt.   0.0
98ds3111  = vcd31.gt.vad31.and.vcd31.gt.vbd31.and.timex.gt.   0.0
98ds3112  = vbd31.lt.vcd31.and.vbd31.lt.vad31.and.timex.gt.   0.0
blank tacs
c ------------------------------------------------------------------------------
c blocco models
c ------------------------------------------------------------------------------
models
input
mm0001 {i(dccodr)}
mm0002 {v(dccodr)}
mm0003 {tacs(gammai)}
mm0004 {i(dcsuve)}
mm0005 {v(dcsuve)}
mm0006 {i(dclucc)}
mm0007 {v(dclucc)}
mm0008 {tacs(gammal)}
output
  uccai
  uccar
  uccal
  iordem
model rad_suv
      --- modello raddrizzatore suvereto --- 
     input  idc,vdc
     output alfa,iorder
     data   vdc_ba{dflt:200.},idc_ba{dflt:1.},kpv{dflt:63},kiv{dflt:5247}
            tc_idc{dflt:0.012},tc_vdc{dflt:0.00531},iordep{dflt:1.0}
			vdcref{dflt:1.0},idcref{dflt:1.1},kpi{dflt:63},kii{dflt:5247}
     var    idc_pu,vdc_pu,idcmed,vdcmed,alfa,errv,errv_int,errv_pro
	        betav,erri,erri_int,erri_pro,betai,beta
			iorder,udinp,vdcoli
     init
 ---- misure
            histdef(idcmed)  :=1.0
            histdef(vdcmed)  :=1.0
 ---- regolatori
            errv             :=0
            errv_pro         :=0
            histdef(errv_int):=163
            erri             :=0.1			
			erri_pro         :=0
			histdef(erri_int):=175
     endinit
  exec
 ---- filtri grandezze dc ---
     idc_pu:=idc/(idc_ba*1000.)
     vdc_pu:=vdc/(vdc_ba*1000.)
     laplace(idcmed/idc_pu):=(1|s0)/(1|s0+tc_idc|s1)
     laplace(vdcmed/vdc_pu):=(1|s0)/(1|s0+tc_vdc|s1)
 ---- vdcol ----
     udinp  := vdcmed {min:0.4,max:0.9}
	 vdcoli := udinp*0.9 + 0.19
     iorder := vdcoli*iordep 
 -------------------------------------------------------------------------
 
 ---- regolatore tensione costante (modalit?? constant-voltage) ----
     errv := (vdcref - vdcmed)
     laplace(errv_int/errv) {dmin:90.,dmax:175.}:=(kiv|s0)/(1|s1)
     errv_pro := kpv*(vdcref - vdcmed)
     betav := (errv_int + errv_pro) {min:90.,max:175.}
	 
 ---- regolatore di corrente (modalit?? cc) ----
     erri := (iorder - idcmed + 0.1)
     laplace(erri_int/erri) {dmin:30.,dmax:175.} := (kii|s0)/(1|s1)	
     erri_pro := kpi*(iorder - idcmed)
     betai := (erri_int	+ erri_pro) {min:30.,max:175}
     
	 beta :=min(betai,betav)	
     alfa := 180. - beta	 
   endexec
endmodel  
model inv_cod
     ---- modello inverter codrongianos ----
     input  idc,vdc,gamma,iorder
     output alfa
     data   vdc_ba{dflt:200.},idc_ba{dflt:1.}
            freq{dflt:50},tc_idc{dflt:0.0012},tc_vdc{dflt:0.02}
            kpi{dflt:36},kii{dflt:3760},kpg{dflt:0.7506},kig{dflt:18.38}
			icodro{dflt:0.75}
     var    idc_pu,vdc_pu,idcmed,vdcmed,alfa
            betaic,erri_int,erri_pro,erri
			gamman,gammin,errg_int,errg_pro,gnlg,gerri,betaig,g[1..12]
			k,inter,betai
     init
 ---- misure
            histdef(idcmed)  :=0.75
            histdef(vdcmed)  :=0.92
 ---- regolatori
            erri             :=0
            erri_pro         :=0
            histdef(erri_int):=34.
			errg_pro         :=0
			histdef(errg_int):=30.
			gamman           :=22.
			gammin           :=180
			k                :=1
			g[1..12]         :=22.
			inter            :=0.0
     endinit
  exec
     --- filtri ---
     idc_pu:=idc/(idc_ba*1000.)
     vdc_pu:=vdc/(vdc_ba*1000.)
     laplace(idcmed/idc_pu):=(1|s0)/(1|s0+tc_idc|s1)
     laplace(vdcmed/vdc_pu):=(1|s0)/(1|s0+tc_vdc|s1)

 -------------------------------------------------------------------------
 --- regolatore corrente (modalita cc) ----
     erri := (iorder*icodro - idcmed)
     laplace(erri_int/erri) {dmin:30.,dmax:90.}:=(kii|s0)/(1|s1)
     erri_pro:=kpi*erri
     betaic:=(erri_int + erri_pro) {min:30.,max:90.}
     
	 
 -------------------------------------------------------------------------
 --- regolatore angolo costante (modalit?? cea) ----
     gerri:=15-gamman
     gnlg:=max(-31,gerri)
     laplace(errg_int/gnlg) {dmin:30.,dmax:90.}:=(kig|s0)/(1|s1)
	 errg_pro:=kpg*gnlg
	 betaig:=(errg_int + errg_pro) {min:30.,max:90.}	 	 
	 if and(t>0.03,inter>30) then
	    g[k]:=gamma
	    k:=k+1
	    if k=13 then k:=1 endif
	    gammin:=min(g[1],g[2],g[3],g[4],g[5],g[6])
	    gamman:=min(gammin,g[7],g[8],g[9],g[10],g[11],g[12])
	    inter:=0.0
	 endif
	 inter:=inter+timestep*360*freq	 
	 
	 betai:=max(betaic,betaig)
	 alfa:=180-betai
   endexec
endmodel
model inv_luc
     ---- modello inverter lucciana ----
     input  idc,vdc,gamma,iorder
     output alfa
     data   vdc_ba{dflt:200.},idc_ba{dflt:1.}
            freq{dflt:50},tc_idc{dflt:0.0012},tc_vdc{dflt:0.02}
            kpi{dflt:36},kii{dflt:3760},kpg{dflt:0.7506},kig{dflt:18.38}
			ilucci{dflt:0.25}
     var    idc_pu,vdc_pu,idcmed,vdcmed,alfa
            betaic,erri_int,erri_pro,erri
			gamman,gammin,errg_int,errg_pro,gnlg,gerri,betaig,g[1..12]
			k,inter,betai
     init
 ---- misure
            histdef(idcmed)  :=0.25
            histdef(vdcmed)  :=0.97
 ---- regolatori
            erri             :=0
            erri_pro         :=0
            histdef(erri_int):=30.
			errg_pro         :=0
			histdef(errg_int):=32.
			gamman           :=15.
			gammin           :=180
			k                :=1
			g[1..12]         :=180.
			inter            :=0.0
     endinit
  exec
     --- filtri ---
     idc_pu:=idc/(idc_ba*1000.)
     vdc_pu:=vdc/(vdc_ba*1000.)
     laplace(idcmed/idc_pu):=(1|s0)/(1|s0+tc_idc|s1)
     laplace(vdcmed/vdc_pu):=(1|s0)/(1|s0+tc_vdc|s1)

 -------------------------------------------------------------------------
 --- regolatore corrente (modalita cc) ----
     erri := (iorder*ilucci - idcmed)
     laplace(erri_int/erri) {dmin:30.,dmax:110.}:=(kii|s0)/(1|s1)
     erri_pro:=kpi*erri
     betaic:=(erri_int + erri_pro) {min:30.,max:110.}
     
	 
 -------------------------------------------------------------------------
 --- regolatore angolo costante (modalit?? cea) ----
     gerri:=15-gamman
     gnlg:=max(-31,gerri)
     laplace(errg_int/gnlg) {dmin:30.,dmax:90.}:=(kig|s0)/(1|s1)
	 errg_pro:=kpg*gnlg
	 betaig:=(errg_int + errg_pro) {min:30.,max:90.}	 	 
	 if and(t>0.03,inter>30) then
	    g[k]:=gamma
	    k:=k+1
	    if k=13 then k:=1 endif
	    gammin:=min(g[1],g[2],g[3],g[4],g[5],g[6])
	    gamman:=min(gammin,g[7],g[8],g[9],g[10],g[11],g[12])
	    inter:=0.0
	 endif
	 inter:=inter+timestep*360*freq	 
	 
	 betai:=max(betaic,betaig)
	 alfa:=180-betai
   endexec
endmodel
use rad_suv as rad_suv
input
  idc:= mm0004
  vdc:= mm0005
data
  vdc_ba:=   200.
  idc_ba:=     1.
  kpv   :=   142.
  kiv   :=  7194.
  kpi   :=    43.
  kii   :=  3547.
  tc_idc:= 0.0012
  tc_vdc:=   0.02
output
  uccar := alfa
  iordem:= iorder
enduse
use inv_cod as inv_cod
input
  idc   := mm0001
  vdc   := mm0002
  gamma := mm0003
  iorder:= iordem
data
  vdc_ba:=   200.
  idc_ba:=     1.
  kpi   :=    50.
  kii   := 476.
  kpg   := 0.4506
  kig   :=  18.38
  tc_idc:= 0.0012
  tc_vdc:=   0.02
  icodro:=   0.75
output
  uccai:=alfa
enduse
use inv_luc as inv_luc
input
  idc   := mm0006
  vdc   := mm0007
  gamma := mm0008
  iorder:= iordem
data
  vdc_ba:=   200.
  idc_ba:=     1.
  kpi   :=   65.
  kii   :=   568.
  kpg   := 0.4506
  kig   :=  18.38
  tc_idc:= 0.0012
  tc_vdc:=   0.02
  ilucci:=   0.25
output
  uccal:=alfa
enduse
record
  inv_cod.gamman as gamman
  inv_cod.betaig as betaig
  inv_cod.betaic as betaic
  inv_cod.erri as errii
  inv_cod.gerri as gerri
  inv_cod.erri_int as erri_i
  inv_cod.erri_pro as erri_p
  inv_cod.errg_int as erri_n
  inv_cod.errg_pro as erri_r
  rad_suv.erri as errir
  rad_suv.alfa as alfar
  rad_suv.idc_pu as idc_pu
  rad_suv.idcmed as idcmed
  rad_suv.vdc_pu as vdc_pu
  rad_suv.vdcmed as vdcmed
  rad_suv.errv as errv
  rad_suv.errv_int as errv_i
  rad_suv.errv_pro as errv_p
  rad_suv.erri_int as erri_i
  rad_suv.erri_pro as erri_p
  rad_suv.betav as betav
  rad_suv.betai as betai
  rad_suv.beta as beta 
endmodels
c        1         2         3         4         5         6         7         8
c 345678901234567890123456789012345678901234567890123456789012345678901234567890
blank models
c ------------------------------------------------------------------------------
c branch cards
c ------------------------------------------------------------------------------
c
c -------------------  impedenze della rete  -----------------------------------
c < n1 >< n2 ><ref1><ref2>< r  >< l  >< c  >
  fontcamisrtc            6.4454205.16                                         1
  fontcbcodrob            6.4454205.16                                         0
  fontcccodroc            6.4454205.16                                         0
  fontsamisrts            2.671285.026                                         1                                                         
  fontsbsuverb            2.671285.026                                         0                                                         
  fontscsuverc            2.671285.026                                         0
  fontlamisrtl            1.366143.483                                         1                                                         
  fontlbluccib            1.366143.483                                         0                                                         
  fontlcluccic            1.366143.483                                         0  
c -----------------  linea di trasmissione  ------------------------------------  
c < n1 >< n2 ><ref1><ref2>< r  >< l  >< c  >
  dcsuvedcsali            2.1904   37.                                         0
  dcsuve                              0.1665                                   0
  dcsali                              0.1665                                   0 
  dcsalidclucc            1.8621 952.2                                         0
  dcsali                               21.16                                   0
  dclucc                               21.16                                   0   
  dcluccdcboni             9.157  156.                                         0
  dclucc                               0.702                                   0
  dcboni                               0.702                                   0   
  dcbonidcsant             0.264  135.                                         0
  dcboni                                  3.                                   0
  dcsant                                  3.                                   0   
  dcsantdccodr             5.032   85.                                         0
  dcsant                              0.3825                                   0
  dccodr                              0.3825                                   0   
c -----------------  trasformatori  --------------------------------------------  
c ----------------- suvereto ---------------------------------------------------
  transformer                         xt001a136.e3                             0
            9999
 1suve00                        271.68  400.
 2vdsuvavdsuvb                  34.963142.85
  transformer xt001a                  xt001b                                   0
 1suverb
 2vdsuvbvdsuvc
  transformer xt001a                  xt001c                                   0
 1suverc
 2vdsuvcvdsuva
  transformer                         xt002a136.e3                             0          
            9999
 1suvera                        271.68  400.
 2vysuvaxt0004                  11.65482.474
  transformer xt002a                  xt002b                                   0
 1suverb
 2vysuvbxt0004
  transformer xt002a                  xt002c                                   0
 1suverc
 2vysuvcxt0004
c -------------------- codrongianos --------------------------------------------
  transformer                         xt005a136.e3                             0                    
            9999
 1codr00                        271.68  400.
 2vycodaxt0006                  11.65484.211
  transformer xt005a                  xt005b                                   0
 1codrob
 2vycodbxt0006
  transformer xt005a                  xt005c                                   0
 1codroc
 2vycodcxt0006
  transformer                         xt007a136.e3                             0
            9999
 1codroa                        271.68  400.
 2vdcodavdcodb                  34.963145.86
  transformer xt007a                  xt007b                                   0
 1codrob
 2vdcodbvdcodc
  transformer xt007a                  xt007c                                   0
 1codroc
 2vdcodcvdcoda
c --------------------- lucciana -----------------------------------------------
  transformer                         xt008a13.4e3                             0                    
            9999
 1lucc00                        26.612   90.
 2vylucaxt0009                   24.09  100.
  transformer xt008a                  xt008b                                   0
 1luccib
 2vylucbxt0009
  transformer xt008a                  xt008c                                   0
 1luccic
 2vyluccxt0009
  transformer                         xt010a13.4e3                             0
            9999
 1luccia                        26.612   90.
 2vdlucavdlucb                  72.281173.21
  transformer xt010a                  xt010b                                   0
 1luccib
 2vdlucbvdlucc
  transformer xt010a                  xt010c                                   0
 1luccic
 2vdluccvdluca
c -----------------------------------------------------------------------------
c                         filtri suvereto 
c -----------------------------------------------------------------------------
c ------------------  passa-alto 2 ordine  ------------------------------------
c < n 1>< n 2><ref1><ref2>< r  >< l  >< c  >
  suveraxf002a                         0.997                                   0
  suverbxf002b                         0.997                                   0
  suvercxf002c                         0.997                                   0
  xf002a                  1639.9                                               0
  xf002b                  1639.9                                               0
  xf002c                  1639.9                                               0
  xf002a                        72.498                                         0
  xf002b                        72.498                                         0
  xf002c                        72.498                                         0
c ------------------  doppio accordo  -----------------------------------------
  suveraxf001a                  35.201 0.999                                   0
  suverbxf001b                  35.201 0.999                                   0
  suvercxf001c                  35.201 0.999                                   0
  xf001a                        17.601                                         0
  xf001b                        17.601                                         0
  xf001c                        17.601                                         0
  xf001a                               1.999                                   0
  xf001b                               1.999                                   0
  xf001c                               1.999                                   0
c -----------------------------------------------------------------------------
c                      filtri codrongianos
c -----------------------------------------------------------------------------
c ------------------  passa-alto 2 ordine  ------------------------------------
c < n 1>< n 2><ref1><ref2>< r  >< l  >< c  >
  codroaxf004a                         0.997                                   1
  codrobxf004b                         0.997                                   0
  codrocxf004c                         0.997                                   0
  xf004a                  1639.9                                               0
  xf004b                  1639.9                                               0
  xf004c                  1639.9                                               0
  xf004a                        72.498                                         0
  xf004b                        72.498                                         0
  xf004c                        72.498                                         0
c ------------------  doppio accordo  -----------------------------------------
  codroaxf003a                  35.201 0.999                                   1
  codrobxf003b                  35.201 0.999                                   0
  codrocxf003c                  35.201 0.999                                   0
  xf003a                        17.601                                         0
  xf003b                        17.601                                         0
  xf003c                        17.601                                         0
  xf003a                               1.999                                   0
  xf003b                               1.999                                   0
  xf003c                               1.999                                   0
c -----------------------------------------------------------------------------
c                       filtri lucciana
c -----------------------------------------------------------------------------
c ------------------  passa-alto 2 ordine  ------------------------------------
c < n 1>< n 2><ref1><ref2>< r  >< l  >< c  >
  lucciaxf006a                        11.315                                   0
  luccibxf006b                        11.315                                   0
  luccicxf006c                        11.315                                   0
  xf006a                  144.57                                               0
  xf006b                  144.57                                               0
  xf006c                  144.57                                               0
  xf006a                         6.391                                         0
  xf006b                         6.391                                         0
  xf006c                         6.391                                         0
c ------------------  doppio accordo  ------------------------------------------
  lucciaxf005a                   3.10311.337                                   0
  luccibxf005b                   3.10311.337                                   0
  luccicxf005c                   3.10311.337                                   0
  xf005a                         1.552                                         0
  xf005b                         1.552                                         0
  xf005c                         1.552                                         0
  xf005a                              22.674                                   0
  xf005b                              22.674                                   0
  xf005c                              22.674                                   0
c ------------------------------------------------------------------------------
c                              reattori 
c ------------------------------------------------------------------------------
  xyz005xyz004            7.9e-4 500.0
  xyz006xyz007            7.9e-4 500.0  
  xyz015xyz016            7.9e-4 500.0  
c ------------------------------------------------------------------------------
c convertitore 12 impulsi suvereto
c ------------------------------------------------------------------------------
  elesuv                    .001                                               0
  uccar uctd11              1.e0                                               0
        uctd11              1.e0                                               0
c -------------------  snubber circuits  ---------------------------------------
  vdsuvacv0711              5.e3        0.05
  vdsuvccv0811              5.e3        0.05
  vdsuvbcv0911              5.e3        0.05
  vdsuvacv1011              5.e3        0.05
  vdsuvccv1111              5.e3        0.05
  vdsuvbcv1211              5.e3        0.05
  vysuvacv0111              5.e3        0.05
  vysuvccv0211              5.e3        0.05
  vysuvbcv0311              5.e3        0.05
  vysuvacv0411              5.e3        0.05
  vysuvccv0511              5.e3        0.05
  vysuvbcv0611              5.e3        0.05
c ---------  resistori di accoppiamento valvole  -------------------------------
  polosucv0711             .0100
  polosucv0911             .0100
  polosucv1111             .0100
  meio11cv0811             .0100
  meio11cv1011             .0100
  meio11cv1211             .0100
  meio11cv0111             .0100
  meio11cv0311             .0100
  meio11cv0511             .0100
  elesuvcv0211             .0100
  elesuvcv0411             .0100
  elesuvcv0611             .0100
c ------------------------------------------------------------------------------
c convertitore 12 impulsi codrongianos
c ------------------------------------------------------------------------------
  uccai uctd21              1.e0                                               0
        uctd21              1.e0                                               0
  elecod                    .001                                               0
c -------------------  snubber circuits  ---------------------------------------
  vdcodacv0721              5.e3        0.05
  vdcodccv0821              5.e3        0.05
  vdcodbcv0921              5.e3        0.05
  vdcodacv1021              5.e3        0.05
  vdcodccv1121              5.e3        0.05
  vdcodbcv1221              5.e3        0.05
  vycodacv0121              5.e3        0.05
  vycodccv0221              5.e3        0.05
  vycodbcv0321              5.e3        0.05
  vycodacv0421              5.e3        0.05
  vycodccv0521              5.e3        0.05
  vycodbcv0621              5.e3        0.05
c ---------  resistori di accoppiamento valvole  -------------------------------
  elecodcv0721             .0100
  elecodcv0921             .0100
  elecodcv1121             .0100
  meio21cv0821             .0100
  meio21cv1021             .0100
  meio21cv1221             .0100
  meio21cv0121             .0100
  meio21cv0321             .0100
  meio21cv0521             .0100
  polococv0221             .0100
  polococv0421             .0100
  polococv0621             .0100
c ------------------------------------------------------------------------------
c convertitore 12 impulsi lucciana 
c ------------------------------------------------------------------------------
  uccal uctd31              1.e0                                               0
        uctd31              1.e0                                               0
  eleluc                    .001                                               0
c -------------------  snubber circuits  ---------------------------------------
  vdlucacv0731              5.e3        0.05
  vdlucccv0831              5.e3        0.05
  vdlucbcv0931              5.e3        0.05
  vdlucacv1031              5.e3        0.05
  vdlucccv1131              5.e3        0.05
  vdlucbcv1231              5.e3        0.05
  vylucacv0131              5.e3        0.05
  vylucccv0231              5.e3        0.05
  vylucbcv0331              5.e3        0.05
  vylucacv0431              5.e3        0.05
  vylucccv0531              5.e3        0.05
  vylucbcv0631              5.e3        0.05
c ---------  resistori di accoppiamento valvole  -------------------------------
  eleluccv0731             .0100
  eleluccv0931             .0100
  eleluccv1131             .0100
  meio31cv0831             .0100
  meio31cv1031             .0100
  meio31cv1231             .0100
  meio31cv0131             .0100
  meio31cv0331             .0100
  meio31cv0531             .0100
  pololucv0231             .0100
  pololucv0431             .0100
  pololucv0631             .0100
blank branch
c ------------------------------------------------------------------------------
c switch cards
c ------------------------------------------------------------------------------
c ---------------------  interruttori per l'inizializzazione  ------------------
c ---------------------  sbarre di commutazione  -------------------------------
  suveraxyz001    -1.000    20.e-3                                             1
  suverbxyz002    -1.000    20.e-3                                             1
  suvercxyz003    -1.000    20.e-3                                             1
  codroaxyz009    -1.000    20.e-3                                             1
  codrobxyz010    -1.000    20.e-3                                             1
  codrocxyz011    -1.000    20.e-3                                             1
  lucciaxyz012    -1.000    20.e-3                                             1
  luccibxyz013    -1.000    20.e-3                                             1
  luccicxyz014    -1.000    20.e-3                                             1
c ---------------------  linea in c.c.  ----------------------------------------  
  xyz004poloco     0.010     50.00      {sw2}
  xyz016pololu     0.010     50.00      {sw2}
  polosuxyz006     0.000     50.00      {sw2}
  xyz006xyz008    -1.000     0.000      {sw4}
c  --------------------  interruttori di misura per la corrente in linea  ------
  dccodrxyz005                                        measuring                1
  xyz007dcsuve                                        measuring                1
  dcluccxyz015                                        measuring                1
c  --------------------  interruttori di misura per la corrente tr  ------------  
  suverasuve00                                        measuring                1
  codroacodr00                                        measuring                1                                         
  luccialucc00                                        measuring                1
  misrtccodroa                                        measuring                1
  misrtssuvera                                        measuring                1
  misrtlluccia                                        measuring                1
c -----------------------  valvole  suvereto  ----------------------------------
11vysuvacv0111                        0.0e-6                    ds1101         0
11vysuvbcv0311                        0.0e-6                    ds1103         0
11vysuvccv0511                        0.0e-6                    ds1105         0
11cv0411vysuva                        0.0e-6                    ds1104         0
11cv0611vysuvb                        0.0e-6                    ds1106         0
11cv0211vysuvc                        0.0e-6                    ds1102         0
11vdsuvacv0711                        0.0e-6                    ds1107         0
11vdsuvbcv0911                        0.0e-6                    ds1109         0
11vdsuvccv1111                        0.0e-6                    ds1111         0
11cv1011vdsuva                        0.0e-6                    ds1110         0
11cv1211vdsuvb                        0.0e-6                    ds1112         0
11cv0811vdsuvc                        0.0e-6                    ds1108         0
c 
c -----------------------  valvole  codrongianos  ------------------------------  
11vycodacv0121                      555.5e-6                    ds2101         0
11vycodbcv0321                      555.5e-6                    ds2103         0
11vycodccv0521                      555.5e-6                    ds2105         0
11cv0421vycoda                      555.5e-6                    ds2104         0
11cv0621vycodb                      555.5e-6                    ds2106         0
11cv0221vycodc                      555.5e-6                    ds2102         0
11vdcodacv0721                      555.5e-6                    ds2107         0
11vdcodbcv0921                      555.5e-6                    ds2109         0
11vdcodccv1121                      555.5e-6                    ds2111         0
11cv1021vdcoda                      555.5e-6                    ds2110         0
11cv1221vdcodb                      555.5e-6                    ds2112         0
11cv0821vdcodc                      555.5e-6                    ds2108         0
c -----------------------  valvole  lucciana  ----------------------------------  
11vylucacv0131                      555.5e-6                    ds3101         0
11vylucbcv0331                      555.5e-6                    ds3103         0
11vylucccv0531                      555.5e-6                    ds3105         0
11cv0431vyluca                      555.5e-6                    ds3104         0
11cv0631vylucb                      555.5e-6                    ds3106         0
11cv0231vylucc                      555.5e-6                    ds3102         0
11vdlucacv0731                      555.5e-6                    ds3107         0
11vdlucbcv0931                      555.5e-6                    ds3109         0
11vdlucccv1131                      555.5e-6                    ds3111         0
11cv1031vdluca                      555.5e-6                    ds3110         0
11cv1231vdlucb                      555.5e-6                    ds3112         0
11cv0831vdlucc                      555.5e-6                    ds3108         0
blank switch
c ------------------------------------------------------------------------------
c source cards
c ------------------------------------------------------------------------------
c ---------------  sorgenti per l'inizializzazione  ----------------------------
c ---------------  sbarre di commutazione  -------------------------------------
14xyz001 1  327.69e3       50.      +22.       0.0                  -1
14xyz002 1  327.69e3       50.      -98.       0.0                  -1
14xyz003 1  327.69e3       50.     -218.       0.0                  -1
14xyz009 1  328.57e3       50.     -20.0       0.0                  -1
14xyz010 1  328.57e3       50.     -140.       0.0                  -1
14xyz011 1  328.57e3       50.     +100.       0.0                  -1
14xyz012 1   74.10e3       50.     -20.0       0.0                  -1
14xyz013 1   74.10e3       50.     -140.       0.0                  -1
14xyz014 1   74.10e3       50.     +100.       0.0                  -1
c ---------------  linea in c.c.  ----------------------------------------------
14xyz004-1     -750.     0.001       0.0       0.0   {i1}           -1     0.015
12poloco-1      750.               0.010       0.0   {i2}          0.0     0.015
14xyz016-1     -250.     0.001       0.0       0.0   {i1}           -1     0.015
12pololu-1      250.               0.010       0.0   {i2}          0.0     0.015
14xyz008 1    200.e3     0.001       0.0       0.0   {s3}           -1     0.015
c ---------------  passaggio angolo di innesco  --------------------------------
60uccar -1                                                                   10.
60uccai -1                                                                   10.
60uccal -1                                                                   10.
c ---------------  sorgenti reti c.a.  -----------------------------------------
c < n 1><>< ampl.  >< freq.  ><phase/t0><   a1   ><   t1   >< tstart >< tstop  >
14fontca 0   323000.       50.      -22.                           -1.       10.
14fontcb 0   323000.       50.     -142.                           -1.       10.
14fontcc 0   323000.       50.       98.                           -1.       10.
14fontsa 0   327690.       50.      22.5                           -1.       10.
14fontsb 0   327690.       50.     -97.5                           -1.       10.
14fontsc 0   327690.       50.     142.5                           -1.       10.
14fontla 0    70266.       50.      22.5                           -1.       10.
14fontlb 0    70266.       50.     -97.5                           -1.       10.
14fontlc 0    70266.       50.     142.5                           -1.       10.
blank source
  suverasuverbsuverccodroacodrobcodrocdcsuvedccodrdcluccluccialuccibluccic