#include "ctrparam.h"
        subroutine ocean4atm
#include "BD2G04.COM"
#include "AGRID.COM"

         do j=1,jm0
           ODATA(1,j,1)=mmsst(j)
           ODATA(1,j,3)=mmicem(j)
           ODATA(1,j,2)=mmfice(j)
           GDATA(1,j,1)=mmsnowm(j)
           GDATA(1,j,3)=mmtice(j)
           GDATA(1,j,7)=0.5*(mmtice2(j)+mmtice1(j))
         enddo

C         print *,'***OCEAN4ATM',mmsst,mmfice,mmtice
CJRS
C           ODATA(1,1,1)=ODATA(1,2,1)
C           ODATA(1,1,3)=ODATA(1,2,3)
C           ODATA(1,1,2)=ODATA(1,2,2)
C           GDATA(1,1,1)=GDATA(1,2,1)
C           GDATA(1,1,3)=GDATA(1,2,3)
C           GDATA(1,1,7)=GDATA(1,2,7)

C           ODATA(1,jm0,1)=ODATA(1,jm0-1,1)
C           ODATA(1,jm0,3)=ODATA(1,jm0-1,3)
C           ODATA(1,jm0,2)=ODATA(1,jm0-1,2)
C           GDATA(1,jm0,1)=GDATA(1,jm0-1,1)
C           GDATA(1,jm0,3)=GDATA(1,jm0-1,3)
C           GDATA(1,jm0,7)=GDATA(1,jm0-1,7)

       return
       end