Voir aussi, comme indiqué dans le `courriel de florent.lyard@legos.obs-mip.fr du 2020-01-20 à 13:24 `_, ``_, téléchargé `ici <8713583_PGS_RM1-E2-2_construction%20calage_vf.pdf>`_. Environment =========== Voir `<~/work/NEA/vars.sh>`_. Marégraphes =========== Origine ------- - Pour commencer : + soit, à partir de ``_, cliquer sur : * `data `_ * `Access to Shom catalog `_ * `Coastal observations `_ + soit aller directement à ``_ - Cliquer sur `Sea level (REFMAR)` - Sélectionner `RONIM tide-gauges` et `Partners tide-gauges` - Pour chaque marégraphe : * cliquer sur le point ou le nom * cliquer sur `Download`, ce qui emmène au formulaire de téléchargement * si il reste des marégraphes à ajouter, cliquer sur `Back to the tide gauges` - Dans le formulaire de téléchargement : * Start date : `01/01/1940` * Type : `TXT` * sélectionner `Validated (delayed mode)` * entrer ses noms et courriels * cliquer sur `Generate archive` (fait à ~ 2020-01-09 11:20) `Courriel du 2020-01-09 12:15 `_ : Une archive restera disponible durant 30 jours au téléchargement à l'adresse suivante : https://services.data.shom.fr/o89z49ygr6c87m2cn967ilcd/telechargement/maregraphie/MAREGRAPHIE_157856510958003472/file/maregraphie.zip D'où :: c $NEArd wget ... c $NEAgd unzip $NEArd/maregraphie.zip list="`ls *_????.txt|sed -re 's/_.*/_/'|sort -u`" for s in $list;do cat $s????.txt>$s.txt done # for REFMAR.mgr, took ~480s detidor --no-date -format REFMAR *_.txt mv -v `ls -t *.mgr |head -n1` REFMAR.mgr # for *.gnu and stations_list, took 32s detidor --no-date -l stations_list -s 2018 -f 2018 -format REFMAR *_.txt sed -i~ -re "s/.*( XYN)/`tail -n+3 stations_list|wc -l`\1/" stations_list rm *_.mgr pocvip -f . ~/work/grd/europe-v04.cst \ settings/ -v arrows -s 0 -v "margins and palette" -s \$noPalette \ . gauges.pocvip $NEAgd/REFMAR.mgr -s "wM2 #000 4 -5" \ -n gauges.png Pour editer le fichier de marégraphes :: pocvip -f . ~/work/grd/europe-v04.cst \ $NEAsd/dump-tides-clamped/M2.sequential.nc \ -v a_eta_LGP1 -s ":p .5#888" -v G_eta_LGP1 -s ,p \ . gauges.pocvip $NEAgd/REFMAR-.mgr -s "wM2 #000 4" sed -re "s/ *STATION.*: *//;t;d" REFMAR.mgr|sort sed -re "; /ST-GILLES-CROIX-DE-VIE/,/LOCALISATION/ s/:.*E /: 46.68722N -1.94685E/; /LES_SABLES_D_OLONNE/,/LOCALISATION/ s/:.*E /: 46.484739N -1.791747E/; /LA_ROCHELLE_LA_PALLICE/,/LOCALISATION/ s/:.*E /: 46.15186N -1.23862E/; /LA_COTINIERE/,/LOCALISATION/ s/:.*E /: 45.909417N -1.328419E/; /BOURCEFRANC_LE_CHAPUS/,/^ *$/ d; /ROYAN/,/LOCALISATION/ s/:.*E /: 45.613599N -1.02609E/; /LAMENA/,/LOCALISATION/ s/:.*E /: 45.33675N -0.79414E/; /PAUILLAC_TROMPELOUP/,/LOCALISATION/ s/:.*E /: 45.21853N -0.74355E/; /ARCACHON_EYRAC/,/LOCALISATION/ s/:.*E /: 44.665251N -1.164062E/; /MIMIZAN/,/^ *$/ d; /AMBES/,/LOCALISATION/ s/:.*E /: 45.040404N -0.603637E/; /BASSENS/,/LOCALISATION/ s/:.*E /: 44.9037157N -0.5371244E/; /BORDEAUX/,/LOCALISATION/ s/:.*E /: 44.859839N -0.552921E/; " REFMAR.mgr >REFMAR-.mgr Pour visualiser les séries supprimées :: for f in 155_.gnu 6144_.gnu;do autoPlot.py $f;done # 155_2007-07-13_to_27.png 6144_2018-08-01_to_10-16.png mogrify -trim -bordercolor "#FFF" -border 1 *.png Forçage NEA =========== Données ------- Dans ``nuwa:/gsa/lyaf/models/nea-2018/simulation-surges/`` : + ``ERA5-detided+FES2014/`` + ``dump-tides/*.sequential.nc`` équivalent de ``ERA5-detided+FES2014/`` clamped puis flather + ``dump-tides-clamped/`` clamped seulement, donc bien meileur comme le montre :: tides-validate -unstructured LGP1 -a dump-tides-clamped/WAVE.sequential.nc -p ./ -v a_eta_LGP1 G_eta_LGP1 -g ../gauges/Residus_tot_V1_TPJ1J2_mgr_corrLoadFES2014_NEA_no_baltic.mgr M2 Copié localement comme l'indique `<~/work/NEA/notes.rst#input>`_ Plot ---- Voir les script de zoom ``_ et ``_. + grille et bathy chacun avec un zoom sur l'entrée Gironde :: pocvip -f . ~/work/grd/europe-v04.cst \ $NEAsd/dump-tides/M2.sequential.nc -v bathymetry -s :p \ . NEA.pocvip -n NEA-bathymetry.png -p "" \ :/ -v p -s "#000 0 140" \ . Gironde.pocvip -n NEA-bathymetry-gironde.png -p "" \ \ 3/ -v bathymetry -s 1 \ 4/ -v bathymetry -s ":p .5#444" \ settings/ -v "margins and palette" -s \$noPalette -v palette -s "" \ . NEA.pocvip -n NEA-grid.png -p "" \ . Gironde.pocvip -n NEA-grid-gironde.png -p "" différence bathy NEA-FES2014 :: pocvip -x -f . ~/work/grd/europe-v04.cst \ $NEAsd/dump-tides/M2.sequential.nc -v bathymetry -s :nea \ ~/work/FES2014/M2.spectral.nc -v bathymetry -s :f14 \ settings/ -v palette -s \$contBlueRed \ :/ -s "p=200*(nea-f14)/(nea+f14)" -v p -s "-1#000 -100 -10 -3 0 3 10 100" \ drawings/ -v title -s "bathymetry relative difference NEA-FES2014 (%)" \ . NEA.pocvip -n NEA-FES2014.png -p "" -q + erreur sur M2 S2 M4 : - avec ../gauges/Residus_tot_V1_TPJ1J2_mgr_corrLoadFES2014_NEA_no_baltic.mgr et aux marégraphes sur région proche Gironde avec la base REFMAR avec les version flather et clamped :: for w in M2 S2 M4;do l= [ $w = M4 ] && l=" 0 0.15" for n in flather clamped;do [ $n = flather ] && e= [ $n = clamped ] && e=-clamped pocvip -x -f . ~/work/grd/europe-v04.cst \ $NEAsd/dump-tides$e/$w.sequential.nc \ -v a_eta_LGP1 -s :p -v G_eta_LGP1 -s ,p -l fi=-1 \ $NEAd/gauges/Residus_tot_V1_TPJ1J2_mgr_corrLoadFES2014_NEA_no_baltic.mgr -s w$w \ $NEAgd/REFMAR-.mgr -s "w$w #444" \ settings/ -v palette -s \$magentaFree \ :/ -v p -s "#000$l" \ drawings/ -v title -s "$n $w" \ . NEA.pocvip -n $n-$w.png -p 728x1063 \ . gauges.pocvip -n $n-$w-gauges.png -p 475x1026 \ . Gironde.pocvip -n $n-$w-gironde.png -p 728x808 -q done done - aux marégraphes sur région proche Gironde avec la base REFMAR avec les version flather :: tides-validate -g $NEAgd/REFMAR-.mgr -p $NEAsd/dump-tides-clamped -a WAVE.sequential.nc -v a_eta_LGP1 G_eta_LGP1 -unstructured LGP1 -d M2 S2 M4 -o NEA-REFMAR + erreur sur surcote detidées extraire les séries temporelles aux marégraphes :: for d in $NEAsd/dump-tides{-clamped,}/extraction-REFMAR;do c $d Nohup.sh - extract-nf-v1 -format nc -start 01/2018 -end 12/2018 -s 900. -i $NEAgd/stations_list -c ../analysis-ROOT.YYYY.MM.nc sed -i~ -re "`head -qn1 $NEAgd/*_.gnu |sed -re "s|# ([^ ]+) .* ./([^.]+).gnu|s/\1.gnu/\2.gnu/|"`" translation done bien regarder le spectre avec :: autoPlot.py sample.PORT-BLOC # applying limits cleans-up around M2 but messes up around S2 #limits="-s 2018-01-08_02:15 -f 2018-12-31_17:45" autoPlot.py -a $limits sample.PORT-BLOC autoPlot.py -a $limits PORT-BLOC.h mettre le résultat dans ``_ appliquer le detiding des séries temporelles et calculer l'erreur sur surcote detidées :: ( for d in $NEAsd/dump-tides{-clamped,}/extraction-REFMAR;do c $d; #Nohup.sh - sample --no-fft -d -s $NEArd/NEA.spectrum -i $NEAgd/stations_list -start 2018-01-02 rm statistics Nohup.sh - analysis-tugo --no-fft -d $NEAgd -e "" -i $NEAgd/stations_list -t translation done ) paste $NEAsd/dump-tides{-clamped,}/extraction-REFMAR/statistics |sed -re "s/ *([^ ]+) .* ([^ ]+)\t.* ([^ ]+)$/\1 \2 \3/" + en regardant PORT-BLOC du 2018-06-26 au 28 avec detiding des obs comme celui du modèle.** :: c $NEAgd detidor --no-date -spectrum $NEArd/NEA.spectrum -s 2018 -f 2018 -format REFMAR *_.txt mv -v `ls -t *.mgr |head -n1` REFMAR-2018-NEAspectrum.mgr Simulations =========== Avec 3 conditions de rivères : étiage, moyenne, crue et M2 S2. Données ------- dans :: shennong:/gsa/lyaf/models/gironde/simulation-tides+DAC+rivers/marest-clamped/ /mnt/srv-ecola/data/sealevel/gironde/SHOM-2018/processing/ atlas copiés localement comme l'indique `<~/work/gironde/notes.rst#simulation>`_ Plot ---- + erreur sur M2 S2 M4 aux marégraphes sur région proche Gironde avec la base REFMAR :: for w in M2 S2 M4;do l= #[ $w = M4 ] && l=" 0 0.15" pocvip -x -f . ~/work/grd/europe-v04.cst \ ~/work/gironde/simulation-tides+DAC+rivers/marest-clamped/$w.sequential.nc \ -v a_eta_LGP1 -s :p -v G_eta_LGP1 -s ,p -l fi=-1 \ $NEAgd/REFMAR-.mgr -s "w$w #444" \ settings/ -v palette -s \$magentaFree \ :/ -v p -s "#000$l" \ drawings/ -v title -s "$w" \ . NEA.pocvip -n gironde-$w.png -p 728x1063 -q done tides-validate -g $NEAgd/REFMAR-.mgr -p $NEAsd/dump-tides-clamped -a WAVE.sequential.nc -v a_eta_LGP1 G_eta_LGP1 -unstructured LGP1 M2 S2 M4 -o Gironde-REFMAR