主页 / user-336457

Maria Karypidou's questions

Martin Hope
Maria Karypidou
Asked: 2019-10-10 02:24:12 +0800 CST
  • 1

我想将分数（在第一个循环中计算分数）提高到小数幂（第二个循环），但是，我总是得到 1 作为结果。我也想将第二个循环的输出存储在一个数组中。任何想法如何解决这个问题？谢谢！

#
# vector of vertical pressure levels
levs=($(seq 200.0 50.0 900.0))
printf "%s\n" "${levs[@]}"
#
# exponent for dry air
rho=$(bc -l <<<'e(l(0.0819)*0.5)')
echo $rho
#
# calculate fraction of P_surf/P_i from Poisson equation for each vertical pressure level
val3=()
#
for i in "${levs[@]}"
 do
   echo $i
   val3+=($(bc -l <<<"1000.0/$i"))
   echo "$val3"
 done
printf "%s\n" "${val3[@]}"
#
# raise fraction of P_surf/P_i to the rho power for dry air (#bc <<< "2 ^ 3")
pow=()
#
for j in "${val3[@]}"
 do
   echo $j
   echo $rho
   pow+=($(bc <<<"$j^rho"))
   #echo $((i*rho))
   echo "$pow"
 done
#
arithmetic
Martin Hope
Maria Karypidou
Asked: 2019-02-13 04:05:22 +0800 CST
  • 0

我有几个文件夹，每个文件夹都包含几个.nc文件。我想遍历所有文件夹，并为每个文件夹使用cdo具有相似模式的文件进行一些计算（合并时间）。到目前为止，我已经写了以下内容：

for dir in /mnt/meteo_a/africa_cordex/historical/0.44/pr/*/     
do 
 dir=${dir%*/}
 echo ${dir##*/} 
 export dir
 cd $dir
 pwd
  for f in `find . -type f -regex /(.*?\day)/*`
  cdo mergetiem io
  done
done

每个文件夹中的文件具有如下名称：

pr_AFR-44_CNRM-CERFACS-CNRM-CM5_historical_r1i1p1_CLMcom-CCLM4-8-17_v1_day_19500101-19501231.nc
pr_AFR-44_CNRM-CERFACS-CNRM-CM5_historical_r1i1p1_CLMcom-CCLM4-8-17_v1_day_19510101-19551231.nc
pr_AFR-44_ICHEC-EC-EARTH_historical_r12i1p1_CLMcom-CCLM4-8-17_v1_day_19491201-19501231.nc
pr_AFR-44_ICHEC-EC-EARTH_historical_r12i1p1_CLMcom-CCLM4-8-17_v1_day_19510101-19551231.nc

我想对之前具有相似模式的所有文件应用 mergetime CDO 命令_day_。

我想合并前两个文件和后两个文件。更具体地说，合并以

"pr_AFR-44_CNRM-CERFACS-CNRM-CM5_historical_r1i1p1_CLMcom-CCLM4-8-17_v1_"

同样合并以

"pr_AFR-44_ICHEC-EC-EARTH_historical_r12i1p1_CLMcom-CCLM4-8-17_v1"

合并时间的输出可能类似于：

 {pattern}_mergetime.nc
bash files