41 |
#ifdef ALLOW_CFC |
#ifdef ALLOW_CFC |
42 |
C == Local variables == |
C == Local variables == |
43 |
_RL SURCFC11(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL SURCFC11(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
44 |
|
_RL xintp(1-OLy:sNy+OLy) |
45 |
INTEGER I,J |
INTEGER I,J |
46 |
INTEGER myYear, lastYear, thisYear |
INTEGER myYear, lastYear, thisYear |
47 |
_RL dtinc, aWght, bWght |
_RL dtinc, aWght, bWght |
48 |
_RL ACFC11north, ACFC11south |
_RL ACFC11north, ACFC11south |
49 |
_RL maxYear |
_RL maxYear |
50 |
_RL a1, a2 |
_RL a1, a2 |
51 |
|
_RL yNorth, ySouth |
52 |
|
|
53 |
|
|
54 |
DO j=1-OLy,sNy+OLy |
DO j=1-OLy,sNy+OLy |
55 |
DO i=1-OLx,sNx+OLx |
DO i=1-OLx,sNx+OLx |
79 |
endif |
endif |
80 |
c print*,'ACFC11north,ACFC11south', ACFC11north,ACFC11south, |
c print*,'ACFC11north,ACFC11south', ACFC11north,ACFC11south, |
81 |
c & lastYear,thisYear |
c & lastYear,thisYear |
82 |
|
c provide gradient between N and S values |
83 |
|
#define OCMIP_GRAD |
84 |
|
#undef STEPH_GRAD |
85 |
|
c STEPH'S INITIAL VERSION |
86 |
|
#ifdef STEPH_GRAD |
87 |
DO j=1-OLy,sNy+OLy |
DO j=1-OLy,sNy+OLy |
88 |
DO i=1-OLx,sNx+OLx |
DO i=1-OLx,sNx+OLx |
89 |
if ((j.gt.int(sNy/2)+3.and.j.le.sNy).or.j.lt.1) then |
if ((j.gt.int(sNy/2)+3.and.j.le.sNy).or.j.lt.1) then |
100 |
endif |
endif |
101 |
ENDDO |
ENDDO |
102 |
ENDDO |
ENDDO |
103 |
|
#endif |
104 |
|
c OCMIP VERSION |
105 |
|
#ifdef OCMIP_GRAD |
106 |
|
yNorth = 10.0 |
107 |
|
ySouth = -10.0 |
108 |
|
DO j=1-OLy,sNy+OLy |
109 |
|
i=1 |
110 |
|
IF(yC(i,j,bi,bj) .GE. yNorth) THEN |
111 |
|
xintp(j) = 1.0 |
112 |
|
ELSE IF(yC(i,j,bi,bj) .LE. ySouth) THEN |
113 |
|
xintp(j) = 0.0 |
114 |
|
ELSE |
115 |
|
xintp(j) = (yC(i,j,bi,bj) - ySouth)/ |
116 |
|
& (yNorth - ySouth) |
117 |
|
ENDIF |
118 |
|
DO i=1-OLx,sNx+OLx |
119 |
|
ATMOSCFC11(i,j,bi,bj)= xintp(j) * ACFC11north |
120 |
|
& + (1.0 - xintp(j))*ACFC11south |
121 |
|
|
122 |
|
ENDDO |
123 |
|
c print*,'QQ cfc11', j, ATMOSCFC11(1,j,bi,bj) |
124 |
|
ENDDO |
125 |
|
#endif |
126 |
c cfc11 air-sea interaction |
c cfc11 air-sea interaction |
127 |
CALL CFC11_SURFFORCING( PTR_CFC11, SURCFC11, |
CALL CFC11_SURFFORCING( PTR_CFC11, SURCFC11, |
128 |
& bi,bj,imin,imax,jmin,jmax, |
& bi,bj,imin,imax,jmin,jmax, |