#[]
#dimensions:    time, level, y1, x1
#units:           The documentation will complain if units have LaTeX characters like ^
#long_name:     
#data:          
#factor:        

# setup for code generator
[VARSET]
# prefix of the generated module
name:     glide
# f90 type containing all necessary data
datatype: glide_global_type   
# module where type is defined
datamod:  glide_types     

[x0]
dimensions:    x0
units:         meter
long_name:     Cartesian x-coordinate, velocity grid
axis:          X
data:          data%general%x0
dimlen:        global_ewn-1

[y0]
dimensions:    y0
units:         meter
long_name:     Cartesian y-coordinate, velocity grid
axis:          Y
data:          data%general%y0
dimlen:        global_nsn-1

[x1]
dimensions:    x1
units:         meter
long_name:     Cartesian x-coordinate
axis: X
data:          data%general%x1
dimlen:        global_ewn
load:1

[y1]
dimensions:    y1
units:         meter
long_name:     Cartesian y-coordinate
axis:          Y
data:          data%general%y1
dimlen:        global_nsn
load:1

[level]
dimensions:    level
units:         1
long_name:     sigma layers
standard_name: land_ice_sigma_coordinate
#formula_terms: sigma: level topo: topg thick: thk
positive:      down
dimlen:        model%general%upn

[staglevel]
dimensions:    staglevel
units:         1
long_name:     stag sigma layers
standard_name: land_ice_stag_sigma_coordinate
positive:      down
dimlen:        model%general%upn-1

[stagwbndlevel]
dimensions:    stagwbndlevel
units:         1
long_name:     stag sigma layers with boundaries
standard_name: land_ice_stag_sigma_coordinate_with_bnd
positive:      down
dimlen:        model%general%upn+1

[lithoz]
dimensions:    lithoz
units:         meter
long_name:     vertical coordinate of lithosphere layer
dimlen:        model%lithot%nlayer

[relx]
dimensions:    time, y1, x1
units:         meter
long_name:     relaxed bedrock topography
data:          data%isostasy%relx
factor:        thk0
load:          1
coordinates:   lon lat

[eus]
dimensions:    time
units:         meter
long_name:     global average sea level
data:          data%climate%eus
factor:        thk0
standard_name: global_average_sea_level_change

[uflx]
dimensions:    time, y0, x0
units:         meter2/year
long_name:     flux in x direction (NOTE: Glide and Glam only)
data:          data%velocity%uflx
factor:        scale_uflx

[vflx]
dimensions:    time, y0, x0
units:         meter2/year
long_name:     flux in x direction (NOTE: Glide and Glam only)
data:          data%velocity%vflx
factor:        scale_uflx

[diffu]
dimensions:    time, y0, x0
units:         meter2/year
long_name:     apparent diffusivity
data:          data%velocity%diffu
factor:        scale_diffu

[btrc]
dimensions:    time, y0, x0
units:         meter/pascal/year
long_name:     basal slip coefficient
data:          data%velocity%btrc
factor:        scale_btrc

[soft]
dimensions:    time, y0, x0
units:         meter/pascal/year
long_name:     bed softness parameter
data:          data%velocity%bed_softness
factor:        scale_btrc
load:          1

[beta]
dimensions:    time, y0, x0
units:         Pa yr/m  
long_name:     higher-order bed stress coefficient
data:          data%velocity%beta
factor:        scale_beta
load:          1

#WHL - added unstagbeta, which lives on the unstaggered ice grid
#      (unlike beta, which lives on the staggered velocity grid)
[unstagbeta]
dimensions:    time, y1, x1
units:         Pa yr/m  
long_name:     higher-order bed stress coefficient on the unstaggered grid (NOTE: this will overwrite beta if both are input)
data:          data%velocity%unstagbeta
factor:        scale_beta
load:          1

#WHL - beta_internal is derived from beta, but may be weighted by f_ground or otherwise adjusted
#      It can be written to history files but should not be read from restart files
[beta_internal]
dimensions:    time, y0, x0
units:         Pa yr/m
long_name:     weighted higher-order bed stress coefficient
data:          data%velocity%beta_internal
factor:        scale_beta

[tauf]
dimensions:    time, y0, x0
units:         Pa
long_name:     higher-order basal yield stress
data:          data%basalproc%mintauf
factor:        scale_tau
load:          1

[btractx]
dimensions:    time, y0, x0
units:         Pa
long_name:     basal traction (x-direction comp)
data:          data%stress%btractx(:,:)
factor:        scale_tau
load:          1

[btracty]
dimensions:    time, y0, x0
units:         Pa
long_name:     basal traction (y-direction comp)
data:          data%stress%btracty(:,:)
factor:        scale_tau
load:          1

[btractx_extend]
dimensions:    time, y1, x1
units:         Pa
long_name:     basal traction (x-direction comp)
data:          data%stress%btractx_extend(:,:)
factor:        scale_tau
load:          1

[btracty_extend]
dimensions:    time, y1, x1
units:         Pa
long_name:     basal traction (y-direction comp)
data:          data%stress%btracty_extend(:,:)
factor:        scale_tau
load:          1

[ubas]
dimensions:    time, y0, x0
units:         meter/year
long_name:     basal slip velocity in x direction
data:          data%velocity%ubas
factor:        scale_uvel
standard_name: land_ice_basal_x_velocity
load:          1
#average:       1

[vbas]
dimensions:    time, y0, x0
units:         meter/year
long_name:     basal slip velocity in y direction
data:          data%velocity%vbas
factor:        scale_uvel
standard_name: land_ice_basal_y_velocity
load:          1
#average:       1

[taux]
dimensions:    time, y0, x0
units:         kilopascal
long_name:     basal shear stress in x direction (NOTE: Glide only)
data:          data%velocity%tau_x
factor:        1e-3*thk0*thk0/len0

[tauy]
dimensions:    time, y0, x0
units:         kilopascal
long_name:     basal shear stress in y direction
data:          data%velocity%tau_y
factor:        1e-3*thk0*thk0/len0

[taudx]
dimensions:    time, y0, x0
units:         Pa
long_name:     driving stress (x-direction comp)
data:          data%stress%taudx(:,:)
factor:        scale_tau
load:          1

[taudy]
dimensions:    time, y0, x0
units:         Pa
long_name:     driving stress (y-direction comp)
data:          data%stress%taudy(:,:)
factor:        scale_tau
load:          1

[thk]
dimensions:    time, y1, x1
units:         meter
long_name:     ice thickness
data:          data%geometry%thck
factor:        thk0
load:          1
standard_name: land_ice_thickness
coordinates:   lon lat

[stagthk]
dimensions:    time, y0, x0
units:         meter
long_name:     staggered ice thickness
data:          data%geomderv%stagthck
factor:        thk0
standard_name: stag_land_ice_thickness
load:          0
coordinates:   lon lat

[calving]
dimensions:    time, y1, x1
units:         meter
long_name:     ice margin calving
data:          data%calving%calving_thck
factor:        thk0
coordinates:   lon lat

[damage]
dimensions:    time, staglevel, y1, x1
units:         unitless [0,1]
long_name:     ice damage
data:          data%calving%damage(up,:,:)
load:          1
coordinates:   lon lat

[damage_column]
dimensions:    time, y1, x1
units:         unitless [0,1]
long_name:     vertically integrated ice damage
data:          data%calving%damage_column
coordinates:   lon lat

[ivol]
dimensions:    time
units:         km3
factor:        thk0*len0*len0*1.e-9
long_name:     ice volume
data:          data%geometry%ivol

[iarea]
dimensions:    time
units:         km2
long_name:     area covered by ice 
factor:        len0*len0*1.e-6
data:          data%geometry%iarea

[iareag]
dimensions:    time
units:         km2
long_name:     area covered by grounded ice 
factor:        len0*len0*1.e-6
data:          data%geometry%iareag

[iareaf]
dimensions:    time
units:         km2
long_name:     area covered by floating ice 
factor:        len0*len0*1.e-6
data:          data%geometry%iareaf

[thkmask]
dimensions:    time, y1, x1
long_name:     mask
units:         1
data:          data%geometry%thkmask
type:          int
coordinates:   lon lat
load:          1

[usurf]
dimensions:    time, y1, x1
units:         meter
long_name:     ice upper surface elevation
data:          data%geometry%usrf
factor:        thk0
load:          1
standard_name: surface_altitude
coordinates:   lon lat

[lsurf]
dimensions:    time, y1, x1
units:         meter
long_name:     ice lower surface elevation
data:          data%geometry%lsrf
factor:        thk0
coordinates:   lon lat

[topg]
dimensions:    time, y1, x1
units:         meter
long_name:     bedrock topography
data:          data%geometry%topg
factor:        thk0
load:          1
standard_name: bedrock_altitude
coordinates:   lon lat

## D. Martin added - fields that need to be passed to POP for ice-ocean coupling
#[floating_mask]
#dimensions:    time, y1, x1
#units:         1
#long_name:     real-valued mask denoting grounded/floating 
#data:          data%geometry%floating_mask
#factor:        1.0
#coordinates:   lon lat

# D. Martin added - fields that need to be passed to POP for ice-ocean coupling
[ice_mask]
dimensions:    time, y1, x1
units:         1
long_name:     real-valued mask denoting ice (1) or no ice (0)
data:          data%geometry%ice_mask
factor:        1.0
coordinates:   lon lat

## D. Martin added - fields that need to be passed to POP for ice-ocean coupling
#[lower_cell_loc]
#dimensions:    time, y1, x1
#units:         meter
#long_name:     location in z of lower cell location
#data:          data%geometry%lower_cell_loc
#factor:        1.0
#coordinates:   lon lat

## D. Martin added - fields that need to be passed to POP for ice-ocean coupling
#[lower_cell_temp]
#dimensions:    time, y1, x1
#units:         degrees K
#long_name:     temperature at lower_cell_loc
#data:          data%geometry%lower_cell_temp
#factor:        1.0
#coordinates:   lon lat

[acab]
dimensions:    time, y1, x1
units:         meter/year
long_name:     accumulation, ablation rate
data:          data%climate%acab
factor:        scale_acab
standard_name: land_ice_surface_specific_mass_balance
coordinates:   lon lat
#average:       1
load:          1

[flux_correction]
dimensions:    time, y1, x1
units:         meter/year
long_name:     flux correction applied in addition to acab
data:          data%climate%flux_correction
factor:        scale_acab
standard_name: land_ice_surface_specific_mass_balance_flux_correction
coordinates:   lon lat
load:          1

[no_advance_mask]
dimensions:    time, y1, x1
units:         1
long_name:     cells where advance is not allowed
data:          data%climate%no_advance_mask
type:          int
standard_name: land_ice_no_advance_mask
coordinates:   lon lat
load:          1

#WHL: scale_bflux = -1, to reverse sign convention from + up to + down
[bheatflx]
dimensions:    time, y1, x1
units:         watt/meter2
long_name:     upward basal heat flux
data:          data%temper%bheatflx
factor:        scale_bflx
load:          1
coordinates:   lon lat

[bmlt_ground]
dimensions:    time, y1, x1
units:         meter/year
long_name:     basal melt rate for grounded ice
data:          data%temper%bmlt_ground
factor:        scale_acab
standard_name: land_ice_basal_melt_rate_grounded
load:          1
coordinates:   lon lat
#average:       1

[bmlt_float]
dimensions:    time, y1, x1
units:         meter/year
long_name:     basal melt rate for floating ice
data:          data%temper%bmlt_float
factor:        scale_acab
standard_name: land_ice_basal_melt_rate_floating
load:          1
coordinates:   lon lat
#average:       1

[bmlt_float_mask]
dimensions:    time, y1, x1
units:         meter/year
long_name:     mask for basal melting of floating ice
data:          data%temper%bmlt_float_mask
standard_name: land_ice_basal_melt_rate_floating_mask
load:          1
coordinates:   lon lat

[bfricflx]
dimensions:    time, y1, x1
units:         watt/meter2
long_name:     basal friction heat flux
data:          data%temper%bfricflx
factor:        1.0
load:          1
coordinates:   lon lat

[bwat]
dimensions:    time, y1, x1
units:         meter
long_name:     basal water depth
data:          data%temper%bwat
factor:        thk0
load:          1
coordinates:   lon lat

[bwatflx]
dimensions:    time, y1, x1
units:         meter3/year
long_name:     basal water flux
data:          data%temper%bwatflx
factor:        thk0
coordinates:   lon lat

[effecpress]
dimensions:    time, y1, x1
units:         Pa
long_name:     effective pressure
data:          data%basal_physics%effecpress
load:          1
coordinates:   lon lat

[C_space_factor]
dimensions:    time, y1, x1
units:         1
long_name:     spatial factor for basal shear stress
data:          data%basal_physics%C_space_factor
load:          1
coordinates:   lon lat

[artm]
dimensions:    time, y1, x1
units:         degree_Celsius
long_name:     annual mean air temperature
data:          data%climate%artm
standard_name: surface_temperature
cell_methods:  time: mean
coordinates:   lon lat
load:          1

[surftemp]
dimensions:    time, y1, x1
units:         degree_Celsius
long_name:     annual mean surface temperature
data:          data%climate%artm
standard_name: surface_temperature
cell_methods:  time: mean
coordinates:   lon lat
load:          1

[btemp]
dimensions:    time, y1, x1
units:         degree_Celsius
long_name:     basal ice temperature
data:          data%temper%temp(data%general%upn,1:data%general%ewn,1:data%general%nsn)
standard_name: land_ice_temperature
coordinates:   lon lat

[dusrfdtm]
dimensions:    time, y1, x1
units:         meter/year
long_name:     rate of upper ice surface elevation change (NOTE: Glide only)
data:          data%geomderv%dusrfdtm
factor:        scale_acab
coordinates:   lon lat

[dthckdtm]
dimensions:    time, y1,x1
units:         meter/year
long_name:     tendency of ice thickness (NOTE: Glide only)
data:          data%geomderv%dthckdtm
factor:        scale_acab
coordinates:   lon lat

[uvel]
dimensions:    time, level, y0, x0
units:         meter/year
long_name:     ice velocity in x direction
data:          data%velocity%uvel(up,:,:)
factor:        scale_uvel
standard_name: land_ice_x_velocity
load:          1

[vvel]
dimensions:    time, level, y0, x0
units:         meter/year
long_name:     ice velocity in y direction
data:          data%velocity%vvel(up,:,:)
factor:        scale_uvel
standard_name: land_ice_y_velocity
load:          1

[uvel_extend]
dimensions:    time, level, y1, x1
units:         meter/year
long_name:     ice velocity in x direction (extended grid)
data:          data%velocity%uvel_extend(up,:,:)
factor:        scale_uvel
standard_name: land_ice_x_velocity
load:          1

[vvel_extend]
dimensions:    time, level, y1, x1
units:         meter/year
long_name:     ice velocity in y direction (extended grid)
data:          data%velocity%vvel_extend(up,:,:)
factor:        scale_uvel
standard_name: land_ice_y_velocity
load:          1

[uvel_2d]
dimensions:    time, y0, x0
units:         meter/year
long_name:     vertically averaged ice velocity in x direction
data:          data%velocity%uvel_2d(:,:)
factor:        scale_uvel
standard_name: land_ice_x_velocity_2d
load:          1

[vvel_2d]
dimensions:    time, y0, x0
units:         meter/year
long_name:     vertically averaged ice velocity in y direction
data:          data%velocity%vvel_2d(:,:)
factor:        scale_uvel
standard_name: land_ice_y_velocity_2d
load:          1

[uvel_2d_extend]
dimensions:    time, y1, x1
units:         meter/year
long_name:     vertically averaged ice velocity in x direction (extended grid)
data:          data%velocity%uvel_2d_extend(:,:)
factor:        scale_uvel
standard_name: land_ice_x_velocity_2d
load:          1

[vvel_2d_extend]
dimensions:    time, y1, x1
units:         meter/year
long_name:     vertically averaged ice velocity in y direction (extended grid)
data:          data%velocity%vvel_2d_extend(:,:)
factor:        scale_uvel
standard_name: land_ice_y_velocity_2d
load:          1

[resid_u]
dimensions:     time, level, y0, x0
units:          Pa/m          
long_name:      u component of residual Ax - b (NOTE: Glam only)
data:           data%velocity%resid_u(up,:,:)
factor:         scale_resid

[resid_v]
dimensions:     time, level, y0, x0
units:          Pa/m          
long_name:      v component of residual Ax - b (NOTE: Glam only)
data:           data%velocity%resid_v(up,:,:)
factor:         scale_resid

[rhs_u]
dimensions:     time, level, y0, x0
units:          Pa/m          
long_name:      u component of b in Ax = b
data:           data%velocity%rhs_u(up,:,:)
factor:         scale_resid

[rhs_v]
dimensions:     time, level, y0, x0
units:          Pa/m          
long_name:      v component of b in Ax = b
data:           data%velocity%rhs_v(up,:,:)
factor:         scale_resid

# used for output of residual
#[ures]
#dimensions:    time, level, y0, x0
#units:         meter/year
#long_name:     ice velocity resid. in x direction
#data:          data%velocity%ures(up,:,:)
#factor:        scale_uvel
#standard_name: land_ice_x_velocity_resid

# used for output of residual
#[vres]
#dimensions:    time, level, y0, x0
#units:         meter/year
#long_name:     ice velocity resid. in y direction
#data:          data%velocity%vres(up,:,:)
#factor:        scale_uvel
#standard_name: land_ice_y_velocity_resid

# used for output of res fields
#[magres]
#dimensions:    time, level, y0, x0
#units:         meter/year
#long_name:     ice velocity resid. magnitude 
#data:          data%velocity%magres(up,:,:)
#factor:        scale_uvel
#standard_name: land_ice_y_velocity_resid

[kinbcmask]
dimensions:   time, y0, x0
units:        1
long_name:    Mask of locations where uvel, vvel value should be held constant
data:         data%velocity%kinbcmask(:,:)
type:         int
load:         1

[dynbcmask]
dimensions:   time, y0, x0
units:        1
long_name:    2d array of higher-order model boundary condition mask values (NOTE: Glam ONLY)
data:         data%velocity%dynbcmask
type:         int
load:         0

[velnorm]
dimensions:    time, level, y0, x0
units:         meter/year
long_name:     Horizontal ice velocity magnitude
data:          data%velocity%velnorm(up,:,:)
factor:        scale_uvel
coordinates:   lon lat

[tau_eff]
dimensions: time, staglevel, y1, x1
units:      Pa
long_name:  effective stress
data:       data%stress%tau%scalar(up,:,:)
factor:     scale_tau

[tau_xz]
dimensions: time, staglevel, y1, x1
units:      Pa
long_name:  X component vertical shear stress 
data:       data%stress%tau%xz(up,:,:)
factor:     scale_tau

[tau_yz]
dimensions: time, staglevel, y1, x1
units:      Pa
long_name:  Y component vertical shear stress 
data:       data%stress%tau%yz(up,:,:)
factor:     scale_tau

[tau_xx]
dimensions: time, staglevel, y1, x1
units:      Pa
long_name:  x component of horiz. normal stress
data:       data%stress%tau%xx(up,:,:)
factor:     scale_tau

[tau_yy]
dimensions: time, staglevel, y1, x1
units:      Pa
long_name:  y component of horiz. normal stress 
data:       data%stress%tau%yy(up,:,:)
factor:     scale_tau

[tau_xy]
dimensions: time, staglevel, y1, x1
units:      Pa
long_name:  horiz. shear stress
data:       data%stress%tau%xy(up,:,:)
factor:     scale_tau

[wvel]
dimensions:    time, level, y1, x1
units:         meter/year
long_name:     vertical ice velocity
data:          data%velocity%wvel(up,:,:)
factor:        scale_wvel
coordinates:   lon lat
load:          1

[wvel_ho]
dimensions:    time, level, y1, x1
units:         meter/year
long_name:     vertical ice velocity relative to ice sheet base from higher-order model (NOTE: Glam only)
data:          data%velocity%wvel_ho(up,:,:)
factor:        scale_wvel
coordinates:   lon lat
load:          1

[wgrd]
dimensions:    time, level, y1, x1
units:         meter/year
long_name:     Vertical grid velocity 
data:          data%velocity%wgrd(up,:,:)
factor:        scale_wvel
coordinates:   lon lat
load:          1

[waterfrac]
dimensions:    time, staglevel, y1, x1
units:         unitless [0,1]
long_name:     internal water fraction
data:          data%temper%waterfrac(up,:,:)
coordinates:   lon lat
load:          1

[enthalpy]
dimensions:    time, stagwbndlevel, y1, x1
units:         J/m^3
long_name:     specific enthalpy
data:          data%temper%enthalpy(up,:,:)
coordinates:   lon lat

[flwa]
dimensions:    time, level, y1, x1
units:         pascal**(-n) year**(-1)
long_name:     Pre-exponential flow law parameter
data:          data%temper%flwa(up,:,:)
factor:        scale_flwa
load:          1
coordinates:   lon lat

[flwastag]
dimensions:    time, staglevel, y1, x1
units:         pascal**(-n) year**(-1)
long_name:     Pre-exponential flow law parameter
data:          data%temper%flwa(up,:,:)
factor:        scale_flwa
# hot=0 because we do not want the stag version to be in the hotvars definition. 
# The conversion to stag happens if needed in check_for_tempstag
# but we do want it loadable in case flwa is listed as hot.
load:          1
coordinates:   lon lat

[efvs]
dimensions:    time, staglevel, y1, x1
units:         Pascal * years
long_name:     effective viscosity               
data:          data%stress%efvs(up,:,:)
factor:        scale_efvs
coordinates:   lon lat
load:          1

[temp]
dimensions:    time, level, y1, x1
units:         degree_Celsius
long_name:     ice temperature
data:          data%temper%temp(up,1:data%general%ewn,1:data%general%nsn)
standard_name: land_ice_temperature
load:          1
coordinates:   lon lat

[tempstag]
dimensions:    time, stagwbndlevel, y1, x1
units:         degree_Celsius
long_name:     ice temperature on staggered vertical levels with boundaries
data:          data%temper%temp(up,1:data%general%ewn,1:data%general%nsn)
standard_name: land_ice_temperature_stag
# hot=0 because we do not want the stag version to be in the hotvars definition.
# The conversion to stag happens if needed in check_for_tempstag.
load:          1
coordinates:   lon lat

[litho_temp]
dimensions:    time, lithoz, y1, x1
units:         degree_Celsius
long_name:     lithosphere temperature
data:          data%lithot%temp
load:          1
coordinates:   lon lat

[ice_age]
dimensions:    time, staglevel, y1, x1
units:         year
long_name:     ice age
data:          data%geometry%ice_age(up,:,:)
standard_name: land_ice_age
factor:        tim0/scyr
load:          1
coordinates:   lon lat

[f_flotation]
dimensions:    time, y1, x1
units:         unitless
long_name:     flotation function
data:          data%geometry%f_flotation
load:          0
coordinates:   lon lat

[f_ground]
dimensions:    time, y0, x0
units:         unitless [0,1]
long_name:     grounded ice fraction
data:          data%geometry%f_ground
standard_name: grounded_fraction
load:          0
coordinates:   lon lat

#TODO - Are these gline fields needed?
#[gl_ew]
#dimensions:    time, y1, x0
#units:         meter
#long_name:     grounding line movement in ew directions
#data:          data%ground%gl_ew
#coordinates:   lon lat

#[gl_ns]
#dimensions:    time, y0, x1
#units:         meter
#long_name:     grounding line movement in ns directions
#data:          data%ground%gl_ns
#coordinates:   lon lat

#[gline_flux]
#dimensions:   time, y1, x1
#units:        meter2/year
#long_name:    grounding line fluxu
#data:         data%ground%gline_flux
#coordinates:  lon f

[rho_ice]
dimensions:    time
units:         kg/meter3
long_name:     ice density
data:          rhoi
factor:        noscale
standard_name: rho_ice

[rho_seawater]
dimensions:    time
units:         kg/meter3
long_name:     seawater density
data:          rhoo
factor:        noscale
standard_name: rho_seawater

[gravity]
dimensions:    time
units:         meter/s/s
long_name:     gravitational acceleration
data:          grav
factor:        noscale
standard_name: gravity

[seconds_per_year]
dimensions:    time
units:         s/yr
long_name:     seconds per year
data:          scyr
factor:        noscale
standard_name: seconds_per_year

[ice_specific_heat]
dimensions:    time
units:         J/kg/K
long_name:     ice specific heat 
data:          shci
factor:        noscale
standard_name: ice_specific_heat

[ice_thermal_conductivity]
dimensions:    time
units:         J/(K kg)
long_name:     ice thermal conductivity 
data:          coni
factor:        noscale
standard_name: ice_thermal_conductivity

[dissip]
dimensions:   time, level, y1, x1
units:        deg C/yr
long_name:    dissipation rate (W m-3) divided by rhoi Ci
data:         data%temper%dissip(up,:,:)
factor:       scyr
load:         1
coordinates:  lon lat

[dissipstag]
dimensions:   time, staglevel, y1, x1
units:        deg C/yr
long_name:    dissipation rate (W m-3) divided by rhoi Ci
data:         data%temper%dissip(up,:,:)
factor:       scyr
load:         1
coordinates:  lon lat

[adv_cfl_dt]
dimensions:    time
units:         years
long_name:     advective CFL maximum time step
data:          data%numerics%adv_cfl_dt

[diff_cfl_dt]
dimensions:    time
units:         years
long_name:     diffusive CFL maximum time step
data:          data%numerics%diff_cfl_dt