#!/bin/bash domain='PERU' #domain='CALIFORNIA' #domain='SENEGAL' #domain='BENGUELA' EUC="./PERU_EUC_uo_seuil_20pct.txt" #EUC="../../PERU_EUC_uo_seuil_20pct.SAV_2019_09_05_OK_01_remake_list_hfds/"${domain}"_"${seuil_uo}".txt" #EUC="../../"${domain}"_"${seuil_uo}".txt.pulsation" #echo ${EUC} line1=3 # premiere ligne des donnees #line2=20 # derniere ligne des donnees pour hfdsO (restricted list) line2=23 # derniere ligne des donnees pour hfdsO (sans les MIROC) #line2=69 # derniere ligne des donnees pour uo seuil_uo="EUC_uo_seuil_20pct" #F_hfds="non" grep -q hfds ${EUC} ; if [ "$?" -eq "0" ]; then F_hfds="oui"; else F_hfds="non"; fi if [ "${F_hfds}" == "oui" ]; then F_Flux=`grep "Heat Flux" PERU_EUC_uo_seuil_20pct.txt | awk -v var=2 -F":" '{print $var}'` fi #------------------------------------------------------------------------------- # Sections PACIFIC / ATLANTIC #------------------------------------------------------------------------------- if [ ${domain} == 'PERU' ] || [ ${domain} == 'CALIFORNIA' ] ; then # PACIFIC #S EUC S EUC S EUC S EUC S EUC S EUC S EUC HFDS (W/m2) T lev1 (SST) T 54m T 87m T 120m T 29-114m T 181m id exp name #98.5W 85W 110.5W 120.25W130W 140.5W 170.5W upA upB upA upB upBas upA upB upBas upA upB upBas upA upB upBas upA upB upBas upA upB upBas #14.3 14.5 14.7 15.5 15.9 16.3 16.8 86.8 108.6 25.3 23.3 24.1 14.1 12.5 13.4 12.9 11.6 12.3 12.3 11.1 11.7 13.4 11.4 13.1 11.2 10.4 10.8 c512 cmip5_GFDL-ESM2G #15.0 14.6 15.0 15.3 15.5 15.7 15.7 31.0 62.0 24.6 20.3 19.5 16.2 15.5 15.2 14.0 14.0 13.9 13.5 13.2 13.1 12.9 7.1 14.9 13.2 12.4 12.2 c509 cmip5_CSIRO-Mk3-6-0 section_list=("98_5W" "85W" "110_5W" "120_25W" "130W" "140_5W" "170_5W") cat ${EUC} | awk -v var=1 -F" " '{print $var}' > S_EUC_98_5W.tmp cat ${EUC} | awk -v var=2 -F" " '{print $var}' > S_EUC_85W.tmp cat ${EUC} | awk -v var=3 -F" " '{print $var}' > S_EUC_110_5W.tmp cat ${EUC} | awk -v var=4 -F" " '{print $var}' > S_EUC_120_25W.tmp cat ${EUC} | awk -v var=5 -F" " '{print $var}' > S_EUC_130W.tmp cat ${EUC} | awk -v var=6 -F" " '{print $var}' > S_EUC_140_5W.tmp cat ${EUC} | awk -v var=7 -F" " '{print $var}' > S_EUC_170_5W.tmp ii=7 else # ATLANTIC #S EUC S EUC S EUC T lev1 (SST) T 54m T 87m T 120m T 29-114m T 181m #30.25W 20.5W 10W upA upB upBas upA upB upBas upA upB upBas upA upB upBas upA upB upBas upA upB upBas #21.5 21.4 21.3 19.8 25.6 25.6 18.1 19.4 19.4 15.8 16.4 16.5 14.0 14.7 15.0 17.2 18.6 18.6 12.1 12.9 13.1 75e1 tr075n_newflx #21.6 21.1 21.4 20.5 25.9 26.0 19.8 23.2 23.1 18.6 19.4 19.6 17.3 17.1 17.4 19.4 21.9 21.9 14.4 14.6 14.8 75e2 tr075n_n02nogm section_list=("30_25W" "20_5W" "10W") cat ${EUC} | awk -v var=1 -F" " '{print $var}' > S_EUC_30_25W.tmp cat ${EUC} | awk -v var=2 -F" " '{print $var}' > S_EUC_20_5W.tmp cat ${EUC} | awk -v var=3 -F" " '{print $var}' > S_EUC_10W.tmp ii=3 fi #------------------------------------------------------------------------------ if [ "${F_hfds}" = "oui" ] ; then ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > HFDS_upA.tmp ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > HFDS_upB.tmp ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > HFDS_upBas.tmp fi ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > T_upA_1m.tmp ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > T_upB_1m.tmp ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > T_upBas_1m.tmp ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > T_upA_54m.tmp ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > T_upB_54m.tmp ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > T_upBas_54m.tmp ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > T_upA_87m.tmp ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > T_upB_87m.tmp ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > T_upBas_87m.tmp ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > T_upA_120m.tmp ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > T_upB_120m.tmp ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > T_upBas_120m.tmp ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > T_upA_29_114m.tmp ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > T_upB_29_114m.tmp ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > T_upBas_29_114m.tmp ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > T_upA_181m.tmp ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > T_upB_181m.tmp ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > T_upBas_181m.tmp ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > id.tmp ii=`expr ${ii} + 1`; cat ${EUC} | awk -v var=${ii} -F" " '{print $var}' > exp.tmp #echo "# ${domain} (${seuil_uo})" for section in ${section_list[*]}; do echo "S_EUC_"${section}"_list=`sed -n -e ${line1},${line2}p S_EUC_${section}.tmp | tr '\n' ',' | sed -e 's/^/[/' | sed -e 's/,$/]/' ` ; " done if [ "${F_hfds}" = "oui" ] ; then for file_root in "HFDS_upA" "HFDS_upB" "HFDS_upBas"; do file="${file_root}" # echo ${file} echo "${file}_list=`sed -n -e ${line1},${line2}p ${file}.tmp | tr '\n' ',' | sed -e 's/^/[/' | sed -e 's/,$/]/' ` ; " done fi for file_root in "T_upA" "T_upB" "T_upBas"; do for depth in "1" "54" "87" "120" "29_114" "181"; do # for depth in "1" "54" "87" "120" "181"; do file="${file_root}_${depth}m" # echo ${file} echo "${file}_list=`sed -n -e ${line1},${line2}p ${file}.tmp | tr '\n' ',' | sed -e 's/^/[/' | sed -e 's/,$/]/' ` ; " done done echo "id_list=`sed -n -e ${line1},${line2}p id.tmp |sed -e "s/^/'/" |sed -e "s/$/'/" | tr '\n' ',' | sed -e 's/^/[/' | sed -e 's/,$/]/'` ; " echo "exp_list=`sed -n -e ${line1},${line2}p exp.tmp |sed -e "s/^/'/" |sed -e "s/$/'/" | tr '\n' ',' | sed -e 's/^/[/' | sed -e 's/,$/]/'` ; " grep 'cmip5_res' ${EUC} grep 'cmip5_2D_interp_method' ${EUC} if [ "${F_hfds}" = "oui" ] ; then echo "F_hfds=True" echo "F_Flux='${F_Flux}'" else echo "F_hfds=False" fi rm -f *.tmp