# Research articles using Fortran **URL:** \ **Category:** Uncategorized\ **Created:** [April 24, 2021, 2:33pm UTC](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101 "2021-04-24T14:33:03Z")\ **Posts on this page:** 20\ **Page:** 3
**Author:** ![vmagnin](https://yyz2.discourse-cdn.com/free1/user_avatar/fortran-lang.discourse.group/vmagnin/32/28_2.png) [@vmagnin](https://fortran-lang.discourse.group/u/vmagnin)\ **Post date:** [October 10, 2021, 8:01pm UTC](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101/41 "2021-10-10T20:01:03Z")
Merci, I did not know that wonderful collection of programming languages resources in many languages: > **[GitHub - EbookFoundation/free-programming-books: Freely available programming...](https://github.com/EbookFoundation/free-programming-books)** > > :books: Freely available programming books. Contribute to EbookFoundation/free-programming-books development by creating an account on GitHub. ---
**Author:** ![nbehrnd](https://avatars.discourse-cdn.com/v4/letter/n/e0b2c6/32.png) [@nbehrnd](https://fortran-lang.discourse.group/u/nbehrnd)\ **Post date:** [October 10, 2021, 8:12pm UTC](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101/42 "2021-10-10T20:12:00Z")
Then have a look on the sibling project (equally accommodating languages _in addition_ to English) [free science books](https://github.com/EbookFoundation/free-science-books/blob/master/free-science-books.md), too. Since Matlab was recently mentioned multiple times in comparison to Fortran, but equally attractive on its own, I speculate Tom O’Haver’s project [_A Pragmatic Introduction to Signal Processing_](https://terpconnect.umd.edu/~toh/spectrum/TOC.html) (filed in the former in section QA/mathematics of the former) may be one of interest for the audience here. Plenty of code ready-to-go (for spreadsheet/Python/Matlab; not [yet] Fortran). ---
**Author:** ![FortranFan](https://avatars.discourse-cdn.com/v4/letter/f/96bed5/32.png) [@FortranFan](https://fortran-lang.discourse.group/u/FortranFan)\ **Post date:** [October 11, 2021, 12:19am UTC](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101/43 "2021-10-11T00:19:20Z")
> [@vmagnin](#): > > Do you know why they use the `.t` extension instead of `.f90` ? … The first time I have seen such an extension being used, to each their own I suppose. Additionally it’s the first I notice a particular nonstandard aspect being employed with a deprecated feature from the standard i.e., the `COMMON` block: > **[src/amrvacdef.t · master · BHAC-release / bhac · GitLab](https://gitlab.itp.uni-frankfurt.de/BHAC-release/bhac/-/blob/master/src/amrvacdef.t#L75)** > > clean, test, debug and document the release version ```Fortran subroutine sub() common, integer :: n end subroutine ``` > C:\Temp\>gfortran -c sub.f90 > sub.f90:2:4: > > ``` > 2 | common, integer :: n > | 1 > > ``` > > Error: Unclassifiable statement at (1) Perhaps further preprocessor “magic” might be at work here, I didn’t bother to investigate. But I personally wouldn’t touch that code with the proverbial ten-foot pole even until there is further simplification of the code using standard facilities. ---
**Author:** ![Beliavsky](https://avatars.discourse-cdn.com/v4/letter/b/ba8739/32.png) [@Beliavsky](https://fortran-lang.discourse.group/u/Beliavsky)\ **Post date:** [October 11, 2021, 12:53am UTC](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101/44 "2021-10-11T00:53:49Z")
> [@FortranFan](#): > > But I personally wouldn’t touch that code with the proverbial ten-foot pole even until there is further simplification of the code using standard facilities. A few scientific journals mandate that papers with computational results be reproducible, with the code made available. Journals should also require that accompanying code be standard-conforming, at least for languages such as C, C++, and Fortran with ISO standards. Considering that a research project may span many years, and that compilers implement standards with a lag, flexibility should be allowed regarding which versions of the standard are acceptable. One could mandate that a Fortran-based paper published today be standard F2003, F2008, or F2018. Journal editors and referees are interested in the scientific content of a paper, not the details of how an underlying program was written. They don’t want to bicker with authors about what features of a programming language can be used. Mandating conformance to a standard outsources the decision of what is acceptable code to a committee that has thought about such things. Once authors are told that their papers would be published if it weren’t for features such as real\*8 declarations sprinkled in their code, they will quickly fix their codes. The results of many papers depend on programming languages without an ISO standard, such as Python or Matlab. I wonder what the guidelines should be for them. Book publishers should also have rules about language standards, since many authors and book editors [are unfamiliar with them](https://fortran-lang.discourse.group/t/numerical-recipes-in-quantum-information-theory-and-quantum-computing-an-adventure-in-fortran-90/1839/2). ---
**Author:** ![vmagnin](https://yyz2.discourse-cdn.com/free1/user_avatar/fortran-lang.discourse.group/vmagnin/32/28_2.png) [@vmagnin](https://fortran-lang.discourse.group/u/vmagnin)\ **Post date:** [October 24, 2021, 4:39pm UTC](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101/45 "2021-10-24T16:39:05Z")
[http://www-personal.umich.edu/~jcaps/code.html](http://www-personal.umich.edu/~jcaps/code.html) > jCODE is a high-performance Fortran-based multiphase/multi-physics flow solver developed and maintained by the Capecelatro Research Group at the University of Michigan. The code is capable of solving the multi-component compressible Navier-Stokes equations on structured curvilinear meshes using a class of high-order energy-stable finite difference operators. It features a range of models and methods including Lagrangian particle tracking, combustion mechanisms, immersed boundaries for complex geometries, shock capturing, and discrete adjoint-based sensitivity. You will find papers here: [http://www-personal.umich.edu/~jcaps/publications.html](http://www-personal.umich.edu/~jcaps/publications.html) ---
**Author:** ![R\_cubed](https://yyz2.discourse-cdn.com/free1/user_avatar/fortran-lang.discourse.group/r_cubed/32/263_2.png) [@R\_cubed](https://fortran-lang.discourse.group/u/R_cubed)\ **Post date:** [October 24, 2021, 11:23pm UTC](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101/46 "2021-10-24T23:23:26Z")
I’d like to add this from an applied stats journal; the editor (Shlomo Sawilowski) is also a Fortran user. > **[JMASM 53: MicceriRD](https://digitalcommons.wayne.edu/jmasm/vol18/iss2/6/)** > > Fortran 77 and 90 modules (REALPOPS.lib) exist for invoking the 8 distributions estimated by Micceri (1989). These respective modules were created by Sawilowsky et al. (1990) and Sawilowsky and Fahoome (2003). The MicceriRD (Micceri’s Real... ---
**Author:** ![Ashok](https://avatars.discourse-cdn.com/v4/letter/a/ed655f/32.png) [@Ashok](https://fortran-lang.discourse.group/u/Ashok)\ **Post date:** [October 26, 2021, 4:08am UTC](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101/47 "2021-10-26T04:08:35Z")
A parallel finite element solver for structural applications from Japan: [https://manual.frontistr.com/en/faq/index.html](https://manual.frontistr.com/en/faq/index.html) Relevant paper can be found here > **[Development of a general‐purpose parallel finite element method for analyzing...](https://onlinelibrary.wiley.com/doi/full/10.1002/eqe.3551)** > > A parallel finite element method (FEM) based on high-fidelity models for solving diverse earthquake engineering problems is presented. Its key feature is a parallel solver that is tuned to solve larg... It is written in Modern Fortran. ---
**Author:** ![vmagnin](https://yyz2.discourse-cdn.com/free1/user_avatar/fortran-lang.discourse.group/vmagnin/32/28_2.png) [@vmagnin](https://fortran-lang.discourse.group/u/vmagnin)\ **Post date:** [October 29, 2021, 6:31pm UTC](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101/48 "2021-10-29T18:31:38Z")
Chen, L.; Ouyang, F. A FORTRAN Program to Model Magnetic Gradient Tensor at High Susceptibility Using Contraction Integral Equation Method. _Minerals,_ 2021, 11, 1129. > **[A FORTRAN Program to Model Magnetic Gradient Tensor at High Susceptibility...](https://www.mdpi.com/2075-163X/11/10/1129)** > > The magnetic gradient tensor provides a powerful tool for detecting magnetic bodies because of its ability to emphasize detailed features of the magnetic anomalies. To interpret field measurements obtained by magnetic gradiometry, the forward... ---
**Author:** ![Beliavsky](https://avatars.discourse-cdn.com/v4/letter/b/ba8739/32.png) [@Beliavsky](https://fortran-lang.discourse.group/u/Beliavsky)\ **Post date:** [October 29, 2021, 7:43pm UTC](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101/49 "2021-10-29T19:43:51Z")
> [@vmagnin](#): > > Chen, L.; Ouyang, F. > A FORTRAN Program to Model Magnetic Gradient Tensor at High Susceptibility Using Contraction > Integral Equation Method. _Minerals,_ 2021, 11, 1129. > [Minerals | Free Full-Text | A FORTRAN Program to Model Magnetic Gradient Tensor at High Susceptibility Using Contraction Integral Equation Method](https://doi.org/10.3390/min11101129) [GitHub - Yonfou/3D-MGTM-HS: 3D Magnetic Gradient Tensor Modeling at High Susceptibility](https://github.com/Yonfou/3D-MGTM-HS) is actually a Fortran program (lower case) using free source form and modules I earlier wrote that journals such require that Fortran (and C and C++) programs be standard-conforming. Before that, journal editors should learn how the language is currently spelled. Looking at code such as that below from [model.f90](https://github.com/Yonfou/3D-MGTM-HS/blob/main/Model.f90), I think the authors should use `parameter`s in their code instead of repeating magic numbers, and otherwise simplify it, but people don’t get tenure for clean code. ``` IF(X(IX)>=550 .AND. X(IX)<=580 .AND. Y(IY)>=350 .AND. Y(IY)<=650 .AND. Z(IZ)>= 50 .AND. Z(IZ)<=150) THEN Mx(IX,IY,IZ)=Mx0 END IF !-------------------------------------------- IF(X(IX)>=250 .AND. X(IX)<=550 .AND. Y(IY)>=640 .AND. Y(IY)<=650 .AND. Z(IZ)>= 50 .AND. Z(IZ)<=70) THEN Mx(IX,IY,IZ)=Mx0 END IF !-------------------------------------------- IF(X(IX)>=250 .AND. X(IX)<=550 .AND. Y(IY)>=630 .AND. Y(IY)<=640 .AND. Z(IZ)>= 60 .AND. Z(IZ)<=80) THEN Mx(IX,IY,IZ)=Mx0 END IF !-------------------------------------------- IF(X(IX)>=250 .AND. X(IX)<=550 .AND. Y(IY)>=620 .AND. Y(IY)<=630 .AND. Z(IZ)>= 70 .AND. Z(IZ)<=90) THEN Mx(IX,IY,IZ)=Mx0 END IF !-------------------------------------------- IF(X(IX)>=250 .AND. X(IX)<=550 .AND. Y(IY)>=610 .AND. Y(IY)<=620 .AND. Z(IZ)>= 80 .AND. Z(IZ)<=100) THEN Mx(IX,IY,IZ)=Mx0 END IF !-------------------------------------------- IF(X(IX)>=250 .AND. X(IX)<=550 .AND. Y(IY)>=600 .AND. Y(IY)<=610 .AND. Z(IZ)>= 90 .AND. Z(IZ)<=110) THEN Mx(IX,IY,IZ)=Mx0 END IF !-------------------------------------------- IF(X(IX)>=250 .AND. X(IX)<=550 .AND. Y(IY)>=590 .AND. Y(IY)<=600 .AND. Z(IZ)>=100 .AND. Z(IZ)<=120) THEN Mx(IX,IY,IZ)=Mx0 END IF !-------------------------------------------- IF(X(IX)>=250 .AND. X(IX)<=550 .AND. Y(IY)>=580 .AND. Y(IY)<=590 .AND. Z(IZ)>=110 .AND. Z(IZ)<=130) THEN Mx(IX,IY,IZ)=Mx0 END IF !-------------------------------------------- IF(X(IX)>=250 .AND. X(IX)<=550 .AND. Y(IY)>=570 .AND. Y(IY)<=580 .AND. Z(IZ)>=120 .AND. Z(IZ)<=140) THEN Mx(IX,IY,IZ)=Mx0 END IF !-------------------------------------------- IF(X(IX)>=250 .AND. X(IX)<=550 .AND. Y(IY)>=560 .AND. Y(IY)<=570 .AND. Z(IZ)>=130 .AND. Z(IZ)<=150) THEN Mx(IX,IY,IZ)=Mx0 END IF ``` ---
**Author:** ![milancurcic](https://yyz2.discourse-cdn.com/free1/user_avatar/fortran-lang.discourse.group/milancurcic/32/2_2.png) [@milancurcic](https://fortran-lang.discourse.group/u/milancurcic)\ **Post date:** [November 9, 2021, 4:32pm UTC](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101/50 "2021-11-09T16:32:36Z")
“Exploring pathways to more accurate machine learning emulation of atmospheric radiative transfer” by Peter Ukkonen > **[ESS Open Archive - Authorea](https://essopenarchive.org/)** > > The ESS Open Archive organization homepage on Authorea. > Abstract > Machine learning (ML) parameterizations of subgrid physics is a growing research area. A key question is whether traditional ML methods such as feed-forward neural networks (FNNs) are better suited for representing only specific processes. Radiation schemes are an interesting example, because they compute radiative flows through the atmosphere using well-established physical equations. The sequential aspect of the problem implies that FNNs may not be well-suited for it. > > This study explores whether emulating the entire radiation scheme is more difficult than its components without vertical dependencies. FNNs were trained to replace a shortwave radiation scheme, its gas optics component, and its reflectance-transmittance computations. In addition, a novel recurrent NN (RNN) method was developed to structurally incorporate the vertical dependence and sequential nature of radiation computations. > > It is found that a bidirectional RNN with an order of magnitude fewer model parameters than FNN is considerably more accurate, while offering a smaller but still significant 4-fold speedup over the original code on CPUs, and a much greater speedup on GPUs. The RNN predicts fluxes with less than 1% error, and heating rates computed from fluxes have a root-mean-square-error of 0.16 K day$^{-1}$ in offline tests using a year of global data. Finally, FNNs emulating gas optics are very accurate while being several times faster. As with RNNs emulating radiative transfer, the smaller dimensionality may be crucial for developing models that are general enough to be used as parameterizations. ---
**Author:** ![vmagnin](https://yyz2.discourse-cdn.com/free1/user_avatar/fortran-lang.discourse.group/vmagnin/32/28_2.png) [@vmagnin](https://fortran-lang.discourse.group/u/vmagnin)\ **Post date:** [November 12, 2021, 7:40pm UTC](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101/51 "2021-11-12T19:40:52Z")
Izmaylova, G.R., Yusupova, L.F. Numerical implementation of the mathematical model of combined effect on the formation using FORTRAN programming language (2021) _Journal of Physics: Conference Series,_ 2032 (1), art. no. 012042, . DOI: 10.1088/1742-6596/2032/1/012042 [https://iopscience.iop.org/article/10.1088/1742-6596/2032/1/012042](https://iopscience.iop.org/article/10.1088/1742-6596/2032/1/012042) A paper about oil production. FORTRAN in uppercase. They use Compaq Visual Fortran (Fortran\<=95) and the screenshot shows it’s fixed form. ---
**Author:** ![certik](https://yyz2.discourse-cdn.com/free1/user_avatar/fortran-lang.discourse.group/certik/32/4_2.png) [@certik](https://fortran-lang.discourse.group/u/certik)\ **Post date:** [November 12, 2021, 10:16pm UTC](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101/52 "2021-11-12T22:16:43Z")
Actually they call it “FORTAN” (in the caption of the figure with the source code). ---
**Author:** ![vmagnin](https://yyz2.discourse-cdn.com/free1/user_avatar/fortran-lang.discourse.group/vmagnin/32/28_2.png) [@vmagnin](https://fortran-lang.discourse.group/u/vmagnin)\ **Post date:** [November 13, 2021, 9:35am UTC](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101/53 "2021-11-13T09:35:24Z")
Another interesting point is that comments in the source code are in Russian which means that Compaq Visual Fortran could manage either [ISO 8859-5](https://en.wikipedia.org/wiki/ISO/IEC_8859-5) (Cyrillic) or UTF. Probably the former as CVF dates from the early 2000’s, and UTF was spread slowly although born in 1991. ---
**Author:** ![nbehrnd](https://avatars.discourse-cdn.com/v4/letter/n/e0b2c6/32.png) [@nbehrnd](https://fortran-lang.discourse.group/u/nbehrnd)\ **Post date:** [November 26, 2021, 4:24pm UTC](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101/54 "2021-11-26T16:24:51Z")
Anton Kasprzhitskii _et al._, «XMHFL: Software for calculating excited and ionized states of molecules by X-ray» in _SoftwareX_ **13** , _2021_, 100647 ([https://doi.org/10.1016/j.softx.2020.100647.](https://doi.org/10.1016/j.softx.2020.100647), open access). [GitHub repository](https://github.com/ElsevierSoftwareX/SOFTX-D-20-00056) of the Fortran 2003 program. Abstract: «XMHFL is a computer program that implements the Single Center Expansion (SCE) method to describe ground, excited and ionized states of a molecular system with one or more vacancies. The simulation mechanism is based on a combination of the three-point difference Numerov scheme and Thomas algorithm. The modular architecture of the XMHFL solvers allows using the generated code as the basis for the development of new methods and algorithms, as well as increasing the possibilities of molecular characterization calculation for X-ray molecular spectroscopy research and training as an educational platform.» ---
**Author:** ![Jcollins](https://yyz2.discourse-cdn.com/free1/user_avatar/fortran-lang.discourse.group/jcollins/32/540_2.png) [@Jcollins](https://fortran-lang.discourse.group/u/Jcollins)\ **Post date:** [November 27, 2021, 7:22pm UTC](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101/55 "2021-11-27T19:22:47Z")
We published several papers on software research in Fortran - see the last section in [SimCon - Software Research](http://simconglobal.com/research.html) . Should they be referenced here? ---
**Author:** ![certik](https://yyz2.discourse-cdn.com/free1/user_avatar/fortran-lang.discourse.group/certik/32/4_2.png) [@certik](https://fortran-lang.discourse.group/u/certik)\ **Post date:** [November 27, 2021, 8:06pm UTC](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101/56 "2021-11-27T20:06:00Z")
> [@Jcollins](#): > > Should they be referenced here? Maybe you can reference a few newest ones. There is obviously tons of Fortran articles in the past 50 years, but maybe we can use this thread to highlight recent work. ---
**Author:** ![nbehrnd](https://avatars.discourse-cdn.com/v4/letter/n/e0b2c6/32.png) [@nbehrnd](https://fortran-lang.discourse.group/u/nbehrnd)\ **Post date:** [November 28, 2021, 6:13pm UTC](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101/57 "2021-11-28T18:13:05Z")
@Beliavsky [In June](https://github.com/Beliavsky/Fortran-tools-on-GitHub), you mentioned your compilation of Fortran code maintained on [GitHub](https://github.com/Beliavsky/Fortran-tools-on-GitHub). Covering multiple fields and managed in sections, e.g., about file I/O, physics, or economics, it offers a showcase what Fortran can do (a bit like the demo pages of e.g., [gnuplot](http://gnuplot.sourceforge.net/demo/)) without aiming to be a package repository like [CTAN](https://ctan.org/). Would it be acceptable for you if - your compilation would contain links of the programs mentioned in this thread (possibly some already are), and - there were an early and prominent link towards your compilation, added to section «other resources/on the web» on fortran-lang’s page dedicated to [learning Fortran](https://fortran-lang.org/learn/) like there is e.g., about the examples by Scivision? ---
**Author:** ![Beliavsky](https://avatars.discourse-cdn.com/v4/letter/b/ba8739/32.png) [@Beliavsky](https://fortran-lang.discourse.group/u/Beliavsky)\ **Post date:** [December 16, 2021, 7:57pm UTC](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101/58 "2021-12-16T19:57:43Z")
The article can be downloaded. Another article from Nature, also in the area of biology and using CUDA Fortran, is also listed below. #### [Computer Physics Communications](https://www.sciencedirect.com/science/journal/00104655) [Volume 273](https://www.sciencedirect.com/science/journal/00104655/273/supp/C), April 2022, 108248 #### [FSEI-GPU: GPU accelerated simulations of the fluid–structure–electrophysiology interaction in the left heart](https://www.sciencedirect.com/science/article/pii/S001046552100360X) Authors: [FrancescoViola](https://www.sciencedirect.com/science/article/pii/S001046552100360X#!) [Vamsi Spandan](https://www.sciencedirect.com/science/article/pii/S001046552100360X#!) [Valentina Meschini](https://www.sciencedirect.com/science/article/pii/S001046552100360X#!) [Joshua Romero](https://www.sciencedirect.com/science/article/pii/S001046552100360X#!) [Massimiliano Fatica](https://www.sciencedirect.com/science/article/pii/S001046552100360X#!) [Marco D. de Tullio](https://www.sciencedirect.com/science/article/pii/S001046552100360X#!) [Roberto Verzicco](https://www.sciencedirect.com/science/article/pii/S001046552100360X#!) ## Abstract The reliability of cardiovascular computational models depends on the accurate solution of the [hemodynamics](https://www.sciencedirect.com/topics/chemistry/hemodynamic), the realistic characterization of the hyperelastic and electric properties of the tissues along with the correct description of their interaction. The resulting fluid–structure–electrophysiology interaction (FSEI) thus requires an immense computational power, usually available in large supercomputing centers, and requires long time to obtain results even if multi–CPU processors are used (MPI acceleration). In recent years, [graphics processing units](https://www.sciencedirect.com/topics/computer-science/graphics-processing-unit) (GPUs) have emerged as a convenient platform for [high performance computing](https://www.sciencedirect.com/topics/computer-science/high-performance-computing), as they allow for considerable reductions of the time–to–solution. This approach is particularly appealing if the tool has to support medical decisions that require solutions within reduced times and possibly obtained by local computational resources. Accordingly, our multi–physics solver [[1]](https://www.sciencedirect.com/science/article/pii/S001046552100360X#br0010) has been ported to GPU architectures using CUDA [Fortran](https://www.sciencedirect.com/topics/computer-science/fortran) to tackle fast and accurate hemodynamics simulations of the human heart without resorting to large–scale [supercomputers](https://www.sciencedirect.com/topics/earth-and-planetary-sciences/supercomputer). This work describes the use of CUDA to accelerate the FSEI on [heterogeneous clusters](https://www.sciencedirect.com/topics/computer-science/heterogeneous-cluster), where both the CPUs and GPUs are used in synergistically with minor modifications of the original source code. The resulting GPU [accelerated code](https://www.sciencedirect.com/topics/computer-science/accelerated-code) solves a single heartbeat within a few hours (from three to ten depending on the grid resolution) running on premises computing facility made of few GPU cards, which can be easily installed in a medical laboratory or in a hospital, thus opening towards a systematic [computational fluid dynamics](https://www.sciencedirect.com/topics/chemistry/computational-fluid-dynamics) (CFD) aided diagnostic. - [**Previous** article](https://www.sciencedirect.com/science/article/pii/S0010465521003647) - [**Next** article](https://www.sciencedirect.com/science/article/pii/S0010465521003659) #### Keywords: Fluid dynamics, Cardiovascular flows, Hemodynamics, Fluid-structure-interaction, Multiphysics model, Computational engineering * * * 2nd paper: GitHub: [Euplectella\_HPC](https://github.com/giacomofalcucci/Euplectella_HPC) - [Published: 21 July 2021](https://www.nature.com/articles/s41586-021-03658-1#article-info) #### [Extreme flow simulations reveal skeletal adaptations of deep-sea sponges](https://www.nature.com/articles/s41586-021-03658-1) - [Giacomo Falcucci](https://www.nature.com/articles/s41586-021-03658-1#auth-Giacomo-Falcucci), - [Giorgio Amati](https://www.nature.com/articles/s41586-021-03658-1#auth-Giorgio-Amati), - [Pierluigi Fanelli](https://www.nature.com/articles/s41586-021-03658-1#auth-Pierluigi-Fanelli), - [Vesselin K. Krastev](https://www.nature.com/articles/s41586-021-03658-1#auth-Vesselin_K_-Krastev), - [Giovanni Polverino](https://www.nature.com/articles/s41586-021-03658-1#auth-Giovanni-Polverino), - [Maurizio Porfiri](https://www.nature.com/articles/s41586-021-03658-1#auth-Maurizio-Porfiri) & - [Sauro Succi](https://www.nature.com/articles/s41586-021-03658-1#auth-Sauro-Succi) [_Nature_](https://www.nature.com/) **volume 595** , pages537–541 (2021) - 6097 Accesses - 2 Citations - 202 Altmetric #### Abstract Since its discovery[1](https://www.nature.com/articles/s41586-021-03658-1#ref-CR1),[2](https://www.nature.com/articles/s41586-021-03658-1#ref-CR2), the deep-sea glass sponge _Euplectella aspergillum_ has attracted interest in its mechanical properties and beauty. Its skeletal system is composed of amorphous hydrated silica and is arranged in a highly regular and hierarchical cylindrical lattice that begets exceptional flexibility and resilience to damage[3](https://www.nature.com/articles/s41586-021-03658-1#ref-CR3),[4](https://www.nature.com/articles/s41586-021-03658-1#ref-CR4),[5](https://www.nature.com/articles/s41586-021-03658-1#ref-CR5),[6](https://www.nature.com/articles/s41586-021-03658-1#ref-CR6). Structural analyses dominate the literature, but hydrodynamic fields that surround and penetrate the sponge have remained largely unexplored. Here we address an unanswered question: whether, besides improving its mechanical properties, the skeletal motifs of _E. aspergillum_ underlie the optimization of the flow physics within and beyond its body cavity. We use extreme flow simulations based on the ‘lattice Boltzmann’ method[7](https://www.nature.com/articles/s41586-021-03658-1#ref-CR7), featuring over fifty billion grid points and spanning four spatial decades. These in silico experiments reproduce the hydrodynamic conditions on the deep-sea floor where _E. aspergillum_ lives[8](https://www.nature.com/articles/s41586-021-03658-1#ref-CR8),[9](https://www.nature.com/articles/s41586-021-03658-1#ref-CR9),[10](https://www.nature.com/articles/s41586-021-03658-1#ref-CR10). Our results indicate that the skeletal motifs reduce the overall hydrodynamic stress and support coherent internal recirculation patterns at low flow velocity. These patterns are arguably beneficial to the organism for selective filter feeding and sexual reproduction[11](https://www.nature.com/articles/s41586-021-03658-1#ref-CR11),[12](https://www.nature.com/articles/s41586-021-03658-1#ref-CR12). The present study reveals mechanisms of extraordinary adaptation to live in the abyss, paving the way towards further studies of this type at the intersection between fluid mechanics, organism biology and functional ecology. ---
**Author:** ![MarDie](https://yyz2.discourse-cdn.com/free1/user_avatar/fortran-lang.discourse.group/mardie/32/54_2.png) [@MarDie](https://fortran-lang.discourse.group/u/MarDie)\ **Post date:** [December 17, 2021, 11:17am UTC](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101/59 "2021-12-17T11:17:53Z")
Saransh Singh, Farangis Ram and Marc De Graef EMsoft: open source software for electron diffraction/image simulations EMsoft is an open source package for the simulation of electron diffraction patterns and images for both SEM and TEM modalities. The package has been under development since the late 1990’s and is currently available as version 3.1 from a GitHub repository. The main component of the package is a source code library of functions and subroutines that cover crystallography, symmetry, electron scattering factors, geometry, rotation representations, basic image filters, defect displacement fields, reciprocal space sampling schemes, and so on. The code is mostly written in fortran-90, with oc- casional routines making use of the C-bindings that are available in fortran-2003. A second library covers routines used to read and write files in the open source HDF5 (Hierarchical Data Format) file format [2], which is the default format for all computational results. At the time of writing of this abstract, more than 200, 000 lines of source code have been written, about half of which are available via the source code repository > **[EMsoft: open source software for electron diffraction/image simulations |...](https://www.cambridge.org/core/journals/microscopy-and-microanalysis/article/emsoft-open-source-software-for-electron-diffractionimage-simulations/E708D64F1183455156010CEF369EB064)** > > EMsoft: open source software for electron diffraction/image simulations - Volume 23 Issue S1 > **[GitHub - EMsoft-org/EMsoft: Public EMsoft repository](https://github.com/EMsoft-org/EMsoft)** > > Public EMsoft repository. Contribute to EMsoft-org/EMsoft development by creating an account on GitHub. ---
**Author:** ![Beliavsky](https://avatars.discourse-cdn.com/v4/letter/b/ba8739/32.png) [@Beliavsky](https://fortran-lang.discourse.group/u/Beliavsky)\ **Post date:** [January 1, 2022, 9:08pm UTC](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101/60 "2022-01-01T21:08:22Z")
[Submitted to arXiv on 23 Dec 2021] ### [A fast computational model for the electrophysiology of the whole human heart](https://arxiv.org/abs/2112.12854) [Giulio Del Corso](https://arxiv.org/search/physics?searchtype=author&query=Del+Corso%2C+G), [Roberto Verzicco](https://arxiv.org/search/physics?searchtype=author&query=Verzicco%2C+R), [Francesco Viola](https://arxiv.org/search/physics?searchtype=author&query=Viola%2C+F) > In this study we present a novel computational model for unprecedented simulations of the whole cardiac electrophysiology. According to the heterogeneous electrophysiologic properties of the heart, the whole cardiac geometry is decomposed into a set of coupled conductive media having different topology and electrical conductivities: (i) a network of slender bundles comprising a fast conduction atrial network, the AV-node and the ventricular bundles; (ii) the Purkinje network; and (iii) the atrial and ventricular myocardium. The propagation of the action potential in these conductive media is governed by the bidomain/monodomain equations, which are discretized in space using an in-house finite volume method and coupled to three different cellular models, the Courtemanche model [1] for the atrial myocytes, the Stewart model [2] for the Purkinje Network and the ten Tusscher-Panfilov model [3] for the ventricular myocytes. The developed numerical model correctly reproduces the cardiac electrophysiology of the whole human heart in healthy and pathologic conditions and it can be tailored to study and optimize resynchronization therapies or invasive surgical procedures. Importantly, the whole solver is GPU-accelerated using CUDA Fortran providing an unprecedented speedup, thus opening the way for systematic parametric studies and uncertainty quantification analyses. Subjects: Medical Physics (physics.med-ph) [Previous page](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101.md?page=2) [Next page](https://fortran-lang.discourse.group/t/research-articles-using-fortran/1101.md?page=4)