Changeset 4748


Ignore:
Timestamp:
Nov 13, 2023, 10:44:33 PM (6 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.

Location:
BOL/LMDZ_Setup
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • BOL/LMDZ_Setup/era2gcm_tuto.sh

    r4615 r4748  
    109109   echo -------- liens de telechargement a actualiser ----
    110110   if [ "$rea" = "ERA5" ] ; then
     111     if [ $an -le 2021 ] ; then
    111112      ANA_DIR='/gpfsstore/rech/psl/rpsl376/ergon/ERA5/NETCDF/GLOBAL_025/4xdaily'
    112       suf="e5.GLOBAL_025"
     113      suf="aphe5.GLOBAL_025"
     114     else
     115      ANA_DIR='/gpfsstore/rech/psl/rpsl376/ergon/ERA5/NETCDF/GLOBAL_025/hourly'
     116      suf="ap1e5.GLOBAL_025"
     117     fi
    113118   elif [ "$rea" = "ERAI" ] ; then
    114119      #ANA_DIR='http://prodn.idris.fr/thredds/dodsC/ipsl_private/rpsl376/ERAI/NETCDF/GLOBAL_1125/4xdaily'
    115120      ANA_DIR='/gpfsstore/rech/psl/rpsl376/ergon/ERAI/NETCDF/GLOBAL_1125/4xdaily'
    116121      #/gpfsstore/rech/psl/rpsl376/ergon/ERAI/NETCDF/GLOBAL_075/4xdaily/AN_PL/2006/u.200601.aphei.GLOBAL_075.nc
    117       suf="ei.GLOBAL_1125"
     122      suf="aphei.GLOBAL_1125"
    118123   else  #OPERA
    119124      #ANA_DIR='http://prodn.idris.fr/thredds/dodsC/ipsl_private/rpsl376/OPERA/NETCDF/4xdaily'
    120125      ANA_DIR='/gpfsstore/rech/psl/rpsl376/ergon/OPERA/NETCDF/4xdaily'
    121       suf=".GLOBAL_1125"
     126      suf="aph.GLOBAL_1125"
    122127   fi
    123128
    … …  
    132137  # variables en format "short" doivent etre transformees en "float" via NCO
    133138  # This is done here with ncap2 ; also possible: "ncpdq --overwrite --unpack fin.nc fout.nc"
    134   fushort="$ANA_DIR/AN_PL/$an/u.$an$mois.aph$suf.nc"
    135   fvshort="$ANA_DIR/AN_PL/$an/v.$an$mois.aph$suf.nc"
    136   ftshort="$ANA_DIR/AN_PL/$an/ta.$an$mois.aph$suf.nc"
    137   fileu="$SCRATCH/u.$an$mois.aph$suf.nc"
    138   filev="$SCRATCH/v.$an$mois.aph$suf.nc"
    139   filet="$SCRATCH/ta.$an$mois.aph$suf.nc"
     139  fushort="$ANA_DIR/AN_PL/$an/u.$an$mois.$suf.nc"
     140  fvshort="$ANA_DIR/AN_PL/$an/v.$an$mois.$suf.nc"
     141  ftshort="$ANA_DIR/AN_PL/$an/ta.$an$mois.$suf.nc"
     142  fileu="$SCRATCH/u.$an$mois.$suf.nc"
     143  filev="$SCRATCH/v.$an$mois.$suf.nc"
     144  filet="$SCRATCH/ta.$an$mois.$suf.nc"
    140145  ncap2 -s 'u=float(u)' $fushort $fileu
    141146  ncap2 -s 'v=float(v)' $fvshort $filev
    142147  ncap2 -s 'ta=float(ta)' $ftshort $filet
    143 else
    144  fileu="$ANA_DIR/AN_PL/$an/u.$an$mois.aph$suf.nc"
    145  filev="$ANA_DIR/AN_PL/$an/v.$an$mois.aph$suf.nc"
    146  filet="$ANA_DIR/AN_PL/$an/ta.$an$mois.aph$suf.nc"
     148  #fileq pas teste ; a faire pour fileq aussi ?
     149  fileq="$ANA_DIR/AN_PL/$an/r.$an$mois.$suf.nc"
     150elif [ "$rea" = "ERA5" -a $an -ge 2022 ] ; then
     151  fu1h="$ANA_DIR/AN_PL/$an/u.$an$mois.$suf.nc"
     152  fv1h="$ANA_DIR/AN_PL/$an/v.$an$mois.$suf.nc"
     153  ft1h="$ANA_DIR/AN_PL/$an/ta.$an$mois.$suf.nc"
     154  fq1h="$ANA_DIR/AN_PL/$an/r.$an$mois.$suf.nc"
     155  fileu="$SCRATCH/u.$an$mois.$suf.nc"
     156  filev="$SCRATCH/v.$an$mois.$suf.nc"
     157  filet="$SCRATCH/ta.$an$mois.$suf.nc"
     158  fileq="$SCRATCH/r.$an$mois.$suf.nc"
     159  cdo selhour,0,6,12,18 $fu1h $fileu
     160  cdo selhour,0,6,12,18 $fv1h $filev
     161  cdo selhour,0,6,12,18 $ft1h $filet
     162  cdo selhour,0,6,12,18 $fq1h $fileq
     163else
     164 fileu="$ANA_DIR/AN_PL/$an/u.$an$mois.$suf.nc"
     165 filev="$ANA_DIR/AN_PL/$an/v.$an$mois.$suf.nc"
     166 filet="$ANA_DIR/AN_PL/$an/ta.$an$mois.$suf.nc"
     167 fileq="$ANA_DIR/AN_PL/$an/r.$an$mois.$suf.nc"
    147168fi
    148 fileq="$ANA_DIR/AN_PL/$an/r.$an$mois.aph$suf.nc"
    149169
    150170# verifier disponibilite des fichiers

  • BOL/LMDZ_Setup/test_ERAfiles_JeanZay.sh

    r4615 r4748  
    1717#
    1818# Choix des "rea"nalyses : ERA5, ERAI, OPERA
    19 rea=ERAI
     19rea=ERA5
    2020
    2121
    … …  
    2525   echo -------- liens de telechargement a actualiser ----
    2626   if [ "$rea" = "ERA5" ] ; then
     27     if [ $an -le 2021 ] ; then
    2728      ANA_DIR='/gpfsstore/rech/psl/rpsl376/ergon/ERA5/NETCDF/GLOBAL_025/4xdaily'
    28       suf="e5.GLOBAL_025"
     29      suf="aphe5.GLOBAL_025"
     30     else
     31      ANA_DIR='/gpfsstore/rech/psl/rpsl376/ergon/ERA5/NETCDF/GLOBAL_025/hourly'
     32      suf="ap1e5.GLOBAL_025"
     33     fi
    2934   elif [ "$rea" = "ERAI" ] ; then
    3035      #ANA_DIR='http://prodn.idris.fr/thredds/dodsC/ipsl_private/rpsl376/ERAI/NETCDF/GLOBAL_1125/4xdaily'
    3136      ANA_DIR='/gpfsstore/rech/psl/rpsl376/ergon/ERAI/NETCDF/GLOBAL_1125/4xdaily'
    3237      #/gpfsstore/rech/psl/rpsl376/ergon/ERAI/NETCDF/GLOBAL_075/4xdaily/AN_PL/2006/u.200601.aphei.GLOBAL_075.nc
    33       suf="ei.GLOBAL_1125"
     38      suf="aphei.GLOBAL_1125"
    3439   else  #OPERA
    3540      #ANA_DIR='http://prodn.idris.fr/thredds/dodsC/ipsl_private/rpsl376/OPERA/NETCDF/4xdaily'
    3641      ANA_DIR='/gpfsstore/rech/psl/rpsl376/ergon/OPERA/NETCDF/4xdaily'
    37       suf=".GLOBAL_1125"
     42      suf="aph.GLOBAL_1125"
    3843   fi
    3944
    … …  
    4651if [ "$rea" = "ERAI" ] ; then
    4752  # variables en format "short" doivent etre transformees en "float" via NCO
    48   fushort="$ANA_DIR/AN_PL/$an/u.$an$mois.aph$suf.nc"
    49   fvshort="$ANA_DIR/AN_PL/$an/v.$an$mois.aph$suf.nc"
    50   fileu="$SCRATCH/u.$an$mois.aph$suf.nc"
    51   filev="$SCRATCH/v.$an$mois.aph$suf.nc"
     53  # This is done here with ncap2 ; also possible: "ncpdq --overwrite --unpack fin.nc fout.nc"
     54  fushort="$ANA_DIR/AN_PL/$an/u.$an$mois.$suf.nc"
     55  fvshort="$ANA_DIR/AN_PL/$an/v.$an$mois.$suf.nc"
     56  ftshort="$ANA_DIR/AN_PL/$an/ta.$an$mois.$suf.nc"
     57  fileu="$SCRATCH/u.$an$mois.$suf.nc"
     58  filev="$SCRATCH/v.$an$mois.$suf.nc"
     59  filet="$SCRATCH/ta.$an$mois.$suf.nc"
    5260  ncap2 -s 'u=float(u)' $fushort $fileu
    5361  ncap2 -s 'v=float(v)' $fvshort $filev
     62  ncap2 -s 'ta=float(ta)' $ftshort $filet
     63  #fileq pas teste ; a faire pour fileq aussi ?
     64  fileq="$ANA_DIR/AN_PL/$an/r.$an$mois.$suf.nc"
     65elif [ "$rea" = "ERA5" -a $an -ge 2022 ] ; then
     66  fu1h="$ANA_DIR/AN_PL/$an/u.$an$mois.$suf.nc"
     67  fv1h="$ANA_DIR/AN_PL/$an/v.$an$mois.$suf.nc"
     68  ft1h="$ANA_DIR/AN_PL/$an/ta.$an$mois.$suf.nc"
     69  fq1h="$ANA_DIR/AN_PL/$an/r.$an$mois.$suf.nc"
     70  fileu="$SCRATCH/u.$an$mois.$suf.nc"
     71  filev="$SCRATCH/v.$an$mois.$suf.nc"
     72  filet="$SCRATCH/ta.$an$mois.$suf.nc"
     73  fileq="$SCRATCH/r.$an$mois.$suf.nc"
     74  cdo selhour,0,6,12,18 $fu1h $fileu
     75  cdo selhour,0,6,12,18 $fv1h $filev
     76  cdo selhour,0,6,12,18 $ft1h $filet
     77  cdo selhour,0,6,12,18 $fq1h $fileq
    5478else
    55  fileu="$ANA_DIR/AN_PL/$an/u.$an$mois.aph$suf.nc"
    56  filev="$ANA_DIR/AN_PL/$an/v.$an$mois.aph$suf.nc"
     79 fileu="$ANA_DIR/AN_PL/$an/u.$an$mois.$suf.nc"
     80 filev="$ANA_DIR/AN_PL/$an/v.$an$mois.$suf.nc"
     81 filet="$ANA_DIR/AN_PL/$an/ta.$an$mois.$suf.nc"
     82 fileq="$ANA_DIR/AN_PL/$an/r.$an$mois.$suf.nc"
    5783fi
    58 filet="$ANA_DIR/AN_PL/$an/ta.$an$mois.aph$suf.nc"
    59 fileq="$ANA_DIR/AN_PL/$an/r.$an$mois.aph$suf.nc"
    6084
    6185# verifier disponibilite des fichiers
Note: See TracChangeset for help on using the changeset viewer.