dimensions:
	nlon = 100 ;
	nlat = 116 ;
	time = UNLIMITED ; // (1 currently)
	d2 = 2 ;
variables:
	float time(time) ;
		time:long_name = "model time" ;
		time:units = "days since 0001-01-01 00:00:00" ;
		time:calendar = "noleap" ;
		time:bounds = "time_bounds" ;
	float time_bounds(time, d2) ;
		time_bounds:long_name = "boundaries for time-averaging interval" ;
		time_bounds:units = "days since 0001-01-01 00:00:00" ;
	float TLON(nlat, nlon) ;
		TLON:long_name = "T grid center longitude" ;
		TLON:units = "degrees_east" ;
	float TLAT(nlat, nlon) ;
		TLAT:long_name = "T grid center latitude" ;
		TLAT:units = "degrees_north" ;
	float ULON(nlat, nlon) ;
		ULON:long_name = "U grid center longitude" ;
		ULON:units = "degrees_east" ;
	float ULAT(nlat, nlon) ;
		ULAT:long_name = "U grid center latitude" ;
		ULAT:units = "degrees_north" ;
		ULAT:comment = "Latitude of NE corner of T grid cell" ;
	float tmask(nlat, nlon) ;
		tmask:long_name = "ocean grid mask" ;
		tmask:coordinates = "TLON TLAT" ;
		tmask:comment = "0 = land, 1 = ocean" ;
	float tarea(nlat, nlon) ;
		tarea:long_name = "area of T grid cells" ;
		tarea:units = "m^2" ;
		tarea:coordinates = "TLON TLAT" ;
	float uarea(nlat, nlon) ;
		uarea:long_name = "area of U grid cells" ;
		uarea:units = "m^2" ;
		uarea:coordinates = "ULON ULAT" ;
	float dxt(nlat, nlon) ;
		dxt:long_name = "T cell width through middle" ;
		dxt:units = "m" ;
		dxt:coordinates = "TLON TLAT" ;
	float dyt(nlat, nlon) ;
		dyt:long_name = "T cell height through middle" ;
		dyt:units = "m" ;
		dyt:coordinates = "TLON TLAT" ;
	float dxu(nlat, nlon) ;
		dxu:long_name = "U cell width through middle" ;
		dxu:units = "m" ;
		dxu:coordinates = "ULON ULAT" ;
	float dyu(nlat, nlon) ;
		dyu:long_name = "U cell height through middle" ;
		dyu:units = "m" ;
		dyu:coordinates = "ULON ULAT" ;
	float HTN(nlat, nlon) ;
		HTN:long_name = "T cell width on North side" ;
		HTN:units = "m" ;
		HTN:coordinates = "TLON TLAT" ;
	float HTE(nlat, nlon) ;
		HTE:long_name = "T cell width on East side" ;
		HTE:units = "m" ;
		HTE:coordinates = "TLON TLAT" ;
	float ANGLET(nlat, nlon) ;
		ANGLET:long_name = "angle grid makes with latitude line on T grid" ;
		ANGLET:units = "radians" ;
		ANGLET:coordinates = "TLON TLAT" ;
	float ANGLE(nlat, nlon) ;
		ANGLE:long_name = "angle grid makes with latitude line on U grid" ;
		ANGLE:units = "radians" ;
		ANGLE:coordinates = "ULON ULAT" ;
	float ice_present(time, nlat, nlon) ;
		ice_present:units = "1" ;
		ice_present:long_name = "fraction of time-averaging interval that any ice is present" ;
		ice_present:coordinates = "TLON TLAT time" ;
	float hi(time, nlat, nlon) ;
		hi:units = "m" ;
		hi:long_name = "grid box mean ice thickness" ;
		hi:coordinates = "TLON TLAT time" ;
		hi:missing_value = 1.e+30f ;
		hi:_FillValue = 1.e+30f ;
		hi:comment = "Ice volume distributed evenly over grid cell" ;
		hi:cell_methods = "time:mean" ;
	float hs(time, nlat, nlon) ;
		hs:units = "m" ;
		hs:long_name = "grid box mean snow thickness" ;
		hs:coordinates = "TLON TLAT time" ;
		hs:missing_value = 1.e+30f ;
		hs:_FillValue = 1.e+30f ;
		hs:comment = "Snow volume distributed evenly over grid cell" ;
		hs:cell_methods = "time:mean" ;
	float Tsfc(time, nlat, nlon) ;
		Tsfc:units = "C" ;
		Tsfc:long_name = "snow/ice surface temperature" ;
		Tsfc:coordinates = "TLON TLAT time" ;
		Tsfc:missing_value = 1.e+30f ;
		Tsfc:_FillValue = 1.e+30f ;
		Tsfc:comment = "Tf is averaged in if no ice is present" ;
		Tsfc:cell_methods = "time:mean" ;
	float aice(time, nlat, nlon) ;
		aice:units = "%" ;
		aice:long_name = "ice area  (aggregate)" ;
		aice:coordinates = "TLON TLAT time" ;
		aice:missing_value = 1.e+30f ;
		aice:_FillValue = 1.e+30f ;
		aice:comment = "none" ;
		aice:cell_methods = "time:mean" ;
	float u(time, nlat, nlon) ;
		u:units = "cm/s" ;
		u:long_name = "ice velocity (x)" ;
		u:coordinates = "ULON ULAT time" ;
		u:missing_value = 1.e+30f ;
		u:_FillValue = 1.e+30f ;
		u:comment = "positive is x direction on U grid" ;
		u:cell_methods = "time:mean" ;
	float v(time, nlat, nlon) ;
		v:units = "cm/s" ;
		v:long_name = "ice velocity (y)" ;
		v:coordinates = "ULON ULAT time" ;
		v:missing_value = 1.e+30f ;
		v:_FillValue = 1.e+30f ;
		v:comment = "positive is y direction on U grid" ;
		v:cell_methods = "time:mean" ;
	float Fswdn(time, nlat, nlon) ;
		Fswdn:units = "W/m^2" ;
		Fswdn:long_name = "down solar flux" ;
		Fswdn:coordinates = "TLON TLAT time" ;
		Fswdn:missing_value = 1.e+30f ;
		Fswdn:_FillValue = 1.e+30f ;
		Fswdn:comment = "positive downward" ;
		Fswdn:cell_methods = "time:mean" ;
	float Flwdn(time, nlat, nlon) ;
		Flwdn:units = "W/m^2" ;
		Flwdn:long_name = "down longwave flux" ;
		Flwdn:coordinates = "TLON TLAT time" ;
		Flwdn:missing_value = 1.e+30f ;
		Flwdn:_FillValue = 1.e+30f ;
		Flwdn:comment = "positive downward" ;
		Flwdn:cell_methods = "time:mean" ;
	float snow(time, nlat, nlon) ;
		snow:units = "cm/day" ;
		snow:long_name = "snowfall rate (cpl)" ;
		snow:coordinates = "TLON TLAT time" ;
		snow:missing_value = 1.e+30f ;
		snow:_FillValue = 1.e+30f ;
		snow:comment = "none" ;
		snow:cell_methods = "time:mean" ;
	float snow_aice(time, nlat, nlon) ;
		snow_aice:units = "cm/day" ;
		snow_aice:long_name = "snowfall rate" ;
		snow_aice:coordinates = "TLON TLAT time" ;
		snow_aice:missing_value = 1.e+30f ;
		snow_aice:_FillValue = 1.e+30f ;
		snow_aice:comment = "weighted by ice area" ;
		snow_aice:cell_methods = "time:mean" ;
	float rain(time, nlat, nlon) ;
		rain:units = "cm/day" ;
		rain:long_name = "rainfall rate (cpl)" ;
		rain:coordinates = "TLON TLAT time" ;
		rain:missing_value = 1.e+30f ;
		rain:_FillValue = 1.e+30f ;
		rain:comment = "none" ;
		rain:cell_methods = "time:mean" ;
	float rain_aice(time, nlat, nlon) ;
		rain_aice:units = "cm/day" ;
		rain_aice:long_name = "rainfall rate" ;
		rain_aice:coordinates = "TLON TLAT time" ;
		rain_aice:missing_value = 1.e+30f ;
		rain_aice:_FillValue = 1.e+30f ;
		rain_aice:comment = "weighted by ice area" ;
		rain_aice:cell_methods = "time:mean" ;
	float Fswabs(time, nlat, nlon) ;
		Fswabs:units = "W/m^2" ;
		Fswabs:long_name = "snow/ice/ocn absorbed solar flux (cpl)" ;
		Fswabs:coordinates = "TLON TLAT time" ;
		Fswabs:missing_value = 1.e+30f ;
		Fswabs:_FillValue = 1.e+30f ;
		Fswabs:comment = "positive downward" ;
		Fswabs:cell_methods = "time:mean" ;
	float Fswabs_aice(time, nlat, nlon) ;
		Fswabs_aice:units = "W/m^2" ;
		Fswabs_aice:long_name = "snow/ice/ocn absorbed solar flux" ;
		Fswabs_aice:coordinates = "TLON TLAT time" ;
		Fswabs_aice:missing_value = 1.e+30f ;
		Fswabs_aice:_FillValue = 1.e+30f ;
		Fswabs_aice:comment = "weighted by ice area" ;
		Fswabs_aice:cell_methods = "time:mean" ;
	float albsni(time, nlat, nlon) ;
		albsni:units = "%" ;
		albsni:long_name = "snw/ice broad band albedo" ;
		albsni:coordinates = "TLON TLAT time" ;
		albsni:missing_value = 1.e+30f ;
		albsni:_FillValue = 1.e+30f ;
		albsni:comment = "none" ;
		albsni:cell_methods = "time:mean" ;
	float Flat(time, nlat, nlon) ;
		Flat:units = "W/m^2" ;
		Flat:long_name = "latent heat flux (cpl)" ;
		Flat:coordinates = "TLON TLAT time" ;
		Flat:missing_value = 1.e+30f ;
		Flat:_FillValue = 1.e+30f ;
		Flat:comment = "positive downward" ;
		Flat:cell_methods = "time:mean" ;
	float Flat_aice(time, nlat, nlon) ;
		Flat_aice:units = "W/m^2" ;
		Flat_aice:long_name = "latent heat flux" ;
		Flat_aice:coordinates = "TLON TLAT time" ;
		Flat_aice:missing_value = 1.e+30f ;
		Flat_aice:_FillValue = 1.e+30f ;
		Flat_aice:comment = "weighted by ice area" ;
		Flat_aice:cell_methods = "time:mean" ;
	float Fsens(time, nlat, nlon) ;
		Fsens:units = "W/m^2" ;
		Fsens:long_name = "sensible heat flux (cpl)" ;
		Fsens:coordinates = "TLON TLAT time" ;
		Fsens:missing_value = 1.e+30f ;
		Fsens:_FillValue = 1.e+30f ;
		Fsens:comment = "positive downward" ;
		Fsens:cell_methods = "time:mean" ;
	float Fsens_aice(time, nlat, nlon) ;
		Fsens_aice:units = "W/m^2" ;
		Fsens_aice:long_name = "sensible heat flux" ;
		Fsens_aice:coordinates = "TLON TLAT time" ;
		Fsens_aice:missing_value = 1.e+30f ;
		Fsens_aice:_FillValue = 1.e+30f ;
		Fsens_aice:comment = "weighted by ice area" ;
		Fsens_aice:cell_methods = "time:mean" ;
	float Flwout(time, nlat, nlon) ;
		Flwout:units = "W/m^2" ;
		Flwout:long_name = "upward longwave flx (cpl)" ;
		Flwout:coordinates = "TLON TLAT time" ;
		Flwout:missing_value = 1.e+30f ;
		Flwout:_FillValue = 1.e+30f ;
		Flwout:comment = "positive downward" ;
		Flwout:cell_methods = "time:mean" ;
	float Flwout_aice(time, nlat, nlon) ;
		Flwout_aice:units = "W/m^2" ;
		Flwout_aice:long_name = "upward longwave flux" ;
		Flwout_aice:coordinates = "TLON TLAT time" ;
		Flwout_aice:missing_value = 1.e+30f ;
		Flwout_aice:_FillValue = 1.e+30f ;
		Flwout_aice:comment = "weighted by ice area" ;
		Flwout_aice:cell_methods = "time:mean" ;
	float evap(time, nlat, nlon) ;
		evap:units = "cm/day" ;
		evap:long_name = "evaporative water flux (cpl)" ;
		evap:coordinates = "TLON TLAT time" ;
		evap:missing_value = 1.e+30f ;
		evap:_FillValue = 1.e+30f ;
		evap:comment = "none" ;
		evap:cell_methods = "time:mean" ;
	float evap_aice(time, nlat, nlon) ;
		evap_aice:units = "cm/day" ;
		evap_aice:long_name = "evaporative water flux" ;
		evap_aice:coordinates = "TLON TLAT time" ;
		evap_aice:missing_value = 1.e+30f ;
		evap_aice:_FillValue = 1.e+30f ;
		evap_aice:comment = "weighted by ice area" ;
		evap_aice:cell_methods = "time:mean" ;
	float Tref(time, nlat, nlon) ;
		Tref:units = "C" ;
		Tref:long_name = "2m reference temperature" ;
		Tref:coordinates = "TLON TLAT time" ;
		Tref:missing_value = 1.e+30f ;
		Tref:_FillValue = 1.e+30f ;
		Tref:comment = "none" ;
		Tref:cell_methods = "time:mean" ;
	float Qref(time, nlat, nlon) ;
		Qref:units = "g/kg" ;
		Qref:long_name = "2m reference specific humidity" ;
		Qref:coordinates = "TLON TLAT time" ;
		Qref:missing_value = 1.e+30f ;
		Qref:_FillValue = 1.e+30f ;
		Qref:comment = "none" ;
		Qref:cell_methods = "time:mean" ;
	float growb(time, nlat, nlon) ;
		growb:units = "cm/day" ;
		growb:long_name = "congelation ice growth" ;
		growb:coordinates = "TLON TLAT time" ;
		growb:missing_value = 1.e+30f ;
		growb:_FillValue = 1.e+30f ;
		growb:comment = "none" ;
		growb:cell_methods = "time:mean" ;
	float frazil(time, nlat, nlon) ;
		frazil:units = "cm/day" ;
		frazil:long_name = "frazil ice growth" ;
		frazil:coordinates = "TLON TLAT time" ;
		frazil:missing_value = 1.e+30f ;
		frazil:_FillValue = 1.e+30f ;
		frazil:comment = "none" ;
		frazil:cell_methods = "time:mean" ;
	float snoice(time, nlat, nlon) ;
		snoice:units = "cm/day" ;
		snoice:long_name = "snow-ice formation" ;
		snoice:coordinates = "TLON TLAT time" ;
		snoice:missing_value = 1.e+30f ;
		snoice:_FillValue = 1.e+30f ;
		snoice:comment = "none" ;
		snoice:cell_methods = "time:mean" ;
	float meltt(time, nlat, nlon) ;
		meltt:units = "cm/day" ;
		meltt:long_name = "top ice melt" ;
		meltt:coordinates = "TLON TLAT time" ;
		meltt:missing_value = 1.e+30f ;
		meltt:_FillValue = 1.e+30f ;
		meltt:comment = "none" ;
		meltt:cell_methods = "time:mean" ;
	float meltb(time, nlat, nlon) ;
		meltb:units = "cm/day" ;
		meltb:long_name = "basal ice melt" ;
		meltb:coordinates = "TLON TLAT time" ;
		meltb:missing_value = 1.e+30f ;
		meltb:_FillValue = 1.e+30f ;
		meltb:comment = "none" ;
		meltb:cell_methods = "time:mean" ;
	float meltl(time, nlat, nlon) ;
		meltl:units = "cm/day" ;
		meltl:long_name = "lateral ice melt" ;
		meltl:coordinates = "TLON TLAT time" ;
		meltl:missing_value = 1.e+30f ;
		meltl:_FillValue = 1.e+30f ;
		meltl:comment = "none" ;
		meltl:cell_methods = "time:mean" ;
	float Fresh(time, nlat, nlon) ;
		Fresh:units = "cm/day" ;
		Fresh:long_name = "frshwtr flx ice to ocn (cpl)" ;
		Fresh:coordinates = "TLON TLAT time" ;
		Fresh:missing_value = 1.e+30f ;
		Fresh:_FillValue = 1.e+30f ;
		Fresh:comment = "positive represents ocean gaining fresh water" ;
		Fresh:cell_methods = "time:mean" ;
	float Fresh_aice(time, nlat, nlon) ;
		Fresh_aice:units = "cm/day" ;
		Fresh_aice:long_name = "frshwtr flx ice to ocn" ;
		Fresh_aice:coordinates = "TLON TLAT time" ;
		Fresh_aice:missing_value = 1.e+30f ;
		Fresh_aice:_FillValue = 1.e+30f ;
		Fresh_aice:comment = "weighted by ice area" ;
		Fresh_aice:cell_methods = "time:mean" ;
	float Fsalt(time, nlat, nlon) ;
		Fsalt:units = "kg/m^2/day" ;
		Fsalt:long_name = "salt flux ice to ocn (cpl)" ;
		Fsalt:coordinates = "TLON TLAT time" ;
		Fsalt:missing_value = 1.e+30f ;
		Fsalt:_FillValue = 1.e+30f ;
		Fsalt:comment = "positive represents ocean gaining salt" ;
		Fsalt:cell_methods = "time:mean" ;
	float Fsalt_aice(time, nlat, nlon) ;
		Fsalt_aice:units = "kg/m^2/day" ;
		Fsalt_aice:long_name = "salt flux ice to ocean" ;
		Fsalt_aice:coordinates = "TLON TLAT time" ;
		Fsalt_aice:missing_value = 1.e+30f ;
		Fsalt_aice:_FillValue = 1.e+30f ;
		Fsalt_aice:comment = "weighted by ice area" ;
		Fsalt_aice:cell_methods = "time:mean" ;
	float Fhnet(time, nlat, nlon) ;
		Fhnet:units = "W/m^2" ;
		Fhnet:long_name = "heat flux ice to ocn (cpl)" ;
		Fhnet:coordinates = "TLON TLAT time" ;
		Fhnet:missing_value = 1.e+30f ;
		Fhnet:_FillValue = 1.e+30f ;
		Fhnet:comment = "positive represents ocean gaining heat" ;
		Fhnet:cell_methods = "time:mean" ;
	float Fhnet_aice(time, nlat, nlon) ;
		Fhnet_aice:units = "W/m^2" ;
		Fhnet_aice:long_name = "heat flux ice to ocean" ;
		Fhnet_aice:coordinates = "TLON TLAT time" ;
		Fhnet_aice:missing_value = 1.e+30f ;
		Fhnet_aice:_FillValue = 1.e+30f ;
		Fhnet_aice:comment = "weighted by ice area" ;
		Fhnet_aice:cell_methods = "time:mean" ;
	float Fswthru(time, nlat, nlon) ;
		Fswthru:units = "W/m^2" ;
		Fswthru:long_name = "SW thru ice to ocean (cpl)" ;
		Fswthru:coordinates = "TLON TLAT time" ;
		Fswthru:missing_value = 1.e+30f ;
		Fswthru:_FillValue = 1.e+30f ;
		Fswthru:comment = "positive represents ice gaining heat" ;
		Fswthru:cell_methods = "time:mean" ;
	float Fswthru_aice(time, nlat, nlon) ;
		Fswthru_aice:units = "W/m^2" ;
		Fswthru_aice:long_name = "SW flux thru ice to ocean" ;
		Fswthru_aice:coordinates = "TLON TLAT time" ;
		Fswthru_aice:missing_value = 1.e+30f ;
		Fswthru_aice:_FillValue = 1.e+30f ;
		Fswthru_aice:comment = "weighted by ice area" ;
		Fswthru_aice:cell_methods = "time:mean" ;
	float strairx(time, nlat, nlon) ;
		strairx:units = "N/m^2" ;
		strairx:long_name = "atm/ice stress (x)" ;
		strairx:coordinates = "ULON ULAT time" ;
		strairx:missing_value = 1.e+30f ;
		strairx:_FillValue = 1.e+30f ;
		strairx:comment = "positive is x direction on U grid" ;
		strairx:cell_methods = "time:mean" ;
	float strairy(time, nlat, nlon) ;
		strairy:units = "N/m^2" ;
		strairy:long_name = "atm/ice stress (y)" ;
		strairy:coordinates = "ULON ULAT time" ;
		strairy:missing_value = 1.e+30f ;
		strairy:_FillValue = 1.e+30f ;
		strairy:comment = "positive is y direction on U grid" ;
		strairy:cell_methods = "time:mean" ;
	float strcorx(time, nlat, nlon) ;
		strcorx:units = "N/m^2" ;
		strcorx:long_name = "coriolis stress (x)" ;
		strcorx:coordinates = "ULON ULAT time" ;
		strcorx:missing_value = 1.e+30f ;
		strcorx:_FillValue = 1.e+30f ;
		strcorx:comment = "positive is x direction on U grid" ;
		strcorx:cell_methods = "time:mean" ;
	float strcory(time, nlat, nlon) ;
		strcory:units = "N/m^2" ;
		strcory:long_name = "coriolis stress (y)" ;
		strcory:coordinates = "ULON ULAT time" ;
		strcory:missing_value = 1.e+30f ;
		strcory:_FillValue = 1.e+30f ;
		strcory:comment = "positive is y direction on U grid" ;
		strcory:cell_methods = "time:mean" ;
	float strocnx(time, nlat, nlon) ;
		strocnx:units = "N/m^2" ;
		strocnx:long_name = "ocean/ice stress (x)" ;
		strocnx:coordinates = "ULON ULAT time" ;
		strocnx:missing_value = 1.e+30f ;
		strocnx:_FillValue = 1.e+30f ;
		strocnx:comment = "positive is x direction on U grid" ;
		strocnx:cell_methods = "time:mean" ;
	float strocny(time, nlat, nlon) ;
		strocny:units = "N/m^2" ;
		strocny:long_name = "ocean/ice stress (y)" ;
		strocny:coordinates = "ULON ULAT time" ;
		strocny:missing_value = 1.e+30f ;
		strocny:_FillValue = 1.e+30f ;
		strocny:comment = "positive is y direction on U grid" ;
		strocny:cell_methods = "time:mean" ;
	float strintx(time, nlat, nlon) ;
		strintx:units = "N/m^2" ;
		strintx:long_name = "internal ice stress (x)" ;
		strintx:coordinates = "ULON ULAT time" ;
		strintx:missing_value = 1.e+30f ;
		strintx:_FillValue = 1.e+30f ;
		strintx:comment = "positive is x direction on U grid" ;
		strintx:cell_methods = "time:mean" ;
	float strinty(time, nlat, nlon) ;
		strinty:units = "N/m^2" ;
		strinty:long_name = "internal ice stress (y)" ;
		strinty:coordinates = "ULON ULAT time" ;
		strinty:missing_value = 1.e+30f ;
		strinty:_FillValue = 1.e+30f ;
		strinty:comment = "positive is y direction on U grid" ;
		strinty:cell_methods = "time:mean" ;
	float strength(time, nlat, nlon) ;
		strength:units = "N/m" ;
		strength:long_name = "compressive ice strength" ;
		strength:coordinates = "TLON TLAT time" ;
		strength:missing_value = 1.e+30f ;
		strength:_FillValue = 1.e+30f ;
		strength:comment = "none" ;
		strength:cell_methods = "time:mean" ;
	float divu(time, nlat, nlon) ;
		divu:units = "%/day" ;
		divu:long_name = "strain rate (divergence)" ;
		divu:coordinates = "TLON TLAT time" ;
		divu:missing_value = 1.e+30f ;
		divu:_FillValue = 1.e+30f ;
		divu:comment = "none" ;
		divu:cell_methods = "time:mean" ;
	float shear(time, nlat, nlon) ;
		shear:units = "%/day" ;
		shear:long_name = "strain rate (shear)" ;
		shear:coordinates = "TLON TLAT time" ;
		shear:missing_value = 1.e+30f ;
		shear:_FillValue = 1.e+30f ;
		shear:comment = "none" ;
		shear:cell_methods = "time:mean" ;
	float sig1(time, nlat, nlon) ;
		sig1:units = "1" ;
		sig1:long_name = "norm. principal stress 1" ;
		sig1:coordinates = "ULON ULAT time" ;
		sig1:missing_value = 1.e+30f ;
		sig1:_FillValue = 1.e+30f ;
		sig1:comment = "sig1 is instantaneous" ;
		sig1:cell_methods = "time:mean" ;
	float sig2(time, nlat, nlon) ;
		sig2:units = "1" ;
		sig2:long_name = "norm. principal stress 2" ;
		sig2:coordinates = "ULON ULAT time" ;
		sig2:missing_value = 1.e+30f ;
		sig2:_FillValue = 1.e+30f ;
		sig2:comment = "sig2 is instantaneous" ;
		sig2:cell_methods = "time:mean" ;
	float dvidtt(time, nlat, nlon) ;
		dvidtt:units = "cm/day" ;
		dvidtt:long_name = "volume tendency thermo" ;
		dvidtt:coordinates = "TLON TLAT time" ;
		dvidtt:missing_value = 1.e+30f ;
		dvidtt:_FillValue = 1.e+30f ;
		dvidtt:comment = "none" ;
		dvidtt:cell_methods = "time:mean" ;
	float dvidtd(time, nlat, nlon) ;
		dvidtd:units = "cm/day" ;
		dvidtd:long_name = "volume tendency dynamics" ;
		dvidtd:coordinates = "TLON TLAT time" ;
		dvidtd:missing_value = 1.e+30f ;
		dvidtd:_FillValue = 1.e+30f ;
		dvidtd:comment = "none" ;
		dvidtd:cell_methods = "time:mean" ;
	float daidtt(time, nlat, nlon) ;
		daidtt:units = "%/day" ;
		daidtt:long_name = "area tendency thermo" ;
		daidtt:coordinates = "TLON TLAT time" ;
		daidtt:missing_value = 1.e+30f ;
		daidtt:_FillValue = 1.e+30f ;
		daidtt:comment = "none" ;
		daidtt:cell_methods = "time:mean" ;
	float daidtd(time, nlat, nlon) ;
		daidtd:units = "%/day" ;
		daidtd:long_name = "area tendency dynamics" ;
		daidtd:coordinates = "TLON TLAT time" ;
		daidtd:missing_value = 1.e+30f ;
		daidtd:_FillValue = 1.e+30f ;
		daidtd:comment = "none" ;
		daidtd:cell_methods = "time:mean" ;
	float opening(time, nlat, nlon) ;
		opening:units = "%/day" ;
		opening:long_name = "lead area opening rate" ;
		opening:coordinates = "TLON TLAT time" ;
		opening:missing_value = 1.e+30f ;
		opening:_FillValue = 1.e+30f ;
		opening:comment = "none" ;
		opening:cell_methods = "time:mean" ;
	float aice1(time, nlat, nlon) ;
		aice1:units = "%" ;
		aice1:long_name = "ice area (category 1)" ;
		aice1:coordinates = "TLON TLAT time" ;
		aice1:missing_value = 1.e+30f ;
		aice1:_FillValue = 1.e+30f ;
		aice1:comment = "Ice range: 0.000   m < hi Cat  1 < 0.645   m" ;
		aice1:cell_methods = "time:mean" ;
	float aice2(time, nlat, nlon) ;
		aice2:units = "%" ;
		aice2:long_name = "ice area (category 2)" ;
		aice2:coordinates = "TLON TLAT time" ;
		aice2:missing_value = 1.e+30f ;
		aice2:_FillValue = 1.e+30f ;
		aice2:comment = "Ice range: 0.645   m < hi Cat  2 < 1.391   m" ;
		aice2:cell_methods = "time:mean" ;
	float aice3(time, nlat, nlon) ;
		aice3:units = "%" ;
		aice3:long_name = "ice area (category 3)" ;
		aice3:coordinates = "TLON TLAT time" ;
		aice3:missing_value = 1.e+30f ;
		aice3:_FillValue = 1.e+30f ;
		aice3:comment = "Ice range: 1.391   m < hi Cat  3 < 2.470   m" ;
		aice3:cell_methods = "time:mean" ;
	float aice4(time, nlat, nlon) ;
		aice4:units = "%" ;
		aice4:long_name = "ice area (category 4)" ;
		aice4:coordinates = "TLON TLAT time" ;
		aice4:missing_value = 1.e+30f ;
		aice4:_FillValue = 1.e+30f ;
		aice4:comment = "Ice range: 2.470   m < hi Cat  4 < 4.567   m" ;
		aice4:cell_methods = "time:mean" ;
	float aice5(time, nlat, nlon) ;
		aice5:units = "%" ;
		aice5:long_name = "ice area (category 5)" ;
		aice5:coordinates = "TLON TLAT time" ;
		aice5:missing_value = 1.e+30f ;
		aice5:_FillValue = 1.e+30f ;
		aice5:comment = "Ice range: 4.567   m < hi Cat  5 < 9.334   m" ;
		aice5:cell_methods = "time:mean" ;
	float alvdr(time, nlat, nlon) ;
		alvdr:units = "%" ;
		alvdr:long_name = "visible direct albedo" ;
		alvdr:coordinates = "TLON TLAT time" ;
		alvdr:missing_value = 1.e+30f ;
		alvdr:_FillValue = 1.e+30f ;
		alvdr:comment = "none" ;
		alvdr:cell_methods = "time:mean" ;
	float alidr(time, nlat, nlon) ;
		alidr:units = "%" ;
		alidr:long_name = "near IR direct albedo" ;
		alidr:coordinates = "TLON TLAT time" ;
		alidr:missing_value = 1.e+30f ;
		alidr:_FillValue = 1.e+30f ;
		alidr:comment = "none" ;
		alidr:cell_methods = "time:mean" ;