QB64.org Forum

Active Forums => QB64 Discussion => Topic started by: bartok on February 20, 2021, 05:56:41 am

Title: Problem to print values on the axis inside and outside WINDOW.
Post by: bartok on February 20, 2021, 05:56:41 am

Hi,
I have a problem and I don't find a solution.
Launching the code, press 1 or 2, and then any key at time.
In the ordinate is drawn a notch every X mc/s.

Given this notch, the target is to print the corresponding value on the axis. The problem is just to print it in the right position.
In the code below, we have to see the subroutine "grafici": lines 302 to 385, but the core of the problem is in lines 365 (inside the command WINDOW) and 374 (outside).

I have thought to print the value with che command _PRINTSTRING. At the moment there is no value, but the string "-". _PRINTSTRING, unfortunately, doesn't use the local coordinates of WINDOW (line 333), but always the general coordinates in pixels of the screen in use (line 331). So, in order to use _PRINTSTRING, the idea was to transform (line 365 and 374) the general coordinates of _PRINTSTRING, into the local coordinates of WINDOW by using fx! and fy!: lines 323-328.

But, as we can see launching the code, there is no way in order to have the strings "-" to fit and match the notches, also introducing a parameter k! (line 362) to do so. In any way, the strings "-" are ot of phase with the notches. The solution is not to use double precision variables for the transformation between coordinates.

So how we can print the values on the axis? I never thought that the problem was this!

Code: QB64: [Select]

1. OPTION BASE 1
2.  
3. mockus: 'dati dell'idrogramma unitario di Mockus
4. DATA 0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2,2.1,2.2,2.3,2.4,2.5,2.6,2.7,2.8,2.9,3,3.1,3.2,3.3,3.4,3.5,3.6,3.7,3.8,3.9,4,4.1,4.2,4.3,4.4,4.5,4.6,4.7,4.8,4.9,5: 't/ta
5. DATA 0.030,0.100,0.190,0.310,0.470,0.660,0.820,0.930,0.990,1.000,0.990,0.930,0.860,0.780,0.680,0.560,0.460,0.390,0.330,0.280,0.244,0.207,0.177,0.147,0.127,0.107,0.092,0.077,0.066,0.055,0.048,0.040,0.035,0.029,0.025,0.021,0.018,0.015,0.013,0.011,0.009,0.008,0.007,0.006,0.006,0.005,0.004,0.003,0.001,0.000: ' t/tp
6.  
7. coefficienti:
8. DATA 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,18,20,22,24,26,28,30,32: 'N volte tc!.
9. DATA 0.1,0.1,0.1,0.2,0.3,0.2,0.3,0.3,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,1,1,1,1,1,1,1,1: 'discretizzazione dei coeff. di Mockus.
10.  
11. TYPE mockus
12.     MOCKUStSUta AS SINGLE
13.     MOCKUSqSUqp AS SINGLE
14. END TYPE
15.  
16. TYPE matrice1
17.     m AS INTEGER '[numero naturale] conta il numero di interazioni fino al raggiungimento di "d".
18.     mdt AS SINGLE '[ore] conta la progessione delle ore fino al raggiungimento di "d".
19.     tminuti AS SINGLE '[minuti] trasforma mdt in minuti.
20.     h1 AS SINGLE '[mm] valore h della curva di possibilità climatica.
21.     i1 AS SINGLE '[mm/ora] intensità di pioggia.
22.     i2 AS SINGLE '[mm/ora] intensità di pioggia secondo lo ietogramma Chicago o costante.
23.     DH AS SINGLE '[mm] pioggia in ogni istante temporale dt.
24.     H2 AS SINGLE '[mm] pioggia cumulata ad ogni istante temporale.
25.     he AS SINGLE '[mm] piggia efficace cumulata ad ogni istante temporale.
26.     Dhe AS SINGLE '[mm] pioggia efficace ad ogni istante temporale.
27.     MOCKUStSUta AS SINGLE
28.     MOCKUSqSUqp AS SINGLE
29.     tmSUta AS SINGLE '[-].
30.     qmSUqp AS SINGLE '[-] il valore di qm/qp corrispondente a tm/ta e' = al valore di q/qp che corrisponde a t/ta=tm/ta.
31.     qm AS SINGLE '[mc/s].
32. END TYPE
33.  
34. TYPE coefficienti
35.     d AS SINGLE
36.     tSUta AS SINGLE
37. END TYPE
38.  
39. TYPE idrogramma
40.     ore AS SINGLE
41.     minuti AS SINGLE
42.     portata AS SINGLE
43. END TYPE
44.  
45. CONST limit = 200 '50
46.  
47. DIM schermo&
48.  
49. DIM SHARED mockus(50) AS mockus
50. DIM SHARED dt!(2, 24)
51. DIM SHARED matrice1(2, 24, 50, 1) AS matrice1
52. DIM SHARED coefficienti(24) AS coefficienti
53. DIM SHARED matrice2!(2, 24, 50, 50)
54. DIM SHARED idrogramma(2, 24, 50, 1) AS idrogramma
55. DIM SHARED massimi1(2, 24, 1) AS idrogramma
56. DIM SHARED massimi2(2) AS idrogramma
57. DIM SHARED i%, n%, z%, ieto%
58. DIM SHARED m%(2, 24)
59. DIM SHARED S2! '[mm] contenuto idrico massimo del terreno.
60. DIM SHARED Ia! '[mm] perdite iniziali.
61. DIM SHARED tc! '[ore] tempo di corrivazione
62. DIM SHARED ta! '[ore] tempo di picco.
63. DIM SHARED qp! '[mc/s] portata al colmo dell'idrogramma unitario.
64. DIM SHARED k! '[-] coeff. della curva di possibilità climatica tado dal sito della Regione.
65. DIM SHARED a2! '[mm/d^n] coeff. della curva di possibilità climatica.
66. DIM SHARED n1! '[-] coeff. della curva di possibilità climatica.
67. DIM SHARED visualizzaieto%
68. DIM A1! '[mq] superficie bacino idrografico.
69. DIM L! '[m] lunghezza asta principale.
70. DIM HmaxL! '[m] punto più alto dell'astra principale.
71. DIM HminL! '[m] punto più basso dell'asta principale.
72. DIM s1! '[-] pendenza media del bacino non espressa in percentuale ma come 0,x.
73. DIM CNII! '[-] tabellato.
74. DIM CNIII '[-] discende da CNII.
75. DIM tl! '[ore] tempo di ritardo.
76. DIM keypress$
77.  
78. schermo& = _NEWIMAGE(_DESKTOPWIDTH, _DESKTOPHEIGHT, 32)
79. SCREEN schermo&
80. _SCREENMOVE 0, 0
81.  
82. A1! = 2641902
83. k! = 2.34
84. a2! = 14.6653
85. n1! = 0.45596
86. L! = 3300
87. HmaxL! = 1981
88. HminL! = 1062
89. s1! = 0.415
90. CNII! = 70
91.  
92. 'INPUT A1!
93. 'INPUT k!
94. 'INPUT a2!
95. 'INPUT n1!
96. 'INPUT L!
97. 'INPUT HmaxL!
98. 'INPUT HminL!
99. 'INPUT s1!
100. 'INPUT CNII!
101. 'INPUT d1!
102.  
103. CNIII! = (23 * CNII!) / (10 + 0.13 * CNII!)
104. tl! = 0.342 * ((L! / 1000) ^ 0.8 / (100 * s1!) ^ 0.5) * (1000 / CNIII! - 9) ^ 0.7 'formula di Mockus.
105. S2! = 25.4 * (1000 / CNIII! - 10)
106. Ia! = 0.1 * S2! 'coeff.=0.03-0.2.
107. tc! = tl! / 0.6
108. ta! = tl! / 0.9
109. qp! = 0.208 * ((A1! / 1000000) / ta!)
110.  
111. RESTORE mockus
112. FOR i% = 1 TO 100
113.     IF i% <= 50 THEN READ mockus(i%).MOCKUStSUta
114.     IF i% > 50 THEN READ mockus(i% - 50).MOCKUSqSUqp
115. NEXT i%
116.  
117. RESTORE coefficienti
118. FOR i% = 1 TO 48
119.     IF i% <= 24 THEN READ coefficienti(i%).d: coefficienti(i%).d = coefficienti(i%).d * tc!
120.     IF i% > 24 THEN READ coefficienti(i% - 24).tSUta
121. NEXT i%
122.  
123. CLS
124. PRINT "Ietogramma Chicago/costante: 1/2> "
125. DO
126.     _LIMIT 60 'limit
127.     keypress$ = INKEY$
128. LOOP UNTIL keypress$ = "1" OR keypress$ = "2"
129. visualizzaieto% = VAL(keypress$)
130.  
131. FOR ieto% = 1 TO 2
132.     CALL creamatrici
133. NEXT ieto%
134.  
135. CLS
136. CALL grafici(massimi1(visualizzaieto%, 24, 1).ore, massimi2(visualizzaieto%).ore, massimi2(visualizzaieto%).portata) 'valori passati per disegnare il grafico.
137. SLEEP
138.  
139. 'CALL risultati
140.  
141. '-----------------------------------------matrice1---------------------------------------------------------------------------------------------------------------------------------
142. SUB creamatrici ()
143.     z% = 0
144.     DO
145.         z% = z% + 1
146.         dt!(ieto%, z%) = coefficienti(z%).tSUta * ta!
147.         m%(ieto%, z%) = coefficienti(z%).d / (coefficienti(z%).tSUta * ta!) 'determina il numero di passi temporali necessari ad arrivare allla durata d considerata con il paso temporale dt!.
148.         FOR i% = 1 TO 50
149.             matrice1(ieto%, z%, i%, 1).MOCKUStSUta = mockus(i%).MOCKUStSUta
150.             matrice1(ieto%, z%, i%, 1).MOCKUSqSUqp = mockus(i%).MOCKUSqSUqp
151.             matrice1(ieto%, z%, i%, 1).m = i%
152.             matrice1(ieto%, z%, i%, 1).mdt = matrice1(ieto%, z%, i%, 1).m * dt!(ieto%, z%)
153.             matrice1(ieto%, z%, i%, 1).tminuti = matrice1(ieto%, z%, i%, 1).mdt * 60
154.             SELECT CASE (_ROUND(10 * (i% * coefficienti(z%).tSUta))) / 10 'crea il vettore dell'idrogramma unitario di Mockus con i multipli del coeff. z% considerato (0.1,0.2,0.3,0.5,1).
155.                 CASE IS <= 5 'valore massimo di t/ta di Mockus
156.                     matrice1(ieto%, z%, i%, 1).tmSUta = (_ROUND(10 * (i% * coefficienti(z%).tSUta))) / 10 '=matrice1(i%).m *  matrice1(i%).mdt / ta!
157.                 CASE ELSE
158.                     matrice1(ieto%, z%, i%, 1).tmSUta = 0
159.             END SELECT
160.         NEXT i%
161.         FOR i% = 1 TO 50
162.             n% = 1
163.             DO 'crea il vettore dell'idrogramma unitario Mockus prendendo i valori corrispondenti di ai coefficienti.
164.                 IF matrice1(ieto%, z%, i%, 1).tmSUta = matrice1(ieto%, z%, n%, 1).MOCKUStSUta THEN
165.                     matrice1(ieto%, z%, i%, 1).qmSUqp = matrice1(ieto%, z%, n%, 1).MOCKUSqSUqp
166.                     EXIT DO
167.                 ELSE
168.                     n% = n% + 1
169.                 END IF
170.             LOOP UNTIL n% = 50
171.             matrice1(ieto%, z%, i%, 1).qm = matrice1(ieto%, z%, i%, 1).qmSUqp * qp!
172.         NEXT i%
173.         FOR i% = 1 TO m%(ieto%, z%)
174.             matrice1(ieto%, z%, i%, 1).h1 = k! * a2! * matrice1(ieto%, z%, i%, 1).mdt ^ n1!
175.             IF i% = 1 THEN
176.                 matrice1(ieto%, z%, i%, 1).i1 = matrice1(ieto%, z%, i%, 1).h1 / dt!(ieto%, z%)
177.             ELSE
178.                 matrice1(ieto%, z%, i%, 1).i1 = (matrice1(ieto%, z%, i%, 1).h1 - matrice1(ieto%, z%, i% - 1, 1).h1) / dt!(ieto%, z%)
179.             END IF
180.         NEXT i%
181.         SELECT CASE ieto%
182.             CASE 1 'Chicago.
183.                 CALL chicago
184.             CASE 2 'costante
185.                 FOR i% = 1 TO m%(ieto%, z%)
186.                     matrice1(ieto%, z%, i%, 1).i2 = matrice1(ieto%, z%, m%(ieto%, z%), 1).h1 / matrice1(ieto%, z%, m%(ieto%, z%), 1).mdt
187.                 NEXT i%
188.         END SELECT
189.         FOR i% = 1 TO m%(ieto%, z%)
190.             matrice1(ieto%, z%, i%, 1).DH = matrice1(ieto%, z%, i%, 1).i2 * dt!(ieto%, z%)
191.             IF i% = 1 THEN
192.                 matrice1(ieto%, z%, i%, 1).H2 = matrice1(ieto%, z%, i%, 1).DH
193.             ELSE
194.                 matrice1(ieto%, z%, i%, 1).H2 = matrice1(ieto%, z%, i%, 1).DH + matrice1(ieto%, z%, i% - 1, 1).H2
195.             END IF
196.             SELECT CASE matrice1(ieto%, z%, i%, 1).H2 - Ia!
197.                 CASE IS < 0
198.                     matrice1(ieto%, z%, i%, 1).he = 0
199.                 CASE IS >= 0
200.                     matrice1(ieto%, z%, i%, 1).he = (matrice1(ieto%, z%, i%, 1).H2 - Ia!) ^ 2 / (matrice1(ieto%, z%, i%, 1).H2 - Ia! + S2!)
201.             END SELECT
202.             IF i% = 1 THEN
203.                 matrice1(ieto%, z%, i%, 1).Dhe = matrice1(ieto%, z%, i%, 1).he
204.             ELSE
205.                 matrice1(ieto%, z%, i%, 1).Dhe = matrice1(ieto%, z%, i%, 1).he - matrice1(ieto%, z%, i% - 1, 1).he
206.             END IF
207.         NEXT i%
208.  
209.         '-------------------------------matrice2-------------------------------------------------------------------------------------------------------------------------
210.  
211.         FOR n% = 1 TO 50 'colonna
212.             FOR i% = n% TO 50 ' riga
213.                 matrice2!(ieto%, z%, i%, n%) = matrice1(ieto%, z%, i% - n% + 1, 1).Dhe * matrice1(ieto%, z%, n%, 1).qm
214.             NEXT i%
215.         NEXT n%
216.         '-------------------------------idrogramma------------------------------------------------------------------------------------------------------------------
217.         FOR i% = 1 TO 50 'riga
218.             FOR n% = 1 TO 50 'colonna
219.                 idrogramma(ieto%, z%, i%, 1).portata = idrogramma(ieto%, z%, i%, 1).portata + matrice2!(ieto%, z%, i%, n%)
220.  
221.             NEXT n%
222.             idrogramma(ieto%, z%, i%, 1).ore = matrice1(ieto%, z%, i%, 1).mdt
223.             idrogramma(ieto%, z%, i%, 1).minuti = matrice1(ieto%, z%, i%, 1).tminuti
224.         NEXT i%
225.         '-------------------------------Massimi-----------------------------------------------------
226.         i% = 1
227.         'massimi1(ieto%, z%, 1).portata = idrogramma(ieto%, z%, 1, 1).portata
228.         DO
229.             IF massimi1(ieto%, z%, 1).portata < idrogramma(ieto%, z%, i%, 1).portata THEN
230.                 massimi1(ieto%, z%, 1).portata = idrogramma(ieto%, z%, i%, 1).portata
231.                 massimi1(ieto%, z%, 1).ore = idrogramma(ieto%, z%, i%, 1).ore
232.                 massimi1(ieto%, z%, 1).minuti = idrogramma(ieto%, z%, i%, 1).minuti
233.             END IF
234.             i% = i% + 1
235.         LOOP UNTIL i% = 50
236.         ' IF massimi2(ieto%).ore < massimi1(visualizzaieto%, z%, 1).ore THEN massimi2(ieto%).ore = massimi1(visualizzaieto%, z%, 1).ore
237.         IF massimi2(ieto%).portata < massimi1(visualizzaieto%, z%, 1).portata THEN
238.             massimi2(ieto%).portata = massimi1(visualizzaieto%, z%, 1).portata
239.             massimi2(ieto%).ore = massimi1(visualizzaieto%, z%, 1).ore
240.         END IF
241.     LOOP UNTIL z% = UBOUND(coefficienti)
242. END SUB
243. '---------------------------------------------------------------------------------------------------------------------------------------------------
244. 'Questa subroutine crea gli ietogrammi chicago per ogni iterazione z%.
245. 'Ad ogni nuovo dt! pone l'ultimo valore di pioggia come ultimo, il penultimo come primo, il terzultimo come penultimo, il quartultimo come secondo e così via.
246. 'Quando vi è un valore di z% che ha un dt! non nuovo, cerca la precedente iterazione z% con lo stesso dt! e copia i suoi valori dello ietogramma chicago in quello dell'iterazione z% in corso. I valori rimanenti sono riempiti
247. 'procedendo come prima.
248. 'Quindi per esempio l'iterazione z%=1 corrisponde a d=1tc con passo temporale dt!. L'iterazione z%=2 corrisponde a d=2tc con il medesimo passo temporale. Per cui lo ietogramma chicago per d=2tc e' dato dallo stesso ietogramma
249. 'chicago che si ha per d=1tc, con l'aggiunta, ai suoi estremi, dei nuovi valori. Il meccanismo si ripete per z%=3, che corrisponde a d=3tc, con il medesimo dt!. Per z%=4, corrispondente a d=4tc, invece, il passo temporale dt!
250. 'cambia.
251. 'Quanto detto e' necessario per evitare piccole discrepanze che potrebbero comportare il fatto che alcuni massimi risultino localmente poco piu' bassi del precedente, anche quando il trend generale e' chiaramente di aumento dei
252. 'massimi all'aumentare della durata della pioggia. In questo modo si fa l'ipotesi che all'aumentare della durata della pioggia, discretizzata con uno stesso passo temporale dt!, ogni ietogramma contenga a sua volta il precedente.
253. 'In tal modo, se si condiera per esempio d=2tc, si assume che l'altezza di pioggia corrispondente a 1tc, segua lo ietogramma relativo a d=1tc, cosi' che l'altezza di pioggia rimamente sia eclusivamente dovuta ai nuovi valori.
254. SUB chicago ()
255.     DIM test%
256.     DIM k%
257.     n% = 1
258.     DO
259.         IF z% = 1 THEN test% = 0: EXIT DO 'permette il calcolo dello ietogramma chicago della prima iterazione. test%=0
260.         IF dt!(ieto%, z%) = dt!(ieto%, z% - n%) THEN ' se trova un'iterazione z% precedente con lo stesso dt! di quella in corso, esegue i comandi successivi
261.             test% = 1
262.             k% = m%(ieto%, z%)
263.             i% = 0
264.             DO 'inserisce i nuovi valori dello ietogramma chicago dell'iterazione z% in corso.
265.                 matrice1(ieto%, z%, m%(ieto%, z%) - i%, 1).i2 = matrice1(ieto%, z%, k%, 1).i1
266.                 k% = k% - 1
267.                 IF k% = m%(ieto%, z% - n%) THEN EXIT DO 'se ha completato l'inserimento dei nuovi valori, passa oltre.
268.                 matrice1(ieto%, z%, i% + 1, 1).i2 = matrice1(ieto%, z%, k%, 1).i1
269.                 k% = k% - 1
270.                 i% = i% + 1
271.             LOOP UNTIL k% = m%(ieto%, z% - n%)
272.             IF matrice1(ieto%, z%, i%, 1).i2 = 0 THEN i% = i% - 1
273.             k% = 1
274.             DO 'inserisce lo ietogramma chicago dell'iterazione z% precedente con lo stesso dt!, nello ietogramma dell'iterazione z% in corso.
275.                 i% = i% + 1
276.                 matrice1(ieto%, z%, i%, 1).i2 = matrice1(ieto%, z% - n%, k%, 1).i2
277.                 k% = k% + 1
278.             LOOP UNTIL k% > m%(ieto%, z% - n%)
279.             EXIT DO 'siccome ha completato lo ietogramma dell'iterazione z% in corso, non serve incrementare n%, quindi passa oltre. test%=1.
280.         ELSE ' se non trova l'iterazione z% precedente con lo stesso dt! di quella in corso, incrementa n% e ripete la ricerca.
281.             test% = 0
282.             n% = n% + 1
283.         END IF
284.     LOOP UNTIL n% = z% - 1 'se non ha trovato iterazioni z% precedenti con lo stesso dt!, allora passa oltre e test%=0.
285.     SELECT CASE test%
286.         CASE IS = 1 'esce dalla subroutine, in quanto lo ietogramma e' definito.
287.             EXIT SUB
288.         CASE IS = 0 'calcola ex-novo lo ietogramma relativo all'iterazione z% in corso in quanto non vi e' un'iterazioni z% precedente con lo stesso dt! dell'iterazione z% in corso.
289.             k% = m%(ieto%, z%)
290.             i% = 0
291.             DO
292.                 matrice1(ieto%, z%, m%(ieto%, z%) - i%, 1).i2 = matrice1(ieto%, z%, k%, 1).i1
293.                 k% = k% - 1
294.                 IF k% = 0 THEN EXIT DO
295.                 matrice1(ieto%, z%, i% + 1, 1).i2 = matrice1(ieto%, z%, k%, 1).i1
296.                 k% = k% - 1
297.                 i% = i% + 1
298.             LOOP UNTIL k% = 0
299.     END SELECT
300. END SUB
301. '-----------------------------------------------------------------------------------------------------------------------------------------------------------------
302. SUB grafici (ore!, oreM!, QM!)
303.     'ore!: ore massime di analisi.
304.     'oreM!: ore corrispondenti alla portata massima.
305.     'QM!: portata massima.
306.  
307.     CONST dxy% = 30
308.     CONST L% = 800
309.     CONST H% = 600
310.  
311.     ' CONST LH! = L% / H% 'costante che, all'interno del comando WINDOW, di portare la scala delle ordinate uguale a quella delle ascisse.
312.     CONST rosso = _RGB32(255, 0, 0)
313.     CONST bianco = _RGB32(255, 255, 255)
314.     CONST giallo = _RGB32(255, 255, 0)
315.  
316.     DIM grafico&
317.     DIM fx! 'all'interno del comando WINDOW e VIEW permette di disegnare in pixel sull'asse x e nelle circonferenze.
318.     DIM fy! 'all'interno del comando WINDOW e VIEW permette di disegnare in pixel sull'asse y.
319.     DIM dx!, dy!
320.     DIM tacca%
321.     DIM k%
322.  
323.     grafico& = _NEWIMAGE(L%, H%, 32)
324.     fx! = L% / (ore! * 1.1) '=168.3349 'pixel/ora
325.     fy! = H% / (QM! * 1.1) '=23.3657 'pixel/mc
326.     dx! = dxy% / fx! '=0.1782162:dx! e' la lunghezza dentro WINDOWS per avere un numero di pixel=dxy%.
327.     dy! = dxy% / fy! '=1.283933
328.     tacca% = QM! \ 5
329.  
330.     _DEST grafico&
331.     SCREEN grafico&
332.     LINE (0, 0)-(L% - 1, H% - 1), rosso, B
333.     WINDOW (0, 0)-(ore! * 1.1 + dx!, QM! * 1.1 + dy!) 'scala ascisse e ordinate in base alle ore e alla portata, estendendole del 10%.
334.     LINE (dx!, QM! * 1.1)-(dx!, dy!): LINE -(ore! * 1.1, dy!), bianco
335.     PSET (dx!, QM! * 1.1!), bianco: DRAW "F20": PSET (dx!, QM! * 1.1), bianco: DRAW "G20" ': LOCATE 7, 65: PRINT "y"
336.     PSET (ore! * 1.1, dy!), bianco: DRAW "G20": PSET (ore! * 1.1, dy!), bianco: DRAW "H20" ': LOCATE 25, 112: PRINT "x"
337.     FOR z% = 1 TO 24
338.         CIRCLE (dx! + massimi1(visualizzaieto%, z%, 1).ore, dy! + massimi1(visualizzaieto%, z%, 1).portata), 3 / fx!, rosso
339.         PAINT (dx! + massimi1(visualizzaieto%, z%, 1).ore + 0.5 / fx!, dy! + massimi1(visualizzaieto%, z%, 1).portata + 0.5 / fx!), rosso
340.         IF z% = 1 THEN
341.             LINE (dx!, dy!)-(dx! + massimi1(visualizzaieto%, z%, 1).ore, dy! + massimi1(visualizzaieto%, z%, 1).portata), rosso
342.         ELSE
343.             LINE (dx! + massimi1(visualizzaieto%, z% - 1, 1).ore, dy! + massimi1(visualizzaieto%, z% - 1, 1).portata)-(dx! + massimi1(visualizzaieto%, z%, 1).ore, dy! + massimi1(visualizzaieto%, z%, 1).portata), rosso
344.         END IF
345.         IF z% = 24 THEN
346.             SELECT CASE visualizzaieto%
347.                 CASE IS = 1
348.                     LINE -(dx!, dy! + massimi1(visualizzaieto%, z%, 1).portata), bianco, , 61440
349.                     LINE (dx! + massimi1(visualizzaieto%, z%, 1).ore, dy! + massimi1(visualizzaieto%, z%, 1).portata)-(dx! + massimi1(visualizzaieto%, z%, 1).ore, dy!), bianco, , 61440
350.                 CASE IS = 2
351.                     LINE (dx! + massimi1(visualizzaieto%, z%, 1).ore, dy! + massimi1(visualizzaieto%, z%, 1).portata)-(dx! + massimi1(visualizzaieto%, z%, 1).ore, dy!), bianco, , 61440
352.             END SELECT
353.         END IF
354.     NEXT z%
355.     IF visualizzaieto% = 2 THEN
356.         LINE (dx!, dy! + QM!)-(dx! + oreM!, dy! + QM!), bianco, , 61440
357.         LINE -(dx! + oreM!, dy!), bianco, , 61440
358.     END IF
359.     ' _PRINTSTRING (L% \ 2, H% \ 2), "O" '_PRINTSTRING non usa le coordinate locali di VIEW, ma sempre quello dello SCREEN in uso. Per tal ragione e' stato necessario creare grafico&.
360.     SLEEP
361.     i% = 1
362.     k! = 0.98
363.     WHILE i% * tacca% <= QM!
364.         LINE (dx! - 5 / fx!, dy! + i% * tacca%)-(dx! + 5 / fx!, dy! + i% * tacca%): SLEEP
365.         _PRINTSTRING (dx! * fx! - 20, H% - 1 - (dy! + i% * tacca%) * fy! * k!), "-": SLEEP
366.         i% = i% + 1
367.     WEND
368.     WINDOW
369.     i% = 1
370.     '  k! = 1
371.     WHILE i% * tacca% <= QM!
372.         '   LINE (dxy% - 1 - 5, H% - 1 - (dy! + i% * tacca%) * fy! * k!)-(dxy% - 1 + 1000, H% - 1 - (dy! + i% * tacca%) * fy! * k!), giallo: SLEEP
373.         '   PRINT (dy! + i% * tacca%) * fy!: SLEEP
374.         _PRINTSTRING (dx! * fx! - 20, H% - 1 - (dy! + i% * tacca%) * fy! * k!), "--": SLEEP
375.         i% = i% + 1
376.     WEND
377.     ' k% = i% * tacca%
378.     'i% = 1
379.     'WHILE i% <= ore!
380.     '    _PRINTSTRING (dxy%, H% - 1 - (dxy% + i% * k%)), "-"
381.     'i% = i% + 1
382.     'WEND
383.     '  _PUTIMAGE ((_DESKTOPWIDTH - L%) \ 2, (_DESKTOPHEIGHT - H%) \ 2), grafico&, schermo&
384.     '   _FREEIMAGE grafico&
385. END SUB
386. '------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
387. SUB risultati ()
388.     FOR z% = 1 TO 24
389.         CLS
390.         PRINT "z%="; z%
391.         PRINT "dt="; dt!(visualizzaieto%, z%); ";"; dt!(visualizzaieto%, z%) * 60
392.         PRINT "visualizzaieto%="; visualizzaieto%
393.         PRINT "m%="; m%(visualizzaieto%, z%)
394.         PRINT "N="; coefficienti(z%).d / tc!
395.         PRINT "n="; coefficienti(z%).tSUta
396.         FOR i% = 1 TO 50
397.             '  _LIMIT limit
398.     PRINT _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).m)), _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).mdt)), _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).tminuti)), _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).h1)),_
399.      _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).i1)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).i2)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).DH)), _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).H2)),_
400.       _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).he)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).Dhe)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).MOCKUStSUta)),_
401.        _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).MOCKUSqSUqp)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).tmSUta)), _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).qmSUqp)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).qm))
402.         NEXT i%
403.         SLEEP
404.         CLS
405.         FOR n% = 1 TO 50
406.             '  _LIMIT limit
407.             LOCATE 1,
408.             FOR i% = 1 TO 50
409.                 LOCATE , n% * 4 - 3
410.                 PRINT matrice2!(visualizzaieto%, z%, i%, n%)
411.             NEXT i%
412.         NEXT n%
413.         SLEEP
414.     NEXT z%
415.     CLS
416.     FOR n% = 1 TO 24
417.         ' _LIMIT limit
418.         SELECT CASE n%
419.             CASE IS <= 16
420.                 FOR i% = 1 TO 50
421.                     LOCATE i%, n% * 15 - 14
422.                     PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).portata))
423.                     'PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).ore))
424.                     'PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).minuti))
425.  
426.                 NEXT i%
427.                 IF n% = 16 THEN SLEEP: CLS
428.             CASE IS > 16
429.                 FOR i% = 1 TO 50
430.                     LOCATE i%, (n% - 16) * 15 - 14
431.                     PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).portata))
432.                     'PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).ore))
433.                     'PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).minuti))
434.                 NEXT i%
435.         END SELECT
436.     NEXT n%
437.     SLEEP
438.     CLS
439.     FOR z% = 1 TO 24
440.         ' _LIMIT limit
441.         PRINT massimi1(visualizzaieto%, z%, 1).ore, massimi1(visualizzaieto%, z%, 1).minuti, massimi1(visualizzaieto%, z%, 1).portata
442.     NEXT z%
443.     SLEEP
444.     PRINT
445.     PRINT massimi2(visualizzaieto%).ore; massimi2(visualizzaieto%).portata
446. END SUB
447.  

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: OldMoses on February 20, 2021, 08:07:07 am

I don't know if I understand this correctly, but if you have a positioning problem between screen absolute coordinates and WINDOW defined coordinates, I've used the map! function to good effect. If you have a value that needs translation from an old  range {minRange!/maxRange!}, it will translate that value to the new range {newMinRange!/newMaxRange!}. The order depends upon which way one is going, from WINDOW to screen or from screen to WINDOW. Hope this helps

Code: QB64: [Select]

1. Function map! (value!, minRange!, maxRange!, newMinRange!, newMaxRange!)
2.  
3.     map! = ((value! - minRange!) / (maxRange! - minRange!)) * (newMaxRange! - newMinRange!) + newMinRange!
4.  
5. End Function 'map!
6.  

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bplus on February 20, 2021, 10:52:47 am

Aren't _PrintString coordinates automatically translated to Window coordinates like the graphics pset, line, circle..., if not, try same trick for translating Mouse Coordinates in a Window use PMap see Wiki.

Use Moses idea if printing outside the Window or returning screen back to normal.

BTW is your graph reversed? Both look flipped: up/down and left/right? Your screen bigger than my laptop, didn't see the axis arrows they were over in the storage shelf left of my laptop on my desk. :)

Ah TempodiBasic and or Kiara are recruiting ;-))

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bartok on February 20, 2021, 11:36:36 am

Quote from: OldMoses on February 20, 2021, 08:07:07 am

I don't know if I understand this correctly, but if you have a positioning problem between screen absolute coordinates and WINDOW defined coordinates, I've used the map! function to good effect. If you have a value that needs translation from an old  range {minRange!/maxRange!}, it will translate that value to the new range {newMinRange!/newMaxRange!}. The order depends upon which way one is going, from WINDOW to screen or from screen to WINDOW. Hope this helps

Code: QB64: [Select]

1. Function map! (value!, minRange!, maxRange!, newMinRange!, newMaxRange!)
2.  
3.     map! = ((value! - minRange!) / (maxRange! - minRange!)) * (newMaxRange! - newMinRange!) + newMinRange!
4.  
5. End Function 'map!
6.  

Hi and thanks for the tip, but it doesn't work with the damned _PRINTSTRING!

As you can see in the semplified code above (press any key to go ahead), outside WINDOW, your function works, because lines drawn in the WHILE cicle of lines 55 to 59, perfetly mach the notches drawn inside WINDOW in the WHILE circle of lines 40 to 44.

But, in the WHILE cicle of lines 62 to 66, no matches! I have also tried to add something like this:

_PRINTSTRING (55, H% - 1 - test!+K)

where K is an empirical and inexplicable value brutally added in order to have the match: no way.

Code: QB64: [Select]

1. CONST ore! = 6.723612
2. CONST oreM! = 6.723612
3. CONST QM! = 29.88063
4.  
5. CONST dxy% = 30
6. CONST L% = 800
7. CONST H% = 600
8.  
9. CONST rosso = _RGB32(255, 0, 0)
10. CONST bianco = _RGB32(255, 255, 255)
11. CONST giallo = _RGB32(255, 255, 0)
12.  
13. DIM grafico&
14. DIM fx! 'all'interno del comando WINDOW e VIEW permette di disegnare in pixel sull'asse x e nelle circonferenze.
15. DIM fy! 'all'interno del comando WINDOW e VIEW permette di disegnare in pixel sull'asse y.
16. DIM dx!, dy!
17. DIM test%
18. DIM tacca%
19. DIM i%
20.  
21. grafico& = _NEWIMAGE(L%, H%, 32)
22. fx! = L% / (ore! * 1.1) '=168.3349 'pixel/ora
23. fy! = H% / (QM! * 1.1) '=23.3657 'pixel/mc
24. dx! = dxy% / fx! '=0.1782162:dx! e' la lunghezza dentro WINDOWS per avere un numero di pixel=dxy%.
25. dy! = dxy% / fy! '=1.283933
26. tacca% = QM! \ 5
27.  
28.  
29.  
30. _DEST grafico&
31. SCREEN grafico&
32. LINE (0, 0)-(L% - 1, H% - 1), rosso, B
33. WINDOW (0, 0)-(ore! * 1.1 + dx!, QM! * 1.1 + dy!) 'scala ascisse e ordinate in base alle ore e alla portata, estendendole del 10%.
34. 'test! = map!(QM! \ 5 + dy!, 0, dy! + QM! * 1.1, 0, H% - 1)
35. LINE (dx!, QM! * 1.1)-(dx!, dy!): LINE -(ore! * 1.1, dy!), bianco
36. PSET (dx!, QM! * 1.1!), bianco: DRAW "F20": PSET (dx!, QM! * 1.1), bianco: DRAW "G20" ': LOCATE 7, 65: PRINT "y"
37. PSET (ore! * 1.1, dy!), bianco: DRAW "G20": PSET (ore! * 1.1, dy!), bianco: DRAW "H20" ': LOCATE 25, 112: PRINT "x"
38. SLEEP
39. i% = 1
40. WHILE i% * tacca% <= QM!
41.     LINE (dx! - 5 / fx!, dy! + i% * tacca%)-(dx! + 5 / fx!, dy! + i% * tacca%): SLEEP
42.     '  _PRINTSTRING (dx! * fx! - 20, H% - 1 - (dy! + i% * tacca%) * fy! * k!), "-": SLEEP
43.     i% = i% + 1
44. WEND
45.  
46. WINDOW
47.  
48. 'i% = 1
49. 'WHILE i% * tacca% <= QM!
50. '_PRINTSTRING (dx! * fx! - 20, H% - 1 - (dy! + i% * tacca%) * fy!), "--": SLEEP
51. 'i% = i% + 1
52. 'WEND
53.  
54. i% = 1
55. WHILE i% * tacca% <= QM!
56.     test! = map!(i% * QM! \ 5 + dy!, 0, dy! + QM! * 1.1, 0, H% - 1)
57.     LINE (35, H% - 1 - test!)-(50, H% - 1 - test!): PRINT test!: SLEEP '------|
58.     i% = i% + 1 '                                                             |
59. WEND '                                                                        |
60. '                                                                             |
61. i% = 1 '                                                                      |---Same (H%-1-test!), different output!
62. WHILE i% * tacca% <= QM! '                                                    |
63.     test! = map!(i% * QM! \ 5 + dy!, 0, dy! + QM! * 1.1, 0, H% - 1) '         |
64.     _PRINTSTRING (55, H% - 1 - test!), "--": PRINT test!: SLEEP '-------------|
65.     i% = i% + 1
66. WEND
67. '---------------------------------------------------------------------------------------
68. FUNCTION map! (value!, minWINDOW!, maxWINDOW!, minPIXEL%, maxPIXEL%)
69.     map! = ((value! - minWINDOW!) / (maxWINDOW! - minWINDOW!)) * (maxPIXEL% - minPIXEL%) + minPIXEL%
70. END FUNCTION
71.  

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bartok on February 20, 2021, 11:42:07 am

Quote from: bplus on February 20, 2021, 10:52:47 am

Your screen bigger than my laptop, didn't see the axis arrows they were over in the storage shelf left of my laptop on my desk. :)

Ah TempodiBasic and or Kiara are recruiting ;-))

I have changed both of the 2 codes with 800x600.

In the second post, I have semplified the code, because the core of my problem is in the "graphs" subroutine of the first code I have posted. However, the program is finalized to calculate the maximum flow of a given drainage basin for a given return period, depending on the duration of the rain.

Quote

Aren't _PrintString coordinates automatically translated to Window coordinates like the graphics pset, line, circle.

no, unfortunately they aren't...

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: SMcNeill on February 20, 2021, 12:04:58 pm

You’re in SCREEN 0, where _PRINTSTRING is nothing more than a shortcut for LOCATE x, y: PRINT.

Here’s the relevant c-snippet for you:

Code: C++: [Select]

1.        void sub__printstring(float x,float y,qbs* text,int32 i,int32 passed){
2.             .... STUFF
3.             if (im->text){
4.                 int oldx = func_pos(0), oldy = func_csrlin();
5.                 qbg_sub_locate(y, x, 0, 0, 0, 3);
6.                 qbs_print(text, 0);
7.                 qbg_sub_locate(oldy, oldx, 0, 0, 0, 3);
8.                 goto printstring_exit;
9.             }
10.             ..... STUFF

If you’re in a text screen, simply use _PRINTSTRING as a LOCATE, PRINT, then restore the previous LOCATE position.

If _PRINTSTRING doesn’t work, just go back to LOCATE: PRINT

Edit: Nevermind...  I see your SCREEN statement down on line 70something...  You’re not trying anything in a text screen after all. 

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bartok on February 20, 2021, 01:06:24 pm

SMcNeill,
yes, in fact I am in graphic mode.

It is not possible it is so complicated to put a text in a graph!:D

(I try to avoid the use of GOTO).

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bplus on February 20, 2021, 01:23:25 pm

Found my old test code.

YES you MUST use PMAP for x, y conversion for both MOUSE and for _PRINTSTRING

demo

Code: QB64: [Select]

1. OPTION _EXPLICIT
2. _TITLE "Test Using Window for Cartesian Coordinate System x = -100 to 100, y = -100 to 100"
3. 'It seems all graphics commands convert x, y to window system BUT...
4. 'Mouse and _PrintString do not convert! BUT...
5. 'Somebody setup a PMAP system that will do these conversions both to grapghics coordinate system and back
6. ' convert mouse qbx, qby to graphics x, y use
7. ' PMAP(qbx, 2), PMAP(qby, 3)
8. ' And to print something at the graphics location use
9. ' PMAP(graphicsX, 0), PMAP(graphicsY, 1) to _PRINTSTRING to a graphics location
10.  
11. CONST xmax = 700, ymax = 700, red = &HFFFF0000, grn = &HFF009900, blu = &HFF0000AA
12. CONST tlx = -100, tly = 100, brx = 100, bry = -100 ' Cartesian Coordinate System corners for WINDOW command
13. SCREEN _NEWIMAGE(xmax, ymax, 32)
14. _SCREENMOVE 300, 40
15.  
16. DIM mx, my, a, r, s$
17.  
18. WINDOW (tlx, tly)-(brx, bry) ' converts QB csreen to Cartesian Coodianate System
19. WHILE _KEYDOWN(27) = 0
20.     CLS
21.     WHILE _MOUSEINPUT: WEND
22.     mx = INT(PMAP(_MOUSEX, 2)): my = INT(PMAP(_MOUSEY, 3)) 'covert mouse to Cart Cood using PMAP
23.     a = atan360(my, mx) 'mouse angle to origin
24.     r = INT((mx * mx + my * my) ^ .5)
25.     LINE (0, 0)-(mx, my), grn 'hypotenuse = radius
26.     LINE (0, 0)-(mx, 0), blu ' r*cos
27.     LINE (mx, 0)-(mx, my), red ' r*sin
28.     arc 0, 0, 10, 0, _R2D(a), &HFFFFFF00 'yellow arc
29.     s$ = "Mouse (" + _TRIM$(STR$(mx)) + ", " + _TRIM$(STR$(my)) + "), Radius: "
30.     s$ = s$ + _TRIM$(STR$(r)) + ", Degrees: " + _TRIM$(STR$(_R2D(a) \ 1))
31.     _PRINTSTRING (PMAP(0, 0), PMAP(-1, 1)), s$
32.     _PRINTSTRING (PMAP(-95, 0), PMAP(95, 1)), "COS = x/r = " + _TRIM$(STR$(mx / r))
33.     _PRINTSTRING (PMAP(-95, 0), PMAP(89, 1)), "SIN = y/r = " + _TRIM$(STR$(my / r))
34.     IF mx <> 0 THEN _PRINTSTRING (PMAP(-95, 0), PMAP(83, 1)), "TAN = y/x = " + _TRIM$(STR$(my / mx)) 'EDIT
35.     _DISPLAY
36.     _LIMIT 60
37. WEND
38.  
39. ' A regular graphic draw still works correctly in a Cartesian Window System!
40. 'note: degrees start due East = 0 and go clockwise: South = 90, West = 180, North = 270
41. SUB arc (x, y, r, degStart, degStop, c~&)
42.  
43.     DIM rs, re, a, ax, ay, lastx, lasty
44.     'x, y origin, r = radius, c = color
45.  
46.     'degStart is first angle clockwise from due East = 0 degrees
47.     ' arc will start drawing there and clockwise until degStop angle reached, unless in Cartesain Window
48.  
49.     IF degStop < degStart THEN
50.         arc x, y, r, degStart, 360, c~&
51.         arc x, y, r, 0, degStop, c~&
52.     ELSE
53.         rs = _D2R(degStart): re = _D2R(degStop)
54.         FOR a = rs TO re STEP 1 / (7 * r)
55.             ax = x + r * COS(a)
56.             ay = y + r * SIN(a)
57.             IF a <> rs THEN LINE (lastx, lasty)-(ax, ay), c~&
58.             lastx = ax: lasty = ay
59.         NEXT
60.     END IF
61. END SUB
62.  
63. 'need to match degrees going around clockwise up to full circle
64. FUNCTION atan360 (y, x)
65.     DIM test
66.     test = _ATAN2(y, x)
67.     IF test < 0 THEN atan360 = _PI(2) + test ELSE atan360 = test
68. END FUNCTION
69.  
70.  

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bartok on February 20, 2021, 01:46:08 pm

bplus,
thank you very much! I did not have time yet to try PMAP after your message, but it was already my purpose to try it. I didn't know the command PMAP before today.

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bartok on February 21, 2021, 04:01:32 am

Thank's you bplus, with PMAP it works better!

But, as you can see in the code above on line 43, I have to add "-7" in order to have the strings "--" fitting. I don't know why! In line 41 PMAP gets the y screen coordinates of WINDOW y coordinates of the line drawn in line 39. So the string "--" printed with _PRINTSTRING in line 43 should match using y calculated by PMAP without this "-7"...

Code: QB64: [Select]

1. CONST ore! = 6.723612
2. CONST oreM! = 6.723612
3. CONST QM! = 29.88063
4.  
5. CONST dxy% = 30
6. CONST rosso = _RGB32(255, 0, 0)
7. CONST bianco = _RGB32(255, 255, 255)
8. CONST giallo = _RGB32(255, 255, 0)
9.  
10. DIM L%, H%
11. DIM grafico&
12. DIM fx! 'all'interno del comando WINDOW e VIEW permette di disegnare in pixel sull'asse x e nelle circonferenze.
13. DIM fy! 'all'interno del comando WINDOW e VIEW permette di disegnare in pixel sull'asse y.
14. DIM dx!, dy!
15. DIM tacca%
16. DIM x%, y%
17.  
18. H% = _DESKTOPHEIGHT * 2 / 3
19. L% = H% * 1.62
20.  
21. grafico& = _NEWIMAGE(L%, H%, 32)
22. fx! = L% / (ore! * 1.1) '=168.3349 'pixel/ora
23. fy! = H% / (QM! * 1.1) '=23.3657 'pixel/mc
24. dx! = dxy% / fx! '=0.1782162:dx! e' la lunghezza dentro WINDOWS per avere un numero di pixel=dxy%.
25. dy! = dxy% / fy! '=1.283933
26. tacca% = QM! \ 5
27.  
28. _DEST grafico&
29. SCREEN grafico&
30. LINE (0, 0)-(L% - 1, H% - 1), rosso, B
31. WINDOW (0, 0)-(ore! * 1.1 + dx!, QM! * 1.1 + dy!) 'scala ascisse e ordinate in base alle ore e alla portata, estendendole del 10%.
32. LINE (dx!, QM! * 1.1)-(dx!, dy!): LINE -(ore! * 1.1, dy!), bianco
33. PSET (dx!, QM! * 1.1!), bianco: DRAW "F20": PSET (dx!, QM! * 1.1), bianco: DRAW "G20" ': LOCATE 7, 65: PRINT "y"
34. PSET (ore! * 1.1, dy!), bianco: DRAW "G20": PSET (ore! * 1.1, dy!), bianco: DRAW "H20" ': LOCATE 25, 112: PRINT "x"
35. SLEEP
36. i% = 1
37. i% = 1
38. WHILE i% * tacca% <= QM!
39.     LINE (dx! - 5 / fx!, dy! + i% * tacca%)-(dx! + 5 / fx!, dy! + i% * tacca%): SLEEP
40.     x% = PMAP(dx! - 5 / fx!, 0)
41.     y% = PMAP(dy! + i% * tacca%, 1)
42.     PRINT x%, y%: SLEEP
43.     _PRINTSTRING (x% + 10, y% - 7), "--": SLEEP '----------------------------------------------------> WHY 7??
44.     i% = i% + 1
45. WEND
46. WINDOW
47.  

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bartok on February 21, 2021, 05:38:33 am

Ok, solved!

PMAP works perfectly with _PRINTSTRING. In this case "--" didn't match the line because the line is drawn at the top and "--" in the middle (_--)-->up side down.
But this is not the solution without using PMAP, because in the preview tries the "--"  changed the distance with the line. So PMAP is necessary.

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: OldMoses on February 21, 2021, 09:20:12 am

I find this interesting as I have not, thus far, investigated PMAP.

It would seem that it is essentially FUNCTION map! that I posted earlier, only optimized for screen ops, and much less typing.

You couldn't use it, like you can map!, for general computations, but it's much easier for graphics. You don't have to do that much range figuring.

Having learned about a new trick, I'll have to do some experiments now...;)

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bplus on February 21, 2021, 01:23:56 pm

Yeah WINDOW is the way to go especially for math apps that you want standard 4 quads upper right #1 and going Counter-clockwise 2, 3, 4 with angle increases.

Anyway with that command (set correctly of course) angles Increase going counter-clockwise and y up is up and x up is right! The way you learned it in school before Basic ruined your mind by flipping it upside down and making angles increase clockwise ;-))

Trouble is I don't know if anyone uses it enough to report bugs when working with all the new extensions installed with QB64 ie how z axis might work with GL 3D?

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bartok on February 21, 2021, 05:21:58 pm

Quote from: bplus on February 21, 2021, 01:23:56 pm

Trouble is I don't know if anyone uses it enough to report bugs when working with all the new extensions installed with QB64 ie how z axis might work with GL 3D?

I am a beginner, so I think that my level is not enough hight to report bugs, but maybe it wouldn't be bad if _PRINTSTRING worked also inside WINDOW. But perhaps there is a reason if it doesn't. For the 3D, maybe some games programmer. I don't think my graphics needs will go beyond 2D.

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bplus on February 21, 2021, 06:45:55 pm

Code: QB64: [Select]

1. OPTION _EXPLICIT
2. _TITLE "But _PrintString Works Fine in a Winodow!"
3. CONST xmax = 512, ymax = 400
4. CONST Top_Left_X = -100, Top_Left_Y = 100, Bottom_Right_X = 100, Bottom_Right_Y = -100 ' Cartesian Coordinate System corners for WINDOW command
5. SCREEN _NEWIMAGE(xmax, ymax, 32)
6. _SCREENMOVE 300, 40
7. WINDOW (Top_Left_X, Top_Left_Y)-(Bottom_Right_X, Bottom_Right_Y) ' converts QB csreen to Cartesian Coodianate System
8. _PRINTSTRING (PMAP(-95, 0), PMAP(60, 1)), "I am a beginner, so I think that my level is not enough high"
9. _PRINTSTRING (PMAP(-95, 0), PMAP(30, 1)), "to report bugs, but maybe it wouldn't be bad if _PRINTSTRING"
10. _PRINTSTRING (PMAP(-95, 0), PMAP(0, 1)), "worked also inside WINDOW. But perhaps there is a reason if "
11. _PRINTSTRING (PMAP(-95, 0), PMAP(-30, 1)), "it doesn't. For the 3D, maybe some games programmer. I don't"
12. _PRINTSTRING (PMAP(-95, 0), PMAP(-60, 1)), " think my graphics needs will go beyond 2D."
13. SLEEP
14.  
15.  
16.  

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bartok on February 22, 2021, 02:38:42 am

You are right. I said that badly. I meant that it would be not bad if _PRINTSTRING, inside WINDOW, worked with the local coordinates of WINDOW, without the convertion with PMAP.

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bplus on February 22, 2021, 03:51:23 am

Quote from: bartok on February 22, 2021, 02:38:42 am

You are right. I said that badly. I meant that it would be not bad if _PRINTSTRING, inside WINDOW, worked with the local coordinates of WINDOW, without the convertion with PMAP.

Oh you mean like pset, circle and line x, y's, yeah that would be nice. But as I was trying to say, to my knowledge, Window is not used so much to keep it updated with newer stuff like mouse and _printstring.

But at least PMap is a resource that works.

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bartok on February 22, 2021, 07:46:28 am

Quote from: bplus on February 22, 2021, 03:51:23 am

Oh you mean like pset, circle and line x, y's, yeah that would be nice. But as I was trying to say, to my knowledge, Window is not used so much to keep it updated with newer stuff like mouse and _printstring.

But at least PMap is a resource that works.

why is not used so much! It's seems fundamental to me.

However, the result for now is not bad.

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bplus on February 22, 2021, 10:41:06 am

@bartok looks good man!

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bartok on February 22, 2021, 10:46:43 am

Quote from: bplus on February 22, 2021, 10:41:06 am

@bartok looks good man!

thank you!

I have an other problem, that it seems inexplicable.

Launch the program.
The purpose is to draw a dashed line from X=3 going up until it find a yellow pixel. As we can see, it works. But, if we deactivate line 403, it doesn't and it goes in loop, as if it doesn't match the yellow pixel. I don't understand, considering that in both cases I am working on grafico& (line 400). I don't want to see grafico&: I just want to work in background on it and then _PUTIMAGE it (line 410) on the visualized screen.

Code: QB64: [Select]

1. OPTION BASE 1
2.  
3. mockus: 'dati dell'idrogramma unitario di Mockus
4. DATA 0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2,2.1,2.2,2.3,2.4,2.5,2.6,2.7,2.8,2.9,3,3.1,3.2,3.3,3.4,3.5,3.6,3.7,3.8,3.9,4,4.1,4.2,4.3,4.4,4.5,4.6,4.7,4.8,4.9,5: 't/ta
5. DATA 0.030,0.100,0.190,0.310,0.470,0.660,0.820,0.930,0.990,1.000,0.990,0.930,0.860,0.780,0.680,0.560,0.460,0.390,0.330,0.280,0.244,0.207,0.177,0.147,0.127,0.107,0.092,0.077,0.066,0.055,0.048,0.040,0.035,0.029,0.025,0.021,0.018,0.015,0.013,0.011,0.009,0.008,0.007,0.006,0.006,0.005,0.004,0.003,0.001,0.000: ' t/tp
6.  
7. coefficienti:
8. DATA 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,18,20,22,24,26,28,30,32: 'N volte tc!.
9. DATA 0.1,0.1,0.1,0.2,0.3,0.2,0.3,0.3,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,1,1,1,1,1,1,1,1: 'discretizzazione dei coeff. di Mockus.
10.  
11. TYPE mockus
12.     MOCKUStSUta AS SINGLE
13.     MOCKUSqSUqp AS SINGLE
14. END TYPE
15.  
16. TYPE matrice1
17.     m AS INTEGER '[numero naturale] conta il numero di interazioni fino al raggiungimento di "d".
18.     mdt AS SINGLE '[ore] conta la progessione delle ore fino al raggiungimento di "d".
19.     tminuti AS SINGLE '[minuti] trasforma mdt in minuti.
20.     h1 AS SINGLE '[mm] valore h della curva di possibilità climatica.
21.     i1 AS SINGLE '[mm/ora] intensità di pioggia.
22.     i2 AS SINGLE '[mm/ora] intensità di pioggia secondo lo ietogramma Chicago o costante.
23.     DH AS SINGLE '[mm] pioggia in ogni istante temporale dt.
24.     H2 AS SINGLE '[mm] pioggia cumulata ad ogni istante temporale.
25.     he AS SINGLE '[mm] piggia efficace cumulata ad ogni istante temporale.
26.     Dhe AS SINGLE '[mm] pioggia efficace ad ogni istante temporale.
27.     MOCKUStSUta AS SINGLE
28.     MOCKUSqSUqp AS SINGLE
29.     tmSUta AS SINGLE '[-].
30.     qmSUqp AS SINGLE '[-] il valore di qm/qp corrispondente a tm/ta e' = al valore di q/qp che corrisponde a t/ta=tm/ta.
31.     qm AS SINGLE '[mc/s].
32. END TYPE
33.  
34. TYPE coefficienti
35.     d AS SINGLE
36.     tSUta AS SINGLE
37. END TYPE
38.  
39. TYPE idrogramma
40.     ore AS SINGLE
41.     minuti AS SINGLE
42.     portata AS SINGLE
43. END TYPE
44.  
45. CONST limit = 200 '50
46.  
47. DIM SHARED schermo&
48.  
49. DIM SHARED mockus(50) AS mockus
50. DIM SHARED dt!(2, 24)
51. DIM SHARED matrice1(2, 24, 50, 1) AS matrice1
52. DIM SHARED coefficienti(24) AS coefficienti
53. DIM SHARED matrice2!(2, 24, 50, 50)
54. DIM SHARED idrogramma(2, 24, 50, 1) AS idrogramma
55. DIM SHARED massimi1(2, 24, 1) AS idrogramma
56. DIM SHARED massimi2(2) AS idrogramma
57. DIM SHARED i%, n%, z%, ieto%
58. DIM SHARED m%(2, 24)
59. DIM SHARED S2! '[mm] contenuto idrico massimo del terreno.
60. DIM SHARED Ia! '[mm] perdite iniziali.
61. DIM SHARED tc! '[ore] tempo di corrivazione
62. DIM SHARED ta! '[ore] tempo di picco.
63. DIM SHARED qp! '[mc/s] portata al colmo dell'idrogramma unitario.
64. DIM SHARED k! '[-] coeff. della curva di possibilità climatica tado dal sito della Regione.
65. DIM SHARED a2! '[mm/d^n] coeff. della curva di possibilità climatica.
66. DIM SHARED n1! '[-] coeff. della curva di possibilità climatica.
67. DIM SHARED visualizzaieto%
68. DIM A1! '[mq] superficie bacino idrografico.
69. DIM L! '[m] lunghezza asta principale.
70. DIM HmaxL! '[m] punto più alto dell'astra principale.
71. DIM HminL! '[m] punto più basso dell'asta principale.
72. DIM s1! '[-] pendenza media del bacino non espressa in percentuale ma come 0,x.
73. DIM CNII! '[-] tabellato.
74. DIM CNIII '[-] discende da CNII.
75. DIM tl! '[ore] tempo di ritardo.
76. DIM keypress$
77.  
78. schermo& = _NEWIMAGE(_DESKTOPWIDTH, _DESKTOPHEIGHT, 32)
79. SCREEN schermo&
80. _SCREENMOVE 0, 0
81.  
82. A1! = 2641902
83. k! = 2.34
84. a2! = 14.6653
85. n1! = 0.45596
86. L! = 3300
87. HmaxL! = 1981
88. HminL! = 1062
89. s1! = 0.415
90. CNII! = 70
91.  
92. 'INPUT A1!
93. 'INPUT k!
94. 'INPUT a2!
95. 'INPUT n1!
96. 'INPUT L!
97. 'INPUT HmaxL!
98. 'INPUT HminL!
99. 'INPUT s1!
100. 'INPUT CNII!
101. 'INPUT d1!
102.  
103. CNIII! = (23 * CNII!) / (10 + 0.13 * CNII!)
104. tl! = 0.342 * ((L! / 1000) ^ 0.8 / (100 * s1!) ^ 0.5) * (1000 / CNIII! - 9) ^ 0.7 'formula di Mockus.
105. S2! = 25.4 * (1000 / CNIII! - 10)
106. Ia! = 0.1 * S2! 'coeff.=0.03-0.2.
107. tc! = tl! / 0.6
108. ta! = tl! / 0.9
109. qp! = 0.208 * ((A1! / 1000000) / ta!)
110.  
111. RESTORE mockus
112. FOR i% = 1 TO 100
113.     IF i% <= 50 THEN READ mockus(i%).MOCKUStSUta
114.     IF i% > 50 THEN READ mockus(i% - 50).MOCKUSqSUqp
115. NEXT i%
116.  
117. RESTORE coefficienti
118. FOR i% = 1 TO 48
119.     IF i% <= 24 THEN READ coefficienti(i%).d: coefficienti(i%).d = coefficienti(i%).d * tc!
120.     IF i% > 24 THEN READ coefficienti(i% - 24).tSUta
121. NEXT i%
122.  
123. CLS
124. 'PRINT "Ietogramma Chicago/costante: 1/2> "
125. 'DO
126. '    _LIMIT 60 'limit
127. 'keypress$ = INKEY$
128. 'LOOP UNTIL keypress$ = "1" OR keypress$ = "2"
129. 'visualizzaieto% = VAL(keypress$)
130. visualizzaieto% = 1
131.  
132. FOR ieto% = 1 TO 2
133.     CALL creamatrici
134. NEXT ieto%
135.  
136. CLS
137. CALL grafici(massimi1(visualizzaieto%, 24, 1).ore, massimi2(visualizzaieto%).ore, massimi2(visualizzaieto%).portata) 'valori passati per disegnare il grafico.
138. SLEEP
139. SYSTEM
140. 'CALL risultati
141.  
142. '-----------------------------------------matrice1---------------------------------------------------------------------------------------------------------------------------------
143. SUB creamatrici ()
144.     z% = 0
145.     DO
146.         z% = z% + 1
147.         dt!(ieto%, z%) = coefficienti(z%).tSUta * ta!
148.         m%(ieto%, z%) = coefficienti(z%).d / (coefficienti(z%).tSUta * ta!) 'determina il numero di passi temporali necessari ad arrivare allla durata d considerata con il paso temporale dt!.
149.         FOR i% = 1 TO 50
150.             matrice1(ieto%, z%, i%, 1).MOCKUStSUta = mockus(i%).MOCKUStSUta
151.             matrice1(ieto%, z%, i%, 1).MOCKUSqSUqp = mockus(i%).MOCKUSqSUqp
152.             matrice1(ieto%, z%, i%, 1).m = i%
153.             matrice1(ieto%, z%, i%, 1).mdt = matrice1(ieto%, z%, i%, 1).m * dt!(ieto%, z%)
154.             matrice1(ieto%, z%, i%, 1).tminuti = matrice1(ieto%, z%, i%, 1).mdt * 60
155.             SELECT CASE (_ROUND(10 * (i% * coefficienti(z%).tSUta))) / 10 'crea il vettore dell'idrogramma unitario di Mockus con i multipli del coeff. z% considerato (0.1,0.2,0.3,0.5,1).
156.                 CASE IS <= 5 'valore massimo di t/ta di Mockus
157.                     matrice1(ieto%, z%, i%, 1).tmSUta = (_ROUND(10 * (i% * coefficienti(z%).tSUta))) / 10 '=matrice1(i%).m *  matrice1(i%).mdt / ta!
158.                 CASE ELSE
159.                     matrice1(ieto%, z%, i%, 1).tmSUta = 0
160.             END SELECT
161.         NEXT i%
162.         FOR i% = 1 TO 50
163.             n% = 1
164.             DO 'crea il vettore dell'idrogramma unitario Mockus prendendo i valori corrispondenti di ai coefficienti.
165.                 IF matrice1(ieto%, z%, i%, 1).tmSUta = matrice1(ieto%, z%, n%, 1).MOCKUStSUta THEN
166.                     matrice1(ieto%, z%, i%, 1).qmSUqp = matrice1(ieto%, z%, n%, 1).MOCKUSqSUqp
167.                     EXIT DO
168.                 ELSE
169.                     n% = n% + 1
170.                 END IF
171.             LOOP UNTIL n% = 50
172.             matrice1(ieto%, z%, i%, 1).qm = matrice1(ieto%, z%, i%, 1).qmSUqp * qp!
173.         NEXT i%
174.         FOR i% = 1 TO m%(ieto%, z%)
175.             matrice1(ieto%, z%, i%, 1).h1 = k! * a2! * matrice1(ieto%, z%, i%, 1).mdt ^ n1!
176.             IF i% = 1 THEN
177.                 matrice1(ieto%, z%, i%, 1).i1 = matrice1(ieto%, z%, i%, 1).h1 / dt!(ieto%, z%)
178.             ELSE
179.                 matrice1(ieto%, z%, i%, 1).i1 = (matrice1(ieto%, z%, i%, 1).h1 - matrice1(ieto%, z%, i% - 1, 1).h1) / dt!(ieto%, z%)
180.             END IF
181.         NEXT i%
182.         SELECT CASE ieto%
183.             CASE 1 'Chicago.
184.                 CALL chicago
185.             CASE 2 'costante
186.                 FOR i% = 1 TO m%(ieto%, z%)
187.                     matrice1(ieto%, z%, i%, 1).i2 = matrice1(ieto%, z%, m%(ieto%, z%), 1).h1 / matrice1(ieto%, z%, m%(ieto%, z%), 1).mdt
188.                 NEXT i%
189.         END SELECT
190.         FOR i% = 1 TO m%(ieto%, z%)
191.             matrice1(ieto%, z%, i%, 1).DH = matrice1(ieto%, z%, i%, 1).i2 * dt!(ieto%, z%)
192.             IF i% = 1 THEN
193.                 matrice1(ieto%, z%, i%, 1).H2 = matrice1(ieto%, z%, i%, 1).DH
194.             ELSE
195.                 matrice1(ieto%, z%, i%, 1).H2 = matrice1(ieto%, z%, i%, 1).DH + matrice1(ieto%, z%, i% - 1, 1).H2
196.             END IF
197.             SELECT CASE matrice1(ieto%, z%, i%, 1).H2 - Ia!
198.                 CASE IS < 0
199.                     matrice1(ieto%, z%, i%, 1).he = 0
200.                 CASE IS >= 0
201.                     matrice1(ieto%, z%, i%, 1).he = (matrice1(ieto%, z%, i%, 1).H2 - Ia!) ^ 2 / (matrice1(ieto%, z%, i%, 1).H2 - Ia! + S2!)
202.             END SELECT
203.             IF i% = 1 THEN
204.                 matrice1(ieto%, z%, i%, 1).Dhe = matrice1(ieto%, z%, i%, 1).he
205.             ELSE
206.                 matrice1(ieto%, z%, i%, 1).Dhe = matrice1(ieto%, z%, i%, 1).he - matrice1(ieto%, z%, i% - 1, 1).he
207.             END IF
208.         NEXT i%
209.  
210.         '-------------------------------matrice2-------------------------------------------------------------------------------------------------------------------------
211.  
212.         FOR n% = 1 TO 50 'colonna
213.             FOR i% = n% TO 50 ' riga
214.                 matrice2!(ieto%, z%, i%, n%) = matrice1(ieto%, z%, i% - n% + 1, 1).Dhe * matrice1(ieto%, z%, n%, 1).qm
215.             NEXT i%
216.         NEXT n%
217.         '-------------------------------idrogramma------------------------------------------------------------------------------------------------------------------
218.         FOR i% = 1 TO 50 'riga
219.             FOR n% = 1 TO 50 'colonna
220.                 idrogramma(ieto%, z%, i%, 1).portata = idrogramma(ieto%, z%, i%, 1).portata + matrice2!(ieto%, z%, i%, n%)
221.  
222.             NEXT n%
223.             idrogramma(ieto%, z%, i%, 1).ore = matrice1(ieto%, z%, i%, 1).mdt
224.             idrogramma(ieto%, z%, i%, 1).minuti = matrice1(ieto%, z%, i%, 1).tminuti
225.         NEXT i%
226.         '-------------------------------Massimi-----------------------------------------------------
227.         i% = 1
228.         'massimi1(ieto%, z%, 1).portata = idrogramma(ieto%, z%, 1, 1).portata
229.         DO
230.             IF massimi1(ieto%, z%, 1).portata < idrogramma(ieto%, z%, i%, 1).portata THEN
231.                 massimi1(ieto%, z%, 1).portata = idrogramma(ieto%, z%, i%, 1).portata
232.                 massimi1(ieto%, z%, 1).ore = idrogramma(ieto%, z%, i%, 1).ore
233.                 massimi1(ieto%, z%, 1).minuti = idrogramma(ieto%, z%, i%, 1).minuti
234.             END IF
235.             i% = i% + 1
236.         LOOP UNTIL i% = 50
237.         ' IF massimi2(ieto%).ore < massimi1(visualizzaieto%, z%, 1).ore THEN massimi2(ieto%).ore = massimi1(visualizzaieto%, z%, 1).ore
238.         IF massimi2(ieto%).portata < massimi1(visualizzaieto%, z%, 1).portata THEN
239.             massimi2(ieto%).portata = massimi1(visualizzaieto%, z%, 1).portata
240.             massimi2(ieto%).ore = massimi1(visualizzaieto%, z%, 1).ore
241.         END IF
242.     LOOP UNTIL z% = UBOUND(coefficienti)
243. END SUB
244. '---------------------------------------------------------------------------------------------------------------------------------------------------
245. 'Questa subroutine crea gli ietogrammi chicago per ogni iterazione z%.
246. 'Ad ogni nuovo dt! pone l'ultimo valore di pioggia come ultimo, il penultimo come primo, il terzultimo come penultimo, il quartultimo come secondo e così via.
247. 'Quando vi è un valore di z% che ha un dt! non nuovo, cerca la precedente iterazione z% con lo stesso dt! e copia i suoi valori dello ietogramma chicago in quello dell'iterazione z% in corso. I valori rimanenti sono riempiti
248. 'procedendo come prima.
249. 'Quindi per esempio l'iterazione z%=1 corrisponde a d=1tc con passo temporale dt!. L'iterazione z%=2 corrisponde a d=2tc con il medesimo passo temporale. Per cui lo ietogramma chicago per d=2tc e' dato dallo stesso ietogramma
250. 'chicago che si ha per d=1tc, con l'aggiunta, ai suoi estremi, dei nuovi valori. Il meccanismo si ripete per z%=3, che corrisponde a d=3tc, con il medesimo dt!. Per z%=4, corrispondente a d=4tc, invece, il passo temporale dt!
251. 'cambia.
252. 'Quanto detto e' necessario per evitare piccole discrepanze che potrebbero comportare il fatto che alcuni massimi risultino localmente poco piu' bassi del precedente, anche quando il trend generale e' chiaramente di aumento dei
253. 'massimi all'aumentare della durata della pioggia. In questo modo si fa l'ipotesi che all'aumentare della durata della pioggia, discretizzata con uno stesso passo temporale dt!, ogni ietogramma contenga a sua volta il precedente.
254. 'In tal modo, se si condiera per esempio d=2tc, si assume che l'altezza di pioggia corrispondente a 1tc, segua lo ietogramma relativo a d=1tc, cosi' che l'altezza di pioggia rimamente sia eclusivamente dovuta ai nuovi valori.
255. SUB chicago ()
256.     DIM test%
257.     DIM k%
258.     n% = 1
259.     DO
260.         IF z% = 1 THEN test% = 0: EXIT DO 'permette il calcolo dello ietogramma chicago della prima iterazione. test%=0
261.         IF dt!(ieto%, z%) = dt!(ieto%, z% - n%) THEN ' se trova un'iterazione z% precedente con lo stesso dt! di quella in corso, esegue i comandi successivi
262.             test% = 1
263.             k% = m%(ieto%, z%)
264.             i% = 0
265.             DO 'inserisce i nuovi valori dello ietogramma chicago dell'iterazione z% in corso.
266.                 matrice1(ieto%, z%, m%(ieto%, z%) - i%, 1).i2 = matrice1(ieto%, z%, k%, 1).i1
267.                 k% = k% - 1
268.                 IF k% = m%(ieto%, z% - n%) THEN EXIT DO 'se ha completato l'inserimento dei nuovi valori, passa oltre.
269.                 matrice1(ieto%, z%, i% + 1, 1).i2 = matrice1(ieto%, z%, k%, 1).i1
270.                 k% = k% - 1
271.                 i% = i% + 1
272.             LOOP UNTIL k% = m%(ieto%, z% - n%)
273.             IF matrice1(ieto%, z%, i%, 1).i2 = 0 THEN i% = i% - 1
274.             k% = 1
275.             DO 'inserisce lo ietogramma chicago dell'iterazione z% precedente con lo stesso dt!, nello ietogramma dell'iterazione z% in corso.
276.                 i% = i% + 1
277.                 matrice1(ieto%, z%, i%, 1).i2 = matrice1(ieto%, z% - n%, k%, 1).i2
278.                 k% = k% + 1
279.             LOOP UNTIL k% > m%(ieto%, z% - n%)
280.             EXIT DO 'siccome ha completato lo ietogramma dell'iterazione z% in corso, non serve incrementare n%, quindi passa oltre. test%=1.
281.         ELSE ' se non trova l'iterazione z% precedente con lo stesso dt! di quella in corso, incrementa n% e ripete la ricerca.
282.             test% = 0
283.             n% = n% + 1
284.         END IF
285.     LOOP UNTIL n% = z% - 1 'se non ha trovato iterazioni z% precedenti con lo stesso dt!, allora passa oltre e test%=0.
286.     SELECT CASE test%
287.         CASE IS = 1 'esce dalla subroutine, in quanto lo ietogramma e' definito.
288.             EXIT SUB
289.         CASE IS = 0 'calcola ex-novo lo ietogramma relativo all'iterazione z% in corso in quanto non vi e' un'iterazioni z% precedente con lo stesso dt! dell'iterazione z% in corso.
290.             k% = m%(ieto%, z%)
291.             i% = 0
292.             DO
293.                 matrice1(ieto%, z%, m%(ieto%, z%) - i%, 1).i2 = matrice1(ieto%, z%, k%, 1).i1
294.                 k% = k% - 1
295.                 IF k% = 0 THEN EXIT DO
296.                 matrice1(ieto%, z%, i% + 1, 1).i2 = matrice1(ieto%, z%, k%, 1).i1
297.                 k% = k% - 1
298.                 i% = i% + 1
299.             LOOP UNTIL k% = 0
300.     END SELECT
301. END SUB
302. '-----------------------------------------------------------------------------------------------------------------------------------------------------------------
303. SUB grafici (ore!, oreM!, QM!)
304.     'ore!: ore corrispondenti al massimo maggiore.
305.     'oreM!: ore corrispondenti alla portata massima.
306.     'QM!: portata massima.
307.  
308.     CONST R& = _RGB32(255, 0, 0)
309.     CONST G& = _RGB32(0, 255, 0)
310.     CONST B& = _RGB32(0, 0, 255)
311.     CONST bianco& = _RGB32(255, 255, 255)
312.     CONST giallo& = _RGB32(255, 255, 0)
313.  
314.     DIM L%, H% 'dimensioni del grafico.
315.     DIM grafico&, colore&
316.  
317.     DIM dx%, dy% 'origine degli assi dentro SCREEN.
318.     DIM dx!, dy! 'posizione dell'origine degli assi dentro WINDOW.
319.     DIM taccaX%, taccaY%
320.     DIM x%, y% 'conversione nelle coordinate di SCREEN delle coordinate dentro WINDOW corrispondenti a determinati valori.
321.  
322.     H% = _DESKTOPHEIGHT * 2 / 3 ' imposta l'altezza dello SCREEN pari ai 2/3 dello schermo visualizzabile, cosi' che su ogni PC sia visualizzabile.
323.     L% = H% * 1.62 'imposta la larghezza dello SCREEN in modo che L%xH% sia un rettangolo aureo.
324.     dx% = 39
325.     dy% = H% - 1 - dx%
326.     taccaX% = ore! \ 5 ': PRINT taccaX%: SLEEP
327.     taccaY% = QM! \ 5
328.     grafico& = _NEWIMAGE(L%, H%, 32)
329.  
330.     _DEST grafico&
331.     ' SCREEN grafico&
332.     LINE (0, 0)-(L% - 1, H% - 1), R, B
333.     WINDOW (0, 0)-(ore! * 1.15, QM! * 1.15) 'scala ascisse e ordinate in base alle ore e alla portata, estendendole del 15%.
334.     dx! = PMAP(dx%, 2)
335.     dy! = PMAP(dy%, 3)
336.     LINE (dx!, QM! * 1.1)-(dx!, dy!): LINE -(ore! * 1.1, dy!), bianco& 'disegna gli assi in modo che si estendano del 10% oltre valori massimi, ma senza arrivare alla massima estensione delle coordinate impostate da WINDOW: cosi' si ha che gli assi non toccano la cornice rossa e i valori massimi che non arrivano al limite degli assi.
337.     PSET (dx!, QM! * 1.1!), bianco&: DRAW "F20": PSET (dx!, QM! * 1.1), bianco&: DRAW "G20"
338.     PSET (ore! * 1.1, dy!), bianco&: DRAW "G20": PSET (ore! * 1.1, dy!), bianco&: DRAW "H20"
339.     _PRINTSTRING (dx%, PMAP(QM! * 1.1, 1) - dx% \ 2), "mc/s"
340.     _PRINTSTRING (PMAP(ore! * 1.1, 0), dy%), "ore"
341.     _PRINTSTRING (dx% \ 4, dy% - 7), "0"
342.     FOR z% = 1 TO 24
343.         SELECT CASE z% MOD 3
344.             CASE IS = 1
345.                 colore& = R&
346.             CASE IS = 2
347.                 colore& = G&
348.             CASE IS = 0
349.                 colore& = B&
350.         END SELECT
351.         CIRCLE (dx! + massimi1(visualizzaieto%, z%, 1).ore, dy! + massimi1(visualizzaieto%, z%, 1).portata), PMAP(3, 2), R&
352.         PAINT (dx! + massimi1(visualizzaieto%, z%, 1).ore + PMAP(0.5, 2), dy! + massimi1(visualizzaieto%, z%, 1).portata + PMAP(0.5, 2)), R&
353.         FOR i% = 1 TO 50
354.             '            _LIMIT limit
355.             IF i% = 1 THEN
356.                 LINE (dx!, dy!)-(dx! + idrogramma(visualizzaieto%, z%, i%, 1).ore, dy! + idrogramma(visualizzaieto%, z%, i%, 1).portata), colore&
357.             ELSE
358.                 LINE (dx! + idrogramma(visualizzaieto%, z%, i% - 1, 1).ore, dy! + idrogramma(visualizzaieto%, z%, i% - 1, 1).portata)-(dx! + idrogramma(visualizzaieto%, z%, i%, 1).ore, dy! + idrogramma(visualizzaieto%, z%, i%, 1).portata), colore&
359.             END IF
360.             '            _PUTIMAGE ((_DESKTOPWIDTH - L%) \ 2, (_DESKTOPHEIGHT - H%) \ 2), grafico&, schermo&
361.         NEXT i%
362.         IF z% = 1 THEN
363.             LINE (dx!, dy!)-(dx! + massimi1(visualizzaieto%, z%, 1).ore, dy! + massimi1(visualizzaieto%, z%, 1).portata), giallo&
364.         ELSE
365.             LINE (dx! + massimi1(visualizzaieto%, z% - 1, 1).ore, dy! + massimi1(visualizzaieto%, z% - 1, 1).portata)-(dx! + massimi1(visualizzaieto%, z%, 1).ore, dy! + massimi1(visualizzaieto%, z%, 1).portata), giallo&
366.         END IF
367.     NEXT z%
368.     LINE (dx!, dy! + QM!)-(dx! + oreM!, dy! + QM!), bianco&, , 61440
369.     LINE -(dx! + oreM!, dy!), bianco&, , 61440
370.     x% = PMAP(dx! + oreM!, 0)
371.     y% = PMAP(dy! + QM!, 1)
372.     _PRINTSTRING (x% + 10, y% - dx%), "(" + _TRIM$(STR$(oreM!)) + ","
373.     _PRINTSTRING (x% + 10, y% - dx% + 16), _TRIM$(STR$(QM!)) + ")"
374.     i% = 1
375.     WHILE i% * taccaX% <= ore!
376.         LINE (dx! + i% * taccaX%, PMAP(dy% + 5, 3))-(dx! + i% * taccaX%, PMAP(dy% - 5, 3))
377.         x% = PMAP(dx! + i% * taccaX%, 0)
378.         _PRINTSTRING (x% - 4 * LEN(_TRIM$(STR$(i% * taccaX%))), H% - 1 - dx% / 1.5), _TRIM$(STR$(i% * taccaX%))
379.         i% = i% + 1
380.     WEND
381.     i% = 1
382.     WHILE i% * taccaY% <= QM!
383.         LINE (PMAP(dx% - 5, 2), dy! + i% * taccaY%)-(PMAP(dx% + 5, 2), dy! + i% * taccaY%)
384.         y% = PMAP(dy! + i% * taccaY%, 1)
385.         _PRINTSTRING (dx% \ 4, y% - 7), _TRIM$(STR$(i% * taccaY%))
386.         i% = i% + 1
387.     WEND
388.     _PUTIMAGE ((_DESKTOPWIDTH - L%) \ 2, (_DESKTOPHEIGHT - H%) \ 2), grafico&, schermo&
389.  
390.     _DEST schermo&
391.     VIEW PRINT 1 TO 10
392.     '  DO
393.     '  CLS 2
394.     ' PRINT "Digitare le ore e i minuti (<= "; _TRIM$(STR$(INT(ore!))); " ore e "; _TRIM$(STR$(INT((ore! - INT(ore!)) * 60))); " minuti) di cui si vuole conoscere la portata di picco."
395.     'INPUT "Ore: ", ore2!
396.     'INPUT "Minuti: ", minuti%
397.     'ore2! = ore2! + minuti% / 60
398.     'LOOP WHILE ore2! > ore!
399.     ore2! = 3
400.     _DEST grafico&
401.     x% = PMAP(dx! + ore2!, 0)
402.     y% = PMAP(dy!, 1)
403.     SCREEN grafico& '------------------------------------------------------------------------------------>?????????????????????????????????????
404.     DO
405.         IF POINT(dx! + ore2!, PMAP(y%, 3)) = giallo& THEN EXIT DO
406.         y% = y% - 1
407.     LOOP
408.     ' _LIMIT limit
409.     LINE (dx! + ore2!, dy!)-(dx! + ore2!, PMAP(y%, 3)), bianco&, , 65280
410.     _PUTIMAGE ((_DESKTOPWIDTH - L%) \ 2, (_DESKTOPHEIGHT - H%) \ 2), grafico&, schermo&
411.     '   _FREEIMAGE grafico&
412. END SUB
413. '------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
414. SUB risultati ()
415.     FOR z% = 1 TO 24
416.         CLS
417.         PRINT "z%="; z%
418.         PRINT "dt="; dt!(visualizzaieto%, z%); ";"; dt!(visualizzaieto%, z%) * 60
419.         PRINT "visualizzaieto%="; visualizzaieto%
420.         PRINT "m%="; m%(visualizzaieto%, z%)
421.         PRINT "N="; coefficienti(z%).d / tc!
422.         PRINT "n="; coefficienti(z%).tSUta
423.         FOR i% = 1 TO 50
424.             '  _LIMIT limit
425.     PRINT _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).m)), _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).mdt)), _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).tminuti)), _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).h1)),_
426.      _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).i1)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).i2)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).DH)), _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).H2)),_
427.       _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).he)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).Dhe)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).MOCKUStSUta)),_
428.        _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).MOCKUSqSUqp)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).tmSUta)), _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).qmSUqp)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).qm))
429.         NEXT i%
430.         SLEEP
431.         CLS
432.         FOR n% = 1 TO 50
433.             '  _LIMIT limit
434.             LOCATE 1,
435.             FOR i% = 1 TO 50
436.                 LOCATE , n% * 4 - 3
437.                 PRINT matrice2!(visualizzaieto%, z%, i%, n%)
438.             NEXT i%
439.         NEXT n%
440.         SLEEP
441.     NEXT z%
442.     CLS
443.     FOR n% = 1 TO 24
444.         ' _LIMIT limit
445.         SELECT CASE n%
446.             CASE IS <= 16
447.                 FOR i% = 1 TO 50
448.                     LOCATE i%, n% * 15 - 14
449.                     PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).portata))
450.                     'PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).ore))
451.                     'PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).minuti))
452.  
453.                 NEXT i%
454.                 IF n% = 16 THEN SLEEP: CLS
455.             CASE IS > 16
456.                 FOR i% = 1 TO 50
457.                     LOCATE i%, (n% - 16) * 15 - 14
458.                     PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).portata))
459.                     'PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).ore))
460.                     'PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).minuti))
461.                 NEXT i%
462.         END SELECT
463.     NEXT n%
464.     SLEEP
465.     CLS
466.     FOR z% = 1 TO 24
467.         ' _LIMIT limit
468.         PRINT massimi1(visualizzaieto%, z%, 1).ore, massimi1(visualizzaieto%, z%, 1).minuti, massimi1(visualizzaieto%, z%, 1).portata
469.     NEXT z%
470.     SLEEP
471.     PRINT
472.     PRINT massimi2(visualizzaieto%).ore; massimi2(visualizzaieto%).portata
473. END SUB
474.  

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bplus on February 22, 2021, 11:21:31 am

Because of math rounding errors in a Window a point might not be exactly where you'd expect down to the pixel.

Use of Point to verify you are on yellow seems smart to me. It's not perfect but we do not live in perfect world.

Maybe I don't understand the problem. Do you want to work on the side and show results on main screen?

Do you know about setting up a separate drawing area and working in that, then use _Putimage of final result for main screen?

You do know that. :)

When using POINT you have to set SOURCE for it. This tells it where you want to read colors from. I bet that might be it.

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bplus on February 22, 2021, 11:39:42 am

@bartok

Yes I just scanned your code and no sign of Source.  So when working in an image on the side set _Source to that just like _Dest and them back when you are done and exiting that part of code. _Source is to let POINT know you are reading point locations in the side thing instead or main screen.

With _Source you can read POINT out of main screen or other images on the side.

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bartok on February 22, 2021, 12:14:06 pm

Perfect. And once again the computer was right!

Thanks!

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bartok on February 25, 2021, 06:41:05 am

Hi, I am in trouble again.
Below, there are 2 codes, that have to do the same thing. Code 1, with a subroutine, and it works. Code 2, splitting the subroutine of Code 1, obteining also a sub-subroutine: it doesn't work.

In the code 1, the line 135 calls the subroutine "grafici".

The subroutine "grafici" starts at line 301 and goes to SYSTEM at line 433 before the end of the subroutine. This code, as said, works fine.

In code 1, I have highlighted a part of the subroutine "grafici" (lines between 396 and 434), because it is the part that in Code 2 will became a sub-subroutine, and that is the core of my problem in the Code 2.

In fact, in the Code 2, I have aptempted to transform the lines between 396 and 434 (of Code 1) in a sub-subroutine, called "idrogrammaprogetto".

Speaking about the Code 2, it is egual to the Code 1 until the line 395, then the sub-subroutine "idrogrammaprogetto" is called.

the sub-subroutine "idrogrammaprogetto" starts at line 444, but the program, this time, goes in loop and never goes to system, as intended in line 399.

With "vwatch.exe", I have seen that, in code 2, the problem is the DO cicle of lines 468 to 472, corresponding to the lines 416 to 419 of code 1.

Similarly to one of my preview problems, this DO cicle searches a yellow pixel and, if it is succeeds, the program have to do something: in this case, it should go ahead and go to SYSTEM, as in Code 1. The other time, the program didn't exit the DO cicle, as none yellow pixel was found, because no _SOURCE was specified. Specyfing the _SOURCE, the problem was solved.

In this case, the problem is apparently the same: the program goes in loop in the DO cicle, because it doesn't match the yellow pixel. But this time, the _SOURCE is specified at line 459.

Futhermore, all the required variables for the sub-subroutine are passed by reference. for example: it is not possible that the DO cicle doesn't find the yellow pixel because the yellow variable is not passed. In fact, it is passed in lines 396 and 444: "giallo&". Similarly, is not possible that the yellow pixel is not found because the image with the graphics is not passed: it is "grafico&". So, I have finished the ideas. To conclude, the problem is not even matter of coordinates.

CODE 1.

Code: QB64: [Select]

1. OPTION BASE 1
2.  
3. mockus: 'dati dell'idrogramma unitario di Mockus
4. DATA 0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2,2.1,2.2,2.3,2.4,2.5,2.6,2.7,2.8,2.9,3,3.1,3.2,3.3,3.4,3.5,3.6,3.7,3.8,3.9,4,4.1,4.2,4.3,4.4,4.5,4.6,4.7,4.8,4.9,5: 't/ta
5. DATA 0.030,0.100,0.190,0.310,0.470,0.660,0.820,0.930,0.990,1.000,0.990,0.930,0.860,0.780,0.680,0.560,0.460,0.390,0.330,0.280,0.244,0.207,0.177,0.147,0.127,0.107,0.092,0.077,0.066,0.055,0.048,0.040,0.035,0.029,0.025,0.021,0.018,0.015,0.013,0.011,0.009,0.008,0.007,0.006,0.006,0.005,0.004,0.003,0.001,0.000: ' t/tp
6.  
7. coefficienti:
8. DATA 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,18,20,22,24,26,28,30,32: 'N volte tc!.
9. DATA 0.1,0.1,0.1,0.2,0.3,0.2,0.3,0.3,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,1,1,1,1,1,1,1,1: 'discretizzazione dei coeff. di Mockus.
10.  
11. TYPE mockus
12.     MOCKUStSUta AS SINGLE
13.     MOCKUSqSUqp AS SINGLE
14. END TYPE
15.  
16. TYPE matrice1
17.     m AS INTEGER '[numero naturale] conta il numero di interazioni fino al raggiungimento di "d".
18.     mdt AS SINGLE '[ore] conta la progessione delle ore fino al raggiungimento di "d".
19.     tminuti AS SINGLE '[minuti] trasforma mdt in minuti.
20.     h1 AS SINGLE '[mm] valore h della curva di possibilità climatica.
21.     i1 AS SINGLE '[mm/ora] intensità di pioggia.
22.     i2 AS SINGLE '[mm/ora] intensità di pioggia secondo lo ietogramma Chicago o costante.
23.     DH AS SINGLE '[mm] pioggia in ogni istante temporale dt.
24.     H2 AS SINGLE '[mm] pioggia cumulata ad ogni istante temporale.
25.     he AS SINGLE '[mm] piggia efficace cumulata ad ogni istante temporale.
26.     Dhe AS SINGLE '[mm] pioggia efficace ad ogni istante temporale.
27.     MOCKUStSUta AS SINGLE
28.     MOCKUSqSUqp AS SINGLE
29.     tmSUta AS SINGLE '[-].
30.     qmSUqp AS SINGLE '[-] il valore di qm/qp corrispondente a tm/ta e' = al valore di q/qp che corrisponde a t/ta=tm/ta.
31.     qm AS SINGLE '[mc/s].
32. END TYPE
33.  
34. TYPE coefficienti
35.     d AS SINGLE
36.     tSUta AS SINGLE
37. END TYPE
38.  
39. TYPE idrogramma
40.     ore AS SINGLE
41.     minuti AS SINGLE
42.     portata AS SINGLE
43. END TYPE
44.  
45. CONST limit = 200 '50
46.  
47. DIM SHARED mockus(50) AS mockus
48. DIM SHARED dt!(2, 24)
49. DIM SHARED matrice1(2, 24, 50, 1) AS matrice1
50. DIM SHARED coefficienti(24) AS coefficienti
51. DIM SHARED matrice2!(2, 24, 50, 50)
52. DIM SHARED idrogramma(2, 24, 50, 1) AS idrogramma
53. DIM SHARED massimi1(2, 24, 1) AS idrogramma
54. DIM SHARED massimi2(2) AS idrogramma
55. DIM SHARED i%, n%, z%, ieto%
56. DIM SHARED m%(2, 24)
57. DIM SHARED S2! '[mm] contenuto idrico massimo del terreno.
58. DIM SHARED Ia! '[mm] perdite iniziali.
59. DIM SHARED tc! '[ore] tempo di corrivazione
60. DIM SHARED ta! '[ore] tempo di picco.
61. DIM SHARED qp! '[mc/s] portata al colmo dell'idrogramma unitario.
62. DIM SHARED k! '[-] coeff. della curva di possibilità climatica tado dal sito della Regione.
63. DIM SHARED a2! '[mm/d^n] coeff. della curva di possibilità climatica.
64. DIM SHARED n1! '[-] coeff. della curva di possibilità climatica.
65. DIM SHARED visualizzaieto%
66. DIM schermo&
67. DIM A1! '[mq] superficie bacino idrografico.
68. DIM L! '[m] lunghezza asta principale.
69. DIM HmaxL! '[m] punto più alto dell'astra principale.
70. DIM HminL! '[m] punto più basso dell'asta principale.
71. DIM s1! '[-] pendenza media del bacino non espressa in percentuale ma come 0,x.
72. DIM CNII! '[-] tabellato.
73. DIM CNIII '[-] discende da CNII.
74. DIM tl! '[ore] tempo di ritardo.
75. DIM keypress$
76.  
77. schermo& = _NEWIMAGE(_DESKTOPWIDTH, _DESKTOPHEIGHT, 32)
78. SCREEN schermo&
79. _SCREENMOVE 0, 0
80.  
81. A1! = 2641902
82. k! = 2.34
83. a2! = 14.6653
84. n1! = 0.45596
85. L! = 3300
86. HmaxL! = 1981
87. HminL! = 1062
88. s1! = 0.415
89. CNII! = 70
90.  
91. 'INPUT A1!
92. 'INPUT k!
93. 'INPUT a2!
94. 'INPUT n1!
95. 'INPUT L!
96. 'INPUT HmaxL!
97. 'INPUT HminL!
98. 'INPUT s1!
99. 'INPUT CNII!
100. 'INPUT d1!
101.  
102. CNIII! = (23 * CNII!) / (10 + 0.13 * CNII!)
103. tl! = 0.342 * ((L! / 1000) ^ 0.8 / (100 * s1!) ^ 0.5) * (1000 / CNIII! - 9) ^ 0.7 'formula di Mockus.
104. S2! = 25.4 * (1000 / CNIII! - 10)
105. Ia! = 0.1 * S2! 'coeff.=0.03-0.2.
106. tc! = tl! / 0.6
107. ta! = tl! / 0.9
108. qp! = 0.208 * ((A1! / 1000000) / ta!)
109.  
110. RESTORE mockus
111. FOR i% = 1 TO 100
112.     IF i% <= 50 THEN READ mockus(i%).MOCKUStSUta
113.     IF i% > 50 THEN READ mockus(i% - 50).MOCKUSqSUqp
114. NEXT i%
115.  
116. RESTORE coefficienti
117. FOR i% = 1 TO 48
118.     IF i% <= 24 THEN READ coefficienti(i%).d: coefficienti(i%).d = coefficienti(i%).d * tc!
119.     IF i% > 24 THEN READ coefficienti(i% - 24).tSUta
120. NEXT i%
121.  
122. CLS
123. PRINT "Ietogramma Chicago/costante: 1/2> "
124. 'DO
125. '_LIMIT 60 'limit
126. 'keypress$ = INKEY$
127. 'LOOP UNTIL keypress$ = "1" OR keypress$ = "2"
128. visualizzaieto% = 1
129.  
130. FOR ieto% = 1 TO 2
131.     CALL creamatrici
132. NEXT ieto%
133.  
134. CLS
135. CALL grafici(massimi1(visualizzaieto%, 24, 1).ore, massimi2(visualizzaieto%).ore, massimi2(visualizzaieto%).portata) 'valori passati per disegnare il grafico.
136. SLEEP
137. SYSTEM
138. 'CALL risultati
139.  
140. '-----------------------------------------matrice1---------------------------------------------------------------------------------------------------------------------------------
141. SUB creamatrici ()
142.     z% = 0
143.     DO
144.         z% = z% + 1
145.         dt!(ieto%, z%) = coefficienti(z%).tSUta * ta!
146.         m%(ieto%, z%) = coefficienti(z%).d / (coefficienti(z%).tSUta * ta!) 'determina il numero di passi temporali necessari ad arrivare allla durata d considerata con il paso temporale dt!.
147.         FOR i% = 1 TO 50
148.             matrice1(ieto%, z%, i%, 1).MOCKUStSUta = mockus(i%).MOCKUStSUta
149.             matrice1(ieto%, z%, i%, 1).MOCKUSqSUqp = mockus(i%).MOCKUSqSUqp
150.             matrice1(ieto%, z%, i%, 1).m = i%
151.             matrice1(ieto%, z%, i%, 1).mdt = matrice1(ieto%, z%, i%, 1).m * dt!(ieto%, z%)
152.             matrice1(ieto%, z%, i%, 1).tminuti = matrice1(ieto%, z%, i%, 1).mdt * 60
153.             SELECT CASE (_ROUND(10 * (i% * coefficienti(z%).tSUta))) / 10 'crea il vettore dell'idrogramma unitario di Mockus con i multipli del coeff. z% considerato (0.1,0.2,0.3,0.5,1).
154.                 CASE IS <= 5 'valore massimo di t/ta di Mockus
155.                     matrice1(ieto%, z%, i%, 1).tmSUta = (_ROUND(10 * (i% * coefficienti(z%).tSUta))) / 10 '=matrice1(i%).m *  matrice1(i%).mdt / ta!
156.                 CASE ELSE
157.                     matrice1(ieto%, z%, i%, 1).tmSUta = 0
158.             END SELECT
159.         NEXT i%
160.         FOR i% = 1 TO 50
161.             n% = 1
162.             DO 'crea il vettore dell'idrogramma unitario Mockus prendendo i valori corrispondenti di ai coefficienti.
163.                 IF matrice1(ieto%, z%, i%, 1).tmSUta = matrice1(ieto%, z%, n%, 1).MOCKUStSUta THEN
164.                     matrice1(ieto%, z%, i%, 1).qmSUqp = matrice1(ieto%, z%, n%, 1).MOCKUSqSUqp
165.                     EXIT DO
166.                 ELSE
167.                     n% = n% + 1
168.                 END IF
169.             LOOP UNTIL n% = 50
170.             matrice1(ieto%, z%, i%, 1).qm = matrice1(ieto%, z%, i%, 1).qmSUqp * qp!
171.         NEXT i%
172.         FOR i% = 1 TO m%(ieto%, z%)
173.             matrice1(ieto%, z%, i%, 1).h1 = k! * a2! * matrice1(ieto%, z%, i%, 1).mdt ^ n1!
174.             IF i% = 1 THEN
175.                 matrice1(ieto%, z%, i%, 1).i1 = matrice1(ieto%, z%, i%, 1).h1 / dt!(ieto%, z%)
176.             ELSE
177.                 matrice1(ieto%, z%, i%, 1).i1 = (matrice1(ieto%, z%, i%, 1).h1 - matrice1(ieto%, z%, i% - 1, 1).h1) / dt!(ieto%, z%)
178.             END IF
179.         NEXT i%
180.         SELECT CASE ieto%
181.             CASE 1 'Chicago.
182.                 CALL chicago
183.             CASE 2 'costante
184.                 FOR i% = 1 TO m%(ieto%, z%)
185.                     matrice1(ieto%, z%, i%, 1).i2 = matrice1(ieto%, z%, m%(ieto%, z%), 1).h1 / matrice1(ieto%, z%, m%(ieto%, z%), 1).mdt
186.                 NEXT i%
187.         END SELECT
188.         FOR i% = 1 TO m%(ieto%, z%)
189.             matrice1(ieto%, z%, i%, 1).DH = matrice1(ieto%, z%, i%, 1).i2 * dt!(ieto%, z%)
190.             IF i% = 1 THEN
191.                 matrice1(ieto%, z%, i%, 1).H2 = matrice1(ieto%, z%, i%, 1).DH
192.             ELSE
193.                 matrice1(ieto%, z%, i%, 1).H2 = matrice1(ieto%, z%, i%, 1).DH + matrice1(ieto%, z%, i% - 1, 1).H2
194.             END IF
195.             SELECT CASE matrice1(ieto%, z%, i%, 1).H2 - Ia!
196.                 CASE IS < 0
197.                     matrice1(ieto%, z%, i%, 1).he = 0
198.                 CASE IS >= 0
199.                     matrice1(ieto%, z%, i%, 1).he = (matrice1(ieto%, z%, i%, 1).H2 - Ia!) ^ 2 / (matrice1(ieto%, z%, i%, 1).H2 - Ia! + S2!)
200.             END SELECT
201.             IF i% = 1 THEN
202.                 matrice1(ieto%, z%, i%, 1).Dhe = matrice1(ieto%, z%, i%, 1).he
203.             ELSE
204.                 matrice1(ieto%, z%, i%, 1).Dhe = matrice1(ieto%, z%, i%, 1).he - matrice1(ieto%, z%, i% - 1, 1).he
205.             END IF
206.         NEXT i%
207.  
208.         '-------------------------------matrice2-------------------------------------------------------------------------------------------------------------------------
209.  
210.         FOR n% = 1 TO 50 'colonna
211.             FOR i% = n% TO 50 ' riga
212.                 matrice2!(ieto%, z%, i%, n%) = matrice1(ieto%, z%, i% - n% + 1, 1).Dhe * matrice1(ieto%, z%, n%, 1).qm
213.             NEXT i%
214.         NEXT n%
215.         '-------------------------------idrogramma------------------------------------------------------------------------------------------------------------------
216.         FOR i% = 1 TO 50 'riga
217.             FOR n% = 1 TO 50 'colonna
218.                 idrogramma(ieto%, z%, i%, 1).portata = idrogramma(ieto%, z%, i%, 1).portata + matrice2!(ieto%, z%, i%, n%)
219.  
220.             NEXT n%
221.             idrogramma(ieto%, z%, i%, 1).ore = matrice1(ieto%, z%, i%, 1).mdt
222.             idrogramma(ieto%, z%, i%, 1).minuti = matrice1(ieto%, z%, i%, 1).tminuti
223.         NEXT i%
224.         '-------------------------------Massimi-----------------------------------------------------
225.         i% = 1
226.         'massimi1(ieto%, z%, 1).portata = idrogramma(ieto%, z%, 1, 1).portata
227.         DO
228.             IF massimi1(ieto%, z%, 1).portata < idrogramma(ieto%, z%, i%, 1).portata THEN
229.                 massimi1(ieto%, z%, 1).portata = idrogramma(ieto%, z%, i%, 1).portata
230.                 massimi1(ieto%, z%, 1).ore = idrogramma(ieto%, z%, i%, 1).ore
231.                 massimi1(ieto%, z%, 1).minuti = idrogramma(ieto%, z%, i%, 1).minuti
232.             END IF
233.             i% = i% + 1
234.         LOOP UNTIL i% = 50
235.         ' IF massimi2(ieto%).ore < massimi1(visualizzaieto%, z%, 1).ore THEN massimi2(ieto%).ore = massimi1(visualizzaieto%, z%, 1).ore
236.         IF massimi2(ieto%).portata < massimi1(visualizzaieto%, z%, 1).portata THEN
237.             massimi2(ieto%).portata = massimi1(visualizzaieto%, z%, 1).portata
238.             massimi2(ieto%).ore = massimi1(visualizzaieto%, z%, 1).ore
239.         END IF
240.     LOOP UNTIL z% = UBOUND(coefficienti)
241. END SUB
242. '---------------------------------------------------------------------------------------------------------------------------------------------------
243. 'Questa subroutine crea gli ietogrammi chicago per ogni iterazione z%.
244. 'Ad ogni nuovo dt! pone l'ultimo valore di pioggia come ultimo, il penultimo come primo, il terzultimo come penultimo, il quartultimo come secondo e così via.
245. 'Quando vi è un valore di z% che ha un dt! non nuovo, cerca la precedente iterazione z% con lo stesso dt! e copia i suoi valori dello ietogramma chicago in quello dell'iterazione z% in corso. I valori rimanenti sono riempiti
246. 'procedendo come prima.
247. 'Quindi per esempio l'iterazione z%=1 corrisponde a d=1tc con passo temporale dt!. L'iterazione z%=2 corrisponde a d=2tc con il medesimo passo temporale. Per cui lo ietogramma chicago per d=2tc e' dato dallo stesso ietogramma
248. 'chicago che si ha per d=1tc, con l'aggiunta, ai suoi estremi, dei nuovi valori. Il meccanismo si ripete per z%=3, che corrisponde a d=3tc, con il medesimo dt!. Per z%=4, corrispondente a d=4tc, invece, il passo temporale dt!
249. 'cambia.
250. 'Quanto detto e' necessario per evitare piccole discrepanze che potrebbero comportare il fatto che alcuni massimi risultino localmente poco piu' bassi del precedente, anche quando il trend generale e' chiaramente di aumento dei
251. 'massimi all'aumentare della durata della pioggia. In questo modo si fa l'ipotesi che all'aumentare della durata della pioggia, discretizzata con uno stesso passo temporale dt!, ogni ietogramma contenga a sua volta il precedente.
252. 'In tal modo, se si condiera per esempio d=2tc, si assume che l'altezza di pioggia corrispondente a 1tc, segua lo ietogramma relativo a d=1tc, cosi' che l'altezza di pioggia rimamente sia eclusivamente dovuta ai nuovi valori.
253. SUB chicago ()
254.     DIM test%
255.     DIM k%
256.     n% = 1
257.     DO
258.         IF z% = 1 THEN test% = 0: EXIT DO 'permette il calcolo dello ietogramma chicago della prima iterazione. test%=0
259.         IF dt!(ieto%, z%) = dt!(ieto%, z% - n%) THEN ' se trova un'iterazione z% precedente con lo stesso dt! di quella in corso, esegue i comandi successivi
260.             test% = 1
261.             k% = m%(ieto%, z%)
262.             i% = 0
263.             DO 'inserisce i nuovi valori dello ietogramma chicago dell'iterazione z% in corso.
264.                 matrice1(ieto%, z%, m%(ieto%, z%) - i%, 1).i2 = matrice1(ieto%, z%, k%, 1).i1
265.                 k% = k% - 1
266.                 IF k% = m%(ieto%, z% - n%) THEN EXIT DO 'se ha completato l'inserimento dei nuovi valori, passa oltre.
267.                 matrice1(ieto%, z%, i% + 1, 1).i2 = matrice1(ieto%, z%, k%, 1).i1
268.                 k% = k% - 1
269.                 i% = i% + 1
270.             LOOP UNTIL k% = m%(ieto%, z% - n%)
271.             IF matrice1(ieto%, z%, i%, 1).i2 = 0 THEN i% = i% - 1
272.             k% = 1
273.             DO 'inserisce lo ietogramma chicago dell'iterazione z% precedente con lo stesso dt!, nello ietogramma dell'iterazione z% in corso.
274.                 i% = i% + 1
275.                 matrice1(ieto%, z%, i%, 1).i2 = matrice1(ieto%, z% - n%, k%, 1).i2
276.                 k% = k% + 1
277.             LOOP UNTIL k% > m%(ieto%, z% - n%)
278.             EXIT DO 'siccome ha completato lo ietogramma dell'iterazione z% in corso, non serve incrementare n%, quindi passa oltre. test%=1.
279.         ELSE ' se non trova l'iterazione z% precedente con lo stesso dt! di quella in corso, incrementa n% e ripete la ricerca.
280.             test% = 0
281.             n% = n% + 1
282.         END IF
283.     LOOP UNTIL n% = z% - 1 'se non ha trovato iterazioni z% precedenti con lo stesso dt!, allora passa oltre e test%=0.
284.     SELECT CASE test%
285.         CASE IS = 1 'esce dalla subroutine, in quanto lo ietogramma e' definito.
286.             EXIT SUB
287.         CASE IS = 0 'calcola ex-novo lo ietogramma relativo all'iterazione z% in corso in quanto non vi e' un'iterazioni z% precedente con lo stesso dt! dell'iterazione z% in corso.
288.             k% = m%(ieto%, z%)
289.             i% = 0
290.             DO
291.                 matrice1(ieto%, z%, m%(ieto%, z%) - i%, 1).i2 = matrice1(ieto%, z%, k%, 1).i1
292.                 k% = k% - 1
293.                 IF k% = 0 THEN EXIT DO
294.                 matrice1(ieto%, z%, i% + 1, 1).i2 = matrice1(ieto%, z%, k%, 1).i1
295.                 k% = k% - 1
296.                 i% = i% + 1
297.             LOOP UNTIL k% = 0
298.     END SELECT
299. END SUB
300. '-----------------------------------------------------------------------------------------------------------------------------------------------------------------
301. SUB grafici (ore!, oreM!, QM!)
302.     'ore!: ore corrispondenti al massimo maggiore.
303.     'oreM!: ore corrispondenti alla portata massima.
304.     'QM!: portata massima.
305.  
306.     CONST R& = _RGB32(255, 0, 0)
307.     CONST G& = _RGB32(0, 255, 0)
308.     CONST B& = _RGB32(0, 0, 255)
309.     CONST bianco& = _RGB32(255, 255, 255)
310.     CONST giallo& = _RGB32(255, 255, 0)
311.     CONST grigio& = _RGB32(127, 127, 127)
312.  
313.     DIM L%, H% 'dimensioni del grafico.
314.     DIM grafico&, graficoNERO&, graficoCOPIA&
315.     DIM colore&
316.  
317.     DIM dx%, dy% 'origine degli assi dentro SCREEN.
318.     DIM dx!, dy! 'posizione dell'origine degli assi dentro WINDOW.
319.     DIM taccaX%, taccaY%
320.     DIM portata!(2)
321.     DIM soglia!
322.     REDIM x%(1), y%(1) 'conversione nelle coordinate di SCREEN delle coordinate dentro WINDOW corrispondenti a determinati valori.
323.  
324.     H% = _DESKTOPHEIGHT * 2 / 3 ' imposta l'altezza dello SCREEN pari ai 2/3 dello schermo visualizzabile, cosi' che su ogni PC sia visualizzabile.
325.     L% = H% * 1.62 'imposta la larghezza dello SCREEN in modo che L%xH% sia un rettangolo aureo.
326.     dx% = 39
327.     dy% = H% - 1 - dx%
328.     taccaX% = ore! \ 5 ': PRINT taccaX%: SLEEP
329.     taccaY% = QM! \ 5
330.     grafico& = _NEWIMAGE(L%, H%, 32)
331.     graficoNERO& = _NEWIMAGE(L%, H%, 32) 'crea un'immagine della stessa dimensione di grafico& che sarà colorata di nero.
332.  
333.     _DEST graficoNERO&
334.     CLS ' colora di nero l'immagine graficoNERO&.
335.     _DEST grafico&
336.     LINE (0, 0)-(L% - 1, H% - 1), R, B
337.     WINDOW (0, 0)-(ore! * 1.15, QM! * 1.15) 'scala ascisse e ordinate in base alle ore e alla portata, estendendole del 15%.
338.     dx! = PMAP(dx%, 2)
339.     dy! = PMAP(dy%, 3)
340.     LINE (dx!, QM! * 1.1)-(dx!, dy!): LINE -(ore! * 1.1, dy!), bianco& 'disegna gli assi in modo che si estendano del 10% oltre valori massimi, ma senza arrivare alla massima estensione delle coordinate impostate da WINDOW: cosi' si ha che gli assi non toccano la cornice rossa e i valori massimi che non arrivano al limite degli assi.
341.     PSET (dx!, QM! * 1.1!), bianco&: DRAW "F20": PSET (dx!, QM! * 1.1), bianco&: DRAW "G20"
342.     PSET (ore! * 1.1, dy!), bianco&: DRAW "G20": PSET (ore! * 1.1, dy!), bianco&: DRAW "H20"
343.     _PRINTSTRING (dx%, PMAP(QM! * 1.1, 1) - dx% \ 2), "mc/s"
344.     _PRINTSTRING (PMAP(ore! * 1.1, 0), dy%), "ore"
345.     _PRINTSTRING (dx% \ 4, dy% - 7), "0"
346.     i% = 1
347.     WHILE i% * taccaX% <= ore!
348.         LINE (dx! + i% * taccaX%, PMAP(dy% + 5, 3))-(dx! + i% * taccaX%, PMAP(dy% - 5, 3))
349.         x%(1) = PMAP(dx! + i% * taccaX%, 0)
350.         _PRINTSTRING (x%(1) - 4 * LEN(_TRIM$(STR$(i% * taccaX%))), H% - 1 - dx% / 1.5), _TRIM$(STR$(i% * taccaX%))
351.         i% = i% + 1
352.     WEND
353.     i% = 1
354.     WHILE i% * taccaY% <= QM!
355.         LINE (PMAP(dx% - 5, 2), dy! + i% * taccaY%)-(PMAP(dx% + 5, 2), dy! + i% * taccaY%)
356.         y%(1) = PMAP(dy! + i% * taccaY%, 1)
357.         _PRINTSTRING (dx% \ 4, y%(1) - 7), _TRIM$(STR$(i% * taccaY%))
358.         i% = i% + 1
359.     WEND
360.     FOR z% = 1 TO 24
361.         SELECT CASE z% MOD 3
362.             CASE IS = 1
363.                 colore& = R&
364.             CASE IS = 2
365.                 colore& = G&
366.             CASE IS = 0
367.                 colore& = B&
368.         END SELECT
369.         CIRCLE (dx! + massimi1(visualizzaieto%, z%, 1).ore, dy! + massimi1(visualizzaieto%, z%, 1).portata), PMAP(3, 2), R&
370.         PAINT (dx! + massimi1(visualizzaieto%, z%, 1).ore + PMAP(0.5, 2), dy! + massimi1(visualizzaieto%, z%, 1).portata + PMAP(0.5, 2)), R&
371.         FOR i% = 1 TO 50
372.             '            _LIMIT limit
373.             IF i% = 1 THEN
374.                 LINE (dx!, dy!)-(dx! + idrogramma(visualizzaieto%, z%, i%, 1).ore, dy! + idrogramma(visualizzaieto%, z%, i%, 1).portata), colore&
375.             ELSE
376.                 LINE (dx! + idrogramma(visualizzaieto%, z%, i% - 1, 1).ore, dy! + idrogramma(visualizzaieto%, z%, i% - 1, 1).portata)-(dx! + idrogramma(visualizzaieto%, z%, i%, 1).ore, dy! + idrogramma(visualizzaieto%, z%, i%, 1).portata), colore&
377.             END IF
378.             '            _PUTIMAGE ((_DESKTOPWIDTH - L%) \ 2, (_DESKTOPHEIGHT - H%) \ 2), grafico&, schermo&
379.         NEXT i%
380.         IF z% = 1 THEN
381.             LINE (dx!, dy!)-(dx! + massimi1(visualizzaieto%, z%, 1).ore, dy! + massimi1(visualizzaieto%, z%, 1).portata), giallo&
382.         ELSE
383.             LINE (dx! + massimi1(visualizzaieto%, z% - 1, 1).ore, dy! + massimi1(visualizzaieto%, z% - 1, 1).portata)-(dx! + massimi1(visualizzaieto%, z%, 1).ore, dy! + massimi1(visualizzaieto%, z%, 1).portata), giallo&
384.         END IF
385.     NEXT z%
386.     PSET (dx!, dy!), giallo&
387.     LINE (dx!, dy! + QM!)-(dx! + oreM!, dy! + QM!), R&, , 65520
388.     LINE -(dx! + oreM!, dy!), R&, , 65520
389.     x%(1) = PMAP(dx! + oreM!, 0)
390.     y%(1) = PMAP(dy! + QM!, 1)
391.     _PRINTSTRING (x%(1) + 10, y%(1) - dx%), "(" + _TRIM$(STR$(oreM!)) + ","
392.     _PRINTSTRING (x%(1) + 10, y%(1) - dx% + 16), _TRIM$(STR$(QM!)) + ")"
393.     _PUTIMAGE ((_DESKTOPWIDTH - L%) \ 2, (_DESKTOPHEIGHT - H%) \ 2), grafico&, schermo&
394.     graficoCOPIA& = _COPYIMAGE(grafico&) 'crea una copia di grafico& prima che venga effettuato il tratteggio grigio: INSERIRE POI NUOVA RICHIESTA.
395.  
396.     '---------------------------------------------------------------------------------------------------------------
397.     _DEST schermo&
398.     VIEW PRINT 1 TO 10
399.     SELECT CASE visualizzaieto%
400.         CASE IS = 1
401.             ' INPUT "Inserire percentuale soglia per portata progetto (numero intero senza simbolo %: ", soglia!
402.             soglia! = 0.1 'soglia! / 100
403.             grafico& = _COPYIMAGE(graficoCOPIA&) ' grafico& torna ad ed essere senza il tratteggio grigio.
404.             _PUTIMAGE ((_DESKTOPWIDTH - L%) \ 2, (_DESKTOPHEIGHT - H%) \ 2), graficoNERO&, schermo& 'incollando grafico& (senza il tratteggio) su schermo&, si vedrebbe comunque il tratteggio in quanto la precedente versione di grafico& (che era con il tratteggio), era stata incollata su schermo&. Quindi prima incollo l'immagine graficoNERO& su schermo&.
405.             _PUTIMAGE ((_DESKTOPWIDTH - L%) \ 2, (_DESKTOPHEIGHT - H%) \ 2), grafico&, schermo& 'copia grafico& (senza tratteggio) su schermo&, che era stato ripulito alla linea precedente.
406.             _DEST grafico&
407.             _SOURCE grafico&
408.             i% = 0
409.             ore2! = 0
410.  
411.             y%(1) = PMAP(dy!, 1)
412.             DO
413.                 i% = i% + 1
414.                 ore2! = ore2! + 1
415.                 REDIM _PRESERVE y%(i%)
416.                 DO
417.                     IF POINT(dx! + ore2!, PMAP(y%(i%), 3)) = giallo& THEN EXIT DO
418.                     y%(i%) = y%(i%) - 1
419.                 LOOP
420.                 IF i% > 1 THEN
421.                     IF PMAP(y%(i%), 3) - dy! < (1 + soglia!) * (PMAP(y%(i% - 1), 3) - dy!) THEN EXIT DO
422.                 END IF
423.             LOOP
424.             ore2! = ore2! - 1
425.             REDIM _PRESERVE y%(1 TO UBOUND(y%) - 1)
426.             portata!(1) = PMAP(y%(UBOUND(y%)), 3) - dy!
427.             _DEST schermo&
428.             PRINT "Portata progetto:"; portata!(1); "mc/s, all'ora n."; ore2!
429.         CASE IS = 2
430.             PRINT "La portata di progetto e' quella visualizzata."
431.     END SELECT
432.     SLEEP
433.     SYSTEM
434.     '--------------------------------------------------------------------------------------------------------------------------------------------------
435.  
436.     DO
437.         '  CLS 2
438.         PRINT "Digitare le ore e i minuti (<= "; _TRIM$(STR$(INT(ore!))); " ore e "; _TRIM$(STR$(INT((ore! - INT(ore!)) * 60))); " minuti) di cui si vuole conoscere la portata di picco."
439.         INPUT "Ore: ", ore2!
440.         INPUT "Minuti: ", minuti%
441.         ore2! = ore2! + minuti% / 60
442.     LOOP WHILE ore2! > ore!
443.     _DEST grafico&
444.     _SOURCE grafico&
445.     y%(1) = PMAP(dy!, 1)
446.     i% = 2
447.     DO
448.         REDIM _PRESERVE y%(i%)
449.         y%(i%) = y%(i% - 1) - 1
450.         IF POINT(dx! + ore2!, PMAP(y%(i% - 1), 3)) = giallo& THEN EXIT DO
451.         i% = i% + 1
452.     LOOP
453.     REDIM _PRESERVE y%(1 TO UBOUND(y%) - 1)
454.     portata!(2) = PMAP(y%(UBOUND(y%)), 3) - dy!
455.     FOR i% = 1 TO UBOUND(y%)
456.         LINE (dx! + ore2!, dy!)-(dx! + ore2!, PMAP(y%(i%), 3)), grigio&, , 65280
457.         _PUTIMAGE ((_DESKTOPWIDTH - L%) \ 2, (_DESKTOPHEIGHT - H%) \ 2), grafico&, schermo&
458.     NEXT i%
459.     x%(1) = PMAP(dx! + ore2!, 0)
460.     i% = 2
461.     DO
462.         REDIM _PRESERVE x%(i%)
463.         x%(i%) = x%(i% - 1) - 1
464.         IF POINT(PMAP(x%(i% - 1), 2), PMAP(y%(UBOUND(y%)), 3)) = bianco& THEN EXIT DO
465.         i% = i% + 1
466.     LOOP
467.     REDIM _PRESERVE x%(1 TO UBOUND(x%) - 1)
468.     FOR i% = 1 TO UBOUND(x%)
469.         LINE (dx! + ore2!, PMAP(y%(UBOUND(y%)), 3))-(PMAP(x%(i%), 2), PMAP(y%(UBOUND(y%)), 3)), grigio&, , 65280
470.         _PUTIMAGE ((_DESKTOPWIDTH - L%) \ 2, (_DESKTOPHEIGHT - H%) \ 2), grafico&, schermo&
471.     NEXT i%
472.     _DEST schermo&
473.     PRINT "Portata:"; portata!(2); "mc/s"
474.     _FREEIMAGE grafico&
475.     _FREEIMAGE graficoNERO&
476.     _FREEIMAGE graficoCOPIA&
477. END SUB
478. '------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
479. SUB risultati ()
480.     FOR z% = 1 TO 24
481.         CLS
482.         PRINT "z%="; z%
483.         PRINT "dt="; dt!(visualizzaieto%, z%); ";"; dt!(visualizzaieto%, z%) * 60
484.         PRINT "visualizzaieto%="; visualizzaieto%
485.         PRINT "m%="; m%(visualizzaieto%, z%)
486.         PRINT "N="; coefficienti(z%).d / tc!
487.         PRINT "n="; coefficienti(z%).tSUta
488.         FOR i% = 1 TO 50
489.             '  _LIMIT limit
490.     PRINT _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).m)), _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).mdt)), _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).tminuti)), _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).h1)),_
491.      _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).i1)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).i2)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).DH)), _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).H2)),_
492.       _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).he)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).Dhe)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).MOCKUStSUta)),_
493.        _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).MOCKUSqSUqp)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).tmSUta)), _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).qmSUqp)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).qm))
494.         NEXT i%
495.         SLEEP
496.         CLS
497.         FOR n% = 1 TO 50
498.             '  _LIMIT limit
499.             LOCATE 1,
500.             FOR i% = 1 TO 50
501.                 LOCATE , n% * 4 - 3
502.                 PRINT matrice2!(visualizzaieto%, z%, i%, n%)
503.             NEXT i%
504.         NEXT n%
505.         SLEEP
506.     NEXT z%
507.     CLS
508.     FOR n% = 1 TO 24
509.         ' _LIMIT limit
510.         SELECT CASE n%
511.             CASE IS <= 16
512.                 FOR i% = 1 TO 50
513.                     LOCATE i%, n% * 15 - 14
514.                     PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).portata))
515.                     'PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).ore))
516.                     'PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).minuti))
517.  
518.                 NEXT i%
519.                 IF n% = 16 THEN SLEEP: CLS
520.             CASE IS > 16
521.                 FOR i% = 1 TO 50
522.                     LOCATE i%, (n% - 16) * 15 - 14
523.                     PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).portata))
524.                     'PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).ore))
525.                     'PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).minuti))
526.                 NEXT i%
527.         END SELECT
528.     NEXT n%
529.     SLEEP
530.     CLS
531.     FOR z% = 1 TO 24
532.         ' _LIMIT limit
533.         PRINT massimi1(visualizzaieto%, z%, 1).ore, massimi1(visualizzaieto%, z%, 1).minuti, massimi1(visualizzaieto%, z%, 1).portata
534.     NEXT z%
535.     SLEEP
536.     PRINT
537.     PRINT massimi2(visualizzaieto%).ore; massimi2(visualizzaieto%).portata
538. END SUB
539.  

CODE 2.

Code: QB64: [Select]

1. OPTION BASE 1
2.  
3. mockus: 'dati dell'idrogramma unitario di Mockus
4. DATA 0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2,2.1,2.2,2.3,2.4,2.5,2.6,2.7,2.8,2.9,3,3.1,3.2,3.3,3.4,3.5,3.6,3.7,3.8,3.9,4,4.1,4.2,4.3,4.4,4.5,4.6,4.7,4.8,4.9,5: 't/ta
5. DATA 0.030,0.100,0.190,0.310,0.470,0.660,0.820,0.930,0.990,1.000,0.990,0.930,0.860,0.780,0.680,0.560,0.460,0.390,0.330,0.280,0.244,0.207,0.177,0.147,0.127,0.107,0.092,0.077,0.066,0.055,0.048,0.040,0.035,0.029,0.025,0.021,0.018,0.015,0.013,0.011,0.009,0.008,0.007,0.006,0.006,0.005,0.004,0.003,0.001,0.000: ' t/tp
6.  
7. coefficienti:
8. DATA 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,18,20,22,24,26,28,30,32: 'N volte tc!.
9. DATA 0.1,0.1,0.1,0.2,0.3,0.2,0.3,0.3,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,1,1,1,1,1,1,1,1: 'discretizzazione dei coeff. di Mockus.
10.  
11. TYPE mockus
12.     MOCKUStSUta AS SINGLE
13.     MOCKUSqSUqp AS SINGLE
14. END TYPE
15.  
16. TYPE matrice1
17.     m AS INTEGER '[numero naturale] conta il numero di interazioni fino al raggiungimento di "d".
18.     mdt AS SINGLE '[ore] conta la progessione delle ore fino al raggiungimento di "d".
19.     tminuti AS SINGLE '[minuti] trasforma mdt in minuti.
20.     h1 AS SINGLE '[mm] valore h della curva di possibilità climatica.
21.     i1 AS SINGLE '[mm/ora] intensità di pioggia.
22.     i2 AS SINGLE '[mm/ora] intensità di pioggia secondo lo ietogramma Chicago o costante.
23.     DH AS SINGLE '[mm] pioggia in ogni istante temporale dt.
24.     H2 AS SINGLE '[mm] pioggia cumulata ad ogni istante temporale.
25.     he AS SINGLE '[mm] piggia efficace cumulata ad ogni istante temporale.
26.     Dhe AS SINGLE '[mm] pioggia efficace ad ogni istante temporale.
27.     MOCKUStSUta AS SINGLE
28.     MOCKUSqSUqp AS SINGLE
29.     tmSUta AS SINGLE '[-].
30.     qmSUqp AS SINGLE '[-] il valore di qm/qp corrispondente a tm/ta e' = al valore di q/qp che corrisponde a t/ta=tm/ta.
31.     qm AS SINGLE '[mc/s].
32. END TYPE
33.  
34. TYPE coefficienti
35.     d AS SINGLE
36.     tSUta AS SINGLE
37. END TYPE
38.  
39. TYPE idrogramma
40.     ore AS SINGLE
41.     minuti AS SINGLE
42.     portata AS SINGLE
43. END TYPE
44.  
45. CONST limit = 200 '50
46.  
47. DIM SHARED mockus(50) AS mockus
48. DIM SHARED dt!(2, 24)
49. DIM SHARED matrice1(2, 24, 50, 1) AS matrice1
50. DIM SHARED coefficienti(24) AS coefficienti
51. DIM SHARED matrice2!(2, 24, 50, 50)
52. DIM SHARED idrogramma(2, 24, 50, 1) AS idrogramma
53. DIM SHARED massimi1(2, 24, 1) AS idrogramma
54. DIM SHARED massimi2(2) AS idrogramma
55. DIM SHARED i%, n%, z%, ieto%
56. DIM SHARED m%(2, 24)
57. DIM SHARED S2! '[mm] contenuto idrico massimo del terreno.
58. DIM SHARED Ia! '[mm] perdite iniziali.
59. DIM SHARED tc! '[ore] tempo di corrivazione
60. DIM SHARED ta! '[ore] tempo di picco.
61. DIM SHARED qp! '[mc/s] portata al colmo dell'idrogramma unitario.
62. DIM SHARED k! '[-] coeff. della curva di possibilità climatica tado dal sito della Regione.
63. DIM SHARED a2! '[mm/d^n] coeff. della curva di possibilità climatica.
64. DIM SHARED n1! '[-] coeff. della curva di possibilità climatica.
65. DIM SHARED visualizzaieto%
66. DIM SHARED schermo&
67. DIM A1! '[mq] superficie bacino idrografico.
68. DIM L! '[m] lunghezza asta principale.
69. DIM HmaxL! '[m] punto più alto dell'astra principale.
70. DIM HminL! '[m] punto più basso dell'asta principale.
71. DIM s1! '[-] pendenza media del bacino non espressa in percentuale ma come 0,x.
72. DIM CNII! '[-] tabellato.
73. DIM CNIII '[-] discende da CNII.
74. DIM tl! '[ore] tempo di ritardo.
75. DIM keypress$
76.  
77. schermo& = _NEWIMAGE(_DESKTOPWIDTH, _DESKTOPHEIGHT, 32)
78. SCREEN schermo&
79. _SCREENMOVE 0, 0
80.  
81. A1! = 2641902
82. k! = 2.34
83. a2! = 14.6653
84. n1! = 0.45596
85. L! = 3300
86. HmaxL! = 1981
87. HminL! = 1062
88. s1! = 0.415
89. CNII! = 70
90.  
91. 'INPUT A1!
92. 'INPUT k!
93. 'INPUT a2!
94. 'INPUT n1!
95. 'INPUT L!
96. 'INPUT HmaxL!
97. 'INPUT HminL!
98. 'INPUT s1!
99. 'INPUT CNII!
100. 'INPUT d1!
101.  
102. CNIII! = (23 * CNII!) / (10 + 0.13 * CNII!)
103. tl! = 0.342 * ((L! / 1000) ^ 0.8 / (100 * s1!) ^ 0.5) * (1000 / CNIII! - 9) ^ 0.7 'formula di Mockus.
104. S2! = 25.4 * (1000 / CNIII! - 10)
105. Ia! = 0.1 * S2! 'coeff.=0.03-0.2.
106. tc! = tl! / 0.6
107. ta! = tl! / 0.9
108. qp! = 0.208 * ((A1! / 1000000) / ta!)
109.  
110. RESTORE mockus
111. FOR i% = 1 TO 100
112.     IF i% <= 50 THEN READ mockus(i%).MOCKUStSUta
113.     IF i% > 50 THEN READ mockus(i% - 50).MOCKUSqSUqp
114. NEXT i%
115.  
116. RESTORE coefficienti
117. FOR i% = 1 TO 48
118.     IF i% <= 24 THEN READ coefficienti(i%).d: coefficienti(i%).d = coefficienti(i%).d * tc!
119.     IF i% > 24 THEN READ coefficienti(i% - 24).tSUta
120. NEXT i%
121.  
122. CLS
123. PRINT "Ietogramma Chicago/costante: 1/2> "
124. 'DO
125. '_LIMIT 60 'limit
126. 'keypress$ = INKEY$
127. 'LOOP UNTIL keypress$ = "1" OR keypress$ = "2"
128. visualizzaieto% = 1 'VAL(keypress$)
129.  
130. FOR ieto% = 1 TO 2
131.     CALL creamatrici
132. NEXT ieto%
133.  
134. CLS
135. CALL grafici(massimi1(visualizzaieto%, 24, 1).ore, massimi2(visualizzaieto%).ore, massimi2(visualizzaieto%).portata) 'valori passati per disegnare il grafico.
136. SLEEP
137. SYSTEM
138. 'CALL risultati
139.  
140. '-----------------------------------------matrice1---------------------------------------------------------------------------------------------------------------------------------
141. SUB creamatrici ()
142.     z% = 0
143.     DO
144.         z% = z% + 1
145.         dt!(ieto%, z%) = coefficienti(z%).tSUta * ta!
146.         m%(ieto%, z%) = coefficienti(z%).d / (coefficienti(z%).tSUta * ta!) 'determina il numero di passi temporali necessari ad arrivare allla durata d considerata con il paso temporale dt!.
147.         FOR i% = 1 TO 50
148.             matrice1(ieto%, z%, i%, 1).MOCKUStSUta = mockus(i%).MOCKUStSUta
149.             matrice1(ieto%, z%, i%, 1).MOCKUSqSUqp = mockus(i%).MOCKUSqSUqp
150.             matrice1(ieto%, z%, i%, 1).m = i%
151.             matrice1(ieto%, z%, i%, 1).mdt = matrice1(ieto%, z%, i%, 1).m * dt!(ieto%, z%)
152.             matrice1(ieto%, z%, i%, 1).tminuti = matrice1(ieto%, z%, i%, 1).mdt * 60
153.             SELECT CASE (_ROUND(10 * (i% * coefficienti(z%).tSUta))) / 10 'crea il vettore dell'idrogramma unitario di Mockus con i multipli del coeff. z% considerato (0.1,0.2,0.3,0.5,1).
154.                 CASE IS <= 5 'valore massimo di t/ta di Mockus
155.                     matrice1(ieto%, z%, i%, 1).tmSUta = (_ROUND(10 * (i% * coefficienti(z%).tSUta))) / 10 '=matrice1(i%).m *  matrice1(i%).mdt / ta!
156.                 CASE ELSE
157.                     matrice1(ieto%, z%, i%, 1).tmSUta = 0
158.             END SELECT
159.         NEXT i%
160.         FOR i% = 1 TO 50
161.             n% = 1
162.             DO 'crea il vettore dell'idrogramma unitario Mockus prendendo i valori corrispondenti di ai coefficienti.
163.                 IF matrice1(ieto%, z%, i%, 1).tmSUta = matrice1(ieto%, z%, n%, 1).MOCKUStSUta THEN
164.                     matrice1(ieto%, z%, i%, 1).qmSUqp = matrice1(ieto%, z%, n%, 1).MOCKUSqSUqp
165.                     EXIT DO
166.                 ELSE
167.                     n% = n% + 1
168.                 END IF
169.             LOOP UNTIL n% = 50
170.             matrice1(ieto%, z%, i%, 1).qm = matrice1(ieto%, z%, i%, 1).qmSUqp * qp!
171.         NEXT i%
172.         FOR i% = 1 TO m%(ieto%, z%)
173.             matrice1(ieto%, z%, i%, 1).h1 = k! * a2! * matrice1(ieto%, z%, i%, 1).mdt ^ n1!
174.             IF i% = 1 THEN
175.                 matrice1(ieto%, z%, i%, 1).i1 = matrice1(ieto%, z%, i%, 1).h1 / dt!(ieto%, z%)
176.             ELSE
177.                 matrice1(ieto%, z%, i%, 1).i1 = (matrice1(ieto%, z%, i%, 1).h1 - matrice1(ieto%, z%, i% - 1, 1).h1) / dt!(ieto%, z%)
178.             END IF
179.         NEXT i%
180.         SELECT CASE ieto%
181.             CASE 1 'Chicago.
182.                 CALL chicago
183.             CASE 2 'costante
184.                 FOR i% = 1 TO m%(ieto%, z%)
185.                     matrice1(ieto%, z%, i%, 1).i2 = matrice1(ieto%, z%, m%(ieto%, z%), 1).h1 / matrice1(ieto%, z%, m%(ieto%, z%), 1).mdt
186.                 NEXT i%
187.         END SELECT
188.         FOR i% = 1 TO m%(ieto%, z%)
189.             matrice1(ieto%, z%, i%, 1).DH = matrice1(ieto%, z%, i%, 1).i2 * dt!(ieto%, z%)
190.             IF i% = 1 THEN
191.                 matrice1(ieto%, z%, i%, 1).H2 = matrice1(ieto%, z%, i%, 1).DH
192.             ELSE
193.                 matrice1(ieto%, z%, i%, 1).H2 = matrice1(ieto%, z%, i%, 1).DH + matrice1(ieto%, z%, i% - 1, 1).H2
194.             END IF
195.             SELECT CASE matrice1(ieto%, z%, i%, 1).H2 - Ia!
196.                 CASE IS < 0
197.                     matrice1(ieto%, z%, i%, 1).he = 0
198.                 CASE IS >= 0
199.                     matrice1(ieto%, z%, i%, 1).he = (matrice1(ieto%, z%, i%, 1).H2 - Ia!) ^ 2 / (matrice1(ieto%, z%, i%, 1).H2 - Ia! + S2!)
200.             END SELECT
201.             IF i% = 1 THEN
202.                 matrice1(ieto%, z%, i%, 1).Dhe = matrice1(ieto%, z%, i%, 1).he
203.             ELSE
204.                 matrice1(ieto%, z%, i%, 1).Dhe = matrice1(ieto%, z%, i%, 1).he - matrice1(ieto%, z%, i% - 1, 1).he
205.             END IF
206.         NEXT i%
207.  
208.         '-------------------------------matrice2-------------------------------------------------------------------------------------------------------------------------
209.  
210.         FOR n% = 1 TO 50 'colonna
211.             FOR i% = n% TO 50 ' riga
212.                 matrice2!(ieto%, z%, i%, n%) = matrice1(ieto%, z%, i% - n% + 1, 1).Dhe * matrice1(ieto%, z%, n%, 1).qm
213.             NEXT i%
214.         NEXT n%
215.         '-------------------------------idrogramma------------------------------------------------------------------------------------------------------------------
216.         FOR i% = 1 TO 50 'riga
217.             FOR n% = 1 TO 50 'colonna
218.                 idrogramma(ieto%, z%, i%, 1).portata = idrogramma(ieto%, z%, i%, 1).portata + matrice2!(ieto%, z%, i%, n%)
219.  
220.             NEXT n%
221.             idrogramma(ieto%, z%, i%, 1).ore = matrice1(ieto%, z%, i%, 1).mdt
222.             idrogramma(ieto%, z%, i%, 1).minuti = matrice1(ieto%, z%, i%, 1).tminuti
223.         NEXT i%
224.         '-------------------------------Massimi-----------------------------------------------------
225.         i% = 1
226.         'massimi1(ieto%, z%, 1).portata = idrogramma(ieto%, z%, 1, 1).portata
227.         DO
228.             IF massimi1(ieto%, z%, 1).portata < idrogramma(ieto%, z%, i%, 1).portata THEN
229.                 massimi1(ieto%, z%, 1).portata = idrogramma(ieto%, z%, i%, 1).portata
230.                 massimi1(ieto%, z%, 1).ore = idrogramma(ieto%, z%, i%, 1).ore
231.                 massimi1(ieto%, z%, 1).minuti = idrogramma(ieto%, z%, i%, 1).minuti
232.             END IF
233.             i% = i% + 1
234.         LOOP UNTIL i% = 50
235.         ' IF massimi2(ieto%).ore < massimi1(visualizzaieto%, z%, 1).ore THEN massimi2(ieto%).ore = massimi1(visualizzaieto%, z%, 1).ore
236.         IF massimi2(ieto%).portata < massimi1(visualizzaieto%, z%, 1).portata THEN
237.             massimi2(ieto%).portata = massimi1(visualizzaieto%, z%, 1).portata
238.             massimi2(ieto%).ore = massimi1(visualizzaieto%, z%, 1).ore
239.         END IF
240.     LOOP UNTIL z% = UBOUND(coefficienti)
241. END SUB
242. '---------------------------------------------------------------------------------------------------------------------------------------------------
243. 'Questa subroutine crea gli ietogrammi chicago per ogni iterazione z%.
244. 'Ad ogni nuovo dt! pone l'ultimo valore di pioggia come ultimo, il penultimo come primo, il terzultimo come penultimo, il quartultimo come secondo e così via.
245. 'Quando vi è un valore di z% che ha un dt! non nuovo, cerca la precedente iterazione z% con lo stesso dt! e copia i suoi valori dello ietogramma chicago in quello dell'iterazione z% in corso. I valori rimanenti sono riempiti
246. 'procedendo come prima.
247. 'Quindi per esempio l'iterazione z%=1 corrisponde a d=1tc con passo temporale dt!. L'iterazione z%=2 corrisponde a d=2tc con il medesimo passo temporale. Per cui lo ietogramma chicago per d=2tc e' dato dallo stesso ietogramma
248. 'chicago che si ha per d=1tc, con l'aggiunta, ai suoi estremi, dei nuovi valori. Il meccanismo si ripete per z%=3, che corrisponde a d=3tc, con il medesimo dt!. Per z%=4, corrispondente a d=4tc, invece, il passo temporale dt!
249. 'cambia.
250. 'Quanto detto e' necessario per evitare piccole discrepanze che potrebbero comportare il fatto che alcuni massimi risultino localmente poco piu' bassi del precedente, anche quando il trend generale e' chiaramente di aumento dei
251. 'massimi all'aumentare della durata della pioggia. In questo modo si fa l'ipotesi che all'aumentare della durata della pioggia, discretizzata con uno stesso passo temporale dt!, ogni ietogramma contenga a sua volta il precedente.
252. 'In tal modo, se si condiera per esempio d=2tc, si assume che l'altezza di pioggia corrispondente a 1tc, segua lo ietogramma relativo a d=1tc, cosi' che l'altezza di pioggia rimamente sia eclusivamente dovuta ai nuovi valori.
253. SUB chicago ()
254.     DIM test%
255.     DIM k%
256.     n% = 1
257.     DO
258.         IF z% = 1 THEN test% = 0: EXIT DO 'permette il calcolo dello ietogramma chicago della prima iterazione. test%=0
259.         IF dt!(ieto%, z%) = dt!(ieto%, z% - n%) THEN ' se trova un'iterazione z% precedente con lo stesso dt! di quella in corso, esegue i comandi successivi
260.             test% = 1
261.             k% = m%(ieto%, z%)
262.             i% = 0
263.             DO 'inserisce i nuovi valori dello ietogramma chicago dell'iterazione z% in corso.
264.                 matrice1(ieto%, z%, m%(ieto%, z%) - i%, 1).i2 = matrice1(ieto%, z%, k%, 1).i1
265.                 k% = k% - 1
266.                 IF k% = m%(ieto%, z% - n%) THEN EXIT DO 'se ha completato l'inserimento dei nuovi valori, passa oltre.
267.                 matrice1(ieto%, z%, i% + 1, 1).i2 = matrice1(ieto%, z%, k%, 1).i1
268.                 k% = k% - 1
269.                 i% = i% + 1
270.             LOOP UNTIL k% = m%(ieto%, z% - n%)
271.             IF matrice1(ieto%, z%, i%, 1).i2 = 0 THEN i% = i% - 1
272.             k% = 1
273.             DO 'inserisce lo ietogramma chicago dell'iterazione z% precedente con lo stesso dt!, nello ietogramma dell'iterazione z% in corso.
274.                 i% = i% + 1
275.                 matrice1(ieto%, z%, i%, 1).i2 = matrice1(ieto%, z% - n%, k%, 1).i2
276.                 k% = k% + 1
277.             LOOP UNTIL k% > m%(ieto%, z% - n%)
278.             EXIT DO 'siccome ha completato lo ietogramma dell'iterazione z% in corso, non serve incrementare n%, quindi passa oltre. test%=1.
279.         ELSE ' se non trova l'iterazione z% precedente con lo stesso dt! di quella in corso, incrementa n% e ripete la ricerca.
280.             test% = 0
281.             n% = n% + 1
282.         END IF
283.     LOOP UNTIL n% = z% - 1 'se non ha trovato iterazioni z% precedenti con lo stesso dt!, allora passa oltre e test%=0.
284.     SELECT CASE test%
285.         CASE IS = 1 'esce dalla subroutine, in quanto lo ietogramma e' definito.
286.             EXIT SUB
287.         CASE IS = 0 'calcola ex-novo lo ietogramma relativo all'iterazione z% in corso in quanto non vi e' un'iterazioni z% precedente con lo stesso dt! dell'iterazione z% in corso.
288.             k% = m%(ieto%, z%)
289.             i% = 0
290.             DO
291.                 matrice1(ieto%, z%, m%(ieto%, z%) - i%, 1).i2 = matrice1(ieto%, z%, k%, 1).i1
292.                 k% = k% - 1
293.                 IF k% = 0 THEN EXIT DO
294.                 matrice1(ieto%, z%, i% + 1, 1).i2 = matrice1(ieto%, z%, k%, 1).i1
295.                 k% = k% - 1
296.                 i% = i% + 1
297.             LOOP UNTIL k% = 0
298.     END SELECT
299. END SUB
300. '-----------------------------------------------------------------------------------------------------------------------------------------------------------------
301. SUB grafici (ore!, oreM!, QM!)
302.     'ore!: ore corrispondenti al massimo maggiore.
303.     'oreM!: ore corrispondenti alla portata massima.
304.     'QM!: portata massima.
305.  
306.     CONST R& = _RGB32(255, 0, 0)
307.     CONST G& = _RGB32(0, 255, 0)
308.     CONST B& = _RGB32(0, 0, 255)
309.     CONST bianco& = _RGB32(255, 255, 255)
310.     CONST giallo& = _RGB32(255, 255, 0)
311.     CONST grigio& = _RGB32(127, 127, 127)
312.  
313.     DIM L%, H% 'dimensioni del grafico.
314.     SHARED grafico&
315.     DIM graficoNERO&, graficoCOPIA&
316.     DIM colore&
317.  
318.     DIM dx%, dy% 'origine degli assi dentro SCREEN.
319.     DIM dx!, dy! 'posizione dell'origine degli assi dentro WINDOW.
320.     DIM taccaX%, taccaY%
321.     DIM portata!(2)
322.     '  DIM soglia!
323.     REDIM x%(1), y%(1) 'conversione nelle coordinate di SCREEN delle coordinate dentro WINDOW corrispondenti a determinati valori.
324.  
325.     H% = _DESKTOPHEIGHT * 2 / 3 ' imposta l'altezza dello SCREEN pari ai 2/3 dello schermo visualizzabile, cosi' che su ogni PC sia visualizzabile.
326.     L% = H% * 1.62 'imposta la larghezza dello SCREEN in modo che L%xH% sia un rettangolo aureo.
327.     dx% = 39
328.     dy% = H% - 1 - dx%
329.     taccaX% = ore! \ 5 ': PRINT taccaX%: SLEEP
330.     taccaY% = QM! \ 5
331.     grafico& = _NEWIMAGE(L%, H%, 32)
332.     graficoNERO& = _NEWIMAGE(L%, H%, 32) 'crea un'immagine della stessa dimensione di grafico& che sarà colorata di nero.
333.  
334.     _DEST graficoNERO&
335.     CLS ' colora di nero l'immagine graficoNERO&.
336.     _DEST grafico&
337.     LINE (0, 0)-(L% - 1, H% - 1), R, B
338.     WINDOW (0, 0)-(ore! * 1.15, QM! * 1.15) 'scala ascisse e ordinate in base alle ore e alla portata, estendendole del 15%.
339.     dx! = PMAP(dx%, 2)
340.     dy! = PMAP(dy%, 3)
341.     LINE (dx!, QM! * 1.1)-(dx!, dy!): LINE -(ore! * 1.1, dy!), bianco& 'disegna gli assi in modo che si estendano del 10% oltre valori massimi, ma senza arrivare alla massima estensione delle coordinate impostate da WINDOW: cosi' si ha che gli assi non toccano la cornice rossa e i valori massimi che non arrivano al limite degli assi.
342.     PSET (dx!, QM! * 1.1!), bianco&: DRAW "F20": PSET (dx!, QM! * 1.1), bianco&: DRAW "G20"
343.     PSET (ore! * 1.1, dy!), bianco&: DRAW "G20": PSET (ore! * 1.1, dy!), bianco&: DRAW "H20"
344.     _PRINTSTRING (dx%, PMAP(QM! * 1.1, 1) - dx% \ 2), "mc/s"
345.     _PRINTSTRING (PMAP(ore! * 1.1, 0), dy%), "ore"
346.     _PRINTSTRING (dx% \ 4, dy% - 7), "0"
347.     i% = 1
348.     WHILE i% * taccaX% <= ore!
349.         LINE (dx! + i% * taccaX%, PMAP(dy% + 5, 3))-(dx! + i% * taccaX%, PMAP(dy% - 5, 3))
350.         x%(1) = PMAP(dx! + i% * taccaX%, 0)
351.         _PRINTSTRING (x%(1) - 4 * LEN(_TRIM$(STR$(i% * taccaX%))), H% - 1 - dx% / 1.5), _TRIM$(STR$(i% * taccaX%))
352.         i% = i% + 1
353.     WEND
354.     i% = 1
355.     WHILE i% * taccaY% <= QM!
356.         LINE (PMAP(dx% - 5, 2), dy! + i% * taccaY%)-(PMAP(dx% + 5, 2), dy! + i% * taccaY%)
357.         y%(1) = PMAP(dy! + i% * taccaY%, 1)
358.         _PRINTSTRING (dx% \ 4, y%(1) - 7), _TRIM$(STR$(i% * taccaY%))
359.         i% = i% + 1
360.     WEND
361.     FOR z% = 1 TO 24
362.         SELECT CASE z% MOD 3
363.             CASE IS = 1
364.                 colore& = R&
365.             CASE IS = 2
366.                 colore& = G&
367.             CASE IS = 0
368.                 colore& = B&
369.         END SELECT
370.         CIRCLE (dx! + massimi1(visualizzaieto%, z%, 1).ore, dy! + massimi1(visualizzaieto%, z%, 1).portata), PMAP(3, 2), R&
371.         PAINT (dx! + massimi1(visualizzaieto%, z%, 1).ore + PMAP(0.5, 2), dy! + massimi1(visualizzaieto%, z%, 1).portata + PMAP(0.5, 2)), R&
372.         FOR i% = 1 TO 50
373.             '            _LIMIT limit
374.             IF i% = 1 THEN
375.                 LINE (dx!, dy!)-(dx! + idrogramma(visualizzaieto%, z%, i%, 1).ore, dy! + idrogramma(visualizzaieto%, z%, i%, 1).portata), colore&
376.             ELSE
377.                 LINE (dx! + idrogramma(visualizzaieto%, z%, i% - 1, 1).ore, dy! + idrogramma(visualizzaieto%, z%, i% - 1, 1).portata)-(dx! + idrogramma(visualizzaieto%, z%, i%, 1).ore, dy! + idrogramma(visualizzaieto%, z%, i%, 1).portata), colore&
378.             END IF
379.             '            _PUTIMAGE ((_DESKTOPWIDTH - L%) \ 2, (_DESKTOPHEIGHT - H%) \ 2), grafico&, schermo&
380.         NEXT i%
381.         IF z% = 1 THEN
382.             LINE (dx!, dy!)-(dx! + massimi1(visualizzaieto%, z%, 1).ore, dy! + massimi1(visualizzaieto%, z%, 1).portata), giallo&
383.         ELSE
384.             LINE (dx! + massimi1(visualizzaieto%, z% - 1, 1).ore, dy! + massimi1(visualizzaieto%, z% - 1, 1).portata)-(dx! + massimi1(visualizzaieto%, z%, 1).ore, dy! + massimi1(visualizzaieto%, z%, 1).portata), giallo&
385.         END IF
386.     NEXT z%
387.     PSET (dx!, dy!), giallo&
388.     LINE (dx!, dy! + QM!)-(dx! + oreM!, dy! + QM!), R&, , 65520
389.     LINE -(dx! + oreM!, dy!), R&, , 65520
390.     x%(1) = PMAP(dx! + oreM!, 0)
391.     y%(1) = PMAP(dy! + QM!, 1)
392.     _PRINTSTRING (x%(1) + 10, y%(1) - dx%), "(" + _TRIM$(STR$(oreM!)) + ","
393.     _PRINTSTRING (x%(1) + 10, y%(1) - dx% + 16), _TRIM$(STR$(QM!)) + ")"
394.     _PUTIMAGE ((_DESKTOPWIDTH - L%) \ 2, (_DESKTOPHEIGHT - H%) \ 2), grafico&, schermo&
395.     graficoCOPIA& = _COPYIMAGE(grafico&) 'crea una copia di grafico& prima che venga effettuato il tratteggio grigio: INSERIRE POI NUOVA RICHIESTA.
396.     CALL idrogrammaprogetto(L%, H%, grafico&, graficoCOPIA&, graficoNERO&, portata!(), x%(), y%(), dx!, dy!, giallo&)
397.     SLEEP
398.  
399.     SYSTEM
400.  
401.     DO
402.         '  CLS 2
403.         PRINT "Digitare le ore e i minuti (<= "; _TRIM$(STR$(INT(ore!))); " ore e "; _TRIM$(STR$(INT((ore! - INT(ore!)) * 60))); " minuti) di cui si vuole conoscere la portata di picco."
404.         INPUT "Ore: ", ore2!
405.         INPUT "Minuti: ", minuti%
406.         ore2! = ore2! + minuti% / 60
407.     LOOP WHILE ore2! > ore!
408.     _DEST grafico&
409.     _SOURCE grafico&
410.     y%(1) = PMAP(dy!, 1)
411.     i% = 2
412.     DO
413.         REDIM _PRESERVE y%(i%)
414.         y%(i%) = y%(i% - 1) - 1
415.         IF POINT(dx! + ore2!, PMAP(y%(i% - 1), 3)) = giallo& THEN EXIT DO
416.         i% = i% + 1
417.     LOOP
418.     REDIM _PRESERVE y%(1 TO UBOUND(y%) - 1)
419.     portata!(2) = PMAP(y%(UBOUND(y%)), 3) - dy!
420.     FOR i% = 1 TO UBOUND(y%)
421.         LINE (dx! + ore2!, dy!)-(dx! + ore2!, PMAP(y%(i%), 3)), grigio&, , 65280
422.         _PUTIMAGE ((_DESKTOPWIDTH - L%) \ 2, (_DESKTOPHEIGHT - H%) \ 2), grafico&, schermo&
423.     NEXT i%
424.     x%(1) = PMAP(dx! + ore2!, 0)
425.     i% = 2
426.     DO
427.         REDIM _PRESERVE x%(i%)
428.         x%(i%) = x%(i% - 1) - 1
429.         IF POINT(PMAP(x%(i% - 1), 2), PMAP(y%(UBOUND(y%)), 3)) = bianco& THEN EXIT DO
430.         i% = i% + 1
431.     LOOP
432.     REDIM _PRESERVE x%(1 TO UBOUND(x%) - 1)
433.     FOR i% = 1 TO UBOUND(x%)
434.         LINE (dx! + ore2!, PMAP(y%(UBOUND(y%)), 3))-(PMAP(x%(i%), 2), PMAP(y%(UBOUND(y%)), 3)), grigio&, , 65280
435.         _PUTIMAGE ((_DESKTOPWIDTH - L%) \ 2, (_DESKTOPHEIGHT - H%) \ 2), grafico&, schermo&
436.     NEXT i%
437.     _DEST schermo&
438.     PRINT "Portata:"; portata!(2); "mc/s"
439.     _FREEIMAGE grafico&
440.     _FREEIMAGE graficoNERO&
441.     _FREEIMAGE graficoCOPIA&
442. END SUB
443. '------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
444. SUB idrogrammaprogetto (L%, H%, grafico&, graficoCOPIA&, graficoNERO&, portata!(), x%(), y%(), dx!, dy!, giallo&)
445.  
446.     DIM ore2!
447.     DIM soglia!
448.  
449.     _DEST schermo&
450.     VIEW PRINT 1 TO 10
451.     SELECT CASE visualizzaieto%
452.         CASE IS = 1
453.             '   INPUT "Inserire percentuale soglia per portata progetto (numero intero senza simbolo %: ", soglia!
454.             soglia! = 0.1
455.             grafico& = _COPYIMAGE(graficoCOPIA&) ' grafico& torna ad ed essere senza il tratteggio grigio.
456.             _PUTIMAGE ((_DESKTOPWIDTH - L%) \ 2, (_DESKTOPHEIGHT - H%) \ 2), graficoNERO&, schermo& 'incollando grafico& (senza il tratteggio) su schermo&, si vedrebbe comunque il tratteggio in quanto la precedente versione di grafico& (che era con il tratteggio), era stata incollata su schermo&. Quindi prima incollo l'immagine graficoNERO& su schermo&.
457.             _PUTIMAGE ((_DESKTOPWIDTH - L%) \ 2, (_DESKTOPHEIGHT - H%) \ 2), grafico&, schermo& 'copia grafico& (senza tratteggio) su schermo&, che era stato ripulito alla linea precedente.
458.             _DEST grafico&
459.             _SOURCE grafico&
460.             i% = 0
461.             ore2! = 0
462.             REDIM x%(1), y%(1)
463.             y%(1) = PMAP(dy!, 1)
464.             DO
465.                 i% = i% + 1
466.                 ore2! = ore2! + 1
467.                 REDIM _PRESERVE y%(i%)
468.                 DO
469.                     IF POINT(dx! + ore2!, PMAP(y%(i%), 3)) = giallo& THEN EXIT DO
470.                     y%(i%) = y%(i%) - 1
471.                     'LINE (dx! + ore2!, PMAP(y%(i%), 3))-(dx! + ore2!, PMAP(y%(i%), 3)): SLEEP
472.                 LOOP
473.                 IF i% > 1 THEN
474.                     IF PMAP(y%(i%), 3) - dy! < (1 + soglia!) * (PMAP(y%(i% - 1), 3) - dy!) THEN EXIT DO
475.                 END IF
476.             LOOP
477.             ore2! = ore2! - 1
478.             REDIM _PRESERVE y%(1 TO UBOUND(y%) - 1)
479.             portata!(1) = PMAP(y%(UBOUND(y%)), 3) - dy!
480.             _DEST schermo&
481.             PRINT "Portata progetto:"; portata!(1); "mc/s, all'ora n."; ore2!
482.         CASE IS = 2
483.             PRINT "La portata di progetto e' quella visualizzata."
484.     END SELECT
485. END SUB
486. '------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
487. SUB risultati ()
488.     FOR z% = 1 TO 24
489.         CLS
490.         PRINT "z%="; z%
491.         PRINT "dt="; dt!(visualizzaieto%, z%); ";"; dt!(visualizzaieto%, z%) * 60
492.         PRINT "visualizzaieto%="; visualizzaieto%
493.         PRINT "m%="; m%(visualizzaieto%, z%)
494.         PRINT "N="; coefficienti(z%).d / tc!
495.         PRINT "n="; coefficienti(z%).tSUta
496.         FOR i% = 1 TO 50
497.             '  _LIMIT limit
498.     PRINT _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).m)), _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).mdt)), _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).tminuti)), _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).h1)),_
499.      _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).i1)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).i2)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).DH)), _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).H2)),_
500.       _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).he)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).Dhe)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).MOCKUStSUta)),_
501.        _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).MOCKUSqSUqp)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).tmSUta)), _TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).qmSUqp)),_TRIM$(STR$(matrice1(visualizzaieto%,z%,i%, 1).qm))
502.         NEXT i%
503.         SLEEP
504.         CLS
505.         FOR n% = 1 TO 50
506.             '  _LIMIT limit
507.             LOCATE 1,
508.             FOR i% = 1 TO 50
509.                 LOCATE , n% * 4 - 3
510.                 PRINT matrice2!(visualizzaieto%, z%, i%, n%)
511.             NEXT i%
512.         NEXT n%
513.         SLEEP
514.     NEXT z%
515.     CLS
516.     FOR n% = 1 TO 24
517.         ' _LIMIT limit
518.         SELECT CASE n%
519.             CASE IS <= 16
520.                 FOR i% = 1 TO 50
521.                     LOCATE i%, n% * 15 - 14
522.                     PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).portata))
523.                     'PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).ore))
524.                     'PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).minuti))
525.  
526.                 NEXT i%
527.                 IF n% = 16 THEN SLEEP: CLS
528.             CASE IS > 16
529.                 FOR i% = 1 TO 50
530.                     LOCATE i%, (n% - 16) * 15 - 14
531.                     PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).portata))
532.                     'PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).ore))
533.                     'PRINT _TRIM$(STR$(idrogramma(visualizzaieto%, n%, i%, 1).minuti))
534.                 NEXT i%
535.         END SELECT
536.     NEXT n%
537.     SLEEP
538.     CLS
539.     FOR z% = 1 TO 24
540.         ' _LIMIT limit
541.         PRINT massimi1(visualizzaieto%, z%, 1).ore, massimi1(visualizzaieto%, z%, 1).minuti, massimi1(visualizzaieto%, z%, 1).portata
542.     NEXT z%
543.     SLEEP
544.     PRINT
545.     PRINT massimi2(visualizzaieto%).ore; massimi2(visualizzaieto%).portata
546. END SUB
547.  

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bplus on February 25, 2021, 11:19:03 am

@bartok

Well I understand your motivation to take working code and refactor it to make it better in some way. Last night I overhauled SB1 to do customized splits and only if the program line needed it, instead of one split for all whether it needs it or not and whether or not the right delimiter :p So clearly I am cutting wasted time, and redundant codes where I have to re-split with correct delimiter.

So what is reason to make a separate sub, load it up with a ton of parameters and do something on the side?
My guess is you are missing one or more values passed to the added sub routine.

Do you know, you can put GOSUBs inside your SUB's?
This would save you all the trouble of loading up another sub with a bunch of values to work out some calculation or do some task.

Doing that, using a GOSUB instead of full-fledged SUB,  all the variables in the SUB are available free of charge in the GOSUB routine and no worries if you forgot something or are in the wrong _Source or _Dest handle. Put the GOSUB at the bottom of your SUB and make sure you exit sub before reaching GOSUB in exec flow.

Code: QB64: [Select]

1. SUB idrogrammaprogetto (L%, H%, grafico&, graficoCOPIA&, graficoNERO&, portata!(), x%(), y%(), dx!, dy!, giallo&)

Is it likely you will need all these parameter / arguments or most of them in other parts of your program present or future?

If no, then go with GoSUB, instead of using another SUB.

If so, then set up a Type for passing all this stuff under a single an UDT (= User Defined Type) it will make passing this stuff off a breeze (maybe not but this > ) or better Shared the Typed variable or array.

You have the problem narrowed down to 4 or so lines in a Do Loop, something is not changing correctly in there. I would suspect a _dest or _source handle missing, but I don't want to dive into detail someone else code and find which variable was mis-typed or missing. :p  This reminds me of Turings problem, "How to get a program to realize it's in an infinite loop and Halt", I suspect it's possible with a map of program and table of variables and a ton of analysis, specially for particular programs if not for "In General" which is what Mathematicians are always after.  I recently saw Internet thing where some lady solved Conway knot that had all the mathematicians tied up, she said, (paraphased) "But it was just one knot" with a tone of what's the Big Deal?  But I digress...

I would set up my own watch by printing the variables in question at the bottom of the loop and then put an input line at bottom with prompt "Ready... press enter"; ready$ 

I label such private watches with ' Debug =============================================
so I know what to remove when problem is discovered.

Sometimes I have to watch a number of places I REDIM SHARED Debug
and turn it on or off at top of main code section
Debug = 0 or -1

Then I start private watches with lines:
If Debug then
   REDIM ready$  ' I just thought of this, private dim of ready$ variable
   print the watch stuff
   input "ready? ... press enter'; ready$
   If ready$ = enoughAlready$ then Debug = 0 ' I just thought of this! turn off debug when you've had enough
end if

You could even output this stuff in a console window so it's out of the way of your program screen, yes you can run two windows in Qb64! But I never tried.

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bartok on February 25, 2021, 01:26:30 pm

bplus,

Quote

So what is reason to make a separate sub, load it up with a ton of parameters and do something on the side?

it is a philosophical question. In general, if I have an idea to do something, I enter a DO cicle until I have found the way to do it. Regardless of the original reason of the idea, the fact is that if I don't solve the problem, I go in loop, and I am too annoyed to go ahead.

Quote

Is it likely you will need all these parameter / arguments or most of them in other parts of your program present or future?

No, I don't. This sub have to be used only at this point of the code. However, the idea to do a separate sub is intended to make a more readable code. In fact, at the moment, I am just building a kind of structure of the program I have in mind, that goes from the start to the end, doing all have to be done, without possibility to navigate into the program or exit it. But finally, I want to implement this capacity. So, for example, at the moment the program flow passes this subroutine only one time and goes ahead, but in the future I want to make the possibility to repeate the procedure if the user want it. So, finally, there will be probably a DO-LOOP cicle or a IF condition. Well, with the simply CALL of the subroutine, the DO-LOOP cicle will be easy to understand what is intended to do. Without the subroutine, the DO-LOOP cicle will be very long and plenty of stuff an less easy to read. So, if I write a part of code that does something very definite, I want that it is a subroutine, so that, at the end, the code will be, as far as possible, a sequence of explicit commands that are evident what are intended to do, without the necessity to read and understand the pure code. I know that the same result can be reached with REM statements, but if the code is self-explaining, is better.

Quote

I would set up my own watch by printing the variables in question at the bottom of the loop and then put an input line at bottom with prompt "Ready... press enter"; ready$ 

I do something similar: I use PRINT (some variable) and SLEEP. I also use LINE or PSET. For example, I know that the problem is not a matter of coordinates passing from the sub-routine to the sub-subroutine, because in the LOOP cicle I have put the command to draw a pixel as the count goes ahead, interspersed by a SLEEP. And in fact the pixels were drawn exactly. Finally, the line of all these pixels met a yellow pixel of the yellow line. Using this method debugging CODE 1 (the one working), when the line met the yellow pixel, nothing happened, because PSET, couloring each pixel at time, never met a yellow pixel. Removing PSET, CODE 1 worked, CODE 2 not, but in this way was clear that the problem wasn't the coordinates.

Quote

Do you know, you can put GOSUBs inside your SUB's?

no, I didn't!

Quote

If no, then go with GoSUB, instead of using another SUB.

Your advice to use GOSUB, is one more time very usefull, because, in this way, I have solved the problem! And the solution is better of the original idea to use a full-fledged SUB for the reasons you said. Thanks you!

Title: Re: Problem to print values on the axis inside and outside WINDOW.
Post by: bplus on February 25, 2021, 10:39:09 pm

@bartok  glad to see :)

Text Only | Text with Attachments