Changeset 4749 for BOL/LMDZ_Setup/test_ERAfiles_JeanZay.sh

Timestamp:

Nov 14, 2023, 2:23:10 PM (8 months ago)

Author:

asima

Message:

Script optimisation :
a/ era2gcm_tuto.sh : preprocessing of ERA files, if needed, is only done for the requested variables (cf options guide_*=y) instead of the full list (u,v,ta,q)
b/ test_ERAfiles_JeanZay.sh : only list the required ERA files, don't do any preprocessing

File:

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

BOL/LMDZ_Setup/test_ERAfiles_JeanZay.sh

@@ -25 +25 @@
    echo -------- liens de telechargement a actualiser ----
    if [ "$rea" = "ERA5" ] ; then
-     if [ $an -le 2021 ] ; then 
+     if [ $an -ge 2022 ] ; then 
+      ANA_DIR='/gpfsstore/rech/psl/rpsl376/ergon/ERA5/NETCDF/GLOBAL_025/hourly'
+      suf="ap1e5.GLOBAL_025"
+     else
       ANA_DIR='/gpfsstore/rech/psl/rpsl376/ergon/ERA5/NETCDF/GLOBAL_025/4xdaily'
       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
@@ -43 +43 @@
    fi
 
-varu=u
-varv=v
-vart=ta # peut etre parfois juste "t"
-varq=q
-#varp=msl
 
-if [ "$rea" = "ERAI" ] ; then
-  # variables en format "short" doivent etre transformees en "float" via NCO 
-  # 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.$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
 
 # verifier disponibilite des fichiers