load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin 
  initial_time = getenv("initial_time")
  output_time = getenv("output_time")
  inout_file = getenv("inout_file")
  tothr = tointeger(getenv("tothr"))
  yyyy = tointeger(getenv("yyyy"))
  mon = tointeger(getenv("mon"))
  dd = tointeger(getenv("dd"))
  hh = tointeger(getenv("hh"))
  hhe=hh+tothr
  units="minutes since 1970-01-01 00:00:00 UTC"
  time_min = tointeger(cd_inv_calendar(yyyy,mon,dd,hh,0,0,units, 0))
  time_max = tointeger(cd_inv_calendar(yyyy,mon,dd,hhe,0,0,units, 0))
  units="minutes since 1970-01-01 00:00:00 UTC"
  reference_time = cd_inv_calendar( yyyy, mon, dd, hh, 00, 00, units, 0 )
  reference_time@units = "minutes since 1970-01-01 00:00:00 UTC" ;
  reference_time@standard_name = "forecast_reference_time" ;
  reference_time@long_name = "model initialization time" ;

  f=addfile(inout_file,"w")
  f@model_initialization_time = initial_time
  f@model_output_valid_time  = output_time
  f@model_output_type="forcing"
  f@model_configuration="analysis_and_assimilation"
  f@model_total_valid_times=tothr
  time=f->time
  time@valid_max=time_max
  time@valid_min=time_min

  f->reference_time=doubletoint(reference_time)

  print("End to replace the precipitation for analysis.")
  exit

end
status_exit(1)