Closed LanPeng-94 closed 1 month ago
Hi,
your field in Fig. 1 does not look lognormally distributed. Try replacing your line gs.transform.normal_to_lognormal(srf)
with srf.transform('normal_to_lognormal')
.
Sorry, I don't know what you mean by "the number of truncations for KL".
Hi,
Below is a new field based on "srf.transform('normal_to_lognormal')".
Dear,
I am utilizing Gstools to create a random field for a slope, and I have retrieved my grid data with coordinates X=data[:,0:1] and Y=data[:,1:2]. Assuming the generation of a random field with a saturated permeability coefficient Ks, and given it follows a lognormal distribution with a mean of (2 \times 10^{-5}) and a standard deviation of (1.6 \times 10^{-5}), it equates to a normal distribution with a mean of -11.067 and a variance of 0.495. I have composed the following code based on these assumptions.
Is this correct? Additionally, how can I manage the number of truncations for KL, and what is the default setting?
(Figure 1)
Below is the result I generated (Figure 1), which doesn't seem right because I have a copy of the reference solution (Figure 2). My results are much steeper. In addition, the reference solution has a truncation of 8. (Figure 2)