load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"   
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"   

begin

  fname = "template_ts.nc"
  fout  = addfile(fname,"c")

  fileAtt = True
  fileAtt@MissingValue  = -999
  fileAtt@FillValue     = -999
  fileAtt@date  = systemfunc("date")
  fileAtt@group = "UCLA"
  fileAtt@author= "Bjorn Stevens"
  fileAtt@case  = "Base Case"
  fileAtt@title = "GCSS 7 Time-series Output"

  fileattdef(fout,fileAtt)
  filedimdef(fout,"time",-1,True)

  filevardef(fout,"time","float","time")
  varAtt = True
  varAtt@units      = "[s]"
  varAtt@longname   = "Time"
  filevarattdef(fout,"time",varAtt)
  delete (varAtt)

  filevardef(fout,"zi_bar","float","time")
  varAtt = True
  varAtt@units      = "[m]"
  varAtt@longname   = "Spatially averaged height of 295 K isentrope"
  filevarattdef(fout,"zi_bar",varAtt)

  filevardef(fout,"zi_var","float","time")
  varAtt = True
  varAtt@units      = "[m:S:2:N:]"
  varAtt@longname   = "Spatial Variance of 295 K isentrope"
  filevarattdef(fout,"zi_var",varAtt)

  filevardef(fout,"zb_bar","float","time")
  varAtt = True
  varAtt@units      = "[m]"
  varAtt@longname   = "Spatially averaged cloud base height"
  filevarattdef(fout,"zb_bar",varAtt)

  filevardef(fout,"zb_var","float","time")
  varAtt = True
  varAtt@units      = "[m:S:2:N:]"
  varAtt@longname   = "Spatial Variance of cloud base"
  filevarattdef(fout,"zb_var",varAtt)

  filevardef(fout,"cfrac","float","time")
  varAtt = True
  varAtt@units      = "[%]"
  varAtt@longname   = "Fraction of columns with cloud present"
  filevarattdef(fout,"cfrac",varAtt)

  filevardef(fout,"lwp_bar","float","time")
  varAtt = True
  varAtt@units      = "[gm:S:-2:N:]"
  varAtt@longname   = "Spatial average of vertically integrated liquid water"
  filevarattdef(fout,"lwp_bar",varAtt)

  filevardef(fout,"lwp_var","float","time")
  varAtt = True
  varAtt@units      = "[g:S:2:N:m:S:-4:N:]"
  varAtt@longname   = "Spatial variance of vertically integrated liquid water"
  filevarattdef(fout,"lwp_var",varAtt)

  filevardef(fout,"tke","float","time")
  varAtt = True
  varAtt@units      = "[m:S:3:N:s:S:2:N:]"
  varAtt@longname   = "Vertically integrated TKE (resolved plus subfilter)"
  filevarattdef(fout,"tke",varAtt)

  filevardef(fout,"w2","float","time")
  varAtt = True
  varAtt@units      = "[m:S:2:N:s:S:-2:N:]"
  varAtt@longname   = "Maximum value of layer averaged vertical velocity variance"
  filevarattdef(fout,"w2",varAtt)

  filevardef(fout,"wstar","float","time")
  varAtt = True
  varAtt@units      = "[ms:S:-1:N:]"
  varAtt@longname   = "Deardorff velocity scale"
  filevarattdef(fout,"wstar",varAtt)

  filevardef(fout,"lhf_bar","float","time")
  varAtt = True
  varAtt@units      = "[Wm:S:-2:N:]"
  varAtt@longname   = "Spatially averaged latent heat flux"
  filevarattdef(fout,"lhf_bar",varAtt)

  filevardef(fout,"shf_bar","float","time")
  varAtt = True
  varAtt@units      = "[Wm:S:-2:N:]"
  varAtt@longname   = "Spatially averaged sensible heat flux"
  filevarattdef(fout,"shf_bar",varAtt)

  filevardef(fout,"ustar","float","time")
  varAtt = True
  varAtt@units      = "[ms:S:-1:N:]"
  varAtt@longname   = "Domain averaged surface friction velocity"
  filevarattdef(fout,"ustar",varAtt)

end