#!/bin/sh

dir=$STORE1/CMIP6/day

mlist="ACCESS-CM2_r4i1p1f1.hssp585 ACCESS-ESM1-5_r34i1p1f1.hssp585 CanESM5_r6i1p2f1.hssp585 EC-Earth3-Veg_r4i1p1f1.hssp585 EC-Earth3-Veg-LR_r3i1p1f1.hssp585 HadGEM3-GC31-LL_r4i1p1f3.hssp585 KACE-1-0-G_r3i1p1f1.hssp585 MPI-ESM1-2-LR_r3i1p1f1.hssp585 MRI-ESM2-0_r5i1p1f1.hssp585 UKESM1-0-LL_r2i1p1f2.hssp585"
mlist="EC-Earth3_r4i1p1f1.hssp585 FGOALS-g3_r4i1p1f1.hssp585"
mlist="ACCESS-CM2_r4i1p1f1.hssp585 ACCESS-ESM1-5_r34i1p1f1.hssp585"
mlist="CanESM5_r6i1p2f1.hssp585 EC-Earth3-Veg_r4i1p1f1.hssp585 EC-Earth3-Veg-LR_r3i1p1f1.hssp585 HadGEM3-GC31-LL_r4i1p1f3.hssp585 KACE-1-0-G_r3i1p1f1.hssp585 MPI-ESM1-2-LR_r3i1p1f1.hssp585 MRI-ESM2-0_r5i1p1f1.hssp585 UKESM1-0-LL_r2i1p1f2.hssp585"
mlist="EC-Earth3_r4i1p1f1.hssp585 FGOALS-g3_r4i1p1f1.hssp585"
mlist="ERA5"

# mlist="ACCESS-CM2_r1i1p1f1.hssp585 ACCESS-ESM1-5_r2i1p1f1.hssp585 CAMS-CSM1-0_r2i1p1f1.hssp585 CanESM5_r1i1p1f1.hssp585 CMCC-ESM2_r1i1p1f1.hssp585 CNRM-CM6-1-HR_r1i1p1f2.hssp585 CNRM-CM6-1_r1i1p1f2.hssp585 CNRM-ESM2-1_r1i1p1f2.hssp585 EC-Earth3_r1i1p1f1.hssp585 EC-Earth3-CC_r1i1p1f1.hssp585 EC-Earth3-Veg_r1i1p1f1.hssp585 EC-Earth3-Veg-LR_r1i1p1f1.hssp585 FGOALS-g3_r1i1p1f1.hssp585 GFDL-CM4_r1i1p1f1.hssp585 GISS-E2-1-G_r1i1p1f2.hssp585 HadGEM3-GC31-LL_r1i1p1f3.hssp585 HadGEM3-GC31-MM_r1i1p1f3.hssp585 INM-CM4-8_r1i1p1f1.hssp585 INM-CM5-0_r1i1p1f1.hssp585 IPSL-CM6A-LR_r1i1p1f1.hssp585 KACE-1-0-G_r1i1p1f1.hssp585 MIROC-ES2L_r1i1p1f2.hssp585 MIROC6_r1i1p1f1.hssp585 MPI-ESM1-2-HR_r1i1p1f1.hssp585 MPI-ESM1-2-LR_r1i1p1f1.hssp585 MRI-ESM2-0_r1i1p1f1.hssp585 NESM3_r1i1p1f1.hssp585 NorESM2-LM_r1i1p1f1.hssp585 NorESM2-MM_r1i1p1f1.hssp585 TaiESM1_r1i1p1f1.hssp585 UKESM1-0-LL_r1i1p1f2.hssp585"

# m=ERA5
# fera=/ccc/store/cont003/gen6877/p25rv/GRIDOBS/ERA5/era5_phi500.fixed.nc
# cdo selseas,3 $fera tmp/phi500.$m.jja.nc
# cdo remapbil,GRID_ANALOGS tmp/phi500.$m.jja.nc tmp/phi500.$m.dom.jja.nc
# cdo timmean tmp/phi500.$m.dom.jja.nc tmp/phi500.$m.dom.jja.mean.nc
# cdo sub tmp/phi500.$m.dom.jja.nc tmp/phi500.$m.dom.jja.mean.nc tmp/phi500_dom.ano.$m.nc

# cdo runmean,30 tmp/phi500.$m.dom.nc tmp/phi500.$m.dom.runmean.nc
# cdo ydaymean tmp/phi500.$m.dom.runmean.nc tmp/phi500.$m.dom.scyc.nc
# cdo sub tmp/phi500.$m.dom.nc tmp/phi500.$m.dom.scyc.nc tmp/phi500_dom.anos.$m.allyear.nc
# cdo selseas,3 tmp/phi500_dom.anos.$m.allyear.nc tmp/phi500_dom.anos.$m.nc

# cdo runmean,30 /ccc/store/cont003/gen6877/p25rv/GRIDOBS/ERA5/tasmax.era5.eur.nc tmp/tasmax.$m.runmean.nc
# cdo ydaymean tmp/tasmax.$m.runmean.nc tmp/tasmax.$m.scyc.nc
# cdo sub /ccc/store/cont003/gen6877/p25rv/GRIDOBS/ERA5/tasmax.era5.eur.nc tmp/tasmax.$m.scyc.nc tmp/tasmax.$m.anos.nc
# cdo selseas,3 tmp/tasmax.$m.anos.nc models/tasmax.$m.eur.jja.anos.nc

for m in $mlist ; do
  echo $m
# cdo selyear,1950/2022 $dir/tasmax/tasmax.$m.nc tmp.nc
# cdo subc,273.15 tmp.nc tmp1.nc
# cdo selseas,3 tmp1.nc tmp2.nc
# cdo remapbil,GRID_TAS tmp2.nc models/tasmax.$m.eur.jja.nc
  cdo selseas,3 $dir/stream500new/stream500.$m.nc tmp/phi500.$m.jja.nc
  cdo remapbil,GRID_ANALOGS tmp/phi500.$m.jja.nc tmp/phi500.$m.dom.jja.nc
  cdo timmean tmp/phi500.$m.dom.jja.nc tmp/phi500.$m.dom.jja.mean.nc
  cdo sub tmp/phi500.$m.dom.jja.nc tmp/phi500.$m.dom.jja.mean.nc tmp/phi500_dom.ano.$m.nc
# gawk '$1>=1950&&$1<=2022{print}' covariates/SGSAT.$m.dat > covariates/SGSAT1950-2022.$m.dat
done