如何从 Julia 的 if 语句内部中断 for 循环？

Question

I am trying to translate a python modelisation (not very effective) into a Julia code (more useful in this context) but I didn't find how to break a for loop and how to print a value in this context.

Python code

for i in np.range(1, nt)
    if t[x == s] == 100:
        s_time = i*dt
        print(s_time)
        break

Julia code I tried :

for i = 1:nt
    if round.(Int, t[x.==s]) == 100
        s_time = i*dt
        println(s_time)
        break
    end
end

The variables here are not important, but dt is a value that allows me to know the time in [y] in the for loop in function of i , and t[x.==s] == 100 allows me to know when I want to break the loop.

The Julia code do not returns me any error but do not break and do not print anything.

Full code (for Jupyter Notebook users) :

# NUMERICS
ti = 30  # surface temperature
td = 1000  # dyke temperature
tel = 30  # right lim temperature
teg = tel  # left lim temperature
x1 = 80  # central position dyke 2
x2 = 110  # central position dyke 1
d1 = 24  # dyke diameter
d2 = 10  # dyke diameter
s = 160  # probe position

yd1 = 0  # dyke 1 intrusion year
yd2 = 15  # dyke 2 intrusion year

d = 40  # rock diffusivity
size = 200  # size of the model
ttot = 30;  # total time in year

nx = 2001  # cells number
dx = size/(nx-1)  # cells size
dt = dx^2/(2.1*d)  # frame condition
nt = round(Int, ttot/dt)  # frame number
x = LinRange(0, size, nx); # size model

# INITIALISATION
indexing_1_l = round(Int, (x1 - d1/2)*nx/size)
indexing_1_r = round(Int, (x1 + d1/2)*nx/size)
indexing_2_l = round(Int, (x2 - d2/2)*nx/size)
indexing_2_r = round(Int, (x2 + d2/2)*nx/size);

t = ti * ones(nx)
t[indexing_1_l:indexing_1_r] .= td;
using Plots
gr() ;

# PLOTTING LOOP
for i = 1:nt 
    
    q = -d*(diff(t)/dx)
    dcdt = -diff(q)/dx
    t[2:end-1] = t[2:end-1] + dcdt*dt
    
    t[1] = teg
    t[end] = tel
    
    if i == round(Int, yd2/dt)
        t[indexing_2_l:indexing_2_r] .= td
    end
    
    # IF STATEMENT THAT DOESN'T WORK
    if round.(Int, t[x.==s]) == 100
        s_time = i*dt
        println("$s_time")
        break
    end

    if mod(i, 1000) == 0
        IJulia.clear_output(true)
        pt = plot(x, t, label = :false, c="red", title=" $(round.(Int, i*dt)) [y] - Dyke intrusion", 
            xlabel="Distance [m]", ylabel="Temperature [°C]", ylims=(0, 1000), size = (800, 600))
        Plots.display(pt)
    end
    
end

Answer 1

This will never return true because an array cannot be equal to a number:

if round.(Int, t[x.==s]) == 100

Look at this

julia> [100] == 100
false

It's not completely clear what you want it to mean, but perhaps this

if all(round.(Int, t[x.==s]) .== 100)

or perhaps any ?

Answer 2

you have a syntax error! println(s_time)