Reaction diffusion equation with Schnakenberg kinetics in 2 dim space
Parallel-computed by MPI using two nodes

<< Graphics`Graphics3D`

SetDirectory["/Users/takasu/home/\:60c5\:5831\:79d1\:5b66\:79d1\:306e\
\:4ed5\:4e8b/\:8b1b\:7fa9/\:5e73\:6210\:ff12\:ff11\:5e74\:5ea6/H21 \:8a08\
\:7b97\:6a5f\:5b9f\:9a132/MPI/testMPI/build/Release/"]

データの読み込み

size  =  51;

    
datau0  =  
        ReadList["RD-2dim-Schnakenberg-u-rank0", {Real, Real,  
            Real}];
datau0  =  
        Partition[Partition[datau0, size], size];

    
datau1  =  
        ReadList["RD-2dim-Schnakenberg-u-rank1", {Real, Real, 
            Real}];
datau1  =  
        Partition[Partition[datau1, size], size];

    
datav0  =  
        ReadList["RD-2dim-Schnakenberg-v-rank0", {Real, Real,  
            Real}];
datav0  =  
        Partition[Partition[datav0, size], size];

    
datav1  =  
        ReadList["RD-2dim-Schnakenberg-v-rank1", {Real, Real, 
            Real}];
datav1  =  
        Partition[Partition[datav1, size], 
          size];


Length[datau0]

2つのファイルを結合

datauAll  =  {};
Do[
AppendTo[ datauAll, 
        Join[datau0[[i]], datau1[[i]]] ], {i, 1, 
        Length[datau0]}
]

    
datavAll  =  {};
Do[
AppendTo[ datavAll, 
        Join[datav0[[i]], datav1[[i]]] ], {i, 1, 
        Length[datav0]}
]

描画

ListSurfacePlot3D[datauAll[[1]], BoxRatios -> {2, 4, 1}, 
      Axes -> True, 
      AxesLabel -> {"x", "y", "u"}]

ListSurfacePlot3D[datavAll[[1]], BoxRatios -> {2, 4, 1}, 
      Axes -> True, 
      AxesLabel -> {"x", "y", "v"}]

ListSurfacePlot3D[Last[datauAll], BoxRatios -> {2, 4, 1}, 
      Axes -> True, PlotRange -> {0, 2}, 
      AxesLabel -> {"x", "y", "u"}]

ListSurfacePlot3D[Last[datavAll], BoxRatios -> {2, 4, 1}, 
      Axes -> True, PlotRange -> {0, 2}, 
      AxesLabel -> {"x", "y", "v"}]

Do[
ListSurfacePlot3D[datauAll[[i]], 
        BoxRatios -> {2, 4, 1}, Axes -> True, 
        PlotRange -> {0, 2}, 
        AxesLabel -> {"x", "y", "u"}], {i, 1, 
        Length[datauAll], 10}
]

Do[
ListSurfacePlot3D[datavAll[[i]], 
        BoxRatios -> {2, 4, 1}, Axes -> True, 
        PlotRange -> {0, 2}, 
        AxesLabel -> {"x", "y", "v"}], {i, 1, 
        Length[datavAll], 10}
]

1