前フリ 卒論の時(もう7年前かー)に書いた、レイリーテイラー不安定による電離圏プラズマバブル発生シミュレーションのプログラムをgithubに上げてみた。Fortran90す。 yumu19/RayleighTaylorInstability · GitHub 一部抜粋 ! -- イオンの数密度 [1/m^3] の計算 -- z1=(dble(z)*dz-zpeak+zbot)/HO(z) ! n1(z)=npeak*cos(chi)*exp(1.0-z1-exp(-z1)/cos(chi)) do i=2,datanumn if (datan(1,i) > height) then n1(z)=1.0e6*(datan(2,i-1) & & +(height-datan(1,i-1)) & & *(datan(2,i)-datan(2,i-1)) & & /(datan(1,i)-data