#!/bin/sh

#####
##### example script that uses CDO to compute the spatial correlation between two fields
#####

file1=CCTdiff_ARI-BASE_1991-2010.nc
file2=RSDSdiff_ARI-BASE_1991-2010.nc

fileland=LANDMASK.nc

fileout=scorr.txt

cdo mul $fileland $file1 kk1
cdo sellonlatbox,-9,53,33.5,66.3 kk1 kk11
cdo mul $fileland $file2 kk2
cdo sellonlatbox,-9,53,33.5,66.3 kk2 kk22
cdo fldcor kk11 kk22 kkc

corr=$(cdo outputf,%g kkc)

echo $corr >> $fileout

rm kk*