! file: THC_4box_with_seaice.h   ( -*- mode: f90 -*- )
character*32 run_mode  &  ! either 'ocn', 'atm' or 'cpl'
     ,FW_parameterization  ! either 'dq_dy_bndry', 'q_target' or 'marotzke'
integer nbegin,nend &
     ,ncoupling_periods,ocn_seaice_steps_per_coupling &
     ,atm_steps_per_coupling_period &
     ,ntime_steps_ocn,ntime_steps_atm,ntime_steps_seaice &
     ,ntime_levels,nboxes_ocn,nboxes_atm,nsurface_boxes_ocn,ntracers_ocn
      
parameter(ntime_levels=3,nboxes_ocn=4,nboxes_atm=2,ntracers_ocn=2 &
     ,nsurface_boxes_ocn=2)
      
real*8 S0,sigma_SB
      
! ocean variables:
real*8 time,tbegin,tend &
     ,robert_filter_coeff &
     ,coupling_period,dt_atm,dt_atm_box,dt_ocn,dt_seaice &
     ,tprint,tprint_short &
     ,g,rho_0,C_u,Cp_ocn &
     ,amplitude_noise_IC(ntracers_ocn) &
     ,amplitude_noise_EmP &
     ,u_thermohaline,u_thermohaline_from_restart &
     ,alpha,beta,Sv,year &
     ,u12,u23,u34,u41 &
     ,kappa_H_srfc,kappa_H_deep,kappa_V_south,kappa_V_north &
     ,Width_ocn,DepthUpper_ocn,DepthLower_ocn &
     ,LengthNorth_ocn,LengthSouth_ocn &
     ,T(ntime_levels,nboxes_ocn,ntracers_ocn), rho(nboxes_ocn) &
     ,T_restart(ntime_levels,nboxes_ocn,ntracers_ocn) &
     ,T_star(nsurface_boxes_ocn,ntracers_ocn) &
     ,delta_rho_shutoff &
     ,air_to_sea_flux(nsurface_boxes_ocn,ntracers_ocn) &
     ,air_to_sea_flux_noise(nsurface_boxes_ocn,ntracers_ocn) &
     ,air_to_sea_flux_from_restart(nsurface_boxes_ocn,ntracers_ocn) &
     ,gamma(ntracers_ocn),gamma_equiv,t_star_equiv(nsurface_boxes_ocn) &
     ,factor_EmP &
     ,L_ocn(nboxes_ocn),D_ocn(nboxes_ocn),W_ocn(nboxes_ocn),V_ocn(nboxes_ocn) &
     ,A_ocn(nboxes_ocn)

! atmospheric variables
integer ny
parameter (ny=25)
real*8 f12_heat(nboxes_atm), f12_water &
     ,incoming_radiation(nboxes_atm),longwave_outgoing(nboxes_atm) &
     ,T_atm(ntime_levels,nboxes_atm) &
     ,L_atm(nboxes_atm),W_atm(nboxes_atm),A_atm(nboxes_atm) &
     ,albedo(nboxes_atm) &
     ,frac_land(nboxes_atm),emissivity(nboxes_atm) &
     ,albedo_clouds(nboxes_atm),albedo_land(nboxes_atm) &
     ,albedo_seaice,albedo_ocean &
     ,solar_incoming_radiation(nboxes_atm) &
     ,K_heat_atm &
     ,K_water_atm_dq_dy_bndry,K_water_atm_q_target,K_water_atm_marotzke &
     ,C_bulk_air_sea_heat_flux &
     ,Cp_atm,convert_to_Pw(nboxes_atm)

! seaice variables:
real*8 T_freezing,ice_fusion_heat,rho_ice,seaice_growth_thickness &
     ,shotrwave_penetration_to_seaice &
     ,V_seaice(ntime_levels,nsurface_boxes_ocn),frac_seaice(nsurface_boxes_ocn) &
     ,ocean_to_seaice_flux(1:nsurface_boxes_ocn,1:ntracers_ocn) &
     ,shortwave_to_seaice_heat_flux(1:nsurface_boxes_ocn) &
     ,thickness_seaice(1:nsurface_boxes_ocn) &
     ,seaice_restoring_time

logical restoring_conditions,use_mixed_bc_and_equiv_t &
     ,write_restart_data,read_restart_data,linear_eqn_state &
     ,stable_stratification_stop_THC,no_reversed_THC &
     ,print_atm_in_Pwatts,add_noise_to_surface_EmP

parameter(Sv=1.e6,year=3600.0*24.0*365.0)

common /Ccoupled/ &
     time,tbegin,tend &
     ,coupling_period,dt_atm,dt_atm_box,dt_ocn,dt_seaice &
     ,tprint,tprint_short &
     ,robert_filter_coeff &
     ! integer:
     ,ncoupling_periods,ocn_seaice_steps_per_coupling &
     ,atm_steps_per_coupling_period &
     ,ntime_steps_ocn,ntime_steps_atm,ntime_steps_seaice

common /Cocn/ &
     ! real:
     T,T_restart,rho,air_to_sea_flux,air_to_sea_flux_noise &
     ,air_to_sea_flux_from_restart &
     ,u12,u23,u34,u41 &
     ,kappa_H_srfc,kappa_H_deep,kappa_V_south,kappa_V_north &
     ,Width_ocn,DepthUpper_ocn,DepthLower_ocn &
     ,LengthNorth_ocn,LengthSouth_ocn &
     ,L_ocn,D_ocn,W_ocn,V_ocn,A_ocn &
     ,g,rho_0,alpha,beta,S0 &
     ,T_star,gamma,factor_EmP,C_u,Cp_ocn &
     ,u_thermohaline,u_thermohaline_from_restart &
     ,delta_rho_shutoff &
     ,gamma_equiv,t_star_equiv &
     ,amplitude_noise_IC,amplitude_noise_EmP &
     ! logicals: 
     ,restoring_conditions,use_mixed_bc_and_equiv_t &
     ,write_restart_data,read_restart_data,linear_eqn_state &
     ,stable_stratification_stop_THC,no_reversed_THC &
     ,add_noise_to_surface_EmP &
     ! character
     ,run_mode

common/Catm/ f12_heat, f12_water &
     ,incoming_radiation,longwave_outgoing &
     ,T_atm &
     ,L_atm,W_atm,A_atm &
     ,sigma_SB,albedo &
     ,frac_land,emissivity &
     ,albedo_clouds,albedo_land &
     ,albedo_seaice,albedo_ocean &
     ,solar_incoming_radiation &
     ,K_heat_atm &
     ,K_water_atm_dq_dy_bndry,K_water_atm_q_target,K_water_atm_marotzke &
     ,C_bulk_air_sea_heat_flux &
     ,Cp_atm,convert_to_Pw &
     ,print_atm_in_Pwatts,FW_parameterization

common/Cseaice/ T_freezing,ice_fusion_heat,rho_ice,seaice_growth_thickness &
     ,frac_seaice,V_seaice &
     ,shortwave_to_seaice_heat_flux,ocean_to_seaice_flux &
     ,thickness_seaice, shotrwave_penetration_to_seaice &
     ,seaice_restoring_time

namelist /nl_coupled/ &
      coupling_period,dt_atm_box,dt_ocn,dt_seaice &
     ,tbegin,tend &
     ,tprint,tprint_short &
     ,write_restart_data,read_restart_data &
     ,restoring_conditions,use_mixed_bc_and_equiv_t &
     ,run_mode,add_noise_to_surface_EmP

namelist /nl_ocn/ &
     Width_ocn,DepthUpper_ocn,DepthLower_ocn &
     ,LengthNorth_ocn,LengthSouth_ocn &
     ,g,rho_0,S0,Cp_ocn,gamma,gamma_equiv,factor_EmP &
     ,kappa_H_srfc,kappa_H_deep,kappa_V_south,kappa_V_north &
     ,linear_eqn_state &
     ,stable_stratification_stop_THC,no_reversed_THC &
     ,delta_rho_shutoff &
     ,robert_filter_coeff &
     ,amplitude_noise_IC,amplitude_noise_EmP,C_u


namelist /nl_atm/ &
     frac_land,emissivity &
     ,albedo_clouds,albedo_land &
     ,albedo_seaice,albedo_ocean &
     ,solar_incoming_radiation &
     ,K_heat_atm &
     ,K_water_atm_dq_dy_bndry,K_water_atm_q_target,K_water_atm_marotzke &
     ,C_bulk_air_sea_heat_flux &
     ,sigma_SB,Cp_atm,FW_parameterization &
     ,print_atm_in_Pwatts

namelist /nl_seaice/ &
     T_freezing,ice_fusion_heat,seaice_growth_thickness,rho_ice &
     ,shotrwave_penetration_to_seaice,seaice_restoring_time