diff --git a/cmor/mpiesm/scripts/pismtest.cmor_ctl b/cmor/mpiesm/scripts/pismtest.cmor_ctl
index 3535befe052a50c182b10d2d48d3e23019ea16ac..d3343b2c46c109400a8b17373f1d8f6a4fbab052 100755
--- a/cmor/mpiesm/scripts/pismtest.cmor_ctl
+++ b/cmor/mpiesm/scripts/pismtest.cmor_ctl
@@ -10,10 +10,10 @@
 n_subsim=1
 #initial and final year of the experiment
 spinup=0
-iniyeararr=( $(( 4000 + spinup )) )
+iniyeararr=( $(( 4000 + $spinup )) )
 finyeararr=( 4999 )
-iniyear=$(( 1 + spinup ))
-finyear=1000
+iniyear=1
+finyear=$(( 1000 - $spinup ))
 #offset compared to startyear
 inioffsetarr=( -3999 )
 
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR.cmor_ctl b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR.cmor_ctl
index 313ce497ef500a5482bd4027cfeb1bbb6e8349dd..117533d641ab3748eafb36c33e283ac91c68dfdd 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR.cmor_ctl
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR.cmor_ctl
@@ -12,10 +12,11 @@ n_subsim=4
 spinup=1000
 iniyeararr=( $(( 3000 + spinup )) 2000 2000 2000 )
 finyeararr=( 7999 8999 8999 8999 )
-iniyear=$(( 1 + spinup ))
-finyear=26000
+iniyearoffset=0 # additional offset for iniyear
+iniyear=$(( $iniyearoffset + 1 ))
+finyear=$(( 26000 + $iniyearoffset - $spinup ))
 #offset compared to startyear
-inioffsetarr=(-2999 3001 10001 17001)
+inioffsetarr=( $(( -2999 + $iniyearoffset )) $(( 3001 + $iniyearoffset )) $(( 10001 + $iniyearoffset )) $(( 17001 + $iniyearoffset )) )
 
 # Time axis (26000-1 BP)
 ##     exp time       total yrs   new exp time   offset
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR.runpp b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR.runpp
index 9d78a7001428cd1fffe00ca86004d6d77b054a9f..996f9f71d703d664d436c544a63673f92cf0842a 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p1f1-CR.runpp
@@ -108,6 +108,7 @@ srfmod_dec_chunk="???0123124"
 srfmod_input_dec_chunk="???1123124"
 ocemod_dec_chunk="???1123124"
 icemod_dec_chunk="???1123124"
+ocemod_fxyear=3000
 
 #InfoTable(s)
 # Define here the "cdocmorinfo"-File or "eum"-Files
@@ -134,10 +135,11 @@ errdir_agg=${errdir}_agg
 #initial and final year of the experiment
 iniyeararr=($((3000+spinup)) 2000 2000 2000) # 26000-1 BP
 finyeararr=(7999 8999 8999 8999)
-iniyear=1
-finyear=$(( 26000 - $spinup ))
+iniyearoffset=0 # additional optional offset
+iniyear=$(( $iniyearoffset + 1 ))
+finyear=$(( 26000 + $iniyearoffset - $spinup ))
 #offset compared to startyear
-inioffsetarr=(-2999 3001 10001 17001) # 26000-1 BP
+inioffsetarr=( $(( -2999 + $iniyearoffset )) $(( 3001 + $iniyearoffset )) $(( 10001 + $iniyearoffset )) $(( 17001 + $iniyearoffset )) ) # 26000-1 BP
 # Time axis (26000-1 BP)
 ##     exp time       total yrs   new exp time   offset
 #  - 26000-21001 BP - 5000 yrs -     1- 5000 -   - 2999
@@ -384,6 +386,7 @@ function run_agg
             mkdir -p ${sdir}/out_aggr
 
             # Create softlinks to raw model output
+            echo "... creating links from '$rawsdir' for $period, $submodel"
             ln -sf $rawsdir/*_${period}* $sdir/ 2>/dev/null 1>&2 || { echo "ERROR creating links for $period ${RAW_EXP_IDS[$y]} $submodel" && exit 1 ; }
 
             # Load and run the aggr ScriptFragment
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR.cmor_ctl b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR.cmor_ctl
index 14f4a23ae3fa9ec60480ab051f17feadb96aadc7..d05bcd8a81c69f4cade52335813e5241d0b82b7b 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR.cmor_ctl
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR.cmor_ctl
@@ -12,10 +12,11 @@ n_subsim=4
 spinup=1000
 iniyeararr=( $(( 4000 + spinup )) 2000 2000 2000 )
 finyeararr=( 8999 8999 8999 8999 )
-iniyear=$(( 1 + spinup ))
-finyear=26000
+iniyearoffset=0 # additional offset for iniyear
+iniyear=$(( $iniyearoffset + 1 ))
+finyear=$(( 26000 + $iniyearoffset - $spinup ))
 #offset compared to startyear
-inioffsetarr=(-3999 3001 10001 17001)
+inioffsetarr=( $(( -3999 + $iniyearoffset )) $(( 3001 + $iniyearoffset )) $(( 10001 + $iniyearoffset )) $(( 17001 + $iniyearoffset )) )
 
 # Time axis (26000-1 BP)
 ##     exp time       total yrs   new exp time   offset
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR.runpp b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR.runpp
index 2ed21f33c7dcef000ffc123bf2f82001699aebf5..4d10d546f99737ecaba9ac1cade42068e9a7590b 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p2f2-CR.runpp
@@ -108,6 +108,7 @@ srfmod_dec_chunk="???0123124"
 srfmod_input_dec_chunk="???1123124"
 ocemod_dec_chunk="???1123124"
 icemod_dec_chunk="???1123124"
+ocemod_fxyear=4000
 
 #InfoTable(s)
 # Define here the "cdocmorinfo"-File or "eum"-Files
@@ -134,10 +135,11 @@ errdir_agg=${errdir}_agg
 #initial and final year of the experiment
 iniyeararr=($((4000+spinup)) 2000 2000 2000) # 26000-1 BP
 finyeararr=(8999 8999 8999 8999)
-iniyear=1
-finyear=$(( 26000 - $spinup ))
+iniyearoffset=0 # additional optional offset
+iniyear=$(( $iniyearoffset + 1 ))
+finyear=$(( 26000 + $iniyearoffset - $spinup ))
 #offset compared to startyear
-inioffsetarr=(-3999 3001 10001 17001) # 26000-1 BP
+inioffsetarr=( $(( -3999 + $iniyearoffset )) $(( 3001 + $iniyearoffset )) $(( 10001 + $iniyearoffset )) $(( 17001 + $iniyearoffset )) ) # 26000-1 BP
 # Time axis (26000-1 BP)
 ##     exp time       total yrs   new exp time   offset
 #a - 26000-21001 BP - 5000 yrs -     1- 5000 -   - 3999
@@ -384,6 +386,7 @@ function run_agg
             mkdir -p ${sdir}/out_aggr
 
             # Create softlinks to raw model output
+            echo "... creating links from '$rawsdir' for $period, $submodel"
             ln -sf $rawsdir/*_${period}* $sdir/ 2>/dev/null 1>&2 || { echo "ERROR creating links for $period ${RAW_EXP_IDS[$y]} $submodel" && exit 1 ; }
 
             # Load and run the aggr ScriptFragment
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR.cmor_ctl b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR.cmor_ctl
index 02b474415eb0f6bbf89e573142c2af414a327e9e..200d2ba74232e4e66ffe5c8c46bc9076996b1c29 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR.cmor_ctl
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR.cmor_ctl
@@ -12,10 +12,11 @@ n_subsim=4
 spinup=1000
 iniyeararr=( $(( 4000 + spinup )) 2000 2000 2000 )
 finyeararr=( 8999 8999 8999 8999 )
-iniyear=$(( 1 + spinup ))
-finyear=26000
+iniyearoffset=0 # additional offset for iniyear
+iniyear=$(( $iniyearoffset + 1 ))
+finyear=$(( 26000 + $iniyearoffset - $spinup ))
 #offset compared to startyear
-inioffsetarr=(-3999 3001 10001 17001)
+inioffsetarr=( $(( -3999 + $iniyearoffset )) $(( 3001 + $iniyearoffset )) $(( 10001 + $iniyearoffset )) $(( 17001 + $iniyearoffset )) )
 
 # Time axis (26000-1 BP)
 ##     exp time       total yrs   new exp time   offset
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR.runpp b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR.runpp
index 6db525f31d3d89c4edf09016bd2bcd88ee90b83f..4f2f0c3e4aa01dcfc323a0e87c3ed3bd2f1a7ec3 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-glac1d_r1i1p3f2-CR.runpp
@@ -108,6 +108,7 @@ srfmod_dec_chunk="???0123124"
 srfmod_input_dec_chunk="???1123124"
 ocemod_dec_chunk="???1123124"
 icemod_dec_chunk="???1123124"
+ocemod_fxyear=4000
 
 #InfoTable(s)
 # Define here the "cdocmorinfo"-File or "eum"-Files
@@ -134,10 +135,11 @@ errdir_agg=${errdir}_agg
 #initial and final year of the experiment
 iniyeararr=($((4000+spinup)) 2000 2000 2000) # 26000-1 BP
 finyeararr=(8999 8999 8999 8999)
-iniyear=1
-finyear=$(( 26000 - $spinup ))
+iniyearoffset=0 # additional optional offset
+iniyear=$(( $iniyearoffset + 1 ))
+finyear=$(( 26000 + $iniyearoffset - $spinup ))
 #offset compared to startyear
-inioffsetarr=(-3999 3001 10001 17001) # 26000-1 BP
+inioffsetarr=( $(( -3999 + $iniyearoffset )) $(( 3001 + $iniyearoffset )) $(( 10001 + $iniyearoffset )) $(( 17001 + $iniyearoffset )) ) # 26000-1 BP
 # Time axis (26000-1 BP)
 ##     exp time       total yrs   new exp time   offset
 #a - 26000-21001 BP - 5000 yrs -     1- 5000 -   - 3999
@@ -384,6 +386,7 @@ function run_agg
             mkdir -p ${sdir}/out_aggr
 
             # Create softlinks to raw model output
+            echo "... creating links from '$rawsdir' for $period, $submodel"
             ln -sf $rawsdir/*_${period}* $sdir/ 2>/dev/null 1>&2 || { echo "ERROR creating links for $period ${RAW_EXP_IDS[$y]} $submodel" && exit 1 ; }
 
             # Load and run the aggr ScriptFragment
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR.cmor_ctl b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR.cmor_ctl
index 352cb7ffcb8653ed0622d657f0d293752922d333..190472bd45329b5b6b8bf4a03df0a7c6903cc757 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR.cmor_ctl
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR.cmor_ctl
@@ -12,10 +12,11 @@ n_subsim=4
 spinup=1000
 iniyeararr=( $(( 3000 + spinup )) 2000 2000 2000 )
 finyeararr=( 7999 8999 8999 8999 )
-iniyear=$(( 1 + spinup ))
-finyear=26000
+iniyearoffset=0 # additional offset for iniyear
+iniyear=$(( $iniyearoffset + 1 ))
+finyear=$(( 26000 + $iniyearoffset - $spinup ))
 #offset compared to startyear
-inioffsetarr=(-2999 3001 10001 17001)
+inioffsetarr=( $(( -2999 + $iniyearoffset )) $(( 3001 + $iniyearoffset )) $(( 10001 + $iniyearoffset )) $(( 17001 + $iniyearoffset )) )
 
 # Time axis (26000-1 BP)
 ##     exp time       total yrs   new exp time   offset
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR.runpp b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR.runpp
index b32ded7fc37306f57814d3421673a7710e728588..0126e5bbbcd3716f8e20f9fdede336978ef6a165 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p1f1-CR.runpp
@@ -108,6 +108,7 @@ srfmod_dec_chunk="???0123124"
 srfmod_input_dec_chunk="???1123124"
 ocemod_dec_chunk="???1123124"
 icemod_dec_chunk="???1123124"
+ocemod_fxyear=3000
 
 #InfoTable(s)
 # Define here the "cdocmorinfo"-File or "eum"-Files
@@ -134,10 +135,11 @@ errdir_agg=${errdir}_agg
 #initial and final year of the experiment
 iniyeararr=($((3000+spinup)) 2000 2000 2000) # 26000-1 BP
 finyeararr=(7999 8999 8999 8999)
-iniyear=1
-finyear=$(( 26000 - $spinup ))
+iniyearoffset=0 # additional optional offset
+iniyear=$(( $iniyearoffset + 1 ))
+finyear=$(( 26000 + $iniyearoffset - $spinup ))
 #offset compared to startyear
-inioffsetarr=(-2999 3001 10001 17001) # 26000-1 BP
+inioffsetarr=( $(( -2999 + $iniyearoffset )) $(( 3001 + $iniyearoffset )) $(( 10001 + $iniyearoffset )) $(( 17001 + $iniyearoffset )) ) # 26000-1 BP
 # Time axis (26000-1 BP)
 ##     exp time       total yrs   new exp time   offset
 #  - 26000-21001 BP - 5000 yrs -     1- 5000 -   - 2999
@@ -384,6 +386,7 @@ function run_agg
             mkdir -p ${sdir}/out_aggr
 
             # Create softlinks to raw model output
+            echo "... creating links from '$rawsdir' for $period, $submodel"
             ln -sf $rawsdir/*_${period}* $sdir/ 2>/dev/null 1>&2 || { echo "ERROR creating links for $period ${RAW_EXP_IDS[$y]} $submodel" && exit 1 ; }
 
             # Load and run the aggr ScriptFragment
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR.cmor_ctl b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR.cmor_ctl
index 6b4923291ba05a6e1477a46a05a42946fdc365e5..b8225136a51c68d3d0299b26f4d6640df9b54997 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR.cmor_ctl
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR.cmor_ctl
@@ -12,10 +12,11 @@ n_subsim=4
 spinup=1000
 iniyeararr=( $(( 4000 + spinup )) 2000 2000 2000 )
 finyeararr=( 8999 8999 8999 8999 )
-iniyear=$(( 1 + spinup ))
-finyear=26000
+iniyearoffset=0 # additional offset for iniyear
+iniyear=$(( $iniyearoffset + 1 ))
+finyear=$(( 26000 + $iniyearoffset - $spinup ))
 #offset compared to startyear
-inioffsetarr=(-3999 3001 10001 17001)
+inioffsetarr=( $(( -3999 + $iniyearoffset )) $(( 3001 + $iniyearoffset )) $(( 10001 + $iniyearoffset )) $(( 17001 + $iniyearoffset )) )
 
 # Time axis (26000-1 BP)
 ##     exp time       total yrs   new exp time   offset
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR.runpp b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR.runpp
index ec23b3e9f56003026268a7a48e99d9cbb78c7731..5a3b0947f1c66fa915d9aa92f0b787105d542dbf 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p2f2-CR.runpp
@@ -108,6 +108,7 @@ srfmod_dec_chunk="???0123124"
 srfmod_input_dec_chunk="???1123124"
 ocemod_dec_chunk="???1123124"
 icemod_dec_chunk="???1123124"
+ocemod_fxyear=4000
 
 #InfoTable(s)
 # Define here the "cdocmorinfo"-File or "eum"-Files
@@ -134,10 +135,11 @@ errdir_agg=${errdir}_agg
 #initial and final year of the experiment
 iniyeararr=($((4000+spinup)) 2000 2000 2000) # 26000-1 BP
 finyeararr=(8999 8999 8999 8999)
-iniyear=1
-finyear=$(( 26000 - $spinup ))
+iniyearoffset=0 # additional optional offset
+iniyear=$(( $iniyearoffset + 1 ))
+finyear=$(( 26000 + $iniyearoffset - $spinup ))
 #offset compared to startyear
-inioffsetarr=(-3999 3001 10001 17001) # 26000-1 BP
+inioffsetarr=( $(( -3999 + $iniyearoffset )) $(( 3001 + $iniyearoffset )) $(( 10001 + $iniyearoffset )) $(( 17001 + $iniyearoffset )) ) # 26000-1 BP
 # Time axis (26000-1 BP)
 ##     exp time       total yrs   new exp time   offset
 #a - 26000-21001 BP - 5000 yrs -     1- 5000 -   - 3999
@@ -384,6 +386,7 @@ function run_agg
             mkdir -p ${sdir}/out_aggr
 
             # Create softlinks to raw model output
+            echo "... creating links from '$rawsdir' for $period, $submodel"
             ln -sf $rawsdir/*_${period}* $sdir/ 2>/dev/null 1>&2 || { echo "ERROR creating links for $period ${RAW_EXP_IDS[$y]} $submodel" && exit 1 ; }
 
             # Load and run the aggr ScriptFragment
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR.cmor_ctl b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR.cmor_ctl
index 758b0daa9f4195149e04f6d9223b871ee8e5dc52..1a151f7fca8bcf5ca03518367db3a40777c78cd9 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR.cmor_ctl
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR.cmor_ctl
@@ -12,10 +12,11 @@ n_subsim=4
 spinup=1000
 iniyeararr=( $(( 4000 + spinup )) 2000 2000 2000 )
 finyeararr=( 8999 8999 8999 8999 )
-iniyear=$(( 1 + spinup ))
-finyear=26000
+iniyearoffset=0 # additional offset for iniyear
+iniyear=$(( $iniyearoffset + 1 ))
+finyear=$(( 26000 + $iniyearoffset - $spinup ))
 #offset compared to startyear
-inioffsetarr=(-3999 3001 10001 17001)
+inioffsetarr=( $(( -3999 + $iniyearoffset )) $(( 3001 + $iniyearoffset )) $(( 10001 + $iniyearoffset )) $(( 17001 + $iniyearoffset )) )
 
 # Time axis (26000-1 BP)
 ##     exp time       total yrs   new exp time   offset
diff --git a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR.runpp b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR.runpp
index 2d9e694f3688e2550f91287312db3534190b8939..f0edf3cf558bf513903d9c469831cc5317ae3027 100755
--- a/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR.runpp
+++ b/cmor/mpiesm/scripts/transient-deglaciation-prescribed-ice6g_r1i1p3f2-CR.runpp
@@ -108,6 +108,7 @@ srfmod_dec_chunk="???0123124"
 srfmod_input_dec_chunk="???1123124"
 ocemod_dec_chunk="???1123124"
 icemod_dec_chunk="???1123124"
+ocemod_fxyear=4000
 
 #InfoTable(s)
 # Define here the "cdocmorinfo"-File or "eum"-Files
@@ -134,10 +135,11 @@ errdir_agg=${errdir}_agg
 #initial and final year of the experiment
 iniyeararr=($((4000+spinup)) 2000 2000 2000) # 26000-1 BP
 finyeararr=(8999 8999 8999 8999)
-iniyear=1
-finyear=$(( 26000 - $spinup ))
+iniyearoffset=0 # additional optional offset
+iniyear=$(( $iniyearoffset + 1 ))
+finyear=$(( 26000 + $iniyearoffset - $spinup ))
 #offset compared to startyear
-inioffsetarr=(-3999 3001 10001 17001) # 26000-1 BP
+inioffsetarr=( $(( -3999 + $iniyearoffset )) $(( 3001 + $iniyearoffset )) $(( 10001 + $iniyearoffset )) $(( 17001 + $iniyearoffset )) ) # 26000-1 BP
 # Time axis (26000-1 BP)
 ##     exp time       total yrs   new exp time   offset
 #a - 26000-21001 BP - 5000 yrs -     1- 5000 -   - 3999
@@ -384,6 +386,7 @@ function run_agg
             mkdir -p ${sdir}/out_aggr
 
             # Create softlinks to raw model output
+            echo "... creating links from '$rawsdir' for $period, $submodel"
             ln -sf $rawsdir/*_${period}* $sdir/ 2>/dev/null 1>&2 || { echo "ERROR creating links for $period ${RAW_EXP_IDS[$y]} $submodel" && exit 1 ; }
 
             # Load and run the aggr ScriptFragment