source: trunk/src/post_plots.sh @ 257

#!/bin/bash

. env.sh


wavemin=1
wavemax=`ls Remaining_space_after_wave_* | cut -d _ -f5 | cut -d. -f1 | sort -n | tail -1`
# wavemax can also be set manually
#wavemax=22
nwave=$(( $wavemax - $wavemin + 1))

##################################################################
# Controle of a sublist for the waves shading on profiles
if [ $wavemax -le 3 ] ; then
  ListWaves=" `seq 1 $wavemax` "
else
  ListWaves=" 1 3 $wavemax "
fi
# could be also set manually
# ListWaves=" 1 3 6 "
# wave_inc=1
# ListWaves=$( seq $wavemin $wave_inc $wavemax )
##################################################################

if [ `python -V   2>&1 | awk ' { print $2 } ' | cut -c1` != 2 ] ; then
        echo You should use python 2 and not 3 ; exit ; fi


###############################################################################
# Some preliminary computations
#------------------------------
#echo $wavemin $wavemax $(seq $wavemin $wavemax) ; exit
for i in $(seq $wavemin $wavemax) ; do
        echo $i
        # Computing scores fore all the SCM simulations if not already done
        if [ ! -f score$i.csv ] ; then echo Runing post_scores.sh first ; pause 10 ; ./post_scores.sh $wavemin $wavemax ; fi
        # Compute the mean and max of the ensembles
        for d in WAVE$i/*/* ; do
            if [ -d $d ] ; then
                cd $d
                for stat in min max avg ; do if [ ! -f ens${stat}.nc ] ; then cdo ens${stat} SCM*nc ens${stat}.nc ; fi ; done
                cd -
           fi
        done
done

##############################################################################
# Ranking the SCM simulations according to scores
#------------------------------------------------

scores=""
for i in $( seq $wavemin $wavemax ) ; do scores="$scores score$i.csv" ; done
cat $scores | awk -F, ' { print  $(NF-1)","$0 } ' | sort >| opt_scores.csv
simss=`head -10 opt_scores.csv | awk -F, ' { print $2 } '`
echo $simss
head -1 WAVE$wavemax/Par*asc >| params_opt.csv
for sim in $simss ; do
        grep $sim */Par*asc >> params_opt.csv
done

iname=0

#############################################################################
# Loop on the best 10 simulations according to the above ranking
#---------------------------------------------------------------

for sims in "$simss" $simss  ; do
echo $sims

set -vx
cat <<eod>| simus.py
#-*- coding:UTF-8 -*-
# files that define the lists of the simulations you want to draw and also define information of the profiles
# time of the profile, color/style of the lines, xmin/xmax, ymin/ymax ranges)
from datetime import datetime

DIRDATA='./'

nomvar='theta'

# couleurs des courbes
#---------------------
# cycle de couleurs
listcoul=['black','red','green','fuchsia','blue','lime','darkviolet','cyan','darkorange','slateblue','brown','gold']

# exemple dictionnaire de couleurs en fonction du nom de fichier (sans le prefixe time_ ou prof_)
#black=reference; blue=bb; violet=turb, orange=pas de temps, Diffusion=vert, resolution=grey, advection=red, resv=brown,domain= microphysique=turquoise
dicocoul={\\
'SCM.nc' : 'black',\\
'LES0.nc' : 'slateblue',\\
'LES1.nc' : 'blue',\\
'LES2.nc' : 'blue',\\
'LES3.nc' : 'blue',\\
'LES4.nc' : 'blue',\\
'LES5.nc' : 'blue',\\
'LES6.nc' : 'blue',\\
'LES7.nc' : 'blue',\\
'LES8.nc' : 'blue',\\
eod

cols=( red green fuchsia lime darkviolet cyan darkorange slateblue brown gold )
ii=0
for sim in $sims ; do
cat <<eod>> simus.py
'$sim.nc' : '${cols[$ii]}',\\
eod
(( ii = $ii + 1 ))
done

cat <<eod>> simus.py
}

dicostyl={\\
'SCM.nc' : '-',\\
'LES0.nc' : '-',\\
'LES1.nc' : '-.',\\
'LES2.nc' : '-.',\\
'LES3.nc' : '-.',\\
'LES4.nc' : '-.',\\
'LES5.nc' : '-.',\\
'LES6.nc' : '-.',\\
'LES7.nc' : '-.',\\
'LES8.nc' : '-.',\\
eod

for sim in $sims ; do
cat <<eod>> simus.py
'$sim.nc' : '-',\\
eod
done

wave_ens=1
cat <<eod>> simus.py
}

# listcoul ou dicocoul ou dicohightune ou ...
#coul1d=listcoul
styl1d=dicostyl
coul1d=dicocoul

# liste des fichiers sélectionnés
listfic=[\\
'LES/'+Case+'/'+SubCase+'/LES0.nc',\\
'LES/'+Case+'/'+SubCase+'/LES1.nc',\\
'LES/'+Case+'/'+SubCase+'/LES2.nc',\\
'LES/'+Case+'/'+SubCase+'/LES3.nc',\\
'LES/'+Case+'/'+SubCase+'/LES4.nc',\\
'LES/'+Case+'/'+SubCase+'/LES6.nc',\\
'LES/'+Case+'/'+SubCase+'/LES7.nc',\\
'LES/'+Case+'/'+SubCase+'/LES8.nc',\\
eod



for sim in $sims ; do
w=`echo $sim | cut -d\- -f2`
echo $sim $w
cat <<eod>> simus.py
'WAVE$w/'+Case+'/'+SubCase+'/$sim.nc',\\
eod
done

cat <<eod>> simus.py
'CTRL/'+Case+'/'+SubCase+'/SCM.nc',\\
]
eod

echo 'list_ens=[\'  >> simus.py
for wave_ens in $ListWaves ; do
cat <<eod>> simus.py
'WAVE$wave_ens/'+Case+'/'+SubCase,\\
eod
done
echo ']'  >> simus.py

./post.sh
mv t.pdf Prof$iname.pdf

mkdir -p PROFILES/BEST$iname
mv profil*  PROFILES/BEST$iname


(( iname = $iname + 1 ))

done