CDO

Hi Guys,

i need to calculate the standard deviation of bandpass filtered sea level pressure data (2.5 - 6 days and 10 -90 days) for winter and summer. Each file contains one year of 6 hourly data.
I would like to make sure that i am doing it correctly. Maybe there is easy way to reach these results.

Thank you very much!

here my code:

#!/bin/bash

i=1
fr1=73 #365/6 days
fr2=146 #365/2.5 days

for d in {1..2}
do
for f in slp.{1979..2008}.nc
do
echo "Processing $f $d"
((i++))
cdo selmon,1,2,12 -bandpass,$fr1,$fr2 del29feb $f rm1.$i.nc
cdo selmon,6,7,8 -bandpass,$fr1,$fr2 $f rm2.$i.nc
done
cdo cat rm1.* djf.nc
cdo cat rm2.* jja.nc

cdo timstd djf.nc std_djf.$d.nc
cdo timstd jja.nc std_jja.$d.nc

rm rm1.* rm2.* djf.nc jja.nc
let i=1
let fr1=4 #365/90 days
let fr2=36.5 #365/10 days
done