commit before TOMACS version
Output includes all replications now
Check Vhat dimension before diag
add GESKExtrpolate function
commit before merge