#!/bin/sh set -x . $HOME/PlotsOfTheDay/LONLAT.def nowdtg=`date -u +"%Y%m%d%H"` nowdate=`$HOME/bin/mandtg -date $nowdtg` nowhour=`$HOME/bin/mandtg -hour $nowdtg` if [ $nowhour -lt 2 ] then basedtg=`$HOME/bin/mandtg $nowdtg + -48` else basedtg=`$HOME/bin/mandtg $nowdtg + -24` fi basedtg=`$HOME/bin/mandtg -date $basedtg`18 baseyyyy=`$HOME/bin/mandtg -year $basedtg` basemm=`$HOME/bin/mandtg -month $basedtg` basedd=`$HOME/bin/mandtg -day $basedtg` basehh=`$HOME/bin/mandtg -hour $basedtg` [ -r /data/hirlam/hl_arc/OPR/$baseyyyy/$basemm/$basedd/$basehh/fc$baseyyyy$basemm${basedd}_${basehh}+0${FL}ve ] || exit # Run not finished yet. ls -1t $HOME/public_html/$basedtg/* $HOME/public_html/${nowdate}??/* | head -1 | grep index && exit # Nothing new to process. echo $0 - DTG of last full run: $basedtg wd=/data/hirlam/WD$$ mkdir $wd || exit cd $wd || exit html=$HOME/public_html todir=$html/$basedtg mkdir $todir # Get the list of regions for which to produce graphs. grep -v '^#' > region.list <Verwachting vanaf $basedtg UTC" >> tmp.html echo "" >> tmp.html echo "" \ "" \ " " \ " " \ " " >> tmp.html echo "
Landgebruikgegevens aan het begin van de verwachting
PNGPSHoe wordt deze kaart gemaaktHoe wordt landgebruik bepaald
" >> tmp.html # Construct the table of animations of the forecast. echo "

Animaties

" >> tmp.html echo "" >> tmp.html echo "" \ "" \ "" \ "" \ "" >> tmp.html echo "" \ "" \ "" \ "" >> tmp.html echo "" \ "" \ "" \ "" >> tmp.html echo "" \ "" \ "" \ "" \ "" >> tmp.html echo "" \ "" \ "" \ "" \ "" >> tmp.html echo "" \ "" \ "" \ "" \ "" >> tmp.html echo "" \ "" \ "" \ "" \ "" >> tmp.html echo "
Lokaal
'Neerslagradar'BewolkingZichtWind
TemperatuurLatente WarmteVoelbare Warmte
Klassiek
Luchtdruk en windStraalstroom en hoogte 500hPa vlak
Bewegings- en potentiële energie
WindsnelheidWindstotenConvective Available Potential Energy
Waterhuishouding
WolkenNeerslag/SoortZichtSneeuwhoogte
Opwarming en afkoeling
TemperatuurLatente warmtestroomVoelbare warmtestroom
'Satelliet' beelden
WolkentoptemperatuurWaterdampWaterdamp met wolkencorrectieWolkenwater fractie
" >> tmp.html # Construct the table of "weercijfers" for the areas predicted by this forecast. echo "

Weercijfers

" >> tmp.html echo "" >> tmp.html echo "" \ "" \ "" \ "" \ "" \ "" \ "" \ >> tmp.html if [ $FL -eq 78 ] then echo "" >> tmp.html fi echo "" >> tmp.html for area in $regionlist do lat=`grep ^$area region.list | awk '{print $2}'` lon=`grep ^$area region.list | awk '{print $3}'` size=`grep ^$area region.list | awk '{print $4}'` tz=`grep ^$area region.list | awk '{print $5}'` echo "" \ "" \ "" \ "" \ "" \ "" \ "" \ >> tmp.html if [ $FL -eq 78 ] then echo "" >> tmp.html fi echo "" >> tmp.html done echo "
Locatie/WaterschapBreedtegraadLengtegraadTijdzoneRoostercellenVandaagMorgenOvermorgen
$area$lat$lon$tz$size x $size`cat $area.WEERCIJFER | awk -F, '{print $1}'``cat $area.WEERCIJFER | awk -F, '{print $2}'``cat $area.WEERCIJFER | awk -F, '{print $3}'`
" >> tmp.html # Construct the table of recent logs of computations and maps of observations used. echo "

Logs van de berekeningen en overzicht van de waarnemingen

" >> tmp.html echo "" >> tmp.html now=`date -u +"%Y%m%d"` for HH in `seq -w 0 6 18` do dtgobs=$now$HH [ -r $html/$dtgobs/observations$HH.png ] || dtgobs=`$HOME/bin/mandtg $dtgobs + -24` DATEobs=`$HOME/bin/mandtg -date $dtgobs` dtglog=$now$HH [ -r $html/$dtglog/HL_Cycle_$dtglog.html ] || dtglog=`$HOME/bin/mandtg $dtglog + -24` DATElog=`$HOME/bin/mandtg -date $dtglog` logfile=${DATElog}${HH}/HL_Cycle_${DATElog}${HH}.html echo "" >> tmp.html echo "" >> tmp.html echo "" >> tmp.html for h in -3 -2 -1 +0 +1 +2 do dtgoff=`$HOME/bin/mandtg $dtgobs + $h` houroff=`$HOME/bin/mandtg -hour $dtgoff` echo "" >> tmp.html done if [ $HH -eq 0 ] then echo "" >> tmp.html else echo "" >> tmp.html fi echo "" >> tmp.html echo "" >> tmp.html echo "" >> tmp.html for h in -3 -2 -1 +0 +1 +2 do dtgoff=`$HOME/bin/mandtg $dtglog + $h` houroff=`$HOME/bin/mandtg -hour $dtgoff` if [ -s $html/$dtglog/use_ascat$houroff.png ] then echo "" >> tmp.html else echo "" >> tmp.html fi done used=`awk '{print $1}' $html/$dtglog/ASCAT_COUNT.txt`; rejected=`awk '{print $2}' $html/$dtglog/ASCAT_COUNT.txt` echo "" >> tmp.html echo "" >> tmp.html echo "" >> tmp.html echo "" >> tmp.html for h in -3 -2 -1 +0 +1 +2 do dtgoff=`$HOME/bin/mandtg $dtglog + $h` houroff=`$HOME/bin/mandtg -hour $dtgoff` if [ -s $html/$dtglog/use_airep$houroff.png ] then echo "" >> tmp.html else echo "" >> tmp.html fi done used=`awk '{print $1}' $html/$dtglog/AIREP_COUNT.txt`; redundant=`awk '{print $2}' $html/$dtglog/AIREP_COUNT.txt`; rejected=`awk '{print $3}' $html/$dtglog/AIREP_COUNT.txt` echo "" >> tmp.html echo "" >> tmp.html echo "" >> tmp.html echo "" >> tmp.html for h in -3 -2 -1 +0 +1 +2 do dtgoff=`$HOME/bin/mandtg $dtglog + $h` houroff=`$HOME/bin/mandtg -hour $dtgoff` if [ -s $html/$dtglog/use_temp$houroff.png ] then echo "" >> tmp.html else echo "" >> tmp.html fi done usedt=`awk '{print $1}' $html/$dtglog/TEMP_COUNT.txt`; rejectedt=`awk '{print $2}' $html/$dtglog/TEMP_COUNT.txt` usedp=`awk '{print $1}' $html/$dtglog/PILOT_COUNT.txt`; rejectedp=`awk '{print $2}' $html/$dtglog/PILOT_COUNT.txt` used=`expr $usedt \+ $usedp`; rejected=`expr $rejectedt \+ $rejectedp` echo "" >> tmp.html echo "" >> tmp.html echo "" >> tmp.html echo "" >> tmp.html for h in -3 -2 -1 +0 +1 +2 do dtgoff=`$HOME/bin/mandtg $dtglog + $h` houroff=`$HOME/bin/mandtg -hour $dtgoff` if [ -s $html/$dtglog/use_other$houroff.png ] then echo "" >> tmp.html else echo "" >> tmp.html fi done used1=`awk '{print $1}' $html/$dtglog/SYNOP_COUNT.txt`; redundant1=`awk '{print $2}' $html/$dtglog/SYNOP_COUNT.txt`; rejected1=`awk '{print $3}' $html/$dtglog/SYNOP_COUNT.txt` used2=`awk '{print $1}' $html/$dtglog/SHIP_COUNT.txt`; redundant2=`awk '{print $2}' $html/$dtglog/SHIP_COUNT.txt`; rejected2=`awk '{print $3}' $html/$dtglog/SHIP_COUNT.txt` used3=`awk '{print $1}' $html/$dtglog/DRIBU_COUNT.txt`; redundant3=`awk '{print $2}' $html/$dtglog/DRIBU_COUNT.txt`; rejected3=`awk '{print $3}' $html/$dtglog/DRIBU_COUNT.txt` used=`expr $used1 \+ $used2 \+ $used3`; redundant=`expr $redundant1 \+ $redundant2 \+ $redundant3`; rejected=`expr $rejected1 \+ $rejected2 \+ $rejected3` echo "" >> tmp.html echo "" >> tmp.html echo "" >> tmp.html echo "" >> tmp.html for h in -3 -2 -1 +0 +1 +2 do dtgoff=`$HOME/bin/mandtg $dtgobs + $h` houroff=`$HOME/bin/mandtg -hour $dtgoff` if [ -s $html/$dtglog/use_span$houroff.png ] then echo "" >> tmp.html else echo "" >> tmp.html fi done used=`awk '{print $1}' $html/$dtglog/SPAN_COUNT.txt`; rejected=`awk '{print $2}' $html/$dtglog/SPAN_COUNT.txt` echo "" >> tmp.html echo "" >> tmp.html done echo "
$HH UTC
`export TZ=UTC; ls -l --time-style=+"%Y%m%d%H%M%S" /data/hirlam/hl_arc/OPR/ob$dtgobs | awk '{print $6}'` UTCIngekomen waarnemingen${h}H PNG " >> tmp.html echo "PSGebruiktRedundantVerworpen
`export TZ=UTC; ls -l --time-style=+"%Y%m%d%H%M%S" $html/$logfile | awk '{print $6}'` LOGVerwerking ASCAT${h}H PNG " >> tmp.html echo "PS$used$rejected
Effect waarnemingenVerwerking AIREP/AMDAR${h}H PNG " >> tmp.html echo "PS$used$redundant$rejected
DruktendensVerwerking TEMP/PILOT${h}H PNG " >> tmp.html echo "PS$used$rejected
Maximum windVerwerking SYNOP/SHIP/DRIBU${h}H PNG " >> tmp.html echo "PS$used$redundant$rejected
Oppervlakte PNG PSVerwerking waarnemingen${h}H PNG " >> tmp.html echo "PS$used$rejected
" >> tmp.html echo "

" \ "De ingekomen waarnemingen zijn gekleurd naar hoogte volgens het cubehelix algoritme. " \ "De verwerkte waarnemingen zijn gekleurd volgens de behandeling (gebruikt, redundant of verworpen)." \ "

" >> tmp.html # Show the evolution of the distance of the model to the observations at the start of data assimilation echo "

Ontwikkeling afstand tot de waarnemingen over de afgelopen 30 dagen.

" >> tmp.html echo "" >> tmp.html echo "" >> tmp.html echo "" >> tmp.html echo "" >> tmp.html echo "" >> tmp.html echo "" >> tmp.html echo "
Afstand tot de waarnemingen aan het begin van de minimalisatie
Beginwaarde afstandHoe dit te bepalen
" >> tmp.html echo "

" \ "Deze afstand geeft aan in hoeverre het model in 6 uur tijd verwijderd is geraakt van de werkelijkheid," \ "en is dus een maat voor de (on)voorspelbaarheid van de atmosfeer in die tijdspanne." \ "

" >> tmp.html # Construct the table of verification data over the past 30 days. echo "

Verificatie over de afgelopen 30 dagen

" >> tmp.html . $HOME/hl_home/OPR/progressPP.log verifdtg=`$HOME/bin/mandtg $DTGPP + -6` echo "" >> tmp.html echo "" >> tmp.html for par in pmsl w10m t2m td2m rh2m cc rr06 rr12 rr24 do echo "" >> tmp.html for fp in `seq -w 6 6 $FL` do file=$verifdtg/verif_time-${par}-${fp}.png if [ -s $html/$file ] then echo "" >> tmp.html else echo "" >> tmp.html fi done echo "" >> tmp.html done echo "" >> tmp.html echo "" >> tmp.html for par in pmsl w10m t2m td2m rh2m cc rr06 rr12 rr24 do echo "" >> tmp.html done echo "" >> tmp.html echo "" >> tmp.html for level in 250 500 700 850 925 1000 do echo "" >> tmp.html echo "" >> tmp.html for par in wind temp dewp geop relh do echo "" >> tmp.html done echo "" >> tmp.html done echo "
Oppervlak (tijdreeks) - Nederland
${par^^}+$fp
Oppervlak (geaggregeerd) - Nederland
${par^^}
Bovenlucht (geaggregeerd) - Europa
$level hPa${par^^}
" >> tmp.html # Loop over the region list, plotting time series of the near surface atmospheric quantities. echo "

Tijdreeksen

" >> tmp.html echo "" >> tmp.html echo "" >> tmp.html echo "" >> tmp.html echo "" >> tmp.html echo "" >> tmp.html echo "" >> tmp.html XTIX="('d0h18' 0" for h in `seq 6 6 $FL` do d=`expr \( $h + 18 \) / 24` hh=`expr \( $h - 6 \) % 24` XTIX=$XTIX",'d${d}h${hh}' $h" done XTIX=$XTIX")" for area in $regionlist do tz=`grep ^$area region.list | awk '{print $5}'` [ -s $todir/$area.WEERBERICHT ] || $HOME/PlotsOfTheDay/Weerbericht $area $basedtg $tz && mv $area.WEERBERICHT $todir echo "" >> tmp.html # Temperature and dew point temperature. plottempobs=""; [ -r $area.TEMPERATUUR-obs ] && plottempobs=",'$area.TEMPERATUUR-obs' using 1:(\$2-273.16):3 title 'T obs' pointtype 5" plotdauwobs=""; [ -r $area.DAUWPUNT-obs ] && plotdauwobs=",'$area.DAUWPUNT-obs' using 1:(\$2-273.16):3 title 'TD obs' pointtype 7" gnuplot << EOF set terminal pngcairo mono size 1440,720 set output '$area-t2m.png' set style data errorbars set key outside set grid linestyle 2 set xtics $XTIX set xrange [-2:`expr $FL + 2`] set xlabel 'Verwachtingstijd (UTC)' set ylabel 'Temperatuur en Dauwpunt (graden C)' set title 'Verwachting $basedtg: $area - 2 meter temperatuur en dauwpunt' plot '$area.TEMPERATUUR' using 1:(\$2-273.16):3 title 'T gem +/- stdev' linetype -1 pointtype 6, \ '$area.TEMPERATUUR' using 1:(\$4-273.16) title 'T min' with lines linetype 0, \ '$area.TEMPERATUUR' using 1:(\$5-273.16) title 'T max' with lines linetype 0, \ '$area.DAUWPUNT' using 1:(\$2-273.16):3 title 'TD gem +/- stdev' linetype -1 pointtype 3, \ '$area.DAUWPUNT' using 1:(\$4-273.16) title 'TD min' with lines linetype 0, \ '$area.DAUWPUNT' using 1:(\$5-273.16) title 'TD max' with lines linetype 0 \ $plottempobs \ $plotdauwobs EOF mv -f $area-t2m.png $todir mv -f $area.TEMPERATUUR $todir mv -f $area.DAUWPUNT $todir echo "" >> tmp.html echo " " >> tmp.html echo " " >> tmp.html # Relative humidity. plotrhobs=""; [ -r $area.RELVOCHT-obs ] && plotrhobs=",'$area.RELVOCHT-obs' using 1:2:3 title 'RH obs' pointtype 5" gnuplot << EOF set terminal pngcairo mono size 1440,720 set output '$area-rh.png' set style data errorbars set key outside set grid linestyle 2 set xtics $XTIX set xrange [-2:`expr $FL + 2`] set xlabel 'Verwachtingstijd (UTC)' set ylabel 'Relatieve vochtigheid (%)' set title 'Verwachting $basedtg: $area - relatieve vochtigheid op 2 meter' plot '$area.RELVOCHT' using 1:(\$2*100):(\$3*100) title 'RH gem +/- stdev' linetype -1 pointtype 6, \ '$area.RELVOCHT' using 1:(\$4*100) title 'RH min' with lines linetype 0, \ '$area.RELVOCHT' using 1:(\$5*100) title 'RH max' with lines linetype 0 \ $plotrhobs EOF mv -f $area-rh.png $todir mv -f $area.RELVOCHT $todir echo "" >> tmp.html echo " " >> tmp.html # Wind direction. plotwdirobs=""; [ -r $area.WINDRICHTING-obs ] && plotwdirobs=",'$area.WINDRICHTING-obs' using 1:2:3 title 'DD obs' pointtype 5" gnuplot << EOF set terminal pngcairo mono size 1440,720 set output '$area-dd10m.png' set style data errorbars set key outside set grid linestyle 2 set xtics $XTIX set xrange [-2:`expr $FL + 2`] set xlabel 'Verwachtingstijd (UTC)' set ytics ('Noord' 0, 'Noordoost' 45, 'Oost' 90, 'Zuidoost' 135, 'Zuid' 180, 'Zuidwest' 225, 'West' 270, 'Noordwest' 315, 'Noord' 360) set ylabel 'Windrichting' set title 'Verwachting $basedtg: $area - 10 meter wind richting' plot '$area.WINDRICHTING' using 1:2:3 title 'DD gem +/- stdev' linetype -1 pointtype 6, \ '$area.WINDRICHTING' using 1:4 title 'DD min' with lines linetype 0, \ '$area.WINDRICHTING' using 1:5 title 'DD max' with lines linetype 0 \ $plotwdirobs EOF mv -f $area-dd10m.png $todir mv -f $area.WINDRICHTING $todir echo "" >> tmp.html echo " " >> tmp.html # Wind speed. plotwspdobs=""; [ -r $area.WINDSNELHEID-obs ] && plotwspdobs=",'$area.WINDSNELHEID-obs' using 1:2:3 title 'FF obs' pointtype 5" plotgustobs=""; [ -r $area.GUST-obs ] && plotgustobs=",'$area.GUST-obs' using 1:2:3 title 'GUST obs' pointtype 7" gnuplot << EOF set terminal pngcairo mono size 1440,720 set output '$area-ff10m.png' set style data errorbars set key outside set grid linestyle 2 set xtics $XTIX set xrange [-2:`expr $FL + 2`] set xlabel 'Verwachtingstijd (UTC)' set yrange [0:] set ylabel 'Wind(stoten) (m/s)' set title 'Verwachting $basedtg: $area - 10 meter wind snelheid en windstoten' plot '$area.WINDSNELHEID' using 1:2:3 title 'FF gem +/- stdev' linetype -1 pointtype 6, \ '$area.WINDSNELHEID' using 1:4 title 'FF min' with lines linetype 0, \ '$area.WINDSNELHEID' using 1:5 title 'FF max' with lines linetype 0, \ '$area.WINDSTOTEN' using 1:2:3 title 'GUST gem +/- stdev' linetype -1 pointtype 3, \ '$area.WINDSTOTEN' using 1:4 title 'GUST min' with lines linetype 0, \ '$area.WINDSTOTEN' using 1:5 title 'GUST max' with lines linetype 0 \ $plotwspdobs \ $plotgustobs EOF mv -f $area-ff10m.png $todir mv -f $area.WINDSNELHEID $todir mv -f $area.WINDSTOTEN $todir echo "" >> tmp.html echo " " >> tmp.html echo " " >> tmp.html # Neerslag. # '$area.CONVNEERSLAG' using 1:2:3 axes x1y1 title 'RRc gem +/- stdev' linetype -1 pointtype 3, # '$area.CONVNEERSLAG' using 1:4 axes x1y1 title 'RRc min' with lines linetype 0, # '$area.CONVNEERSLAG' using 1:5 axes x1y1 title 'RRc max' with lines linetype 0, awk -f $HOME/PlotsOfTheDay/neso.awk $area.NEERSLAGSRT > $area.neso awk 'BEGIN {prevmin=0; prevmax=0}; {if ($1 > 0) {min=$4-prevmin; max=$5-prevmax; if (min > max) max=min; prevmin=$4; prevmax=$5; print $1, min, max, 0}}' $area.FRNTNEERSLAG > $area.RR1 plotprecobs=""; [ -r $area.NEERSLAG-obs ] && plotprecobs=",'$area.NEERSLAG-obs' using 1:2:3 axes x1y1 title 'RR1 obs' pointtype 5 pointsize 1" gnuplot << EOF set terminal pngcairo mono size 1440,720 set output '$area-rr.png' set style data errorbars set key outside set grid linestyle 2 set xtics $XTIX set xrange [-2:`expr $FL + 2`] set xlabel 'Verwachtingstijd (UTC)' set ytics nomirror autofreq set y2range [0:1.05] set y2tics 0,0.1,1 set ylabel 'Neerslag (kg/m^2)' set y2label ' Neerslagsoort (%)' set title 'Verwachting $basedtg: $area - neerslag en neerslagsoort' plot '$area.FRNTNEERSLAG' using 1:2:3 axes x1y1 title 'RRf gem +/- stdev' linetype -1 pointtype 6, \ '$area.FRNTNEERSLAG' using 1:4 axes x1y1 title 'RRf min' with lines linetype 0, \ '$area.FRNTNEERSLAG' using 1:5 axes x1y1 title 'RRf max' with lines linetype 0, \ '$area.RR1' using 1:2:4 axes x1y1 title 'RR1 min' with points pointtype 8 pointsize 2, \ '$area.RR1' using 1:3:4 axes x1y1 title 'RR1 max' with points pointtype 12 pointsize 2, \ '$area.neso' using 1 axes x1y2 title 'Regen' with points pointtype 7 pointsize 2, \ '$area.neso' using 2 axes x1y2 title 'Sneeuw' with points pointtype 9 pointsize 2, \ '$area.neso' using 3 axes x1y2 title 'Sneeuw' with points pointtype 9 pointsize 2, \ '$area.neso' using 4 axes x1y2 title 'Ijzel' with points pointtype 4 pointsize 2, \ '$area.neso' using 5 axes x1y2 title 'Ijsregen' with points pointtype 11 pointsize 2, \ '$area.neso' using 6 axes x1y2 title 'Ijsregen/natte sneeuw' with points pointtype 13 pointsize 2, \ '$area.neso' using 7 axes x1y2 title 'Natte sneeuw' with points pointtype 15 pointsize 2 \ $plotprecobs EOF mv -f $area-rr.png $todir # mv -f $area.FRNTNEERSLAG $todir # Later ... mv -f $area.CONVNEERSLAG $todir echo "" >> tmp.html echo " " >> tmp.html # echo " " >> tmp.html # Cloud cover. plotccobs=""; [ -r $area.BEWOLKING-obs ] && plotccobs=",'$area.BEWOLKING-obs' using 1:(\$2*0.01):3 title 'CC obs' pointtype 5" gnuplot << EOF set terminal pngcairo mono size 1440,720 set output '$area-cc.png' set style data errorbars set key outside set grid linestyle 2 set xtics $XTIX set xrange [-2:`expr $FL + 2`] set xlabel 'Verwachtingstijd (UTC)' set ytics ('0' 0, '1' 0.125, '2' 0.25, '3' 0.375, '4' 0.5, '5' 0.625, '6' 0.75, '7' 0.875, '8' 1.0) set yrange [-0.8:1.1] set y2tics autofreq set y2range [:600] set ylabel 'Bewolkingsgraad (octa)' set y2label 'Latente en Voelbare warmtestroom (W/m^2)' set title 'Verwachting $basedtg: $area - bewolking; latente en voelbare warmtestroom' plot '$area.BEWOLKING' using 1:2:3 axes x1y1 title 'CC gem +/- stdev' linetype -1 pointtype 6, \ '$area.BEWOLKING' using 1:4 axes x1y1 title 'CC min' with lines linetype 0, \ '$area.BEWOLKING' using 1:5 axes x1y1 title 'CC max' with lines linetype 0, \ '$area.LATHFLUX' using 1:2:3 axes x1y2 title 'LATHFLX gem +/- stdev' linetype -1 pointtype 3, \ '$area.LATHFLUX' using 1:4 axes x1y2 title 'LATHFLX min' with lines linetype 0, \ '$area.LATHFLUX' using 1:5 axes x1y2 title 'LATHFLX max' with lines linetype 0, \ '$area.SENSHFLUX' using 1:2:3 axes x1y2 title 'SNSHFLX gem +/- stdev' linetype -1 pointtype 4, \ '$area.SENSHFLUX' using 1:4 axes x1y2 title 'SNSHFLX min' with lines linetype 0, \ '$area.SENSHFLUX' using 1:5 axes x1y2 title 'SNSHFLX max' with lines linetype 0 \ $plotccobs EOF mv -f $area-cc.png $todir mv -f $area.BEWOLKING $todir mv -f $area.LATHFLUX $todir mv -f $area.SENSHFLUX $todir echo "" >> tmp.html echo " " >> tmp.html echo " " >> tmp.html echo " " >> tmp.html # Visibility. plotvisobs=""; [ -r $area.ZICHT-obs ] && plotvisobs=",'$area.ZICHT-obs' using 1:2:3 title 'VIS obs' pointtype 5" gnuplot << EOF set terminal pngcairo mono size 1440,720 set output '$area-vis.png' set style data errorbars set key outside set grid linestyle 2 set xtics $XTIX set xrange [-2:`expr $FL + 2`] set xlabel 'Verwachtingstijd (UTC)' set logscale y set yrange [10:50000] set ylabel 'Zicht (m)' set title 'Verwachting $basedtg: $area - zicht' plot '$area.ZICHT' using 1:2:3 title 'VIS gem +/- stdev' linetype -1 pointtype 6, \ '$area.ZICHT' using 1:4 title 'VIS min' with lines linetype 0, \ '$area.ZICHT' using 1:5 title 'VIS max' with lines linetype 0 \ $plotvisobs EOF mv -f $area-vis.png $todir mv -f $area.ZICHT $todir echo "" >> tmp.html echo " " >> tmp.html # Snow heigth. plotsnhobs=""; [ -r $area.SNOWDEPTH-obs ] && plotsnhobs=",'$area.SNOWDEPTH-obs' using 1:2:3 title 'SNH obs' pointtype 5" gnuplot << EOF set terminal pngcairo mono size 1440,720 set output '$area-sn.png' set style data errorbars set key outside set grid linestyle 2 set xtics $XTIX set xrange [-2:`expr $FL + 2`] set xlabel 'Verwachtingstijd (UTC)' set ylabel 'Sneeuwhoogte (m)' set title 'Verwachting $basedtg: $area - sneeuwhoogte' plot '$area.SNEEUWHOOGTE' using 1:2:3 title 'SNH gem +/- stdev' linetype -1 pointtype 6, \ '$area.SNEEUWHOOGTE' using 1:4 title 'SNH min' with lines linetype 0, \ '$area.SNEEUWHOOGTE' using 1:5 title 'SNH max' with lines linetype 0 \ $plotsnhobs EOF mv -f $area-sn.png $todir mv -f $area.SNEEUWHOOGTE $todir echo "" >> tmp.html echo " " >> tmp.html # Convective Available Potential Energy. gnuplot << EOF set terminal pngcairo mono size 1440,720 set output '$area-cape.png' set style data errorbars set key outside set grid linestyle 2 set xtics $XTIX set xrange [-2:`expr $FL + 2`] set xlabel 'Verwachtingstijd (UTC)' set ytics nomirror autofreq set ylabel 'CAPE en CIN (J/kg)' set y2tics autofreq set y2label 'Convectieve neerslag (kg/m^2)' set title 'Verwachting $basedtg: $area - Convective Available Potential Energy' plot '$area.CAPE' using 1:2:3 axes x1y1 title 'CAPE gem +/- stdev' linetype -1 pointtype 6, \ '$area.CAPE' using 1:4 axes x1y1 title 'CAPE min' with lines linetype 0, \ '$area.CAPE' using 1:5 axes x1y1 title 'CAPE max' with lines linetype 0, \ '$area.CIN' using 1:2:3 axes x1y1 title 'CIN gem +/- stdev' linetype -1 pointtype 4, \ '$area.CIN' using 1:4 axes x1y1 title 'CIN min' with lines linetype 0, \ '$area.CIN' using 1:5 axes x1y1 title 'CIN max' with lines linetype 0, \ '$area.FRNTNEERSLAG' using 1:2:3 axes x1y2 title 'RR gem +/- stdev' linetype -1 pointtype 3, \ '$area.FRNTNEERSLAG' using 1:4 axes x1y2 title 'RR min' with lines linetype 0, \ '$area.FRNTNEERSLAG' using 1:5 axes x1y2 title 'RR max' with lines linetype 0 EOF mv -f $area-cape.png $todir mv -f $area.CAPE $todir mv -f $area.CIN $todir mv -f $area.FRNTNEERSLAG $todir echo "" >> tmp.html echo " " >> tmp.html echo " " >> tmp.html # Incoming and outgoing radiation. gnuplot << EOF set terminal pngcairo mono size 1440,720 set output '$area-rad.png' set style data errorbars set key outside set grid linestyle 2 set xtics $XTIX set xrange [-2:`expr $FL + 2`] set xlabel 'Verwachtingstijd (UTC)' set ylabel 'Zonlicht (J/m^2)' set y2tics autofreq set y2label 'Infrarood (J/m^2)' set title 'Verwachting $basedtg: $area - Zonne-instraling en infrarood-uitstraling' plot '$area.ZONLICHT' using 1:2:3 axes x1y1 title 'ZL gem +/- stdev' linetype -1 pointtype 6, \ '$area.ZONLICHT' using 1:4 axes x1y1 title 'ZL min' with lines linetype 0, \ '$area.ZONLICHT' using 1:5 axes x1y1 title 'ZL max' with lines linetype 0, \ '$area.INFRAROOD' using 1:2:3 axes x1y2 title 'IR gem +/- stdev' linetype -1 pointtype 3, \ '$area.INFRAROOD' using 1:4 axes x1y2 title 'IR min' with lines linetype 0, \ '$area.INFRAROOD' using 1:5 axes x1y2 title 'IR max' with lines linetype 0 EOF mv -f $area-rad.png $todir mv -f $area.ZONLICHT $todir mv -f $area.INFRAROOD $todir echo "" >> tmp.html echo " " >> tmp.html echo " " >> tmp.html done echo "
Locatie/WaterschapGroothedenTabellen (verwachtingstijd gemiddelde stddev min max ...)
$area
Temperatuur en dauwpuntTemperatuurDauwpunt
Relatieve vochtigheidRelatieve vochtigheid
WindrichtingWindrichting
Windsnelheid en windstotenWindsnelheidWindstoten
Neerslag; neerslagsoortNeerslag
Convectieve neerslag
Bewolking; latente en voelbare warmtestroomBewolkingLatente warmtestroomVoelbare warmtestroom
ZichtZicht
SneeuwhoogteSneeuwhoogte
Convective Available Potential Energy; InhibitionConvective Available Potential EnergyConvective Inhibition
Zonlicht en infraroodZonlichtInfrarood
" >> tmp.html echo "" >> tmp.html echo "" >> tmp.html mv -f tmp.html $todir/index.html rm -f $html/index.html ln -s $todir/index.html $html/index.html cd rm -rf $wd exit 0