Remapping algorithms. - ECCO2

Hybridcoordinates&KernelsCommonformula+on• Levelmodels(z‐coordinate)– Eulerianalgorithm– FVapproachin3D– Con+nuousinterpreta+onofstructure• Layeredmodels– Lagrangianalgorithm– FVapproachinhorizontal• Ver+caliscontainedinformula+on– Piecewiseconstantinterpreta+onofvariablesCommonsoOware• Samenumericalschemesinmanymodels– Duplicityofeffort/soOware• Differentstaggeringofvariables– B‐grid,C‐grid,…• Differentchoiceofindexing– NE[u(i,j)istorightofT(i,j)]– SW[u(i,j)istoleOofT(i,j)]• SharedsoOware?• Programmer:NikiZadeh• Generalcoordinateapproach?• Post‐doc:LaurentWhite

Kernels• RepsforMITgcm,HyCOM,POP,ROMS,MOM,GOLD• Libraryofrou+nes– Proto‐typingwithhorizontaladvec+veflux– Extendwithver+caladvec+onandre‐gridding/re‐mappingnext• Finegrained– Ini+ally• Use“soO”conven+ons– EnablesNESWtransla+on• Dopevectors– Worksformul+plememorylayouts(fordifferentmodels)• BLAS‐orNAG‐like• Unittests– Validatessolu+ons– Testsfeatures• code.google.com– searchforhome‐kernel• Opensource• Openfordebate• Openforcontribu+on(willbe)• Butfornow,verymuchunderourcontrol!• Proofofconcept– Reallywanttoextendtoparameteriza+onsOpenissues:• i,j,kalterna+ves• Language(!)• Adjoint?(loopbounds)• Performance• Quiteafewothers