Add GIF to README

README.md

@@ -1,7 +1,5 @@
 <!-- Title -->
-<h1 align="center">
- Subzero.jl
-</h1>
+![Subzero.jl](https://github.com/Caltech-OCTO/Subzero.jl/blob/main/docs/src/assets/title.gif)
 
 <!-- description -->
 <p align="center">
@@ -27,10 +25,10 @@
 Subzero.jl is a native [Julia](https://julialang.org/) discrete-element model (DEM) for exploring fine-scale sea ice dynamics, 
 reimplementing MATLAB model [SubZero by Manucharyan and Montemuro](https://doi.org/10.1029/2022MS003247).
 
-Subzero.jl is **fast** and **flexible**.
+**Subzero.jl is fast and flexible**.
 
 - 🚀 Runs over **35 times faster** that original MATLAB model for title simulation!
-- 🧩 Modular simulation model makes it easy to customize simulations within a **single_ run file**!
+- 🧩 Modular simulation model makes it easy to customize simulations within a **single file**!
  - Enable and disable physical processes such as fracturing, ridging, and welding
  - Choose algorithms for key processes (or add your own!)
 

docs/src/assets/title.jl

@@ -7,7 +7,7 @@ const Δgrid = 2e3
 const hmean = 0.25
 const Δh = 0.0
 const Δt = 20
-const nΔt = 5000
+const nΔt = 15000
 
 # Model instantiation
 grid = RegRectilinearGrid((0.0, Lx), (0.0, Ly), Δgrid, Δgrid)
@@ -61,14 +61,12 @@ function plot_logo(floe_fn, Lx, Ly, dir)
  time_polys = @lift(sim_polys[$time])
 
  # Set up figure
- fig = Figure(; size = (1200, 300))
+ fig = Figure(; size = (800, 200))
 
- floe_rgb = (165, 222, 242)
+ floe_rgb = (255, 255, 255)
  floe_color = RGBf((floe_rgb ./ 255)...)
- ocean_rgb = (4, 31, 74)
+ ocean_rgb = (0, 157, 196)
  ocean_color = RGBf((ocean_rgb ./ 255)...)
- text_outline_rgb = (219, 76, 0)
- text_outline_color = RGBf((text_outline_rgb ./ 255)...)
 
  ax = Axis(fig[1, 1]; limits = (0.0, Lx, 0.0, Ly), backgroundcolor = ocean_color)
  hidedecorations!(ax)
@@ -80,8 +78,8 @@ function plot_logo(floe_fn, Lx, Ly, dir)
  align = (:center, :center),
  text = "Subzero.jl",
  font = :bold,
- fontsize = 175,
- strokecolor = text_outline_color,
+ fontsize = 130,
+ strokecolor = :white,
  strokewidth = 5,
  )
 

src/physical_processes/simplification.jl

@@ -26,7 +26,7 @@ function dissolve_floe!(floe, grid::RegRectilinearGrid, domain, dissolved)
  yidx = shift_cell_idx(yidx, grid.Ny + 1, domain.north)
  # If centroid is within bounds after wrapping, add mass to dissolved
  if 0 < xidx <= grid.Nx && 0 < yidx <= grid.Ny
- dissolved[xidx, yidx] += floe.mass
+ dissolved[yidx, xidx] += floe.mass
  end
  return
 end