Equation with constants and two variables that are ranges of values

2 views (last 30 days)
Erin Summerlin-Donofrio
Erin Summerlin-Donofrio on 11 Sep 2023
Commented: Star Strider on 20 Sep 2023
My final equation is dependent on variables I set and define, as well as variables that are dependent on a range of values, Beta and T. I would like the equation to run for example, with Beta=0 and T=473.15 to T=1473.15 in 100 degree steps. Basically I would like to calculate a drf value for each Beta value from 0.0 to 1.0 in 0.1 increments, and for each T from 473.15 to 1473.15. I've been doing this with naming each: Beta1 and T1 individually and calculating drf1, while naming all of the other variables with 1's as well. There must be a more straightforward way...?
Beta=(0.0:0.10:1.0);
T=(473.15:100:1473.15);
BigDeltaPlgH2O = (2.91-0.76*Beta) * (10^6/T^2) - 3.41 - 0.41*Beta;
BigDelta=BigDeltaPlgH2O;
d18Oin=0;
d18Or=5.5;
%Q=Mwater/Mo
Q=(0.01:0.01:10);
fw=0.2;
%Conc=Cro/Cwo
Conc=0.5;
Qo=(fw+(1-fw)*Conc);
drf=(BigDelta+d18Oin+(d18Or-BigDelta-d18Oin)*exp(-Q/Qo));

Sign in to answer this question.

Accepted Answer

Star Strider
Star Strider on 11 Sep 2023
Ran in:
Try this —
Beta=(0.0:0.10:1.0);
T=(473.15:100:1473.15);
BigDeltaPlgH2O = (2.91-0.76*Beta) .* (10^6./T.^2) - 3.41 - 0.41*Beta;
BigDelta=BigDeltaPlgH2O;
d18Oin=0;
d18Or=5.5;
%Q=Mwater/Mo
Q=(0.01:0.01:10);
VectorLengths = [numel(Beta); numel(T); numel(Q)] % Potential Problem With Vector Lengths
VectorLengths = 3×1
11 11 1000
fw=0.2;
%Conc=Cro/Cwo
Conc=0.5;
Qo=(fw+(1-fw)*Conc);
% drf=(BigDelta+d18Oin+(d18Or-BigDelta-d18Oin)*exp(-Q/Qo));
drf=(BigDelta+d18Oin+(d18Or-BigDelta-d18Oin)).*exp(-Q/Qo).'; % <— LIGHTLY EDITED
figure
surfc(T, Q, drf)
colormap(turbo)
colorbar
xlabel('T (°K)')
ylabel('Q')
.

More Answers (1)

Erin Summerlin-Donofrio
Erin Summerlin-Donofrio on 20 Sep 2023
%Oxygen
BetaA=0.0;%plgAn#
TA = 373.15;%100C
BigDeltaPlgH20A = (2.91-0.76*BetaA) * (10^6/TA^2) - 3.41 - 0.41*BetaA;
BigDeltaA=BigDeltaPlgH20A;
d18OinA=0;
d18OrA=6;
%Q=Mwater/Mo
Q=(0.001:0.001:10);
fw=0.2;
%Conc=Cro/Cwo
ConcA=0.5;
QoA=(fw+(1-fw)*ConcA);
drfA=(BigDeltaA+d18OinA+(d18OrA-BigDeltaA-d18OinA)*exp(-Q/QoA));
%Hydrogen
TA = 373.15;%100C
BigDeltaEpH20A = 9.3 * (10^6/TA^2) - 61.9;
DBigDeltaA=BigDeltaEpH20A;
dDinA=0;
dDrA=-75;
%Q=Mwater/Mo
Q=(0.001:0.001:500);
fw=0.2;
%Conc=Cro/Cwo=0.5/11.11
DConcA=0.0045;
DQoA=(fw+(1-fw)*DConcA);
DdrfA=(DBigDeltaA+dDinA+(dDrA-DBigDeltaA-dDinA)*exp(-Q/DQoA));
How can I modify the code to vary (T, Beta, d18Or, and Q) for the first equation and (T, dDr, and Q) for the second equation. My ultimate goal is to chart how moving the d18Or and dDr values changes the trajectory of the curves when plotting drf vs dDrf. Beta will vary between 0.0 and 1.0 so that's 11 values as above. My other variables I wish to vary need to be an array of 11 as well yes? Thank you!
  1 Comment
Star Strider
Star Strider on 20 Sep 2023
Ran in:
I am a bit lost.
To get them all as 11-element arrays, use the linspace function. So for example:
N = 11;
Tv = linspace(min(T), max(T), N);
Betav = linspace(min(Beta), max(Beta), N);
and so for the rest of them. Then, if you want to, use the ndgrid function to create matrices from the vectors. You can then use those matrices as arguments to the appropriate calculations, avoiding nested loops. The problem is that the result will have the same dimensions — (11 x 11 x 11) in this instance — and that could be a challenge to plot. It might be necessary to use hold and then loop through the ‘pages’of the result matrix, plotting each ‘page separately, in layers in the surf plot.
A simplified version of that mnight go something like this —
N = 11;
Tv = linspace(153.15, 453.15, N);
Betav = linspace(0, 1, N);
d18Orv = linspace(0, 5, N);
Qv = linspace(0.001, 10, N);
[Tm,Bm,Om,Qm] = ndgrid(Tv, Betav, d18Orv, Qv);
TA = Tm;
BetaA = Bm;
d18OrA = Om;
Q = Qm;
%Oxygen
BetaA=0.0;%plgAn#
TA = 373.15;%100C
BigDeltaPlgH20A = (2.91-0.76*BetaA) * (10^6./TA.^2) - 3.41 - 0.41*BetaA;
BigDeltaA=BigDeltaPlgH20A;
d18OinA=0;
% d18OrA=6;
%Q=Mwater/Mo
% % Q=(0.001:0.001:10);
fw=0.2;
%Conc=Cro/Cwo
ConcA=0.5;
QoA=(fw+(1-fw)*ConcA);
drfA=(BigDeltaA+d18OinA+(d18OrA-BigDeltaA-d18OinA).*exp(-Q/QoA))
drfA =
drfA(:,:,1,1) = 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 0.0291 drfA(:,:,2,1) = 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 0.5283 drfA(:,:,3,1) = 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 1.0275 drfA(:,:,4,1) = 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 1.5266 drfA(:,:,5,1) = 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 2.0258 drfA(:,:,6,1) = 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 2.5250 drfA(:,:,7,1) = 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 3.0241 drfA(:,:,8,1) = 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 3.5233 drfA(:,:,9,1) = 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 4.0225 drfA(:,:,10,1) = 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 4.5216 drfA(:,:,11,1) = 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 5.0208 drfA(:,:,1,2) = 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 14.1907 drfA(:,:,2,2) = 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 14.2850 drfA(:,:,3,2) = 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 14.3793 drfA(:,:,4,2) = 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 14.4736 drfA(:,:,5,2) = 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 14.5679 drfA(:,:,6,2) = 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 14.6622 drfA(:,:,7,2) = 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 14.7565 drfA(:,:,8,2) = 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 14.8508 drfA(:,:,9,2) = 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 14.9451 drfA(:,:,10,2) = 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 15.0394 drfA(:,:,11,2) = 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 15.1337 drfA(:,:,1,3) = 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 16.8660 drfA(:,:,2,3) = 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 16.8838 drfA(:,:,3,3) = 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 16.9016 drfA(:,:,4,3) = 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 16.9194 drfA(:,:,5,3) = 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 16.9372 drfA(:,:,6,3) = 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 16.9550 drfA(:,:,7,3) = 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 16.9728 drfA(:,:,8,3) = 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 16.9906 drfA(:,:,9,3) = 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 17.0085 drfA(:,:,10,3) = 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 17.0263 drfA(:,:,11,3) = 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 17.0441 drfA(:,:,1,4) = 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 17.3713 drfA(:,:,2,4) = 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 17.3747 drfA(:,:,3,4) = 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 17.3781 drfA(:,:,4,4) = 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 17.3814 drfA(:,:,5,4) = 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 17.3848 drfA(:,:,6,4) = 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 17.3882 drfA(:,:,7,4) = 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 17.3915 drfA(:,:,8,4) = 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 17.3949 drfA(:,:,9,4) = 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 17.3982 drfA(:,:,10,4) = 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 17.4016 drfA(:,:,11,4) = 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 17.4050 drfA(:,:,1,5) = 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 17.4668 drfA(:,:,2,5) = 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 17.4674 drfA(:,:,3,5) = 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 17.4681 drfA(:,:,4,5) = 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 17.4687 drfA(:,:,5,5) = 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693 17.4693
figure
hold on
for k = 1:numel(d18Orv)
surf(drfA(:,:,k,1)) % It Is Only Possible To Plot 2D Matrices With 'surf', So This Is A Compromise, Choose Any Two Dimensions To Plot
end
hold off
colormap(turbo)
view(30,30)
I am not certain what ranges your variables have or what you want to do with them. This is simply an illustration. I must necessarily leave the rest to you, since I have no idea what you are doing.
.

Sign in to comment.

Categories

Find more on Install Products in Help Center and File Exchange

Tags

Products


Release

R2023a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!