ppth=/data/public/HYCOM/shian/hindcast/GLBv0.08 ncks=/cluster/gcc-4.8.5/tools/bin/ncks rm -rf record.txt for yy in {2003..2004..01} do for mm in {01..12..01} do mkdir -p ${ppth}/${yy}/sst/${mm} mkdir -p ${ppth}/${yy}/sss/${mm} mkdir -p ${ppth}/${yy}/ssuv/${mm} mkdir -p ${ppth}/${yy}/ssh/${mm} mkdir -p ${ppth}/${yy}/others/${mm} #get file year=${yy} path_ssh=${ppth}/${yy}/ssh/${mm} path_sst=${ppth}/${yy}/sst/${mm} path_sss=${ppth}/${yy}/sss/${mm} path_ssuv=${ppth}/${yy}/ssuv/${mm} path_others=${ppth}/${yy}/others/${mm} vv=536 case ${yy}${mm} in 200301|200302|200303|200304|200305|200306) vv=532 ;; *) esac case ${yy}${mm} in 200307|200308|200309|200310|200311|200312) vv=533 ;; *) esac case ${yy} in 2004|2005) vv=533 ;; *) esac case ${yy} in 2006) vv=534 ;; *) esac case ${yy}${mm} in 200507|200508|200509|200510|200511|200512|200701|200702|200703|200704|200705|200706) vv=534 ;; *) esac case ${yy}${mm} in 200707|200708|200709|200710|200711|200712) vv=535 ;; *) esac case ${yy} in 2008) vv=535 ;; *) esac case ${yy}${mm} in 200901|200902|200903|200904|200905|200906) vv=535 ;; *) esac case ${yy}${mm} in 201107|201108|201109|201110|201111|201112) vv=537 ;; *) esac case ${yy} in 2012) vv=537 ;; *) esac case ${yy}${mm} in 201301|201302|201303|201304|201305|201306) vv=537 ;; *) esac case ${yy}${mm} in 201307|201308|201309|201310|201311|201312|201401|201402|201403|201404|201405|201406) vv=538 ;; *) esac case ${yy}${mm} in 201407|201408|201409|201410|201411|201412|201601|201602|201603|201604) vv=563 ;; *) esac case ${yy} in 2015) vv=563 ;; *) esac case ${yy} in 2017) vv=928 ;; *) esac case ${yy}${mm} in 201605|201606|201607|201608|201609|201610|201611|201612|201701) vv=572 ;; *) esac case ${yy}${mm} in 201706|201707|201708|201709) vv=577 ;; *) esac case ${yy}${mm} in 201710|201711|201712) vv=929 ;; *) esac case ${yy} in 2018|2019|2020) vv=930 ;; *) esac bb=GLBv0.08 nc0=.nc nc1=.nc nc2=.nc case ${yy} in 2017|2018|2019|2020) bb=glbv nc0=_ssh.nc nc1=_ts3z.nc nc2=_uv3z.nc ;; *) esac case ${yy}${mm} in 201701|201706|201707|201708|201709) bb=GLBv0.08 nc0=.nc nc1=.nc nc2=.nc ;; *) esac rm -rf ${path}/${mm}/record.txt for j in {01..31..01} do for k in {00..18..06} do echo "do ${yy}${mm}${j}12_t0${k}" ######## #get ssh file=hycom_ssh_${yy}${mm}${j}12_t0${k}.nc panosize=`ls -l ${path_ssh}/${file} |cut -d ' ' -f 5` if [[ ! -f ${path_ssh}/${file} ]] || [[ ${panosize} -le 29000000 ]];then ${ncks} -v surf_el ${ppth}/${yy}/${mm}/hycom_${bb}_${vv}_${yy}${mm}${j}12_t0${k}${nc0} ${path_ssh}/hycom_ssh_${yy}${mm}${j}12_t0${k}.nc fi if [ -f ${path_ssh}/${file} ]; then echo "ssh_${yy}${mm}${j}_t0${k} ok" >> record.txt else echo "ssh_${yy}${mm}${j}_t0${k} false" >> record.txt fi ######## #get sst file=hycom_sst_${yy}${mm}${j}12_t0${k}.nc panosize=`ls -l ${path_sst}/${file} |cut -d ' ' -f 5` if [[ ! -f ${path_sst}/${file} ]] || [[ ${panosize} -le 29000000 ]];then ${ncks} -v water_temp -d depth,0 ${ppth}/${yy}/${mm}/hycom_${bb}_${vv}_${yy}${mm}${j}12_t0${k}${nc1} ${path_sst}/hycom_sst_${yy}${mm}${j}12_t0${k}.nc fi if [ -f ${path_sst}/${file} ]; then echo "sst_${yy}${mm}${j}_t0${k} ok" >> record.txt else echo "sst_${yy}${mm}${j}_t0${k} false" >> record.txt fi ######## #get sss file=hycom_sss_${yy}${mm}${j}12_t0${k}.nc panosize=`ls -l ${path_sss}/${file} |cut -d ' ' -f 5` if [[ ! -f ${path_sss}/${file} ]] || [[ ${panosize} -le 29000000 ]];then ${ncks} -v salinity -d depth,0 ${ppth}/${yy}/${mm}/hycom_${bb}_${vv}_${yy}${mm}${j}12_t0${k}${nc1} ${path_sss}/hycom_sss_${yy}${mm}${j}12_t0${k}.nc fi if [ -f ${path_sss}/${file} ]; then echo "sss_${yy}${mm}${j}_t0${k} ok" >> record.txt else echo "sss_${yy}${mm}${j}_t0${k} false" >> record.txt fi ######## #get ssu & ssv file=hycom_ssv_${yy}${mm}${j}12_t0${k}.nc panosize=`ls -l ${path_ssuv}/${file} |cut -d ' ' -f 5` if [[ ! -f ${path_ssuv}/${file} ]] || [[ ${panosize} -le 29000000 ]];then ${ncks} -v water_u -d depth,0 ${ppth}/${yy}/${mm}/hycom_${bb}_${vv}_${yy}${mm}${j}12_t0${k}${nc2} ${path_ssuv}/hycom_ssu_${yy}${mm}${j}12_t0${k}.nc ${ncks} -v water_v -d depth,0 ${ppth}/${yy}/${mm}/hycom_${bb}_${vv}_${yy}${mm}${j}12_t0${k}${nc2} ${path_ssuv}/hycom_ssv_${yy}${mm}${j}12_t0${k}.nc fi if [ -f ${path_ssuv}/${file} ]; then echo "ssuv_${yy}${mm}${j}_t0${k} ok" >> record.txt else echo "ssuv_${yy}${mm}${j}_t0${k} false" >> record.txt fi ####### #get M20 M15 MLD file=hycom_MLD_${yy}${mm}${j}12_t0${k}.nc panosize=`ls -l ${path_others}/${file} |cut -d ' ' -f 5` if [[ ! -f ${path_others}/${file} ]] || [[ ${panosize} -le 29000000 ]];then export hfile=${ppth}/${yy}/${mm}/hycom_${bb}_${vv}_${yy}${mm}${j}12_t0${k}${nc1} export ohfile1=${path_others}/hycom_M15_${yy}${mm}${j}12_t0${k}.nc export ohfile2=${path_others}/hycom_M20_${yy}${mm}${j}12_t0${k}.nc export ohfile3=${path_others}/hycom_MLD_${yy}${mm}${j}12_t0${k}.nc echo "${hfile}" python get2.py fi if [ -f ${path_others}/${file} ]; then echo "M15M20MLD_${yy}${mm}${j}_t0${k} ok" >> record.txt else echo "M15M20MLD_${yy}${mm}${j}_t0${k} false" >> record.txt fi echo "done ${yy}${mm}${j}12_t0${k}" done done done done