igsm/src/ocean4atm.F


#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
1
2	#include "ctrparam.h"
3	subroutine ocean4atm
4	#include "BD2G04.COM"
5	#include "AGRID.COM"
6
7	do j=1,jm0
8	ODATA(1,j,1)=mmsst(j)
9	ODATA(1,j,3)=mmicem(j)
10	ODATA(1,j,2)=mmfice(j)
11	GDATA(1,j,1)=mmsnowm(j)
12	GDATA(1,j,3)=mmtice(j)
13	GDATA(1,j,7)=0.5*(mmtice2(j)+mmtice1(j))
14	enddo
15
16	C print ,'**OCEAN4ATM',mmsst,mmfice,mmtice
17	CJRS
18	C ODATA(1,1,1)=ODATA(1,2,1)
19	C ODATA(1,1,3)=ODATA(1,2,3)
20	C ODATA(1,1,2)=ODATA(1,2,2)
21	C GDATA(1,1,1)=GDATA(1,2,1)
22	C GDATA(1,1,3)=GDATA(1,2,3)
23	C GDATA(1,1,7)=GDATA(1,2,7)
24
25	C ODATA(1,jm0,1)=ODATA(1,jm0-1,1)
26	C ODATA(1,jm0,3)=ODATA(1,jm0-1,3)
27	C ODATA(1,jm0,2)=ODATA(1,jm0-1,2)
28	C GDATA(1,jm0,1)=GDATA(1,jm0-1,1)
29	C GDATA(1,jm0,3)=GDATA(1,jm0-1,3)
30	C GDATA(1,jm0,7)=GDATA(1,jm0-1,7)
31
32	return
33	end