CAMx模擬結果之壓縮

$ diff shkavg4.4.f90 shkavgcb6.f90
9,10c9,10
<    integer,parameter::nlst=32,nvlst=36  
<    integer,parameter::nout=8  
---
>    integer,parameter::nlst=32,nvlst=62 
>    integer,parameter::nout=9  

diff shkavg4.4.f90 shkavgPar.f90
130,131c131,132
<    lsvname(17)  = 'ISPD'
<    lsvname(18)  = 'ISPX' !
---
>    lsvname(17)  = 'ISP'  !ISP = ISOP
>    lsvname(18)  = 'ISPD' !
145c146
<    lsvname(32)  = 'TO2' !TO2 =  Peroxy radical from OH addition to TOL 7 173.1
>    lsvname(32)  = 'TOLA' !TOLA = TOL or ARO1
149c150
<    lsvname(36)  = 'XLO2' !Peroxy radical from OH addition to XYL 8 187.1
---
>    lsvname(36)  = 'XYLA' !XYLA = XYL or ARO2

kuang@master ~/bin
$ cat shk_camx.cs
if [ $HOSTNAME == '114-32-164-198.HINET-IP.hinet.net' ];then NCO='/opt/anaconda3/bin'
elif [ $HOSTNAME == 'master' ];then NCO='/cluster/netcdf/bin'
elif [ $HOSTNAME == 'node01' ];then NCO='/cluster/netcdf/bin'
elif [ $HOSTNAME == 'centos8' ];then NCO='/opt/anaconda3/envs/py37/bin'
elif [ $HOSTNAME == 'node03' ];then NCO='/opt/miniconda3/bin'
else NCO='/usr/bin'
fi
NCKS='pncgen -f uamiv'
NCRCAT=${NCO}/ncrcat

VAR='TFLAG,CO,NO2,SO2,O3,PAR,PSO4,PNO3,PNH4,FPRM,FCRS'

$NCKS -O --slice LAY,0,0 -v $VAR $1 $2
#$NCKS -O -v $VAR -d LAY,0,0 $1 $2

kuang@master ~/bin
$ cat shk.cs
nc=$1
if [ $HOSTNAME == '114-32-164-198.HINET-IP.hinet.net' ];then NCO='/opt/anaconda3/bin'
elif [ $HOSTNAME == 'master' ];then NCO='/cluster/netcdf/bin'
elif [ $HOSTNAME == 'node01' ];then NCO='/cluster/netcdf/bin'
elif [ $HOSTNAME == 'node02' ];then NCO='/cluster/netcdf/bin'
elif [ $HOSTNAME == 'centos8' ];then NCO='/opt/anaconda3/envs/py37/bin'
elif [ $HOSTNAME == 'node03' ];then NCO='/opt/miniconda3/bin'
else NCO='/usr/bin'
fi
NCKS=${NCO}/ncks
NCRCAT=${NCO}/ncrcat
NCDUMP=${NCO}/ncdump
v=$($NCDUMP -h $nc|grep PM25|wc -l)
a=$($NCDUMP -h $nc|grep AVERAGE|wc -l)
if [ $v != 0 ];then
  if  [ $a == 0 ];then
    #cmaq combined files
    VAR='TFLAG,CO,NO2,SO2,O3,PM25_NO3,PM25_SO4,PM25_TOT,PM10,VOC'
  else
  #version 700 nc directly from CAMx
    VAR='TFLAG,CO,NO2,SO2,O3,PNO3,PSO4,PNH4,PAR,CCRS,FCRS,CPRM,FPRM'
  fi
#camx>400, nc generated by pncgen from uamiv file
else
  VAR='TFLAG,CO,NO2,SO2,O3,PNO3,PSO4,PNH4,PAR,CCRS,FCRS,CPRM,FPRM'
fi
if ! [ -e $2 ];then
  touch $2
  res=${VAR//[^,]}
  nvars=$(echo ${#res})
  $NCKS -O -v $VAR -d LAY,0,0 -d VAR,1,$nvars $1 $2
  ncatted -a NVARS,global,o,i,$nvars $2
fi

kuang@master /nas1/cmaqruns/2019base/Annual
$ cat ./do_shk
#note the ncks or ncrcat can not run in parallel, will conflick in memory
ls -r /nas1/cmaqruns/2016base/data/out*/POST/COMBINE_ACONC*_sCh*10.nc > fnames.txt
sort fnames.txt>a;mv a fnames.txt
for i in $(cat fnames.txt);do ymd=$(echo $i|cut -d_ -f11);shk.cs $i $ymd.nc;done &

kuang@master ~/bin
$ cat shk_Days_DM.cs
#$1:first mmdd
#$2:last  mmdd
#$3:domain d01~d04

if [ $HOSTNAME == '114-32-164-198.HINET-IP.hinet.net' ];then NCO='/opt/anaconda3/bin'
elif [ $HOSTNAME == 'master' ];then NCO='/cluster/netcdf/bin'
elif [ $HOSTNAME == 'centos8' ];then NCO='/opt/anaconda3/envs/py37/bin'
elif [ $HOSTNAME == 'node03' ];then NCO='/opt/miniconda3/bin'
else NCO='/usr/bin'
fi
NCKS=${NCO}/ncks
NCRCAT=${NCO}/ncrcat

VAR='TFLAG,CO,NO2,SO2,O3,PM25_NO3,PM25_SO4,PM25_TOT,PM10,VOC'

CASE=10
ROOT=$3
test $ROOT == 'd04' && GRID='TWN_3X3'
test $ROOT == 'd01' && GRID='EAsia_81K'
test $ROOT == 'd02' && GRID='sChina_27k'

yymm=$(echo $PWD|cut -d'_' -f5|cut -d'/' -f1)
bdate=`echo $(ls -rt COMBINE_ACONC*${GRID}_${CASE}.nc|head -n1)|cut -d'_' -f7`
bjul=$(date -d "$bdate" +%Y%j)
fjul=$(date -d "2016$1" +%Y%j)
ljul=$(date -d "2016$2" +%Y%j)
bi=$(( $fjul - $bjul + 1))
ei=$(( $ljul - $bjul + 1))
for ((i=$bi;i<=$ei;i+=1));do
  j=$(( 10#$i - 1 ))
  r=$( echo $j/4+5|bc )
  cdate=$(date -d "$bdate +${j} day" +%Y%m%d)
  file=COMBINE_ACONC_v53_gcc_${yymm}_run${r}_${cdate}_${GRID}_${CASE}.nc
  echo $file
  $NCKS -O -v $VAR -d LAY,0,0 $file ${cdate}_$ROOT.nc
done
$NCRCAT -O 20*_$ROOT.nc $ROOT.nc