地面濃度模擬結果之整併與切割

module purge
module load rcec/tools pkg/Anaconda3

conda create -n py310 python==3.10
conda activate py310
conda install -c conda-forge wrf-python
conda install matplotlib cartopy

#sinotec2@lgn302 /work/sinotec2/cmaq_recommend/work/2023-06/grid03/cctm.ContEms/daily
#$ cat post.cs
#!/usr/bin/bash
module purge
module load rcec/tools pkg/Anaconda3
work=/work/sinotec2/cmaq_recommend/work
ym=2023-06
nd=4
GRD=( 'grid45'     'grid09'  'grid03' )
DOM=(grid45 grid09 Taiwan)
ITs=( 1 1 1 )
BEGD=2023-06-14
dates=();datep=()
for id in {0..11};do dates=( ${dates[@]} $(date -d "$BEGD +${id}days" +%Y-%m-%d) ); datep=( ${datep[@]} $(date -d "$BEGD +${id}days" +%Y%m%d) ); done

BIN=/work/sinotec2/opt/cmaq_recommend/bin


i=2
if [[ $1 == "comb" ]];then
  idb=0
  for ((id=$idb;id <= $nd;d+=1));do nc=$work/$ym/${GRD[$i]}/cctm.ContEms/daily/CCTM_ACONC_v532_intel_${DOM[$i]}_${datep[$id]}.nc;if [[ -e $nc ]];then sbatch $work/combine.sh $nc;fi;done
  exit 0
fi

it=${ITs[$i]}
dt=$(( 10#$it * 10000 ))
  for fn in $(ls PMs20*|tail -n${nd});do
    nt=$(( $($BIN/pr_tflag.py $fn|wc -l) - 1 ))
    j=$(echo $fn|cut -c4-11)
    for s in PM25_TOT PM10;do
      p=$s;test $p == "PM25_TOT" && p="PM2.5"
      ncks -O -d VAR,0 -d TSTEP,0,${nt},$it -v TFLAG,$s $fn ${p}_$j.nc
      ncatted -a TSTEP,global,o,i,$dt ${p}_$j.nc;done;done
  for s in SO2 CO O3 NO2;do for fn in $(ls CCTM_ACONC*|tail -n${nd});do
    j=$(echo $fn|cut -d'_' -f6)
    ncks -O -d LAY,0 -d VAR,0 -d TSTEP,0,${nt},$it -v TFLAG,$s $fn ${s}_$j
    ncatted -a TSTEP,global,o,i,$dt ${s}_$j;done;done
  for s in PM2.5 PM10 SO2 CO O3 NO2;do
    ncrcat -O ${s}_2*.nc ${s}.nc;rm ${s}_2*.nc;done
  for s in PM2.5 PM10 SO2 CO O3 NO2;do
    $BIN/m3nc2gif.py $s.nc;done
exit 0

#$ cat $work/combine.sh
#!/bin/sh
#SBATCH -A *********             # Account name/project number
#SBATCH -J comb                  # Job name
#SBATCH -p ct56                  # Partiotion name
#SBATCH -n 10                    # Number of MPI tasks (i.e. processes)
#SBATCH -c 1                     # Number of cores per MPI task
#SBATCH -N 1                     # Maximum number of nodes to be allocated
#SBATCH --ntasks-per-node=10    # Maximum number of tasks on each node
#SBATCH -o rsl.out.%j            # Path to the standard output file
#SBATCH -e rsl.error.%j          # Path to the standard error ouput file
module load compiler/intel/2021 IntelMPI/2021 hdf5/1.12 netcdf/4.7.4 pnetcdf/1.12.2

nc=$1
ii=$(echo $nc|cut -d'/' -f7|cut -c5-6)
ymd=$(echo $nc|cut -d'_' -f7|cut -d'.' -f1)
ym=$(echo $nc|cut -d'/' -f6)

export BASE=/work/sinotec2/cmaq_recommend/work
export BLD=/work/sinotec2/opt/cmaq_recommend/POST/combine/scripts/BLD_combine_v532_intel
export EXEC=$BLD/combine_v532.exe
export m3input=${BASE}/$ym/grid$ii/mcip
export cctmout=${BASE}/$ym/grid$ii/cctm.ContEms/daily

# user define
#> File [1]: CMAQ conc/aconc file
#> File [2]: MCIP METCRO3D file
#> File [3]: CMAQ APMDIAG file
#> File [4]: MCIP METCRO2D file
export INFILE2=${m3input}/METCRO3D_Taiwan.nc
export INFILE4=${m3input}/METCRO2D_Taiwan.nc

# programs
export LC_ALL=C
export LANG=C
export GENSPEC=N
export SPECIES_DEF=${BASE}/SpecDef_cb6r3_ae7_aq.txt
export INFILE1=$nc
export INFILE3=${nc/ACONC/APMDIAG}
export OUTFILE=${cctmout}/out${ymd}.conc.nc
if [ -e ${OUTFILE} ]; then rm ${OUTFILE};fi

time mpirun -bootstrap slurm -n 10 ${EXEC} >& ${BASE}/cmb.out
if [ -e ${cctmout}/PMs$ymd.nc ];then rm ${cctmout}/PMs$ymd.nc;fi
${BASE}/shk.cs $OUTFILE ${cctmout}/PMs$ymd.nc
if [ -e ${OUTFILE} ]; then rm ${OUTFILE};fi