2d:null 2e:["$","title",null,{"children":"Microclima Bahia | Vinícius Mello | Observable"}] 33:T63c,{ let display = html`

`; let map = html`

`; let mesh = buildMeshes(topo); let mapSvg = svg``; map.appendChild(mapSvg); let mapSvgSel = d3.select(mapSvg); mapSvgSel.append("path").datum(mesh.outerBoundary) .attr("class", "out-boundary") .attr("d", mesh.path); mapSvgSel.append("path").datum(mesh.divisionBoundaries) .attr("class", "in-boundary") .attr("d", mesh.path); //mapSvgSel.selectAll("circle").data(projectedSamples).enter().append("circle") // .attr("r",2).attr("cx", (d)=>d[0]).attr("cy", (d)=>d[1]).style("fill","red"); let overlay = html`

`; display.appendChild(map); overlay.appendChild(context.canvas); display.appendChild(overlay); return display; } // viewof canvas = context.canvas34:T481,context = { let context = DOM.context2d(width, height); context.canvas.width = width; context.canvas.height = height; window.t=0; let xMin = bounds.x; let xMax = bounds.x+bounds.width; let yMin = bounds.y; let yMax = bounds.y+bounds.height; function draw() { let prev = context.globalCompositeOperation; context.globalCompositeOperation = "destination-in"; context.fillRect(0, 0, width, height); context.globalCompositeOperation = prev; context.beginPath(); R.forEach((particle) => { let x = particle.x; let y = particle.y; context.moveTo(x, y); x = particle.xt; y = particle.yt; context.lineTo(x, y); particle.x = particle.xt; particle.y = particle.yt; }, particles); context.stroke(); }; context.lineWidth = 1.0; context.fillStyle = settings.fadeFillStyle; function animate() { requestAnimationFrame(animate); evolve(window.t); window.t = window.t + 0.001; if (window.t>=Time-2) window.t=0; draw(); d3.select("#time").text(`t=${window.t.toFixed(2)}`); } animate(); context.canvas.id = "canvas"; return context; } 35:[] 31:{"id":"dbf193d3486d9a4f","slug":"microclima-bahia","title":"Microclima Bahia","thumbnail":"89779e8b5f21306fba87f7da9ec5cf51d4adfcefa93c345f6b9ebcd556a17dd1","access_level":"public","create_time":"$D2018-10-11T19:19:20.748Z","update_time":"$D2019-01-08T18:41:08.329Z","view_time":null,"first_public_time":"$D2018-10-14T15:29:37.069Z","has_unpublished_edits":null,"fork":false,"stars":4,"starred":false,"owner":{"type":"team","id":"7a4ecd430a93cd6f","login":"vinicius-mello","name":"Vinícius Mello","avatar_url":"https://avatars-next.observableusercontent.com/b5d4671aa488d8bcb984f575f985332403cde5f63942ea883a6a1a630941ec05","tier":"starter_2024"},"version":555,"publish_version":555,"role":null,"authored":false,"authors":[{"type":"individual","id":"ecc4854ab7ad95b6","login":"vinicius-mello","name":"Vinícius Mello","avatar_url":"https://avatars-next.observableusercontent.com/b5d4671aa488d8bcb984f575f985332403cde5f63942ea883a6a1a630941ec05","tier":null}],"members":[{"type":"individual","id":"ecc4854ab7ad95b6","login":"vinicius-mello","name":"Vinícius Mello","avatar_url":"https://avatars-next.observableusercontent.com/b5d4671aa488d8bcb984f575f985332403cde5f63942ea883a6a1a630941ec05","tier":null,"role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"lapsed":null,"subscribed":null,"forks":2,"collection":null,"type":"notebook","body":{"title":"Microclima Bahia","cells":[{"id":0,"value":"md`# Microclima Bahia`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":22,"value":"$33","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":516,"value":"{ \n // window.t = window.t +1;\n let canvas = context.canvas;\n context.clearRect(0,0,canvas.width, canvas.height);\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":19,"value":"height = width","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":49,"value":"particles = R.map(randomPoint, \n R.times( (t) => ({age: rand(0,settings.maxParticleAge)}),settings.particleCount))\n","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":17,"value":"$34","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":54,"value":"function evolve(t) {\n R.forEach((particle) => {\n if (particle.age > settings.maxParticleAge) {\n randomPoint(particle).age = 0;\n };\n let x = particle.x;\n let y = particle.y;\n let v = vector(x, y, t);\n let xt = x + v.x;\n let yt = y + v.y;\n particle.xt = xt;\n particle.yt = yt;\n particle.age = particle.age + 1;\n }, particles);\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":15,"value":"settings = ({\n particleCount: 800,\n maxParticleAge: 100,\n animate: true,\n frameRate: 40,\n velocityScale: 0.5,\n fadeFillStyle: \"rgba(255, 255, 255, 0.97)\"})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":402,"value":"function bilinearInterpolation(array, dim_x, dim_y, x, y) {\n if(x<0) x=0;\n if(y<0) y=0;\n let i0 = Math.floor(x*dim_x);\n let j0 = Math.floor(y*dim_y);\n i0 = R.clamp(0,dim_x-2,i0);\n j0 = R.clamp(0,dim_y-2,j0);\n let u = x*dim_x-i0; \n let v = y*dim_y-j0;\n if(u>1) u=1;\n if(v>1) v=1;\n let a00 = array[dim_y*j0+i0];\n let a10 = array[dim_y*j0+i0+1];\n let a01 = array[dim_y*(j0+1)+i0];\n let a11 = array[dim_y*(j0+1)+i0+1];\n return (1-v)*(a00*(1-u)+a10*u)+v*(a01*(1-u)+a11*u);\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":25,"value":"vector = {\n let size = west_east*south_north;\n return function(x, y, t) {\n t = R.clamp(0,Time-0.000001,t);\n let i = Math.floor(t);\n let s0 = size*i;\n let s1 = s0+size;\n let array0U = wind.U10.subarray(s0,s1);\n let array0V = wind.V10.subarray(s0,s1);\n let array1U = wind.U10.subarray(s0+size,s1+size);\n let array1V = wind.V10.subarray(s0+size,s1+size);\n x = (x-bounds.x)/bounds.width;\n y = (y-bounds.y)/bounds.height;\n let vx0 = bilinearInterpolation(array0U,west_east,south_north,x,y)*settings.velocityScale;\n let vy0 = bilinearInterpolation(array0V,west_east,south_north,x,y)*settings.velocityScale;\n let vx1 = bilinearInterpolation(array1U,west_east,south_north,x,y)*settings.velocityScale;\n let vy1 = bilinearInterpolation(array1V,west_east,south_north,x,y)*settings.velocityScale;\n t = t - i;\n return {x: vx0*(1-t)+vx1*t, y: vy0*(1-t)+vy1*t};\n }\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":117,"value":"function createProjection(lng0, lat0, lng1, lat1) {\n let projection = d3.geoMercator() \n //let projection = d3.geoAzimuthalEqualArea()\n .rotate([-((lng0 + lng1) / 2), 0])\n .center([0, (lat0 + lat1) / 2])\n .scale(1)\n .translate([0, 0]);\n let p0 = projection([lng0, lat0]);\n let p1 = projection([lng1, lat1]);\n let s = 1 / Math.max((p1[0] - p0[0]) / width, (p0[1] - p1[1]) / height) /* 0.95*/;\n let t = [width / 2, height / 2];\n return projection.scale(s).translate(t);\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":121,"value":"function createDisplayBounds(lng0, lat0, lng1, lat1, projection) {\n let upperLeft = projection([lng0, lat1]).map(Math.floor);\n let lowerRight = projection([lng1, lat0]).map(Math.ceil);\n return {\n x: upperLeft[0],\n y: upperLeft[1],\n width: lowerRight[0] - upperLeft[0] + 1,\n height: lowerRight[1] - upperLeft[1] + 1,\n }\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":125,"value":"function buildMeshes(topo) {\n let path = d3.geoPath().projection(projection);\n let outerBoundary = topojson.mesh(topo, topo.objects.main, (a, b) => a === b);\n let divisionBoundaries = topojson.mesh(topo, topo.objects.main, (a, b) => a!== b);\n return {\n path: path,\n outerBoundary: outerBoundary,\n divisionBoundaries: divisionBoundaries\n }\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":6,"value":"rand = (min, max) => min + (max - min)*Math.random()","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":13,"value":"randomPoint = (o) => {\n o.x = rand(bounds.x, bounds.x+bounds.width);\n o.y = rand(bounds.y, bounds.y+bounds.height);\n return o;\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":500,"value":"wind = {\n const size = west_east*south_north;\n let U10 = new Float32Array(Time*size);\n let V10 = new Float32Array(Time*size);\n for(let tt=0;tt {\n c.U10 = new Float32Array(c.U10);\n c.V10 = new Float32Array(c.V10);\n return c;\n });\n return c;\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":10,"value":"bounds = createDisplayBounds(bbox[0], bbox[1], bbox[2], bbox[3],projection)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":276,"value":"quadtree = d3.quadtree(projectedSamples)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":279,"value":"function quadtreeSearch(quadtree, cx, cy, radius) {\n const inCircle = (x, y, cx, cy, radius) => (x - cx) * (x - cx) + (y - cy) * (y- cy) < radius * radius;\n const rectCircleIntersect = (rx, ry, rwidth, rheight, cx, cy, radius) => {\n let dx = cx - Math.max(rx, Math.min(cx, rx + rwidth));\n let dy = cy - Math.max(ry, Math.min(cy, ry + rheight));\n return (dx * dx + dy * dy) < (radius * radius);\n }\n let internal = [];\n quadtree.visit(function(node, x1, y1, x2, y2) {\n if (!node.length) {\n do {\n var d = node.data;\n if(inCircle(d[0], d[1], cx, cy, radius)) internal.push(d);\n } while (node = node.next);\n }\n return !rectCircleIntersect(x1, y1, x2 - x1, y2 - y1, cx, cy, radius);\n });\n return internal;\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":285,"value":"function createInterpolator(quadtree, radius, array) {\n return function(x,y) {\n let internal = quadtreeSearch(quadtree, x, y, radius);\n let weight = R.map( (p) => [1.0/(0.000001+(x-p[0])*(x-p[0])+(y-p[1])*(y-p[1])), array[p[2]]], internal);\n let sum = R.reduce( (s,w) => [s[0]+w[0],s[1]+w[0]*w[1]], [0,0], weight);\n return sum[1]/sum[0];\n }\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":235,"value":"projectedSamples = R.addIndex(R.map)( (p,i) => [...projection(p),i] ,samples)\n","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":240,"value":"samples = R.zipWith( (long,lat)=> [long,lat], coords.XLONG, coords.XLAT) ","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":268,"value":"projection = createProjection(bbox[0], bbox[1], bbox[2], bbox[3])\n","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":109,"value":"bbox = [R.reduce(R.min,Infinity,coords.XLONG),\n R.reduce(R.min,Infinity,coords.XLAT),\n R.reduce(R.max,-Infinity,coords.XLONG),\n R.reduce(R.max,-Infinity,coords.XLAT)]","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":218,"value":"coords = {\n let c = d3.json(\"https://raw.githubusercontent.com/vinicius-mello/climabahia/master/coord.json\").then( (c)=> {\n c.XLAT = c.XLAT.slice(0,west_east*south_north);\n c.XLONG = c.XLONG.slice(0,west_east*south_north);\n return c;\n });\n return c;\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":390,"value":"south_north = 99","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":392,"value":"west_east = 99","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":400,"value":"Time = 25","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":66,"value":"topo = d3.json(\"https://raw.githubusercontent.com/vinicius-mello/climabahia/master/bahia-topo4.json\")","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":123,"value":"topojson = require(\"topojson\")","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":68,"value":"d3 = require(\"d3\")","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":27,"value":"R =require(\"ramda@0.25\")","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":435,"value":"import {slider} from \"@jashkenas/inputs\"","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":555,"value":"","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"}],"files":"$W35","stdlib":"1","titling":"auto"}} 32:"$31:body"