ppth=/data/public/HYCOM/shian/hindcast/GLBy0.08
ncks=/cluster/gcc-4.8.5/tools/bin/ncks
rm -rf record.txt
for yy in {2023..2023..01}
do
  for mm in {11..11..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=930
      bb=glby
      nc0=_ssh.nc
      nc1=_ts3z.nc
      nc2=_uv3z.nc

      rm -rf ${path}/${mm}/record.txt


    for j in {16..30..01}
    do
      for k in {00..00..00}
      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