« previous next »
Pages: 1 [2]   Go Down

Author Topic: Problem to print values on the axis inside and outside WINDOW.  (Read 5123 times)

0 Members and 1 Guest are viewing this topic.

Offline bartok

  • Newbie
  • Posts: 80
    • View Profile
Re: Problem to print values on the axis inside and outside WINDOW.
« Reply #15 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.
Logged

Offline bplus

  • Global Moderator
  • Forum Resident
  • Posts: 8053
  • b = b + ...
    • View Profile
Re: Problem to print values on the axis inside and outside WINDOW.
« Reply #16 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.
Logged

Offline bartok

  • Newbie
  • Posts: 80
    • View Profile
Re: Problem to print values on the axis inside and outside WINDOW.
« Reply #17 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.

« Last Edit: February 22, 2021, 07:48:09 am by bartok »
Logged

Offline bplus

  • Global Moderator
  • Forum Resident
  • Posts: 8053
  • b = b + ...
    • View Profile
Re: Problem to print values on the axis inside and outside WINDOW.
« Reply #18 on: February 22, 2021, 10:41:06 am »
@bartok looks good man!
Logged

Offline bartok

  • Newbie
  • Posts: 80
    • View Profile
Re: Problem to print values on the axis inside and outside WINDOW.
« Reply #19 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]
  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
  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].
  33.  
  34. TYPE coefficienti
  35.     d AS SINGLE
  36.     tSUta AS SINGLE
  38.  
  39. TYPE idrogramma
  40.     ore AS SINGLE
  41.     minuti AS SINGLE
  42.     portata AS SINGLE
  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.  
  79. SCREEN schermo&
  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.  
  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.  
  137. CALL grafici(massimi1(visualizzaieto%, 24, 1).ore, massimi2(visualizzaieto%).ore, massimi2(visualizzaieto%).portata) 'valori passati per disegnare il grafico.
  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)
  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
  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&
  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
  474.  
« Last Edit: February 22, 2021, 10:51:02 am by bartok »
Logged

Offline bplus

  • Global Moderator
  • Forum Resident
  • Posts: 8053
  • b = b + ...
    • View Profile
Re: Problem to print values on the axis inside and outside WINDOW.
« Reply #20 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.
« Last Edit: February 22, 2021, 11:35:42 am by bplus »
Logged

Offline bplus

  • Global Moderator
  • Forum Resident
  • Posts: 8053
  • b = b + ...
    • View Profile
Re: Problem to print values on the axis inside and outside WINDOW.
« Reply #21 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.
« Last Edit: February 22, 2021, 11:42:52 am by bplus »
Logged

Offline bartok

  • Newbie
  • Posts: 80
    • View Profile
Re: Problem to print values on the axis inside and outside WINDOW.
« Reply #22 on: February 22, 2021, 12:14:06 pm »
Perfect. And once again the computer was right!

Thanks!
Logged

Offline bartok

  • Newbie
  • Posts: 80
    • View Profile
Re: Problem to print values on the axis inside and outside WINDOW.
« Reply #23 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]
  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
  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].
  33.  
  34. TYPE coefficienti
  35.     d AS SINGLE
  36.     tSUta AS SINGLE
  38.  
  39. TYPE idrogramma
  40.     ore AS SINGLE
  41.     minuti AS SINGLE
  42.     portata AS SINGLE
  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.  
  78. SCREEN schermo&
  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.  
  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.  
  135. CALL grafici(massimi1(visualizzaieto%, 24, 1).ore, massimi2(visualizzaieto%).ore, massimi2(visualizzaieto%).portata) 'valori passati per disegnare il grafico.
  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)
  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
  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&
  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
  539.  

CODE 2.
Code: QB64: [Select]
  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
  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].
  33.  
  34. TYPE coefficienti
  35.     d AS SINGLE
  36.     tSUta AS SINGLE
  38.  
  39. TYPE idrogramma
  40.     ore AS SINGLE
  41.     minuti AS SINGLE
  42.     portata AS SINGLE
  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.  
  78. SCREEN schermo&
  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.  
  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.  
  135. CALL grafici(massimi1(visualizzaieto%, 24, 1).ore, massimi2(visualizzaieto%).ore, massimi2(visualizzaieto%).portata) 'valori passati per disegnare il grafico.
  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)
  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
  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&
  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
  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
  547.  

« Last Edit: February 25, 2021, 08:48:41 am by bartok »
Logged

Offline bplus

  • Global Moderator
  • Forum Resident
  • Posts: 8053
  • b = b + ...
    • View Profile
Re: Problem to print values on the axis inside and outside WINDOW.
« Reply #24 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.
« Last Edit: February 25, 2021, 11:48:13 am by bplus »
Logged

Offline bartok

  • Newbie
  • Posts: 80
    • View Profile
Re: Problem to print values on the axis inside and outside WINDOW.
« Reply #25 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!
« Last Edit: February 25, 2021, 01:30:39 pm by bartok »
Logged

Offline bplus

  • Global Moderator
  • Forum Resident
  • Posts: 8053
  • b = b + ...
    • View Profile
Re: Problem to print values on the axis inside and outside WINDOW.
« Reply #26 on: February 25, 2021, 10:39:09 pm »
@bartok  glad to see :)
Logged
Pages: 1 [2]   Go Up
« previous next »