DO j_atm=2,jm0-1

        IF (ocnArea(j_atm).GT.1. _d -32) THEN

     &      cfice(j_atm)= cfice(j_atm)/ocnArea(j_atm)

          IF (cfice(j_atm).GT.1. _d -32) THEN

            ctice(j_atm)= ctice(j_atm)/ocnArea(j_atm)/cfice(j_atm)

            csAlb(j_atm)= csAlb(j_atm)/ocnArea(j_atm)/cfice(j_atm)

            csAlbNIR(j_atm)= csAlbNIR(j_atm)/ocnArea(j_atm)/cfice(j_atm)

          IF ((1. _d 0-cfice(j_atm).GT.1. _d -32).AND.doAll)

     &        ctocn(j_atm)= ctocn(j_atm)/ocnArea(j_atm)

     &                        /(1. _d 0-cfice(j_atm))

          mmsst(j_atm)= ctocn(j_atm)

          mmfice(j_atm)= cfice(j_atm)

          mmco2flux(j_atm)= cco2flux(j_atm)

        mmtice(j_atm)= ctice(j_atm)

        mmsAlb(j_atm)= csAlb(j_atm)

        mmsAlbNIR(j_atm)= csAlbNIR(j_atm)

      INTEGER j_atm  ! loop counter for the atm grid

      DO j_atm=1,jm0

          ctocn(j_atm)=0. _d 0

          cfice(j_atm)=0. _d 0

          cco2flux(j_atm)=0. _d 0

        ctice(j_atm)=0. _d 0

        csAlb(j_atm)=0. _d 0

        csAlbNIR(j_atm)=0. _d 0

     &         (yG(1,j,1,1).LT.atm_yG(jj+1))) j_atm=jj

        atm_oc_ind(j)=j_atm

        enda1 = sin(atm_yG(j_atm) *deg2rad)

        enda2 = sin(atm_yG(j_atm+1) *deg2rad)

        IF ( yG(1,j+1,1,1) .GT. atm_yG(j_atm+1) ) THEN

           enda3 = sin(atm_yG(j_atm+2) *deg2rad)

     &        ( yG(1,j,1,1) .LT. atm_yC(jj+1))) j_atm=jj

        tauv_jpt(j)= j_atm

        tauv_jwght(j)= 1. _d 0 - (yG(1,j,1,1) - atm_yC(j_atm)) /

     &                 (atm_yC(j_atm+1) - atm_yC(j_atm))

      DO j_atm=1,jm0

        cflan(j_atm)=0. _d 0

        ocnArea(j_atm)=0. _d 0

      DO j_atm=3,jm0-2

        cflan(j_atm)=1. _d 0 - ocnArea(j_atm) /

     &    (sin(atm_yG(j_atm+1)*deg2rad) - sin(atm_yG(j_atm)*deg2rad)))

        if (cflan(j_atm).LT.1. _d -14) cflan(j_atm)=0. _d 0

      DO j_atm=1,jm0

          atau(j_atm,mn)=0. _d 0

          atav(j_atm,mn)=0. _d 0

          awind(j_atm,mn)=0. _d 0

      DO j_atm=1,jm0

          sum_tauu_ta(j_atm,mn)= 0. _d 0

          sum_tauv_ta(j_atm,mn)= 0. _d 0

          sum_wsocean_ta(j_atm,mn)= 0. _d 0

          sum_ps4ocean_ta(j_atm,mn)= 0. _d 0

      INTEGER j_atm, mn

      DO j_atm=2,jm0

        atm_yG(j_atm)=atm_yG(j_atm-1)+atm_dyG(j_atm-1)

        atm_yC(j_atm-1)=(atm_yG(j_atm-1)+atm_yG(j_atm))/2.0

      INTEGER j_atm

      DO j_atm=1,jm0

        sum_tauu_ta(j_atm,inMonth)= sum_tauu_ta(j_atm,inMonth) +

     &                           tauu(j_atm)*dtatmo

        sum_tauv_ta(j_atm,inMonth)= sum_tauv_ta(j_atm,inMonth) +

     &                           tauv(j_atm)*dtatmo

        sum_wsocean_ta(j_atm,inMonth)= sum_wsocean_ta(j_atm,inMonth) +

     &                           wsocean(j_atm)*dtatmo

        sum_ps4ocean_ta(j_atm,inMonth)= sum_ps4ocean_ta(j_atm,inMonth) +

     &                           ps4ocean(j_atm)*dtatmo

      DO j_atm=1,jm0

        atm_tauu(j_atm) = tauu(j_atm)

        atm_tauv(j_atm) = tauv(j_atm)

        atm_tair(j_atm) = tempr(j_atm)

        atm_precip(j_atm)    = -precip(j_atm)/secDay1000

        atm_runoff(j_atm)    = -arunoff(j_atm)/secDay1000

        atm_evap_ice(j_atm)  = -evai(j_atm)/secDay1000

        atm_evap_ocn(j_atm)  = -evao(j_atm)/secDay1000

        atm_qnet_ice(j_atm)  = -hfluxi(j_atm)

        atm_qnet_ocn(j_atm)  = -hfluxo(j_atm)

        atm_dFdt_ice(j_atm)  = -dhfidtg(j_atm)

        atm_dFdt_ocn(j_atm)  = -dhfodtg(j_atm)

        atm_dLdt_ice(j_atm)  = 0. _d 0  ! -devidtg(j_atm)/secDay1000

        atm_dLdt_ocn(j_atm)  = -devodtg(j_atm)/secDay1000

        atm_dFdt_iceq(j_atm) = -dhfidtgeq(j_atm)

        atm_dFdt_ocnq(j_atm) = -dhfodtgeq(j_atm)

        atm_dLdt_iceq(j_atm) = 0. _d 0  ! -devidtgeq(j_atm)/secDay1000

        atm_dLdt_ocnq(j_atm) = -devodtgeq(j_atm)/secDay1000

        atm_solarinc(j_atm)  = solarinc_ice(j_atm)

        atm_solar_ocn(j_atm) = solarnet_ocean(j_atm)

        atm_solar_ice(j_atm) = solarnet_ice(j_atm)

        atm_windspeed(j_atm) = wsocean(j_atm)

        atm_slp(j_atm) = ps4ocean(j_atm)*1013.25/984.0 - 1013.25

        atm_pco2(j_atm) = co24ocean(j_atm)

      INTEGER j_atm

      DO j_atm=2,jm0-1

        atm_tauu(j_atm)= atm_tauu(j_atm) -

     &       (wght0*atau(j_atm,intime0) + wght1*atau(j_atm,intime1))

        atm_tauv(j_atm)= atm_tauv(j_atm) -

     &       (wght0*atav(j_atm,intime0) + wght1*atav(j_atm,intime1))

        atm_windspeed(j_atm)= atm_windspeed(j_atm) -

     &       (wght0*awind(j_atm,intime0) + wght1*awind(j_atm,intime1))

        DO j_atm=1,jm0

          sum_tauu_ta(j_atm,mn) = sum_tauu_ta(j_atm,mn) / norm_factor

          sum_tauv_ta(j_atm,mn) = sum_tauv_ta(j_atm,mn) / norm_factor

          sum_wsocean_ta(j_atm,mn) = sum_wsocean_ta(j_atm,mn) /

          sum_ps4ocean_ta(j_atm,mn) = sum_ps4ocean_ta(j_atm,mn) /

        DO j_atm=1,jm0

          sum_tauu_ta(j_atm,mn) = 0. _d 0

          sum_tauv_ta(j_atm,mn) = 0. _d 0

          sum_wsocean_ta(j_atm,mn) = 0. _d 0

          sum_ps4ocean_ta(j_atm,mn) = 0. _d 0

      INTEGER i,j,mn,j_atm