source: trunk/INPUT/ANT-LBq-15km/moyennes_par_fleuves-bis.sh @ 278

# utilise la carte de fleuves pour faire des moyennes

ntime=14
grdmath Netcdf-Resu/AntI2S21_class01_01.nc?Hdot[$ntime] 1 MUL = file1.grd
grdmath Netcdf-Resu/AntI2S18_class01_01.nc?Hdot[$ntime] 1 MUL = file2.grd
grdmath Netcdf-Resu/AntI2S17_class01_01.nc?Hdot[$ntime] 1 MUL = file3.grd
grdmath Netcdf-Resu/AntI2S14_class01_01.nc?Hdot[$ntime] 1 MUL = file4.grd
grdmath Netcdf-Resu/AntI2S12_class01_01.nc?Hdot[$ntime] 1 MUL = file5.grd
grdmath Netcdf-Resu/AntI2S13_class01_01.nc?Hdot[$ntime] 1 MUL = file6.grd
grdmath Netcdf-Resu/AntI2S15_class01_01.nc?Hdot[$ntime] 1 MUL = file7.grd
grdmath Netcdf-Resu/AntI2S16_class01_01.nc?Hdot[$ntime] 1 MUL = file8.grd
grdmath Netcdf-Resu/AntI2S19_class01_01.nc?Hdot[$ntime] 1 MUL = file9.grd
grdmath Netcdf-Resu/AntI2S20_class01_01.nc?Hdot[$ntime] 1 MUL = file10.grd

rm Resu_total

# tableau coef de départ 1 sur grounded, NaN ailleurs
cp masque_grounded_15km_Nan.grd coef_gen.grd

for ((i = 1; i <= 210 ; i += 1))
do
  echo $i
  echo $i > resultat_fleuve

# cree un fichier par fleuve 
# fleuves-2byrd.grd est le fichier qui contient les fleuves
# avec Byrd coupe en basin (4) et glacier (210)

  grdmath fleuves-2byrd.grd $i EQ = toto.grd     # trie un fleuve donne : numero i
  grdfilter toto.grd -D0 -Fu30000 -Gfleuve_$i.grd   # elargit le fleuve d'un rang

# nombre de points dans ce fleuve = MEAN * NX*NY
    grdmath fleuve_$i.grd MEAN = nb_fleuve_$i.grd
    
nbpoints=`grd2xyz -ZBL fleuve_$i.grd  |awk '{if ($1 == 1) sum=$1+sum; print sum}' |tail -1`
 echo $nbpoints >> resultat_fleuve

# Pour chaque fichier donnee
  for ((k = 1; k <= 10 ;k += 1))
    do
    name="file"$k.grd
    echo $name
    cp $name file.grd
    grdmath fleuve_$i.grd file.grd MUL = Hdot_fleuve_$i.grd


# moyenne du hdot sur ce fleuve
# pour faire une carte
#    grdmath Hdot_fleuve_$i.grd MEAN nb_fleuve_$i.grd DIV fleuve_$i.grd MUL = Hdot_mean_fleuve_$i.grd

    htot=`grd2xyz -ZBL Hdot_fleuve_$i.grd  |awk '{ sum=$1+sum; print sum}' |tail -1`

    echo $htot $nbpoints 

    echo  $htot $nbpoints | awk '{print $1/$2}' >> resultat_fleuve

    rm   Hdot_fleuve_$i.grd 
  done   # fin de la boucle sur les runs

# calcule le beta
  beta=`cherche_beta`

# donne le coeff de la region et l'additionne a coef_gen
  grdmath fleuve_$i.grd 0 NAN = masque_fleuve.grd
  grdmath fleuve_$i.grd $beta MUL masque_fleuve.grd MUL coef_gen.grd AND = coef_gen.grd
 

  cat resultat_fleuve  | tr '\n'  ' ' > titi
cp  resultat_fleuve  recapitulatif_fleuves/resultat_fleuve_$i
  cat resultat_fleuve  | tr '\n'  ' ' > titi
  awk '{print $0, '\n'}' titi >>  Resu_total

# rangements


mv fleuve_$i.grd   Travail_sur_fleuves
rm toto.grd nb_fleuve_$i.grd
rm  resultat_fleuve



done  # fin de la boucle fleuve


# moyenne du delta hdot sur ce fleuve (5-6)
#grdmath  Delta_Hdot_fleuve_$i.grd  MEAN nb_fleuve_$i.grd DIV fleuve_$i.grd MUL =  Delta_Hdot_mean_fleuve_$i.grd

# rapport entre les deux
#grdmath Hdot_mean_fleuve_$i.grd Delta_Hdot_mean_fleuve_$i.grd DIV = rapport_fleuve_$i.grd


gzip -f Travail_sur_fleuves/*grd