#!/bin/bash

## Just run code in recipes, use it carefully
## example: use cdo to make cache file. Since cdo is quicker than ncl.
#;DIAG_NORCPM; RUNTHESECODES: echo '    no codes run here.'

#!/bin/bash
obsfn="HadISST_with_model_grid_ano.nc"
outfn="rhisto_NPac.txt"
lonlatbox="120,200,10,30"
modfns=$(ls model_sst_ano_mem*.nc)
operator="-ensrkhisttime"

cmd="cdo -s -outputkey,value ${operator}"
cmd="${cmd} -sellonlatbox,${lonlatbox} ${obsfn}"
for i in $modfns; do
    cmd="${cmd} -sellonlatbox,${lonlatbox} ${i}"
done
${cmd} > ${outfn}