REAL*8 hfluxo(jm0)    ! ocean surf. heat flux (W/m2) (+=downward)

       REAL*8 hfluxi(jm0)    ! seaice surf. heat flux (W/m2) (+=downward)

       REAL*8 dhfodtg(jm0)   ! dF/dT ocean, consistent with above units (-)

       REAL*8 dhfodtgeq(jm0) ! dF/dT ocean, alternative formulation

       REAL*8 dhfidtg(jm0)   ! dF/dT seaice

       REAL*8 dhfidtgeq(jm0) ! dF/dT seaice alternative formulation

       REAL*8 precip(jm0)    ! precip (mm/day) (+=into ocean)

       REAL*8 evao(jm0)      ! ocean evap. (mm/day) (+=out of ocean)

       REAL*8 evai(jm0)      ! seaice evap. (mm/day) (+=out of ocean)

       REAL*8 devodtg(jm0)   ! dL/dT ocean, consistent with above units (+)

       REAL*8 devodtgeq(jm0) ! dL/dT ocean, alternative formulation

       REAL*8 devidtg(jm0)   ! dL/dT ice

       REAL*8 devidtgeq(jm0) ! dL/dT ice, alternative formulation

       REAL*8 tauu(jm0)      ! zonal mom flux at lower boundary (N/m2)

       REAL*8 tauv(jm0)      ! merid. mom flux at lower boundary (N/m2)

       REAL*8 arunoff(jm0)   ! runoff (kg/day) (+=into ocean)

       REAL*8 solarinc_ice(jm0)    ! solar incoming (+=into ocean)

       REAL*8 solarnet_ice(jm0)    ! net solar radation to ice (not used)

       REAL*8 solarinc_ocean(jm0)  ! solar incoming (+=into ocean)(unused)

       REAL*8 solarnet_ocean(jm0)  ! net radiation solar to ocean

       REAL*8 wsocean(jm0)   ! windspeed at ocean surface (m/s)

       REAL*8 co24ocean(jm0) ! atmospheric CO2  (ppm)

       REAL*8 ps4ocean(jm0)  ! surf pres (mb), to be normalized to sea level

       REAL*8 tempr(jm0)     ! precipitation temp (used over seaice) (C)

       REAL*8 qyz4ocean(jm0,lm0)  ! atmos mixing ratio structure (kg/kg)

       REAL*8 tyz4ocean(jm0,lm0)  ! atmos temp structure (K)

       REAL*8 mmsst(jm0)    ! zonal mean SST (C)

       REAL*8 mmfice(jm0)   ! fraction of ocean area with seaice cover

       REAL*8 mmtice(jm0)   ! zonal mean seaice TSurf (C)

       REAL*8 mmco2flux(jm0)! total ocean->atmos CO2 flux (mol/s for latitude band)

       REAL*8 cflan(jm0)    ! fraction of land in latitude circle

       REAL*8 mmsAlb(jm0)   ! zonal mean visual seaice albedo

       REAL*8 mmsAlbNIR(jm0)   ! zonal mean near IR seaice albedo

       REAL*8 mmicem(jm0)

       REAL*8 mmtice1(jm0)

       REAL*8 mmtice2(jm0)

       REAL*8 mmsnowm(jm0)

      common/OCEAN2ATM/osst(jm0),aoice(jm0),foice(jm0),

       REAL*4 atm_yC(jm0)    ! center points of atmos grid

     & snowice(jm0),tice1(jm0),tice2(jm0)

      common/mltemp/osst2(jm0),osst3(jm0),rseaice(jm0)

     & ,CLAND4OCEAN(io0,jm0)

       REAL*4 atm_yG(jm0)    ! southern boundary of atmos grid cells

       _RL atm_dFdt_ocnq(jm0) ! dF/dT ocean, alternate form (W/m2/K) (+)

       _RL atm_solarinc(jm0)  ! solar incoming (W/m2 +=into ocean)

       _RL atm_solar_ice(jm0) ! net solar rad to ice (unused)(+=into ocean)

       _RL atm_solar_ocn(jm0) ! solar incoming to ocean (+=into ocean)

       _RL atm_windspeed(jm0) ! windspeed at ocean surface (m/s)

       _RL atm_slp(jm0)       ! SLP anomaly (mb)

       _RL atm_pco2(jm0)      ! atmospheric pCO2 (ppmv)

       _RL ctocn(jm0)    ! zonal mean ice-free ocean temp

       _RL ctice(jm0)    ! zonal mean seaice surface temp

       _RL cfice(jm0)    ! zonal mean ice fraction

       _RL csAlb(jm0)    ! zonal mean seaice albedo

       _RL csAlbNIR(jm0)    ! zonal mean seaice albedo

       _RL ocnArea(jm0)  ! ocean area of latitude strip on atm grid (m2)

       _RL cco2flux(jm0) ! zonally integrated flux of CO2 from ocean->atm (mol/m2/s)

       _RL atau(jm0,nForcingPer) ! zonal wind stress

       _RL atav(jm0,nForcingPer) ! meridional wind stress

       _RL awind(jm0,nForcingPer)! (total) wind speed

       _RL sum_tauu_ta(jm0,nForcingPer)

       _RL sum_tauv_ta(jm0,nForcingPer)

       _RL sum_wsocean_ta(jm0,nForcingPer)

       _RL sum_ps4ocean_ta(jm0,nForcingPer)

       _RL atm_tauu(jm0)      ! zonal mom flux at lower boundary (N/m2)

       _RL atm_tauv(jm0)      ! merid. mom flux at lower boundary (N/m2)

       _RL atm_Tair(jm0)      ! precipitation temp (used over seaice) (C)

       _RL atm_precip(jm0)    ! precip (m/s) (+=out of ocean, - definite)

       _RL atm_runoff(jm0)    ! runoff (m3/s) (+=out of ocean, normally - )

       _RL atm_evap_ice(jm0)  ! seaice evap. (m/s) (-=out of ocean, normally -)

       _RL atm_evap_ocn(jm0)  ! ocean evap. (m/s) (-=out of ocean, normally -)

       _RL atm_qnet_ice(jm0)  ! seaice surf. heat flux (W/m2) (+=upward)

       _RL atm_qnet_ocn(jm0)  ! ocean surf. heat flux (W/m2) (+=upward)

       _RL atm_dLdt_ice(jm0)  ! dL/dT ice (m/s/K) (-)

       _RL atm_dLdt_ocn(jm0)  ! dL/dT ocean (m/s/K) (-)

       _RL atm_dFdt_ice(jm0)  ! dF/dT ice (W/m2/K) (+)

       _RL atm_dFdt_ocn(jm0)  ! dF/dT ocean (W/m2/K) (+)

       _RL atm_dLdt_iceq(jm0) ! dL/dT ice, alternate form (m/s/K) (-)

       _RL atm_dLdt_ocnq(jm0) ! dL/dT ocean, alternate form (m/s/K) (-)

       _RL atm_dFdt_iceq(jm0) ! dF/dT ice, alternate form (W/m2/K) (+)

      INTEGER jm0,lm0

      PARAMETER (jm0=N_LAT,lm0=N_LEV)

       REAL*8 pC_atm(jm0)

       REAL*8 fluxco2(jm0)

      DO j_atm=2,jm0-1

          mmsst(jm0)= ctocn(jm0-1)

          mmfice(jm0)= cfice(jm0-1)

          mmco2flux(jm0)= 0. _d 0 ! converted to mol/s; pole point contribution in jm0-1

        mmtice(jm0)= ctice(jm0-1)

        mmsAlb(jm0)= csAlb(jm0-1)

        mmsAlbNIR(jm0)= csAlbNIR(jm0-1)

      DO j_atm=1,jm0

        DO j=1,jm0

        DO j=1,jm0

        DO j=1,jm0

        DO j=1,jm0

        DO j=1,jm0

        DO j=1,jm0

        DO j=1,jm0

      atm_yG(jm0)=90. _d 0

        DO jj=2,jm0-1

        DO jj=1,jm0-1

      DO j_atm=1,jm0

      DO j_atm=3,jm0-2

      cflan(jm0-1)= 1.d0 - ocnArea(jm0-1) /

     &             (1. _d 0-sin(atm_yG(jm0-1)*deg2rad)))

      IF (cflan(jm0-1).LT.1. _d -14) cflan(jm0-1)=0. _d 0

      cflan(jm0)=cflan(jm0-1)

      DO j_atm=1,jm0

     &        ACCESS='direct', RECL=8*jm0*nForcingPer,

     &        ACCESS='direct', RECL=8*jm0*nForcingPer,

     &        ACCESS='direct', RECL=8*jm0*nForcingPer,

      DO j_atm=1,jm0

      _RS atm_dyG(jm0)  ! southern point/(boundary) of atmos grid

      DO j_atm=2,jm0

      atm_yC(jm0)=atm_yG(jm0)+atm_dyG(jm0)/2.0

      atm_oc_ind(sNy)=jm0-1

     &      (sin((atm_yG(jm0)+atm_dyG(jm0))*deg2rad)-

     &      sin(atm_yG(jm0-1)*deg2rad))

      DO j_atm=1,jm0

      DO j_atm=1,jm0

      IF (cflan(jm0-1).NE.1. _d 0)

     &          CALL COMBINE_ENDS(endwgt1,endwgt2,jm0,jm0-1,rsumwgt)

      DO j_atm=2,jm0-1

        DO j_atm=1,jm0

     &     ACCESS='direct', RECL=8*jm0*nForcingPer, FORM='unformatted')

     &     ACCESS='direct', RECL=8*jm0*nForcingPer, FORM='unformatted')

     &     ACCESS='direct', RECL=8*jm0*nForcingPer, FORM='unformatted')

     &     ACCESS='direct', RECL=8*jm0*nForcingPer, FORM='unformatted')

        DO j_atm=1,jm0