The documentation for blupf90 is available separately as a PDF document. The following lines can be added at the end of parameter file for blupf90: Available options: OPTION conv_crit 1e-12 Set convergence criteria (deault 1e-12). OPTION maxrounds 10000 Set maximum number of rounds (default 1000). OPTION solv_method FSPAK Selection solutions by FSPAK, SOR or PCG (default PCG). OPTION r_factor 1.6 Set relaxation factor for SOR (default 1.4). OPTION sol se Store solutions and s.e. OPTION missing -999 Missing observation (default 0). OPTION residual y-hat and residual will be shown in yhat_residual. OPTION blksize 3 Set block size for preconditioner (default 1).