移流拡散のfortranプログラミング
土壌中のイオンの深度別濃度分布を移流拡散して、新たな濃度量を求めたいのですが、プログラミングがわからなくて困っています。
濃度分布ですが、
縦軸は深度で、単位はkg/m2,横軸は濃度でmol/kgです。移流の単位がkg/m2/yr,拡散の単位がkg2/m4/yrです。
求めたいものは毎年の深度10kg/m2間隔での、イオンの量(濃度×深度)(mol/m2)です。
しかし、私は初めの年の濃度分布を、10kg/m2間隔で層を作り、イオンの量(mol/m2)を格納しています。
というのも、次の年では上から新たなイオン量が降ってくるので、その降ってきたものを一番表層の0-10kg/m2の層に足して、層ごとにイオン量を格納したいのです。なので、次の年の、移流拡散が起こった時の各層のイオン量が知りたいのですが、どう計算していいかわかりません。
濃度と深度の計算式はy=10exp(-0.3x)(y=濃度、x=深度)です。
普通であれば、この式に移流拡散を入れて、積分すれば、10kg/m2間隔でのイオンの量(mol/m2)は求められると思うのですが、イオンの量をすでに層ごとに格納している場合はどうすればよいか、どなたかアドバイスいただけないでしょうか?
恐縮ですが、よろしくお願い致します。
補足
何度もご回答ありがとうございます。 極端な例で,標高100mと標高200mの柱状図でも,各データ伴に深度0から考えて良いということですよね・・?解釈が間違っていたら申し訳ありません。