User Tools

Site Tools

Trace:
initial_conditions:oceanic

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision 		Previous revision
initial_conditions:oceanic [2016/01/08 14:37]
eexarcho [[[how_to_generate_new_ocean_initial_conditions_|Generation of new ocean initial condition]]] 		initial_conditions:oceanic [2022/08/23 18:08] (current)
etourign
Line 1: Line 1:
  
-==== Brief description ====+==== Shared in EC-Earth consortium ====
  
 +Details are here https://dev.ec-earth.org/projects/ecearth3/wiki/Ocean_initial_conditions_for_climate_predictions
 +
 +==== Brief description ====
  
 Please note that Ec-Earth2 includes NEMO2.2 and Ec-Earth3 includes NEMO3.3 Please note that Ec-Earth2 includes NEMO2.2 and Ec-Earth3 includes NEMO3.3
  
 ^ ID^ OCEAN MODEL^ GRID^ FIRST^ LAST^ INTERVAL (MONTH)^ MEMBERS^ Comments^ References | ^ ID^ OCEAN MODEL^ GRID^ FIRST^ LAST^ INTERVAL (MONTH)^ MEMBERS^ Comments^ References |
-| s4| NEMO2.2 / NEMO3.2 / NEMO3.3| ORCA1 and ORCAL46| 19580101| 20141201| 02, 05, 08, 11| fc0 - f4| ocean reanalysis - mask corrected for ORCA1 - vertically interpolated and extrapolated for ORCA1L46| Mogensen et al(2011), Balmaseda et al (2012) |+| s4| NEMO3.6 | ORCA1L46| 19600101| 20160701| all months | fc0 - f4| ocean reanalysis - mask corrected for ORCA1 - vertically interpolated and extrapolated for ORCA1L75| Mogensen et al(2011), Balmaseda et al (2012) | 
 +| s4| NEMO2.2 / NEMO3.2 / NEMO3.3| ORCA1 and ORCAL46| 19580101| 20141201| 02, 05, 08, 11| fc0 - f4| ocean reanalysis - mask corrected for ORCA1 - vertically interpolated and extrapolated for ORCA1L46 and ORCA1L75| Mogensen et al(2011), Balmaseda et al (2012) |
 | GLORYS2V1| NEMO3.3| ORCA1L46| 19921101| 20091101| 05,11| fc0| Interpolated initial conditions from an ocean reanalysis| Ferry et al (2010) | | GLORYS2V1| NEMO3.3| ORCA1L46| 19921101| 20091101| 05,11| fc0| Interpolated initial conditions from an ocean reanalysis| Ferry et al (2010) |
 | GLORYS2V1| NEMO3.3| ORCA025L46| 19921101| 20091101| 05,11| fc0| Interpolated initial conditions from an ocean reanalysis| Ferry et al (2010) | | GLORYS2V1| NEMO3.3| ORCA025L46| 19921101| 20091101| 05,11| fc0| Interpolated initial conditions from an ocean reanalysis| Ferry et al (2010) |
Line 17: Line 21:
 | m01w| NEMO3.3| ORCA025L46| 19921101| 20091101| 02,05,08,11| fc0| Initial conditions produced by nudging the ocean component in EC-Earth3.0.1 toward the monthly T and S from GLORYS2V1 ocean reanalysis| no reference yet | | m01w| NEMO3.3| ORCA025L46| 19921101| 20091101| 02,05,08,11| fc0| Initial conditions produced by nudging the ocean component in EC-Earth3.0.1 toward the monthly T and S from GLORYS2V1 ocean reanalysis| no reference yet |
 | m01u| NEMO3.3| ORCA025L75| 19921101| 20091101| 02,05,08,11| fc0| Initial conditions produced by nudging the ocean component in EC-Earth3.0.1 toward the monthly T and S from GLORYS2V1 ocean reanalysis| no reference yet | | m01u| NEMO3.3| ORCA025L75| 19921101| 20091101| 02,05,08,11| fc0| Initial conditions produced by nudging the ocean component in EC-Earth3.0.1 toward the monthly T and S from GLORYS2V1 ocean reanalysis| no reference yet |
- |+| official_repro| NEMO3.6| ORCA01L75| 19500101 | 19500101 | --- | fc00-fc04| Initial conditions for the official EC-Earth reproducibility test| https://github.com/plesager/ece3-postproc#reproducibility-test | 
 +| a1e1/a1ue| NEMO3.6| ORCA1| 19600101| 20151201 | all |fc00-fc04|5-category sea ice product using DFS5.2 surface forcing fields, surface restoring and subsurface nudging to ORAS4 with EC-Earth3.2 (NEMO-LIM only sims). Member fc00 should be taken from a1ue and members fc01 to fc04 should be taken from a1e1. For 2015 only member fc00 is available, so use a22b instead for 2015. | no reference | 
 +| a1n1| NEMO3.6| ORCA025| 19600101| 20151201 | all |fc04-fc00|5-category sea ice product using DFS5.2 surface forcing fields, surface restoring and subsurface nudging to ORAS5 with EC-Earth3.2 (NEMO-LIM only sims). | In progress | 
 +| a22b| NEMO3.6| ORCA1| 20100101| 20181201 | all |fc00-fc04|5-category sea ice product using ERA-Interim surface forcing fields, surface restoring and subsurface nudging to ORAS4 with EC-Earth3.2 (NEMO-LIM only sims). | no reference | 
 +| a2fd| NEMO3.6| ORCA1| 20180101| 20191201 | all |fc01-fc05|5-category sea ice product using ERA5 surface forcing fields, surface restoring and subsurface nudging to ORAS5 with EC-Earth3.2 (NEMO-LIM only sims). | no reference |
  
  
Line 45: Line 53:
  
 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. 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.
-==== How_to_generate_new_ocean_initial_conditions_|Generation of new ocean initial condition ====+==== How_to_generate_new_ocean_initial_conditions ====
  
-== ORAP5, ORCA025L75 == +== ORAS4, ORCA1L42 and ORCA1L75 == 
  
 1. Update the ecaccess-certificate-create tool, which gives 7 days of access to ECMWF repositories.  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**)+2. clone the project https://earth.bsc.es/gitlab/ec-earthCPWG/Oce_IC.git
  
-  #!/bin/bash +3. run  export PYTHONPATH=$PYTHONPATH:$path_to_the_Oce_IC
-  #SBATCH -n 1 +
-  #SBATCH -t 10-12:00:00 +
-  #SBATCH -J fetchORAP5 +
-  #SBATCH --mem=60000 +
-  #SBATCH -o fetch-orap5-%j.log +
- +
  
-  set -evx +4. go into the path_to_the_Oce_IC and run fetch_ORAS4.sh after filling the start dates you want to download from EMCWF and interpolate.
-  TMPDIR=/scratch/$(whoami)/$$ +
-  mkdir -p $TMPDIR +
-  cd $TMPDIR+
  
-  module load ECAC+== ORAP5, ORCA025L75 == 
  
-  DATASOURCE="ec:/RDX/prepIFS/g0tx" # Root of the path where to find ORAS4 on ECFS +1. Update the ecaccess-certificate-create tool, which gives 7 days of access to ECMWF repositories. 
-  YEARI=2000 # First eyar to download   +
-  YEARF=2013 # Last year to download +
-  MEMBER_LST="0" # List of members to download +
-  YEAR=$YEARI +
  
-  workdir=`pwd` +2. git clone https://earth.bsc.es/gitlab/es/cpg_tools.git -b develop-eleftheria 
-  +
  
-  echo $YEAR +2. Fetch the data locally (for example when in fat nodes) 
-  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+
  
 +  cd cpg_tools/prep_restarts/generate_restarts_orap5_HR
 +  sbatch dl_restarts.sh
 +  
 +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**.       
  
-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:**+- Also, we check the size in the original data, because some dates **have half size**, which means that they do not have all the variables we need for interpolation. For example, the restart of 19930430 is not suitable for using as initial condition:
  
   ll /esnas/reconstructions/ecmwf/orap5/restarts/1993/*199304*   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 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 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 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 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 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 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+  -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: +After the manual checking, we run the following script:   
 +  cd cpg_tools/prep_restarts/generate_restarts_orap5_HR 
 +  sbatch copy_orap5.sh 
 +   
 +   
 +  
 +== ORAP5, ORCA1L46 == 
  
 +We use a software called **SOSIE** [[http://sosie.sourceforge.net/]] to interpolate the data. We do the following steps (on fat nodes):
  
 +  cd cpg_tools/prep_restarts/generate_restarts_orap5_LR
 +  sbatch launchscript.sh 
  
-dir=/esnas/reconstructions/ecmwf/orap5/restarts/  +  
-outdir=/esnas/releases/ic/ocean/ORCA025L75/ORAP5 +== GLOSEA5, ORCA025L75 == 
-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 +
  
 +We need to extrapolate the files to the mask we use locally. We do the following separately per year (here as example for 19940501, in the fat nodes)
  
 +   cd cpg_tools/prep_restarts/generate_restarts_gs5_HR/19940501
 +   sbatch 19940501.sh 
  
-# November start dates  +== GLOSEA5, ORCA1L46 == 
- 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    +
  
 +Similarly, as in ORAP5 ORCA1L46, we need to interpolate by using a software called **SOSIE** [[http://sosie.sourceforge.net/]] to interpolate the data. We do the following steps on the fat nodes: 
  
 +   sbatch launchscript.sh
initial_conditions/oceanic.1452263823.txt.gz · Last modified: 2016/01/08 14:37 by eexarcho

Page Tools

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki