Please note that Ec-Earth2 includes NEMO2.2 and Ec-Earth3 includes NEMO3.3

|

These are ranked by priority :

1) ORAS4 ORCA1L75 configuration for 1958 to present

2) GLORYS2v1 ORCA1L75 configuration for 1993 to 2009

3) ORAS4 all configurations whole period

4) GLORYS2V3 all configurations whole period

Balmaseda M.A., Mogensen K.S., Weaver A.T., 2012, Evaluation of the ECMWF Ocean Reanalysis ORAS4. Quart J Roy Meteor Soc DOI 10.1002/qj.2063.

Ferry N et al., 2010, Mercator global eddy permitting ocean reanalysis GLORYS1v1: description and results. Mercator Ocean Quart. Newsl, 36, 15-27.

Mogensen K.S., Balmaseda M.A., Weaver A., 2011, The NEMOVAR ocean data assimilation as implemented in the ECMWF ocean analysis for system 4. ECMWF Technical Memorandum, 668.

Mogensen K.S., Balmaseda M.A., Weaver A., Martin M.J., Vidard A., 2009, NEMOVAR: A variational data assimilation system for the NEMO ocean model. ECMWF newsletter.

Hao Z., Balmaseda M. A., Mogensen K., 2015, The new eddy-permitting ORAP5 ocean reanalysis: description, evaluation and uncertainties in climate signals. Climate Dynamics (2015): 1-21.

1. Update the ecaccess-certificate-create tool, which gives 7 days of access to ECMWF repositories.

2. Fetch the data locally by running the script.sh below (in moore or amdhahl with sbatch script.sh)

#!/bin/bash
#SBATCH -n 1
#SBATCH -t 10-12:00:00
#SBATCH -J fetchORAP5
#SBATCH --mem=60000
#SBATCH -o fetch-orap5-%j.log

set -evx
TMPDIR=/scratch/$(whoami)/$$
mkdir -p $TMPDIR
cd $TMPDIR

module load ECAC

DATASOURCE="ec:/RDX/prepIFS/g0tx" # Root of the path where to find ORAS4 on ECFS
YEARI=2000 # First eyar to download  
YEARF=2013 # Last year to download
MEMBER_LST="0" # List of members to download
YEAR=$YEARI 

workdir=`pwd`

echo $YEAR
for year in `seq $YEARI $YEARF`; do
mkdir -p /esnas/reconstructions/ecmwf/orap5/restarts/$YEAR
for mem in $MEMBER_LST ; do
  for file in $( ecaccess-file-dir ${DATASOURCE}/opa${mem}/restart/${YEAR}/g0tx_${YEAR}*_*_restart.tar.gz )
  do
    if [[ "$file" != "ECFS" ]];then
      echo $file
      if [[ ! -f /esnas/reconstructions/ecmwf/orap5/restarts/${YEAR}/$file ]]; then
        while [[ ! -f /esnas/reconstructions/ecmwf/orap5/restarts/${YEAR}/$file ]]
        do
          ecaccess-file-get -progress ${DATASOURCE}/opa${mem}/restart/${YEAR}/$file $file
          mv $file /esnas/reconstructions/ecmwf/orap5/restarts/${YEAR}/.
        done
      fi
    fi
  done
done
YEAR=$((YEAR+1))
done

3. Copy the data using the correct naming, in the repository of the initial conditions. Pay attention that for initial conditions for 1st May and 1st November we use the closest available dates in the original data, which are not always the same for each year. Also, we check the size in the original data, because some dates have half size, so they do not have all the variables we need for interpolation, for example the restart for 19930430 is not suitable for using as initial condition:

ll /esnas/reconstructions/ecmwf/orap5/restarts/1993/*199304*

-rw-r--r-- 1 pbretonn Earth 9184196770 Oct  7 04:31 /esnas/reconstructions/ecmwf/orap5/restarts/1993/g0tx_19930401_00021960_restart.tar.gz

-rw-r–r– 1 pbretonn Earth 9181022277 Oct 7 04:54 /esnas/reconstructions/ecmwf/orap5/restarts/1993/g0tx_19930406_00022320_restart.tar.gz -rw-r–r– 1 pbretonn Earth 9183757715 Oct 7 05:17 /esnas/reconstructions/ecmwf/orap5/restarts/1993/g0tx_19930411_00022680_restart.tar.gz -rw-r–r– 1 pbretonn Earth 9181268805 Oct 7 05:45 /esnas/reconstructions/ecmwf/orap5/restarts/1993/g0tx_19930416_00023040_restart.tar.gz -rw-r–r– 1 pbretonn Earth 9178049931 Oct 7 06:13 /esnas/reconstructions/ecmwf/orap5/restarts/1993/g0tx_19930421_00023400_restart.tar.gz -rw-r–r– 1 pbretonn Earth 9179559914 Oct 7 06:35 /esnas/reconstructions/ecmwf/orap5/restarts/1993/g0tx_19930426_00023760_restart.tar.gz -rw-r–r– 1 pbretonn Earth 4927336766 Oct 7 06:50 /esnas/reconstructions/ecmwf/orap5/restarts/1993/g0tx_19930430_00024048_restart.tar.gz

After the manual checking, we run the following script:

dir=/esnas/reconstructions/ecmwf/orap5/restarts/ outdir=/esnas/releases/ic/ocean/ORCA025L75/ORAP5 curdir=`pwd`

for year in {2005..2007} ; do ##May start dates

echo $year 
dirin=${dir}/${year} 
if [ -e $dirin  ] ; then   
	cd $dirin 
              for file in $( ls   g0tx_${year}0428* ) ; do 
		inname=$(echo $file | cut -f1 -d"." )
		ofile=ORAP5_fc0_${year}0430_restart.nc.gz
		if [ ! -e  ${outdir}/${ofile}  ]  ; then   
		cd $curdir 
		gunzip -c   ${dirin}/$file  >  ${inname}.tar  
		tar xfv ${inname}.tar  
		gzip -c ${inname}.nc >  ${outdir}/${ofile} 
		rm -rf  ${curdir}/g0tx*tar  ${curdir}/g0t*nc 
	fi 
done 
fi 

# November start dates

echo $year 
dirin=${dir}/${year} 
if [ -e $dirin  ] ; then   
	cd $dirin 
              for file in $( ls   g0tx_${year}1030* ) ; do 
		inname=$(echo $file | cut -f1 -d"." )
		ofile=ORAP5_fc0_${year}1031_restart.nc.gz
		if [ ! -e  ${outdir}/${ofile}  ]  ; then   
		cd $curdir 
		gunzip -c   ${dirin}/$file  >  ${inname}.tar  
		tar xfv ${inname}.tar  
		gzip -c ${inname}.nc >  ${outdir}/${ofile} 
		rm -rf  ${curdir}/g0tx*tar  ${curdir}/g0t*nc 
	fi 
done 
fi 

done