if [ $# -ne 1 ] ; then
  echo "Usage $0 <input.nc>"
fi
inp=$1
out=${inp%.nc}_dephy.nc

ncks $inp -A $out \
  -v MEANWIND,MEAN_PRE,MEAN_RC,MEAN_RT,MEAN_RV,.MEAN_RR,.MEAN_RI,.MEAN_RS,.MEAN_RG
ncks $inp -A $out \
  -v MEAN_RHO,MEAN_REHU,MEAN_CF,MEAN_TH,MEAN_THV,MEAN_THL,MEAN_U,MEAN_V,.MEAN_KE
ncks $inp -A $out \
  -v .SBG_RT2,.SBG_THL2,.SBG_U2,.SBG_V2,.SBG_W2,.SBG_WRT,.SBG_WTHL,.SBG_TKE,.SBG_WU,.SBG_WV
ncks $inp -A $out \
  -v .RES_RT2,.RES_THL2,.RES_U2,.RES_V2,.RES_W2,.RES_WRT,.RES_WTHL,.RES_TKE,.RES_WU,.RES_WV
ncks $inp -A $out \
  -v .FRCUP_MF,.RCUP_MF,.RTUP_MF,.THLUP_MF,.THVUP_MF
ncks $inp -A $out \
  -v Q0,E0,.SWU,.SWD,.LWU,.LWD
ncks $inp -A $out \
  -v level_les,time_les

ncrename $out \
  -v MEANWIND,WND \
  -v MEAN_PRE,pf \
  -v MEAN_RC,ql \
  -v MEAN_RT,qt \
  -v MEAN_RV,qv \
  -v MEAN_RR,qr \
  -v MEAN_RHO,rho \
  -v MEAN_REHU,rh \
  -v MEAN_CF,rneb \
  -v MEAN_TH,theta \
  -v MEAN_THV,thetav \
  -v MEAN_THL,thl \
  -v MEAN_U,u \
  -v MEAN_V,v \
  -v .SBG_RT2,rt2_sbg \
  -v .SBG_THL2,thl2_sbg \
  -v .SBG_U2,uu_sbg \
  -v .SBG_V2,vv_sbg \
  -v .SBG_W2,ww_sbg \
  -v .SBG_WRT,wrt_sbg \
  -v .SBG_WTHL,wthl_sbg \
  -v .SBG_TKE,tke_sbg \
  -v .SBG_WU,uw_sbg \
  -v .SBG_WV,vw_sbg \
  -v .RES_RT2,rt2_res \
  -v .RES_THL2,thl2_res \
  -v .RES_U2,uu_res \
  -v .RES_V2,vv_res \
  -v .RES_W2,ww_res \
  -v .RES_WRT,wrt_res \
  -v .RES_WTHL,wthl_res \
  -v .RES_TKE,tke_res \
  -v .RES_WU,uw_res \
  -v .RES_WV,vw_res \
  -v .FRCUP_MF,alphau_shcon \
  -v .RCUP_MF,qlu_shcon \
  -v .RTUP_MF,qtu_shcon \
  -v .THLUP_MF,thlu_shcon \
  -v .THVUP_MF,thvu_shcon \
  -v .SWU,swu -v .SWD,swd -v .LWU,lwu -v .LWD,lwd \
  -v level_les,zf \
  -v time_les,time

ncap2 $out -O $out \
  -s "Q0=Q0*1004" \
  -s "E0=E0*2500000"\
  -s "temp=theta*(pf/100000.0)^0.286"

ncap2 $out -O $out \
  -s "rt2=rt2_sbg +  rt2_res" \
  -s "thl2=thl2_sbg + thl2_res"  \
  -s "uu=uu_sbg + uu_res"    \
  -s "vv=vv_sbg + vv_res"    \
  -s "ww=ww_sbg + ww_res"    \
  -s "wrt=wrt_sbg + wrt_res" \
  -s "wthl=wthl_sbg + wthl_res"  \
  -s "tke=tke_sbg + tke_res" \
  -s "uw=uw_sbg + uw_res" \
  -s "vw=vw_sbg + vw_res"