XIMUL package « readme file » (version 2001) Author : R.ABABOU

This is a brief "readme" file for using XIMUL (version 2001) _________________________________________________________ Fortran source files of XIMUL: CSI-MUL.FCSI-TURN99.FCSI-GAUSSIAN123.F

or CSI-EXPONENTIAL123.F (or any other covariance/spectrum subroutine)

Example compilation of XIMUL: xlf -o ximul CSI-MUL.F CSI-TURN99.F CSI-GAUSSIAN123.F Using the program file CSI_NOUS.F instead of CSI-MUL.F gives the output-data in the correct format for xmgr (at least in the case of 1D fields). Thus: xlf -o ximul CSI_NOUS.F CSI-TURN99.F CSI-GAUSSIAN123.F

Random number generator in subroutine CSI_TURN99.F: In the subroutine CSI_TURN99.F, the random number generator (RNG) of the LPT3D code is implemented. LPT3D is a 3D Lagrangian Particle Tracking program (R.Ababou, A.Fadili, M.Spiller 1997-2001). This RNG may lead to better results than the previous version. Specifying or changing the RNG parameters in the input-file "csinput" are not necessary, as all important constants are direcly fixed inside the code.

XIMUL package « readme file » (version 2001) Author : R.ABABOU

Input file "csiname" The file CSINAME contains the names of input files, e.g., for unconditional random field generation, only the "csinput" file name needs to be given. In the example below, we have "csinput"="csin_big". =======================================================================

"csiname" file example : ----------------------------------------------------------------------This file, csiname, is the first input file read by the csimul program: 1rst line LINE No. 11 CONTAINS SEQUENCE # FOR SIMULTANEOUS RUNS (Ldummy = 1 to 5). LINES No. 12,13,14,15 CONTAIN THE NAMES OF INPUT FILES AS FOLLOWS: Line 12 : Basic input file ............................... (mandatory) Line 13 : Data file, i.e. measurement values and locations (optional) Line 14 : Grid file, i.e. node numbers and locations ..... (optional) Line 15: Turning Band field as input .................... (optional) *************** 123456789012345 Use no more than 15 characters for the name: *************** csifield.outcsifield.out 1 csin_big / basic input file defining the problem and options 'csindata' --> measurement values and nodal locations (only if M > 0 ) 'csingrid' --> location and ordering of nodes (only for irregular grid) 'csinturn' --> separate input field from turning band (rarely used) OUTPUTS: 'csiout.out' --> short output file containing basic options, date, etc. 'csiturn.out' --> output field(s) from turning band module (only on option) 'csifield.out' --> output field(s) from csimul (conditional realizations) 'csistatTB.out' --> output statistics for unconditional fields (Turning Band) 'csiblumean.out' --> output theoretical conditional mean field (B.L.U.E.) 'csiblucovar.out' --> output theoretical posterior covariance matrix 'csibluvar.out' --> output theoretical posterior variance = diag(covar.matrix) 'csimean.out' --> output computed conditional mean over MC realizations 'csivar.out' --> output computed conditional variance over MC realizations

XIMUL package « readme file » (version 2001) Author : R.ABABOU

Input file "csinput" (or any other specified name): =========================================================================== "csinput" file example (the name of this file was actually "csin_big"): -----------------------------------------------------------------------/********************************************************** ⇐ 1RST LINE /BASIC INPUT FILE FOR CONDITIONAL SIMULATION CODE "csimul": /********************************************************** / /---------------------------------------------------------------/Global data and options for Conditional Simulation (Csimul): /---------------------------------------------------------------/ 050799,01/ IDATE,INUM (date,identification number) 1.,1.,1.,1.,1./GMEAN,SIGMY,CORLL1,CORLL2,CORLL3 0.1,0.1,0.25/ DX1,DX2,DX3 1000,1000,1/ N1 ,N2 ,N3 (must be compatible with LDIM) 1000000/ NTOTAL (must be compatible with Ni ) 0/ MDATA (number of data points, may be zero) 1/ MCARLO (number of Monte-Carlo realizations) / 2/ LDIM (dimensionality of space: 1,2,3) 1/ LTURN (1 for Turning Band, 0 for reading existing r.f.) 1/ LSTATB (TB stat output:-1/0/1/2=No/Last(short)/Last(long)/All(long) 0/ LFORM (i/o format: 0=formatted ascii, 1=unformatted binary) 0/ LGRID (0=regular grid, 1=irregular grid) 1/ LOGNOR (1=generate log-normal field, 0=generate normal field) 1/ LISOTR (1= isotropic input field, 0=non-isotropic input field) 0/ LSTATC (Conditional stat outputs: -1/0/1 = None/LastOne/All) 0/ LFIELD (Conditional field outputs:-1/0/1 = None/LastOne/All) 1/ LCOVAR (Posterior covar matrix output:-1/0/1/2 = No/Diag/Full/-) / /-------------------------------------------------------------------/Additional inputs for the Turning Band module (Turnip): /-------------------------------------------------------------------/ 10000/ NLINE (number of T.B. lines to be generated in m-dim space) 0.05,1000000/ DELZET,NZTEST (mesh for line processes, estimated No.points) 62.83,0.04/ BIGK,DK (wnumber cutoff and increment, normalized values) 117899/ ISEED (initial seed for random number generator: caution) 1048576,1027 /MRAND,LRAND (R.N.G modulus & multiplier: 0,0 for default) 1,0,0,0/ LIRAND,LINEW,LICENT,LIFFT (some options for T.B. algorithms) / /---------------------------------------------------------------------/Additional inputs for statistical analysis module (Rfstat): /---------------------------------------------------------------------/ -3,+3,100/ YLOWER,YUPPER,NPDF (p.d.f. histogram: bounds, No. intervals) / /---------------------------------------------------------------------/End of file.