I recently started using GSLIB and I am trying to generate a 2D surface using a
histogram and the power law variogram with SASIM.
I am getting the following error and I am guessing that I am not doing
something correct.
run-time error M6203: MATH
- **: OVERFLOW error
I am also unclear on how to specify the 'w' parameter (i.e., shape) in the
parameter file for the power law. My goal is to simulate a non-stationary 2D
surface for a given histogram.
Example of my PAR file:
Below is an example of my PAR file:
Parameters for SASIM
********************
START OF PARAMETERS:
1 1 0 0 0 \components: hist,varg,ivar,corr,cpdf
1 4 0 0 0 \weight: hist,varg,ivar,corr,cpdf
0 \0=no transform, 1=log transform
1 \number of realizations
200 15.0 30.0 \grid definition: nx,xmn,xsiz
200 15.0 30.0 \ ny,ymn,ysiz
1 15.0 30.0 \ nz,zmn,zsiz
69069 \random number seed
3 \debugging level
C:\Data_314\Proj_Ras_Accuracy\Applications\GSLIB_Jobs\Data\sasim.dbg \file
for debugging output
C:\Data_314\Proj_Ras_Accuracy\Applications\GSLIB_Jobs\Data\sasim.out \file
for simulation output
1 \schedule (0=automatic,1=set below)
0 0.00001 4000000 400000 5 0.0000001 \ schedule:
t0,redfac,ka,k,num,Omin
1000 10 \ maximum perturbations, reporting
50 \ maximum number without a change
0 \conditioning data:(0=no, 1=yes)
../xx \ file with data
1 2 0 3 \ columns: x,y,z,attribute
-1.0e21 1.0e21 \ trimming limits
1 \file with histogram:(0=no, 1=yes)
C:\Data_314\Proj_Ras_Accuracy\Applications\GSLIB_Jobs\Data\test_noCoords.dat
\file with histogram
1 1 \ column for value and weight
1 \ number of quantiles for obj. func.
0 \number of indicator variograms
10.0 20.0 30.0 40.0 50.0 \ indicator thresholds
../xx \file with gridded secondary data
1 \ column number
1 \ vertical average (0=no, 1=yes)
0.7 \correlation coefficient
../xx \file with paired data
2 1 0 \ columns for primary, secondary, wt
-0.5 100.0 \ minimum and maximum
5 \ number of primary thresholds
5 \ number of secondary thresholds
6400 \Variograms: number of lags ([(2K+1)^2]-1);
k=range/cell size
1 \ standardize sill (0=no,1=yes)
1 0.05 \ nst, nugget effect
4 .95 90.0 0.0 0.0 0.1 1 \ it,cc,ang1,ang2,ang3,anis1, anis2
1500.0 1400.0 0.0 \ a_hmax, a_hmin, a_vert 1500.0 800.0 0.0
Thank you for your assistance,
Michael
+
+ To post a message to the list, send it to ai-geost...@jrc.ec.europa.eu
+ To unsubscribe, send email to majordomo@ jrc.ec.europa.eu with no subject and
"unsubscribe ai-geostats" in the message body. DO NOT SEND
Subscribe/Unsubscribe requests to the list
+ As a general service to list users, please remember to post a summary of any
useful responses to your questions.
+ Support to the forum can be found at http://www.ai-geostats.org/
