Changeset 4748 for BOL/LMDZ_Setup/test_ERAfiles_JeanZay.sh

Timestamp:

Nov 13, 2023, 10:44:33 PM (11 months ago)

Author:

asima

Message:

For ERA5 : starting 2022, 6-hourly files are gradually not saved anymore in the IPSL repository, only the hourly ones are available.
So the script starts by extracting the hours 0,6,12 and 18 from the hourly files.

File:

• BOL/LMDZ_Setup/test_ERAfiles_JeanZay.sh (modified) (3 diffs)

BOL/LMDZ_Setup/test_ERAfiles_JeanZay.sh

@@ -17 +17 @@
 #
 # Choix des "rea"nalyses : ERA5, ERAI, OPERA
-rea=ERAI
+rea=ERA5
 
 
@@ -25 +25 @@
    echo -------- liens de telechargement a actualiser ----
    if [ "$rea" = "ERA5" ] ; then
+     if [ $an -le 2021 ] ; then 
       ANA_DIR='/gpfsstore/rech/psl/rpsl376/ergon/ERA5/NETCDF/GLOBAL_025/4xdaily'
-      suf="e5.GLOBAL_025"
+      suf="aphe5.GLOBAL_025"
+     else
+      ANA_DIR='/gpfsstore/rech/psl/rpsl376/ergon/ERA5/NETCDF/GLOBAL_025/hourly'
+      suf="ap1e5.GLOBAL_025"
+     fi 
    elif [ "$rea" = "ERAI" ] ; then
       #ANA_DIR='http://prodn.idris.fr/thredds/dodsC/ipsl_private/rpsl376/ERAI/NETCDF/GLOBAL_1125/4xdaily'
       ANA_DIR='/gpfsstore/rech/psl/rpsl376/ergon/ERAI/NETCDF/GLOBAL_1125/4xdaily'
       #/gpfsstore/rech/psl/rpsl376/ergon/ERAI/NETCDF/GLOBAL_075/4xdaily/AN_PL/2006/u.200601.aphei.GLOBAL_075.nc 
-      suf="ei.GLOBAL_1125"
+      suf="aphei.GLOBAL_1125"
    else  #OPERA
       #ANA_DIR='http://prodn.idris.fr/thredds/dodsC/ipsl_private/rpsl376/OPERA/NETCDF/4xdaily'
       ANA_DIR='/gpfsstore/rech/psl/rpsl376/ergon/OPERA/NETCDF/4xdaily'
-      suf=".GLOBAL_1125"
+      suf="aph.GLOBAL_1125"
    fi
 
@@ -46 +51 @@
 if [ "$rea" = "ERAI" ] ; then
   # variables en format "short" doivent etre transformees en "float" via NCO 
-  fushort="$ANA_DIR/AN_PL/$an/u.$an$mois.aph$suf.nc"
-  fvshort="$ANA_DIR/AN_PL/$an/v.$an$mois.aph$suf.nc"
-  fileu="$SCRATCH/u.$an$mois.aph$suf.nc"
-  filev="$SCRATCH/v.$an$mois.aph$suf.nc"
+  # This is done here with ncap2 ; also possible: "ncpdq --overwrite --unpack fin.nc fout.nc"
+  fushort="$ANA_DIR/AN_PL/$an/u.$an$mois.$suf.nc"
+  fvshort="$ANA_DIR/AN_PL/$an/v.$an$mois.$suf.nc"
+  ftshort="$ANA_DIR/AN_PL/$an/ta.$an$mois.$suf.nc"
+  fileu="$SCRATCH/u.$an$mois.$suf.nc"
+  filev="$SCRATCH/v.$an$mois.$suf.nc"
+  filet="$SCRATCH/ta.$an$mois.$suf.nc"
   ncap2 -s 'u=float(u)' $fushort $fileu
   ncap2 -s 'v=float(v)' $fvshort $filev
+  ncap2 -s 'ta=float(ta)' $ftshort $filet
+  #fileq pas teste ; a faire pour fileq aussi ?
+  fileq="$ANA_DIR/AN_PL/$an/r.$an$mois.$suf.nc"
+elif [ "$rea" = "ERA5" -a $an -ge 2022 ] ; then
+  fu1h="$ANA_DIR/AN_PL/$an/u.$an$mois.$suf.nc"
+  fv1h="$ANA_DIR/AN_PL/$an/v.$an$mois.$suf.nc"
+  ft1h="$ANA_DIR/AN_PL/$an/ta.$an$mois.$suf.nc"
+  fq1h="$ANA_DIR/AN_PL/$an/r.$an$mois.$suf.nc"
+  fileu="$SCRATCH/u.$an$mois.$suf.nc"
+  filev="$SCRATCH/v.$an$mois.$suf.nc"
+  filet="$SCRATCH/ta.$an$mois.$suf.nc"
+  fileq="$SCRATCH/r.$an$mois.$suf.nc"
+  cdo selhour,0,6,12,18 $fu1h $fileu
+  cdo selhour,0,6,12,18 $fv1h $filev
+  cdo selhour,0,6,12,18 $ft1h $filet
+  cdo selhour,0,6,12,18 $fq1h $fileq
 else
- fileu="$ANA_DIR/AN_PL/$an/u.$an$mois.aph$suf.nc"
- filev="$ANA_DIR/AN_PL/$an/v.$an$mois.aph$suf.nc"
+ fileu="$ANA_DIR/AN_PL/$an/u.$an$mois.$suf.nc"
+ filev="$ANA_DIR/AN_PL/$an/v.$an$mois.$suf.nc"
+ filet="$ANA_DIR/AN_PL/$an/ta.$an$mois.$suf.nc"
+ fileq="$ANA_DIR/AN_PL/$an/r.$an$mois.$suf.nc"
 fi
-filet="$ANA_DIR/AN_PL/$an/ta.$an$mois.aph$suf.nc"
-fileq="$ANA_DIR/AN_PL/$an/r.$an$mois.aph$suf.nc"
 
 # verifier disponibilite des fichiers