| 1 | # Script pour tester la presence de fichiers de guidage ERA* sur Jean-Zay |
|---|
| 2 | # 2021/04/15 |
|---|
| 3 | # A. Sima (adriana.sima@lmd.ipsl.fr) |
|---|
| 4 | |
|---|
| 5 | |
|---|
| 6 | ##################################################################### |
|---|
| 7 | # A choisir par l'utilisateur : |
|---|
| 8 | ##################################################################### |
|---|
| 9 | |
|---|
| 10 | # Guidage en vent(u&v) et/ou temperature, humidite |
|---|
| 11 | guide_uv=y |
|---|
| 12 | guide_t=y |
|---|
| 13 | guide_q=y |
|---|
| 14 | |
|---|
| 15 | an=2000 |
|---|
| 16 | mois=01 |
|---|
| 17 | # |
|---|
| 18 | # Choix des "rea"nalyses : ERA5, ERAI, OPERA |
|---|
| 19 | rea=ERAI |
|---|
| 20 | |
|---|
| 21 | |
|---|
| 22 | ##################################################################### |
|---|
| 23 | # Lien avec les fichiers netcdf contenant les d0 ECMWF |
|---|
| 24 | ##################################################################### |
|---|
| 25 | echo -------- liens de telechargement a actualiser ---- |
|---|
| 26 | if [ "$rea" = "ERA5" ] ; then |
|---|
| 27 | ANA_DIR='/gpfsstore/rech/psl/rpsl376/ergon/ERA5/NETCDF/GLOBAL_025/4xdaily' |
|---|
| 28 | suf="e5.GLOBAL_025" |
|---|
| 29 | elif [ "$rea" = "ERAI" ] ; then |
|---|
| 30 | #ANA_DIR='http://prodn.idris.fr/thredds/dodsC/ipsl_private/rpsl376/ERAI/NETCDF/GLOBAL_1125/4xdaily' |
|---|
| 31 | ANA_DIR='/gpfsstore/rech/psl/rpsl376/ergon/ERAI/NETCDF/GLOBAL_1125/4xdaily' |
|---|
| 32 | #/gpfsstore/rech/psl/rpsl376/ergon/ERAI/NETCDF/GLOBAL_075/4xdaily/AN_PL/2006/u.200601.aphei.GLOBAL_075.nc |
|---|
| 33 | suf="ei.GLOBAL_1125" |
|---|
| 34 | else #OPERA |
|---|
| 35 | #ANA_DIR='http://prodn.idris.fr/thredds/dodsC/ipsl_private/rpsl376/OPERA/NETCDF/4xdaily' |
|---|
| 36 | ANA_DIR='/gpfsstore/rech/psl/rpsl376/ergon/OPERA/NETCDF/4xdaily' |
|---|
| 37 | suf=".GLOBAL_1125" |
|---|
| 38 | fi |
|---|
| 39 | |
|---|
| 40 | varu=u |
|---|
| 41 | varv=v |
|---|
| 42 | vart=ta # peut etre parfois juste "t" |
|---|
| 43 | varq=q |
|---|
| 44 | #varp=msl |
|---|
| 45 | |
|---|
| 46 | if [ "$rea" = "ERAI" ] ; then |
|---|
| 47 | # 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" |
|---|
| 52 | ncap2 -s 'u=float(u)' $fushort $fileu |
|---|
| 53 | ncap2 -s 'v=float(v)' $fvshort $filev |
|---|
| 54 | else |
|---|
| 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" |
|---|
| 57 | fi |
|---|
| 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" |
|---|
| 60 | |
|---|
| 61 | # verifier disponibilite des fichiers |
|---|
| 62 | if [ "$guide_uv" = "y" ] ; then ls $fileu ; ls $filev ; fi |
|---|
| 63 | if [ "$guide_t" = "y" ] ; then ls $filet ; fi |
|---|
| 64 | if [ "$guide_q" = "y" ] ; then ls $fileq ; fi |
|---|
