From fdcd2c714a3ab45b13f3b208cf22263c285cc420 Mon Sep 17 00:00:00 2001 From: Michael Denes Date: Tue, 5 Aug 2025 15:20:47 +0200 Subject: [PATCH 1/6] Draft blog post on PolyTEOS10_bsq kernel issue. --- .../posts/polyteos10-kernel/depth_vs_r0.png | Bin 0 -> 20633 bytes src/posts/polyteos10-kernel/index.md | 30 ++++++++++++++++++ 2 files changed, 30 insertions(+) create mode 100644 public/posts/polyteos10-kernel/depth_vs_r0.png create mode 100644 src/posts/polyteos10-kernel/index.md diff --git a/public/posts/polyteos10-kernel/depth_vs_r0.png b/public/posts/polyteos10-kernel/depth_vs_r0.png new file mode 100644 index 0000000000000000000000000000000000000000..6f3c5c9ac50821bac266e558ce4d778bba8130d2 GIT binary patch literal 20633 zcmbTegY9SwycGsA_5b5kAi?0YK4IGcnEkAffoY8lRQL(C*UXI<7au0|M$0u z(Rok)zxVl%2Y$-^k&J*KiXbcTR^1(OFB>@vNAhoIaInF0!uN!D5}Y`Syu5+}d*fe- zFG4pwESFc*g`RABWF7J0a+v&6mQZ?&^CS-oM>!mVMcYSf`Q{i}|McT5m!M$1V7}%0 z;;~*u1#|0}`h}o|U9o9jEV+neDoU6?N>iGs1O5=*sDvnzM5nl@r`$mR_;?mo2KKEj zgeDifKTTRd#Kb^s#E#2x^m-{P4Bl@swmqSufY{PgK-bQR61>6tF*Q$QaE>`8QP06W z9-1rgz6joZkFM`~wN@;jy3lBs>!!rJ;QnqPf!1wzdg-peV`QpIFZ2ar-TP50FLv(; zTxN?4#;lucw*|`Vvzwf4S)o~_Yu9!*W#Be%{`z>j_FXI)Gy<>Q3vct+z`($-U+09C zN(*ZzDwjgBB>bu+3WN_^uML`$5H9ANU}xKtccG_O+;J@4Qc=zDvdW|+9Fc@d-7iec z%sVBSwgMYjP73Vqs)@b@G|SXelsjm=0bPN}C$@LDCnM)k0xsK0DJdaiT~4!KKf291 zNbmjqvi0oE```X4A_vu|t+w+5%fff3<1pbXn%IPdg!aqj=Y5&imn*T@X&slD+Uo~i z2W@}8Ic?N=Xl1_nxZr)FxFTXcm?Y2kEok>L=#rM+#&@zR%&z0^%~ss!wRdi5`s*XP zerw5EW8oD#^>i>NCnxDBg130Za>?C^(jE8L72z5kzCJAtQtD(rhh511y}aShP_J|S z8smUu`G@wej0ZJbNU$+4BBOcl9ZdKG8& zr<#-#@`^ASl@(&(AUJtfiEO>4Hru5Qp zA^s>v{&9Uc-7?36ju114t*T7o(^VzQt(A%@nx`wy#;50OX;81PvWW?Qm*V*DV?2M1 zUuH7g^QW;-E#~@OZ~e9T{1m<5_@a5Ix=_##3iZCMVdVS#7Tn(7Mmt!=qUSbK-+`fx zmMk8uyV9_YX#Bg+#k?i1^-P_4Yu+p9DTXjZS+mQ|F;Uz3UmNGhOCk@uOLRr!lkGya zUxEj%4|fNCrD;8p_{19B_+qdCoADAwB-hr(SmA319esx$+?Bpq`l;uyTm55YomTtf zzlXfccre;O>WyOfNq5UatA~?7 zy%c7ize48dsZCU=^AkkJ*^iAr=QEb%ezy*(YzE2->>QB8vjn1uCzRlE?R;*(jtCj( z!8cMM@N6-Q=zc0Ql)`FdozBo7h~8~^x;YxualKu^C|Jw|Q&HWdYaqs;h=PVyroJ}% zrPe+#(C9|uSZ?RV#TOv^?O*MH<+x0&_0B#V18iA{DfclQ}@L4>Uq zl>XkDjqG%f?1sD~W>LxU+UL*7Xq;-jZt$#ZLRxW4Kw!XOMnk{3X>;Av)is>n+Cyxe z_m1YC>(7|D{K<<58C7e&)3o+Bmhij3+F&J5!u%MDOy=HsNzhd}@@}O!nqcCf{l=F0 z7DS>|s)5r_p5DF&&P_%=#7*kN9%*tPvnv9~=howhqqne8vl@c2{A&GgoW?)Z8n7Y` zQjL6T3-r*{O^%j(WPctn_4&67KioPE_3h`8?Nz!5#n3%X{NiSwucK@Y&vqq3Peczc zewh7YqvaGCL-7_TA~~{+?1)C1?W92zU6X17oyeyFd-QD6^rXRMXDZdeHQh}H|GGcs zmpoQ@N*afEimX&mpsx|}MK#kRd-EUfx58`bfR-D+vXE{~RQp1YlU=c_1k zG`jE$-wmITqCV%1Y`3v6aicgEa{jLOMwc&k#}j=M*J5wEV`^%$a72iU#fd*&pKfUq zv; z-ni}S69sA}A&|=GQ6PJ%GSFqTMHBOu6nyK?I4{%~)h445FU;vKnAyXeUr!;0it#v3 zYG{pq$}8&@9JJ$c!er-9Z*}my-hwY!qqJF|N;md6A}{A1iLvLdYYQ`7Q+NwGKe@`@ zK*l4lq<`C1=IRlbYh-+&GNcBBPQptq$fQ6vA%#`FCiJ_eE-Ne3L?~dd4ct#ydq?6b z;oP?i+Y;2H{S*Sj38FU_FuGrwP3+lYVBE$UnQw8^>)N%C>X;8TzNPE(K3HhO8IQmX zqsBl#cmAF4oK^Pz{derU*%mju=g$qaYJTIcX-L!;=H>dF>zvG5!G657-6}(wgxuxr zeD1i9=3VF{l8g;mG2)0l&-E3cBMxO~{h9L0E^Q6{M*!g`<7xBv-2B~K*CD}? zvZ^IsJ!-=0(rFO8y<-ZIMgo21m*-(IdOMg)4&!Z{dwfHS`!(=b5fx=|?WSE5`@qqq zl@mOoyu+)oq%tDwA4d}?TZi0@QX{&t?BAEZ+Q_yX_U6hpcCY-Ch7?C_{oN-F6J$V- zdiDFh0ADt8fF|ap???_*Q>u3%bRZThJp`9^3YS?^8*y$WE+<0Y;P;u=*`(IhKmjy3 zeeRmo91T}kI}hUL5RF#jOVBnX8i{&?bMB}x5d&Xmu4BgBy@+9dJ5MD-H_fW|)z}z~ zhN@%s^eq8K!^Q5cG|ltj&Qmi&^Aqy14YljhM`+#p#G(HWMRahwW|ZzF@f&BM!2`S} z*mrIfJAF-va&WJUec|tk)O3E`FNnDN?`O%)Wh6uTsKAz8d0(ReT{AJEKLMHHp zD?E?KPcZtt^gXo22`}t4$*n4Q2reFq-jL~z>+ZEXUXtU!j5;h!ph+mS{KZ0*A+YZ~ zTavDC8*UO&6~eQGns3I=d28l&-YN*shRA#`O+wH|%^Ow-zIR6`ybdRaJzEL04I$!q zGrmT1QIK?1S^vw_;(Jx&f`HwEkg%3-v6Zm*#k`v{)5GnCh!;HS1J2cf--BSUuhh}+ zU$;3AAbqcI%nJ$$zbIAHzGhJ%>MR%IjXuZ3+XQf!v~@3BI1L?$IfvABA(d^4jd2t6 z*;6nO1WDQ)ZrSCeVC8dosbH z6!a?alz@&Mnm3fDZ5I8ipgt`_0v?W^v`?Q)qM9xdhTHF+%i}KX?lkJM+$Q}ePjoae zZ9qkDJsxuyGx)9Nhgom=+kfMTKYaV55V3R68Mm@zE><51Qq2e^6h;g9lCno~P|4T% z8QJ}|&9p}zpn+!ep#x3kCt!LjaV-a~uHBUivud`V z;kBw-LkvrxOvcLeKUFvV(;@d7QY$KY*REYkJHP-5R}|eV`%&$eY=~qU01jZ8P?az< z8)oLzo8d?g8z79~R4fkoa0}3jxIqwvsa64zaR^2BsoJe)byIcK` z*iHIJxtFKa5fB;j-1NP*k9%okmsM5I(~ak&s0p~#WyxhxzVCVqBxoYkQLF_&$L|hB zP)8<~XJ4B{qL6kyFe%|PY1)Ob=W-oFbZ2@hKC=WN}Ww|!XC$hJN@xo4dKiQa_oh}p7@s*@U z^fhsY0><8E9wOG~JZDBgw+Kk2cs^z5cjsbJT<<@%F~}oU@O1hjM)Dhj zoN`_j$}k~HNmGig!~N-*K3qt?~#u>#5A^~ zZAatATom{H=KlKQQgox{Seo`ZExyIamxAqjs6^22L1FLPgRwpcenK4;Q7BYUv(|+A z9ZxLovkF2*+?1ZBcX$~!QzM=ct`q{f&7hC?U1(+`K=ig3o9F|dBv(YS_vLCM%Ya-L z_ZL3S@5d$sm@|~mgjtxRIPxzm;id9;7d0rOmjpE|vbB)3J?N-yZ>R>}`5DyxnVu-~ z8lP(^w1*i#7hg|e%rSqYXC0}*^Sz=raZzx2=Tr_NO~0e&eTG!xgTXc0@Li@|dp@!G zVtmLs&d%kigy!6B-B$zlln*)MMAamnfxV zX8+=_8`W+2IwO8JDr**~T=(_UXhXwX+FluEVH69@PTK$JX$(QCa9v@a`Ij6qw~&JX zo?gG3*~Lw>Ba@oG$O|@kNtRZhaau09_g^zO&Qrx=7i$}{E_xlqtQRBMIx0zNO z5RweRqIQ6ZShVBJH(I@{WK36;ZWE?Aj)m)g|iJ=+A}|jiG+dc5%S&(6~mW5 zD?w=s^t%28<38J|WTeEa$cuNjy<$N&2qF` zshTu_(YqUDS!~zAL*;NAaN%C;GYc!Id7dmAkb+OMXWcKS6LX7Fx|m9yK7B3?t(6Zi zF1l2@jk>sEVR{G4QZR=WybeK~X+Arsaosv0eDQMAGE}7C_ZGF_%_b!^!Zf_qC1B)q zFQ?)Cv~uuNaWeHCXA<8>=(@s;`1*#$2_eh%E~VvDRZ8kQRG<{d(Tuxl zr^lkvT{*(9<)S-qe#=E4kIjp}P$yHbbsc5qh2?i?7B_ekUcSS6)~LzC?@Nu*_YMiw zklQs%!syWHXHM=eztyvx<9DVx$1XaEo2L#nxV<0F+r%|E{;fe5EJc(C2|})HisKg9 zeAW|72u-AOL(`id(Np2to2hKE(nOE&%5pK6711^4v8UT>MOhGa<2LD_mZ!q`j0 zCQ=y7@<1~q_fov`MIc^*zox_TWd0qNuiA`bxM)8ryGcY_&J`j9~f zF=4=hgDnY5ziL-D7O5}Q3PHj^#I(oYv;T*{ih7bQK<1Fi;5q1gp{zQl zw(vT@TcVQ|>6mNiaucmt!Ze%`0!|b1NM0fR|0J(Y(I4rtY;cNM8=2#9e-s(m7AY74 z%5y}gQVE%alpHivDOxHB)jyh=4?6)Vu|6(m74_M#K6Rw#pyi~Sf^Q><8qeOpldgB< ztaEyTq9+=+xmlLZW$Bnj==HN^*M9j&cuLfqK(%GA&y=S%|ANN`sdUY|z8D?Y+T*Q0 z3D`pxGBS^zP?f}|6Y$(VT>4HzHh|d@1Nl&IKq3~&*?c_6OlmUBI!iX)ZM*V9Cbw&q z5kiH*3K2Z_I#t&V#_-M@o%vI&wXhcb+2UtL3O)RnHB_1r-Ik1 z-^|^Jhm6oMF#=F9@)o7pIOcMg;xLi*U^znRTy-bUJpj#u{>@^iL)hnO(>`*(kT)Ny0Spq z>wgK@0P<^r68RNb>%(e~QqBV6Dm56fFl6?eoc4%d{{lDc*~$dyc+zqw%bXU>vEIsN zJMKr!b!O7a*F|3IiBiz5a#0YV2(P!Y={7L1va-f?qCJj- z5Yo%1(XJj|$4E&~`TQ+bFvAWd0w)!^ZMRiZshtc#DD1U zy|eyP9c@MO0NBN`6PvLj)PSvhOrIZNm@$XVxQLixQ~`T8(Cqd&?h+Y#0Z(hd*9`Ns zV<4lGgcbTd+^6+AcmvLiJHEQ(-t%PQy2W8dY6qd)=jMFgp^p@>UY0DM4O(1>iS{PS z5^+@@pHe9{yVPk$KD&}iiYY^IDLx|v^7~knDJ6K*Z`^V#a z+y4dMeaG+3VQDSa2!E6xyomMh>D2cax5_Sk_0$D0!_WivJQ7EX;RWJG-@pM#V*ez( zm~awD2b91Bru{d8_DD)m7s358Y*1VhpI|&hOolj73)lT>@%*e*3aAE!)6>Veh!MSh zJEDX-Alt*7s)){JA)H?fKSrds14Hh{;;Su`55_q+V){VT5lIg{6C`QI1m_V5sJamn zq@5^NZr+ql$0>q~Nvdvv$jPFpamj@*YPA%cCK;kWMld+TiFyylF5icU((@#}3SEFz z!&rQy2WnnS7B`%i-$J%FZ8=j%auk_@W`r+RT5DO8;%IR;SPp`!iBH!^B~3{L^ehXFOnyXMsr z7@Ll`oIIe(PjAM=cnL|10O$P~6!2oQoegm7&|hIO_G4=Km6XM_5GrFl$_i*P?^|$m z8*B@Fl~eZZPIly=4%Q1rAqef39_)@%8;Zu zFp8~H4c62IRVv0HZg6cBHQ8{{0&J8%V{(OJIM!W{>Xjr#5K=eCbg2878f7IlEiEKm ziGn21hdnt*?i62s4L<+_YdXX2g=1Bm(y>JA|F8!v{YJoDlSp$oEK58e2mI@&SQnT5 zxf#z@d+*axY0z{l$jUAq1fgP+z+9L?P~bxJ&dDejG4_Sb5=KKu<(pmyz7>wxY&xyZ1ozFG<;_YO&dgK4Be{jjjZZ&hMj z%$fhQLE?SlbM`NW9;qXp-0{{^m0KtAYy94qhkysQKm6BqdFa*~MR*xSDmaBk5emM& zib^sxC@Zqb|LmJmaa=a{J2}kk^Ca?5<8at#5R-sYMSejLmy@&4MYC6y1NG94WozCn zn5nnaYfPc0re5*itGj(fRTN~zcM`8Nv(#%9`d$0E9Gm4~qr;vDZ5LqBnINR^0ll(4 zT8P}-OS}o@%4VX}P01)*&YYbi^v&Jo>CG|o z6HWS)Qb4$o3OJVSolsua-%klyqOezm z1via{mhO+Kxw=|eS#6e0R%p#aqlurSv>{4f{9b6S>wdSrTxa{=(3t~%dIb7V@+wi* zw_N+UhPAjIEHqeQ{~g@7nlos2P8V<~UwdTo0RgDn=m?@nBDk}pwM+QghjIieuvZ4i zj4GY_UL1@0W|v4j*77X3mxP?FgUp{Lc&os8!J~SlY7x;?IV>KSorwuBFQb>X)!(4%o@=kD8 zN_zV-VchF;b&jatl=y#wQB9>SH|>C!0Ejp?92ix&XeGAsvHFR0M-fgWSxb#+`l@l{ z%ENP)&pu4kj!doK0pXy2nAu1VH)Jc82p*%%5t6JMQ&XwbF=vGwiT^W$e{L9&6_aFv zpQ>XbGBy7f^_NkbG;Gk~R0WV=kdTzV!o}<)r1`eWFec&>z}lH^=&LW3=cDsTs066Q z4sZ+3+V!1Isc@jR#Z}I?&wlx@u=`be` zo?~m%GDE_}D0nX}ds_!K&&+e>k_5E)B`Ky*>@AF>Ty0g4OuvZxe8`m)+hcqlOag`+ zOj%>n*vF?)VOo|TyEI0>B-hB^qd@|0i-4wPNdo>_zn3vi2G;Y>i!xraQxQo?AqQ=? z2x*WBx>kag%0h(^_&9CdJRWE?;Px&JWy~<`u>+z#&dUxJkBchD)NIS`2dE$%GXIS; z#fRw`GmI)ylnP_okJ|=-vz#vMyiB6hE?EW1G~S=1K#MhXxP)yf(de%IuvTM`3gZFk zO^}CqE+qpL_rWefjsdbnZSM#VVkd%Q7GU)(x|{|Nm9z<8x{|k4xqrlSX;cUD+I%T}+C{95m|-!(N&wUz$`iA(1-`ECcl!igLbKaG4rh zHwTnDs&J+3H@0)ypADfy{J3c@m-XB4K zZ?McF#Ej>w#U<>tzv!!Aq0&?npW!?^+o^70raVdX6uuFX(Wj}2BfL=FL0$QJi`?nEU7->HRP9a+aPi1G}_Md zE_$!jbF78wwYne0kP2RZfEBMA{I9A)!3{v#yC9?u9!m$JZ8DPvB;Vuq@<=igueiAQ zDVDa$V_ue_F!f{-*ki2x2gn+Z*7pMW@0ql!hqKBms+%KZQ6K4cC8b{g74w$02-k z4A7OsCW-KcCJsLsE0t&`0o+=HK|_F_ac%~;a1m!;Xf95e^>q@_n={c5hPD{moMHY$ zsJC|a>u*=^b0?byz5I6domKAq@@+{Dm_5^zKLN;VLam&)qN5G+sfqn00}F0>liZH- zgFYZp>a+%y1IayA_ih(wN=6CF^dD7QGMurLqlDt_|5d-H!U#ojbc)UWt+tfAiYmNS z#^)WXPJ@ZU8}Jc1#@uWpMS)4P#dOPlL-Kk|4b#7pMS@}hrE+#lCqVzRV>*^>{BO$m ze(HsCnKjISAo|0d@UC_lg3$n0X!4<57DAZlA<|Tl7ac8kH5nyTBkgLlvilty->jWx zg=Vx1YgeiY2bcoW%+Br8K%6|;qEmH@Jx08>M>)09)Wp;+@YzeT_?~0z=KY15-K}g@6?@7myJfe-?-Ps?KxBi5UiSbtl!&kgvnplUKGNOFxCQ%4$lSiC{k zIV18;$Nkx7EhC>ZAXmRv&`C{73QVw-q+mm>$$pOPtI3$DUzOEoL0ZfF9xCY;iyZ|D zJxmVF$?hmv4>3TY@rvfGQx2CIELklngcSJw*kWyd;nXSIjm+x=i-5dlAs|OIb zZ`;;0YVS(Ekx<;+~fJ46tOR%D}aLEyaM6Y182jpj5z+}~8;JoXSk5ClEM z6U9Dr+kE4zAwVoCqd=J5-05f=!lb;IP@rdkT^~POnn)K*Vvl2BVEi;=8MnMpV-o15 zFDCFQy(uORv<}3eC-|=Rh#$ceafx1Y)9`|^gXi7CJ-X?DRw*hju1(wk?T0Oo2AYb6 zh{~=am*=N*Juj9PKi{rrugDaDPyoz2b<(MNAaWQW#IE)6mbPUHJv&;A#OgodkkG|o zZFlFhg08!lU^vX0mHlJspb#+x{t}R!V%{k_=J#~lE9rX;h4BcSyCx=cQa(YI8o{(l za?u(6ue(U`w^FFhmRCJu*j*zI$|*2zFa^oGw{#9TvwQ)B0-SOZa|C9i8MeD`8#%Le z=JGYcjiGkaHI!)Yw!d1=oL(?>jklh& zMFFu@U=55v;t$n<=f42t_h;Tfpa&EYRcGGYK$0%+1u_BW_rM^uj=Aon6!00CC^jV6 zr!PkQ?i0n6zu6}-s;!3*+p>2}D?AGh8wj5TW+wcf{mTD}Qw^^N=yMe^-|V#c7cM&83D9@Msz}&kA+s0gWNsm>qq&)9&QW!%$ zOSt;2xEuVZ#}X~#SIXpV??!&vhsR-wBhAVB<8qPaX7om+7mKqoWvX5(rR9!YJ%|PV za$M>TevZdt{7xu?t8@euT2A(O(HUBX9#3(l9wLKTe2Ye@0f?kg*LQldD!5^;Y|MS z+_&*cITJ(NMAY5bRgi4LomH70)dmx8C`4TntVx-ma*+*$|E{hkKql&vCJ#6YyvUhuSevt$`<}{#qlZB}87N||< zo#AXm`2v5GtAg{B=;J!CL00ffa-@`xN0n*e1(?}$0D;K?EvLvk|Kr(43Hm;8{tz-|SjXme-!ifyTH{8v*pEkBn&K zn?35IH!jS7-AY~_{d^pu>GqG zvy@hv=o!vXEa@js>Hw?FUXkYJDath4$540^>_?%YA zMifrRwTvKglRz#DC zXTTP-Kghi_s7*!V((z9w)kpH5{IU$R(^F1V#15TT;8YEp84m9tm}ac%QKc^oF{j7mW58QCy_!cYd0 zCo1(R?+n5LZdi&O6qdP7!d)Bey>zu_WvA*wZMsa6C_z!3tuMV~D}3(cDCKKWP-S7b zh&crjqFGqRnL*}yLd+IL@a2XO4HnEW9|!@w&MEefjAEHuSp%^gb{6TpOjdI87tnuF zOUKk0eX@Ba>iR@~fC9oHXM8`Plq;XI|7`R?(0LOhlrPI`{~Ioo`upE`h#O>u@a;dJ zM{Y-eUy)I2N8A*oId^QnGB}f$NJIoI&od&B__yM)-~p(7Pk!OE>NYs84FZ#^(Ea5K ziO~1>eg(v}7Z;^TZ_@$88x$vMyYEm47ijt(rZoD0RWB>HlHeXUa*R^8fI}9!pmz)r zPXzP?v_Q)!7!$1J_8^6wGkhQvF#UnjGIF{#J|cJ^)bh${?Mr?sY1kZ?rhIuVIkr*~I11%HAl+KJd5VM0~HuXgno%a&%Pq0m`2-g%%JjJ|*HdzXSwx z_TeM!j}0~h5+WCMzT!L@Z6ajP4wou5uK=kWUMUg!irn`n&!3w1kz>^S@6`^VwDYUQ zUF6@$TH*3>#a>Hu&y9VQKENnWrdr;wCo#(WzD{@<)&lZZ{vKY~@bGZl?5J#F{nt&E zb>6wo=I(}u2H@vd9C_wTNJz*+RL3T7{8(zJm|t_* z7CbMegIUK_MLHVwXsix{$BOz79;+U62}=P8x~@hklWWo7XE)6zN6uAVKwluZ%@;DJ zE?Nz1UwlKO4>S!IHN{G6apiq_;BAeog)8zN+Ow}yf{;>9#A~bnKJxT@*6tnM!*}v~ z44FC-cX1oGxVT5zM419`UW1aojArP%FzDHYn)y%)Pcct9lW4Fi1qHVMJ+)^!P&#Pk zZs}(Kbp_6N`#f4qO*-SP01^F#^Q-F%4ShtqYGZH@z9`07*DmlrzR5C^s2xW2khZK+ zCKk;R=tb6iJ%g>McNSrnOfD3Y#dPW|Rw?{#;IsU%a;a?5Uo(U}96>|17qIrIPLKqA zv$?>imHasoVV>|jRV>Z)skOAU;zWhkNC(Dd45(3VUSU_k@2i|g6XdQU{3}p|giHo* zL_W0M?S{mjTM?Rc-TLJR{F49e!G&4<#+)7!59mm~H8yGi65o2ntI&cIi`$b?HOw#< z$QE);6;E&3vwcn_#5cdsac&e$g8nwUFmu5(-h6Df`NLK(*bRD8KEM0hgN_HoP}81k z8>TDltj;+XvK1UII8##i@mAjgF6Mi$!58*T6fUs;C}ZWmbR>;jWl4=a2(B`2yJ*Q z!TVzoY_7lV3T9-f-zud3%;JRLw?lh+-9IGo0(+%}MV%E6dP>c|baIp{fGJHWox>V} z@m1o1a^?w<6RxOgh|gb-bBzeq1A+@cicSI)ln_8d0JO#|EBRx{n11T4oXk-&YNTVE zcdR?yW4!jv+BN?C#8$k8VkW2Qk@hm7ic8NHq_?|r;MF36Gq~7;gWPssOx9bD9Ip*= zDF*y}0|>W|>v{F0lOP4;00%C6zIxq%-|Xylm*#XTPm#A!6-%pQFsNl-_Gu?Of>KzqyTiM3b_TeXAg(b-%6o#KQS zM1iWwZMzJpZvV7NRJ$b@+)yTVDe zW<$OOO;dcqAKo^?hu>;^D0Ku4s{I?z4MzJ?axWz0V5tC*W%b33uLrYteP*awzj#_R zriQJuN`V4Ki9D3@wP9;o2T4+m@13WoAIZrTz#lI-8##EYso0mEJ?J#8d6F)?Lxc@E z<+*MKfNC#M2Rf;QiLC6mZ!xr)D67)`gs2wZhObzySMhsa-oM>p4E_D$7NB+NccaI~ zVxRtSr5PH^Syz!F&)??ylvkd(fpq?7aQp}G%eRY%FQgmStd1HLB29 z02S_U72C_VcC1AH1DHxSzNpDQ_7^r|uRi`=nV}9K%y}yrIHEj|a0*QEG0gpKxavTjdi{g@CnB>OZ;evvYq7SxnAA0I}8i7y&3uu@O!XGS~!uQZOpC)9kWNqg6Z zK4+84nt5Yb0q*pjA!<%1uHvz{D>%5z%~uAyA9z0gPhZ=YKbG+-d zH@pRVKQsVp6rk3qUmrB*JeqR6?G; zEcHCwen~EzHH*ipJDkZIi%knQ$vl)Zgz{0cn=IkUIqJ;*1c&N=vF7V{a>1P*Bk858DW1cueG; z6mwSV_Ego5d%4~txr2>wTRdUOTc+98>g7gwUXAaz#fbj50J$$@=dt|&`W7aztS+A6 z5s9~~4J_a^?Fmm60*3cM#AE9p>hnkvRC)hc^JS{>tqNQ{41LSoy2E5vJ-lE^-rvU@ z=#^UBWh^77r_Ktx8s+M>-Yn8`3+D5Kln_*N0QV*NQ2x&EW`ku0Z#E2NJB+5M4r5q{LWE@_NeHu1k!F5O4#lN`6m!(I|`hka8S z7P@@Ywwkw@qW%G8O4`)Z{_gsWt6nOw$TptykzDeYkU1V;oHZMoCNFYs^8nHUKWy|X~uivIE9sV~?%aIxt3 z0O(=@Lhqv+QGBPrllXhl`f`TFi4z7iDw&Ae6XcN1&EJyE28*HrNggfc8%+?i;= z6ct%Lt>xY#nZ|S~+T9<{XBb?Zz89*o#brVN=HoQ|Y{n!_m@QGQSc>}Jg!oj^Uen`6#Q8-P~K zq3eL9hXkQL{26gdip%g5wB?LJVWujX_f6q2i>b zZ}NeB49TbJBfpIcIlOl)pVN0!9(EN_C{2aog=F#292v|>n;3Do)+$v? z?<0&Q4_kX8weTp95v+UpM!-2-eBhoOZOXk}hMt3n5{* zLckW1x;dt1D6>ilsg?Ebop|mMUFK`B{U;26D5(r>t8Chvmi7yZJ7!vg`%VJu%@7U! zyob`Qc8LOGF}&4UkL>Ycz*Zo1*Mo*6EvvzpiWNwENMCoFBqHER2dNeN_fyBFabF}w zGk?P07Ty~%954hhm5(_X2Qd{(M#O$MxpmU})tNOJa>Z3EP$2*JBnbj6J9Is0xZoe> zBT^g6aNiXm(^a3U8)xCjAvQqy9a)TWek`JPZWPw~v^|+bIFIoI2q-`eaB~<$ zi{j}6Z{FL{ia-p1PFrxQo?pY504UE+;`*s-SrPRw`p1Fzyl$ZqCh=mPOf1^Tr6>~q zogwp~Mz@|kgwW&A6x<;iX^Ou0=1Xi)0FE_47)p6<)DF-RY=oH^OlJ0R1dN)%o&C)h z(_qW>Fg8pX5(xjdzK1vW8&SJ+pdM&z_~@N6QuWt4H75fgxzYV`D%cGN8Xq1(!4}Jr z9MEliHFM$nZ91!D2*fZ4p~)&>n+Z@^th#f!%c&cX{a(!%y|{bH_5oE3_h9ICN5iVe z&0eFG#?M@l4zS-K1wP*ewtF!D^;yLWK%%U1yy@C#Ki-Y(BruvEG{RddCH-t?upiJ)oOLmI%oTVHRCaQo9cz&lpUeD1j%-l+itbR-WG;L<- zLOWms9R9)G#p1(iA}mb|wAZd^;%YU{TVs5HI&B5*2p|ep0FCYbs3d){PEJqfw-_b> zZAVv%{*oj`2+Ha@S)Xf+PPyaEYO3vfEzW^xX#U|5SRh3Fjx>@`uxqml!X>zrqG!95 zlCDb3Ys-EAf8B?&$$MrNoi?wFuLS@Vj->xBj%)hVItxl4x4AmqG@#`e?FTjMTT8$4 zAiZoH3kF^%1_&h$yTek-O)o1|*G9Filg7;_=@hXvmmdko55aVrI~=CJfesNA(HW&7 z_;B9QQ`8&MeNVUn%2$>VC zc6a^Wzg*9l+MR6RfF+8E{BN_cJ?7)u=&+mXey8&-n@HVvXDVh*^)%SiS{g_HG)687 ziJpDC#pNOG_POAbtUr~1!-z+e-9|jg{|}}oP-=>Ch^$YuUdq@5Ob#&|E0_3G6?GH~ zuu_B^iW^XaJ>C8=ZM0zIuM?5Q_dBf(W7(H{$v9mO_EcU0JkX>^1m4N1eeOm&$YNCf zS4&)gSN@Efbk9NG{nWIUNfmC@)R!8A21owZBjn$b5_ zYWqiO)5rg&>WLQoOw<&)UqzIhsw+3)Xd`{1s4-t0<+Z53`v;G%45TD-LMLl&(plMF z3zQLid-CXG1?9qTY;eKNFN7UR+D){XQWK2u8s4FQhCcJRG0MSuw0z3B zRcPa~69qqB#c(+B`ZijmViwsHGPz~XURL&Mm_;zj{gi>>2InkCbfQ9b+(_@;mr@)= zKF&FGIw85$n!@#`3wm=$K2B5{(sD2fAa0cvaH*q(%TwI_s~y3!yD>YO>d8DOMAI0e zHJYv6Kg-^pEK{d*_!iW>3e0qL<)_Qp2yoYj=99xeLvG$}LOu&2P>*M&Y6%hef)x%L zzXHECY5YW5q~;?zdkw`!FPGy&N3D5&7`t;`^Wkz9cWdW_&N@S^ZfCvtqxId=scuFs zVY@Cc=#uLdlbj&Ak`HP(SS`+dY1ycp|CzqULhuHE1Ix`;Js5QjnxUojzy_HF?F3@~ zAe99Q317+==b`&t-V?jUWyPCx1athmiE?ul*gTQe-c&;do%P>=pjlQZu5Gv6VjB`V z*^p_$_%>;%<{n?u;0j;U#|Mtkrs(qIz3Dv1R4EOCT-{U%8X z$WNzF2s7WhBh?`wyt;V&TL2eiGff*mIut!(fKRZsr`(|SbjyXQp-eBngG64VJt@ZP zCiZ$C8OfbLS7fS0&5LoBnutJc&?s`GC77AH$t0!8^oM8CVAazy3Ggwkw2+omeCJ_9AR7*MyaQot;`Gv9LIkN5<5_)Y*WFS`6O zm@L<50e$9x=XW7=Aj!kK*5ID6XviEcCG=7FyK?U8v*axv|Bd> z)HVzJqVLOY}eQ#6HQTM`UgexAA6qNJQ7GbH7WmFchPUZ>4D~<}_U= z6N$O*|HYq^ojGLvg(Q}GBy^T!<=^KG6~i=aTNq9(%rI$gB@3Rgbd{PI?=sd^)t4RE zamLNux4MPNT^hM1k1L;Asa;}L%sIdvUw-oX3(k5GQ5j=S$C`X`3R`|M^IOVjd*cVA zQqTYamG@+0#BVb{ltl{cuFnEn+Ffi&N@;FcMhutiEPJEBY3U5(ru4HGDuo|2%cdk;i z%xnS96T+G;N{aw@#=OQ@f@Ut40*faad}S$g30f*3!7$U)!=<8$@O0;ZKe^F1`U$yo z+V1X$<&xs!8Syc+N2zljF3PLh3;+0~GYs>2XLqfsNBXEhXS|zN_fWG-5T)1xHB9KF zKf;TX7P}MamUTG%ut$kerm^jf%1^A-WVi^AOkQ+SSm3shwXH$AFHs<{v?jGw9 zr^3x2!0|QW-)jjd`XHJYlpF=Gh`1SCcSI8ap?69DuiBTq&d_V`3ovo|E{D^nZOniv z0FK>O-@EHer^sKN!B#-Ry?Y~sN(p*Uj)H|S8J$ZhG{^mGevOe;w)Yo`?t^q$vy){MXES*#FLL12NQw{1)F7o z$gfl{MN_D}xDsN4Mq`F@m|X|Uax)C)9O2^SbpsHyHk2w-8ymnIo+!NOGBE=zW!-Is zhCet3bm}($pEk}tp6UIM<63i)ThU5sh1gtTh(Z`K=2E#tsvT@};hd6~T*jzVXfCBO zw-YfX9rx4_$E8zfkwsWp<$7FmJwi+8y>ouQ$M2ute;yuR->+?ty+7~w>-~JazFXW+ zKzB7C`+`f)WuY7xBbet2U-fJ;4m3t{Z_zRJS>3IetpVm|uH86cdnMhxU&Xh!RbR(! z|K!<~s3!%l-j7*_T>rY?eHDbWJ)2SR67L3Tmcze+l28L=R!x#$YwEZ(VG4aApQhXI zW2?;$R8X5;{P!*)sHgMlj9uolo9tCY+c{gyZ+Onz_+i+C1;V=#B+sh{f#E%z!!;$g zYVhh}m|!wt#()0z;)sdR{ZJ^&%2y(aT<;k#orS=J2HO3dphRlO-O5>&5msLs`V7^# z0qoCG3b-i{Gy#DwW~$PqQ1o7rnWPxlj3ZfrpO3=lOP7W8fp7@30R7|PSb}^VvbfAg z(h(9LoGiXA&JDF}$hA!MEqA*I#uHennon1CeFz&$(r_G7r?K9S?L+>PI6>N$G}@kR zvDj(dbKx`^3XR_b*mFKo{O!`v3v~UK%t&?y7jQiT|2Jw!r5X9s*E_Oo*usZ$QDC8m z;Ra6~bIxNj9R8S^ntB$uiMc;|+Hq_Q!c10G87it4I73DZOnx)`M}-bBH09DiBm1N-V6RKwNG_z2(1Hk2LwBuO8(=g7qH(s; ze712t#8iZklD3Bz6w*C(mpr|fFI0S7!~NK3WLOx!(F^DxVxsUoC0VH|U+Pnonyg(DeyO%7@nq&|?HqWNdiGszOOA80ctU3|E5@F^jYrq_M!QW&d?5u?tjcHA-?qH0%x`6@(D2`E9Lp5lG zKM7(x2U@npuAQo-Y6{GM@G;I_%TK^bQhw(-ymSg@V_cnAR#S3suE0n%&SSR#sD^gS zN|Bx-(iNQ(`;?`l#|`#%qV1ayulG)o5J;Ruw>7};w>TvLp#!U)C-b@hNvAZcq$!cM zXB(P*7cOM5^e>ZjAGqFgFh<4)X3<|JsErMZOKt034NH>iw+Aw2tj7es`3NXk?Ah*w7=CXAecjM9mXRnH0Yfw-I! z66e0jG1we`+#UDs499#f{{$y}*fXy~XLbCfH~~cq_dRpCt567Eq5};Yytth)Nc>3Q zrV2k^zL*cuB!01bM(a~)-DId2rJCXFI4;JdfSnz1ua7pHs76Bnl2xTK{b z4YxyT7bq?oIUcMEFR!*4_j@Tz;ApJ-BlB+E*>*+d2J42ZicA@xeCy-+&CMjUvTon{ zy4xRra}W$r4*RX{#3B-_O($w#bA~O0(oFt7$XQJ^0NL8tHUr+EO3_`F1{AQ3Q<N z)=6Vfof|_$UK7tWJLP{?=i{0iHjB28)yYs_CLa$ z@(VmuemD`My-qvYQ4Qn^ysYYt|F~t&mN1 z9qqpd8wLXFqmu?Yso+An+AaVO#|J>E0((lOPNK~tsL&8a3x||`NlU7K=5pazG1Pg& zb>06CAcX}51#nO^x47~hCqW!u1LeT^1dvj)etqOIXemu#kLj3>Yj=tL{`KQ`Z7;(t zW7C)f)nVg`!4%?d6+I{ayA>+E(Q6AzU55>J83vb&Kb5}CyQN}XOqfcBU*}>OfAu~s zj3uL4O)k5zd-nB1Xb@FvulqxS@2v$>W0%1m(#4VZQRt)O`f~85y#O}!d{zy;CV{MW zf}_=vj_g~RdE=-XVrGK@)0~u3_i-SNMU;d$K%mo{kGAY1xp)B*56l|WKH|{N&ly{N zUnCBQUFt^OG|p2n-$!l~d*}CK6;~oND~k~U!pB=UH4+rF_MKs<5(CClGc|UPerrTHQBdwBIoU0M=t`;C%03-8Y5n%{u5X76>oRz>nl;Lo9GmkPb)QD z7=(?=r(;f+!koZL()WVnR^%!w+i(6v-3JxwLw~-igtc{BKVHP~ox16)aa`Lfs4^jh z%Rh;H5>UU(&oC7wqd!A@&<&Vv`#h&G(Mp^2MEe{&vrxlSQ3?uk1tK^=507U)Zf=#6 zMZ8lox0>!DaAVWcj=$Dz5L!*H=hnJMcGR8kExi-+z9#tG)l3p#9iACX<}-@faWO8% z;uL6VNn0X{iLxyK(fU{Of(h;>yDud~Nxt>}!g0lVD?0bAiyTEN(`>1pH^PQnu>5F@jU+ngap;WK zUQKjzY6*#LIlW_=|0r%2W}QdGbW;R0JbT1U1}c}KZDVPl41#1b&UwCk$scWRsAZ9^ z8K$A7BR8MgGqLyepU-VGEG{OlJ1FVqT44RNEa^nd#)DUNAT!K5n@5$67Q~8TBC0p4 z6IOv(KxZ?~@+*8dV1uz5`nq}GSkSm5T85YVrIE?a9Z4mU?@h2e&%{lQh}XW2J4FwU zEv0)!+@R;qo%^8?3V#eHzAhr}1lX!L|1IIXa$hyXqA=XJaVmdk&R-$PgRNuIVZD9p7O1r>yJ8CgZ`?vk4XMXRV!xfA+cfa@;uEAIseYOgpZM%~Hr|1_|5 z-!=GjM6SB$((AD2#P|O8R$b;mezU@znAIqoa}Q9k+HQ8Voc!9%^(cZy8U;pRH{K2{ zacIM4p^KB$^}dG#_JIl|@_&~{Xw^lgftSdV2SsY@ojg~O<4hO&$Ju%J2`H#GkeKIA mjfE`_LAef||9@}39POih$MdM&%fCdifc7>=38hx#tA7I!`r;e_ literal 0 HcmV?d00001 diff --git a/src/posts/polyteos10-kernel/index.md b/src/posts/polyteos10-kernel/index.md new file mode 100644 index 00000000..a26074d4 --- /dev/null +++ b/src/posts/polyteos10-kernel/index.md @@ -0,0 +1,30 @@ +--- +title: 'PolyTEOS10_bsq Kernel Bug Identified' +date: '2025-08-05' +authors: + - name: Michael Denes + github: michaeldenes +summary: 'An issue with the `PolyTEOS10_bsq` kernel has been identified, where incorrect values of seawater density were being computed. In this post, we explain the issue and how it may impact your simulations.' +--- + +**Update to `PolyTEOS10_bsq` kernel** + +In recent days we’ve come across a bug in an application kernel of `parcels`. The `PolyTEOS10_bsq` kernel is used to calculate the density of seawater from the temperature and salinity fields. The kernel is based on equation (13) from +[Roquet et al. (2014)](https://doi.org/10.1016/j.ocemod.2015.04.002) where the density $\rho$ is determined from the summation of a vertical reference profile $r_0$ and a residual function or density anomaly $r$. That is, $\rho(S,T,z) = r_0(z) + r(S,T,z)$. + + +Up until now the kernel only computed the density anomaly, and was missing the vertical reference profile. At the ocean surface this isn’t a problem, as $r_0(z=0) = 0$ and so the computed density was correct. However, $r_0(z)>0$ below the ocean surface, and at relatively deep depths, this value can be significant. For example, at 100 m depth this value is 0.46443 kg/m$^3$ and at 500 m depth this value is 2.31175 kg/m$^3$, see the figure below. + +![Vertical reference profile as a function of depth](/posts/polyteos10-kernel/depth_vs_r0.png) + +By using the old kernel implementation, the seawater density at depths below the ocean surface were being underestimated, as the actual seawater density would have otherwise been denser. Simulations that computed vertical velocities (such as sinking/rise velocities) using density differences between a particle and its surrounding seawater may be impacted. + +Below we’ve compiled a list of several papers that may be impacted, but we urge the community to check their own projects for use of this kernel at depth. This list may be updated as we become aware of other potentially impacted papers. + +**Impacted papers** + +[Global mass of buoyant marine plastics dominated by large long-lived debris, *Nature Geoscience* (2023)](https://doi.org/10.1038/s41561-023-01216-0) +[Modeling carbon export mediated by biofouled microplastics in the Mediterranean Sea, *Limnology and Oceanography* (2023)](https://doi.org/10.1002/lno.12330) +[Modelling submerged biofouled microplastics and their vertical trajectories, *Biogeosciences* (2022)](https://doi.org/10.5194/bg-19-2211-2022) +[Influence of Particle Size and Fragmentation on Large-Scale Microplastic Transport in the Mediterranean Sea, *Environmental Science & Technology* (2022)](https://doi.org/10.1021/acs.est.2c03363) +[Global Modeled Sinking Characteristics of Biofouled Microplastic, *JGR Oceans* (2021)](https://doi.org/10.1029/2020JC017098) From 7f15a4a37d5f375b0e3bfe2da2d1137a27f07708 Mon Sep 17 00:00:00 2001 From: "pre-commit-ci[bot]" <66853113+pre-commit-ci[bot]@users.noreply.github.com> Date: Tue, 5 Aug 2025 13:43:55 +0000 Subject: [PATCH 2/6] [pre-commit.ci] auto fixes from pre-commit.com hooks for more information, see https://pre-commit.ci --- src/posts/polyteos10-kernel/index.md | 15 +++++++-------- 1 file changed, 7 insertions(+), 8 deletions(-) diff --git a/src/posts/polyteos10-kernel/index.md b/src/posts/polyteos10-kernel/index.md index a26074d4..8f22d072 100644 --- a/src/posts/polyteos10-kernel/index.md +++ b/src/posts/polyteos10-kernel/index.md @@ -9,11 +9,10 @@ summary: 'An issue with the `PolyTEOS10_bsq` kernel has been identified, where i **Update to `PolyTEOS10_bsq` kernel** -In recent days we’ve come across a bug in an application kernel of `parcels`. The `PolyTEOS10_bsq` kernel is used to calculate the density of seawater from the temperature and salinity fields. The kernel is based on equation (13) from +In recent days we’ve come across a bug in an application kernel of `parcels`. The `PolyTEOS10_bsq` kernel is used to calculate the density of seawater from the temperature and salinity fields. The kernel is based on equation (13) from [Roquet et al. (2014)](https://doi.org/10.1016/j.ocemod.2015.04.002) where the density $\rho$ is determined from the summation of a vertical reference profile $r_0$ and a residual function or density anomaly $r$. That is, $\rho(S,T,z) = r_0(z) + r(S,T,z)$. - -Up until now the kernel only computed the density anomaly, and was missing the vertical reference profile. At the ocean surface this isn’t a problem, as $r_0(z=0) = 0$ and so the computed density was correct. However, $r_0(z)>0$ below the ocean surface, and at relatively deep depths, this value can be significant. For example, at 100 m depth this value is 0.46443 kg/m$^3$ and at 500 m depth this value is 2.31175 kg/m$^3$, see the figure below. +Up until now the kernel only computed the density anomaly, and was missing the vertical reference profile. At the ocean surface this isn’t a problem, as $r_0(z=0) = 0$ and so the computed density was correct. However, $r_0(z)>0$ below the ocean surface, and at relatively deep depths, this value can be significant. For example, at 100 m depth this value is 0.46443 kg/m$^3$ and at 500 m depth this value is 2.31175 kg/m$^3$, see the figure below. ![Vertical reference profile as a function of depth](/posts/polyteos10-kernel/depth_vs_r0.png) @@ -23,8 +22,8 @@ Below we’ve compiled a list of several papers that may be impacted, but we urg **Impacted papers** -[Global mass of buoyant marine plastics dominated by large long-lived debris, *Nature Geoscience* (2023)](https://doi.org/10.1038/s41561-023-01216-0) -[Modeling carbon export mediated by biofouled microplastics in the Mediterranean Sea, *Limnology and Oceanography* (2023)](https://doi.org/10.1002/lno.12330) -[Modelling submerged biofouled microplastics and their vertical trajectories, *Biogeosciences* (2022)](https://doi.org/10.5194/bg-19-2211-2022) -[Influence of Particle Size and Fragmentation on Large-Scale Microplastic Transport in the Mediterranean Sea, *Environmental Science & Technology* (2022)](https://doi.org/10.1021/acs.est.2c03363) -[Global Modeled Sinking Characteristics of Biofouled Microplastic, *JGR Oceans* (2021)](https://doi.org/10.1029/2020JC017098) +[Global mass of buoyant marine plastics dominated by large long-lived debris, _Nature Geoscience_ (2023)](https://doi.org/10.1038/s41561-023-01216-0) +[Modeling carbon export mediated by biofouled microplastics in the Mediterranean Sea, _Limnology and Oceanography_ (2023)](https://doi.org/10.1002/lno.12330) +[Modelling submerged biofouled microplastics and their vertical trajectories, _Biogeosciences_ (2022)](https://doi.org/10.5194/bg-19-2211-2022) +[Influence of Particle Size and Fragmentation on Large-Scale Microplastic Transport in the Mediterranean Sea, _Environmental Science & Technology_ (2022)](https://doi.org/10.1021/acs.est.2c03363) +[Global Modeled Sinking Characteristics of Biofouled Microplastic, _JGR Oceans_ (2021)](https://doi.org/10.1029/2020JC017098) From 095abb71b321078f16692601f76c22716b1387dc Mon Sep 17 00:00:00 2001 From: Vecko <36369090+VeckoTheGecko@users.noreply.github.com> Date: Tue, 5 Aug 2025 15:54:05 +0200 Subject: [PATCH 3/6] Make bullet list --- src/posts/polyteos10-kernel/index.md | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/src/posts/polyteos10-kernel/index.md b/src/posts/polyteos10-kernel/index.md index 8f22d072..228b074c 100644 --- a/src/posts/polyteos10-kernel/index.md +++ b/src/posts/polyteos10-kernel/index.md @@ -22,8 +22,8 @@ Below we’ve compiled a list of several papers that may be impacted, but we urg **Impacted papers** -[Global mass of buoyant marine plastics dominated by large long-lived debris, _Nature Geoscience_ (2023)](https://doi.org/10.1038/s41561-023-01216-0) -[Modeling carbon export mediated by biofouled microplastics in the Mediterranean Sea, _Limnology and Oceanography_ (2023)](https://doi.org/10.1002/lno.12330) -[Modelling submerged biofouled microplastics and their vertical trajectories, _Biogeosciences_ (2022)](https://doi.org/10.5194/bg-19-2211-2022) -[Influence of Particle Size and Fragmentation on Large-Scale Microplastic Transport in the Mediterranean Sea, _Environmental Science & Technology_ (2022)](https://doi.org/10.1021/acs.est.2c03363) -[Global Modeled Sinking Characteristics of Biofouled Microplastic, _JGR Oceans_ (2021)](https://doi.org/10.1029/2020JC017098) +- [Global mass of buoyant marine plastics dominated by large long-lived debris, _Nature Geoscience_ (2023)](https://doi.org/10.1038/s41561-023-01216-0) +- [Modeling carbon export mediated by biofouled microplastics in the Mediterranean Sea, _Limnology and Oceanography_ (2023)](https://doi.org/10.1002/lno.12330) +- [Modelling submerged biofouled microplastics and their vertical trajectories, _Biogeosciences_ (2022)](https://doi.org/10.5194/bg-19-2211-2022) +- [Influence of Particle Size and Fragmentation on Large-Scale Microplastic Transport in the Mediterranean Sea, _Environmental Science & Technology_ (2022)](https://doi.org/10.1021/acs.est.2c03363) +- [Global Modeled Sinking Characteristics of Biofouled Microplastic, _JGR Oceans_ (2021)](https://doi.org/10.1029/2020JC017098) From 4be917b6d2ab3231cf86baf71b07353775329190 Mon Sep 17 00:00:00 2001 From: Michael Denes Date: Tue, 5 Aug 2025 16:25:10 +0200 Subject: [PATCH 4/6] Update src/posts/polyteos10-kernel/index.md Co-authored-by: Erik van Sebille --- src/posts/polyteos10-kernel/index.md | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/posts/polyteos10-kernel/index.md b/src/posts/polyteos10-kernel/index.md index 228b074c..c47418d7 100644 --- a/src/posts/polyteos10-kernel/index.md +++ b/src/posts/polyteos10-kernel/index.md @@ -16,7 +16,7 @@ Up until now the kernel only computed the density anomaly, and was missing the v ![Vertical reference profile as a function of depth](/posts/polyteos10-kernel/depth_vs_r0.png) -By using the old kernel implementation, the seawater density at depths below the ocean surface were being underestimated, as the actual seawater density would have otherwise been denser. Simulations that computed vertical velocities (such as sinking/rise velocities) using density differences between a particle and its surrounding seawater may be impacted. +By using the kernel implementation before the fix in PR [#2133](https://github.com/OceanParcels/Parcels/pull/2133), the seawater density at depths below the ocean surface were being underestimated, as the actual seawater density would have otherwise been higher. Simulations that computed vertical velocities (such as sinking/rise velocities) using density differences between a particle and its surrounding seawater may have been impacted. Below we’ve compiled a list of several papers that may be impacted, but we urge the community to check their own projects for use of this kernel at depth. This list may be updated as we become aware of other potentially impacted papers. From ee1306c89c2c18dcd2fd36f6cd9cdfc90c221319 Mon Sep 17 00:00:00 2001 From: Michael Denes Date: Tue, 5 Aug 2025 16:25:27 +0200 Subject: [PATCH 5/6] Update src/posts/polyteos10-kernel/index.md Co-authored-by: Erik van Sebille --- src/posts/polyteos10-kernel/index.md | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/posts/polyteos10-kernel/index.md b/src/posts/polyteos10-kernel/index.md index c47418d7..016d9d08 100644 --- a/src/posts/polyteos10-kernel/index.md +++ b/src/posts/polyteos10-kernel/index.md @@ -12,7 +12,7 @@ summary: 'An issue with the `PolyTEOS10_bsq` kernel has been identified, where i In recent days we’ve come across a bug in an application kernel of `parcels`. The `PolyTEOS10_bsq` kernel is used to calculate the density of seawater from the temperature and salinity fields. The kernel is based on equation (13) from [Roquet et al. (2014)](https://doi.org/10.1016/j.ocemod.2015.04.002) where the density $\rho$ is determined from the summation of a vertical reference profile $r_0$ and a residual function or density anomaly $r$. That is, $\rho(S,T,z) = r_0(z) + r(S,T,z)$. -Up until now the kernel only computed the density anomaly, and was missing the vertical reference profile. At the ocean surface this isn’t a problem, as $r_0(z=0) = 0$ and so the computed density was correct. However, $r_0(z)>0$ below the ocean surface, and at relatively deep depths, this value can be significant. For example, at 100 m depth this value is 0.46443 kg/m$^3$ and at 500 m depth this value is 2.31175 kg/m$^3$, see the figure below. +Up until the fix in [PR #2133](https://github.com/OceanParcels/Parcels/pull/2133), the kernel only computed the density anomaly $r$, and did not include the vertical reference profile. At the ocean surface this isn’t a problem, as $r_0(z=0) = 0$ and so the computed density was correct. However, $r_0(z)>0$ below the ocean surface, and at relatively deep depths, this value can be significant. For example, at 100 m depth this value is 0.46443 kg/m$^3$ and at 500 m depth this value is 2.31175 kg/m$^3$, see the figure below. ![Vertical reference profile as a function of depth](/posts/polyteos10-kernel/depth_vs_r0.png) From 1cf97dba49e7ff1c2c66f208bb50315a2be8505a Mon Sep 17 00:00:00 2001 From: Michael Denes Date: Tue, 5 Aug 2025 16:27:13 +0200 Subject: [PATCH 6/6] Slight rewording of papers list --- src/posts/polyteos10-kernel/index.md | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/posts/polyteos10-kernel/index.md b/src/posts/polyteos10-kernel/index.md index 016d9d08..f83e91c8 100644 --- a/src/posts/polyteos10-kernel/index.md +++ b/src/posts/polyteos10-kernel/index.md @@ -20,7 +20,7 @@ By using the kernel implementation before the fix in PR [#2133](https://github.c Below we’ve compiled a list of several papers that may be impacted, but we urge the community to check their own projects for use of this kernel at depth. This list may be updated as we become aware of other potentially impacted papers. -**Impacted papers** +**Potentially impacted papers** - [Global mass of buoyant marine plastics dominated by large long-lived debris, _Nature Geoscience_ (2023)](https://doi.org/10.1038/s41561-023-01216-0) - [Modeling carbon export mediated by biofouled microplastics in the Mediterranean Sea, _Limnology and Oceanography_ (2023)](https://doi.org/10.1002/lno.12330)