1 |
jscott |
1.1 |
|
2 |
|
|
#include "ctrparam.h" |
3 |
|
|
subroutine ocean4atm |
4 |
|
|
#include "BD2G04.COM" |
5 |
jscott |
1.2 |
#include "AGRID.h" |
6 |
jscott |
1.1 |
|
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 |