/[MITgcm]/MITgcm/pkg/cfc/cfc11_forcing.F
ViewVC logotype

Diff of /MITgcm/pkg/cfc/cfc11_forcing.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by stephd, Mon Oct 6 20:03:06 2003 UTC revision 1.2 by stephd, Tue May 3 17:02:26 2005 UTC
# Line 41  C     == Routine arguments == Line 41  C     == Routine arguments ==
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
# Line 76  c       &write(0,*) 'myYear = ',myYear,l Line 79  c       &write(0,*) 'myYear = ',myYear,l
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
# Line 92  c    & Line 100  c    &
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,
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22