all-moc.sh

#!/bin/bash

echo "Start create AMOC of ALL"

data="/gpfsES/geo/zywang/CESM/Archive/MASK/ocn/hist/months"
path="/gpfsES/geo/the/MocArchieve/ALL"

# Create ALL path
if [ -d "$path" ]; then
	rm -r "$path"
	mkdir "$path"
else
	mkdir "$path"
fi
echo "Create ALL path complete!"

# Get AMOC from archive data of ALL
mkdir "$path/month"
for ((y=1; y<=2000; y++))
do
	ystr=$y
	if [ $y -lt 1000 ]; then
		if [ $y -lt 10 ]; then
			ystr="000"$ystr
		elif [ $y -lt 100 ]; then
			ystr="00"$ystr
		else
			ystr="0"$ystr
		fi
	fi
	for ((m=1; m<=12; m++))
	do
		mstr=$m
		if [ $m -lt 10 ]; then
			mstr="0"$m	
		fi
		ncks -v MOC $data/MASK.pop.h.$ystr-$mstr.nc $path/month/Moc.ALL.$ystr-$mstr.nc
	done
done
echo "Get AMOC data of ALL complete!"

# Create .nc to save monthly data and annual mean data of ALL
mkdir "$path/monthly"
mkdir "$path/annual"
echo "Start month2year_ALL.ncl"
ncl month2year_ALL.ncl
echo "Create .nc to save monthly data and annual mean data of ALL complete!"

# Draw maximum series and do wavelet transform of ALL
echo "Start drawseries_ALL.ncl"
ncl drawseries_ALL.ncl
echo "Draw maximum series and do wavelet transform of ALL complete!"