/[MITgcm]/MITgcm/model/src/ini_forcing.F
ViewVC logotype

Diff of /MITgcm/model/src/ini_forcing.F

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

revision 1.5 by cnh, Mon Jun 8 21:43:01 1998 UTC revision 1.6 by cnh, Tue Jun 9 18:57:38 1998 UTC
# Line 36  C     tauX   - Local zonal wind stress Line 36  C     tauX   - Local zonal wind stress
36        _RL lY        _RL lY
37        _RL tauX, tauMax        _RL tauX, tauMax
38    
39          _BARRIER
40    
41  C--   Initialise surface bc arrays  C--   Initialise surface bc arrays
42  C     In cartesian yc, delY and ly are meters.        IF ( zonalWindFile .EQ. ' ' ) THEN
43  C     In spherical polar yc, delY and ly are degrees  C      In cartesian yc, delY and ly are meters.
44        tauMax = 0.1 _d 0  C      In spherical polar yc, delY and ly are degrees
45        lY = 0. _d 0         tauMax = 0.1 _d 0
46        DO j=1,nY-1         lY = 0. _d 0
47         lY = lY + delY(j)         DO j=1,nY-1
48        ENDDO          lY = lY + delY(j)
49        DO bj = myByLo(myThid), myByHi(myThid)         ENDDO
50         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
51          DO j=1-OLy,sNy+OLy          DO bi = myBxLo(myThid), myBxHi(myThid)
52           DO i=1-OLx,sNx+OLx           DO j=1-OLy,sNy+OLy
53            fu(i,j,bi,bj) = 0. _d 0            DO i=1-OLx,sNx+OLx
54            fv(i,j,bi,bj) = 0. _d 0             fu(i,j,bi,bj) = 0. _d 0
55               fv(i,j,bi,bj) = 0. _d 0
56              ENDDO
57           ENDDO           ENDDO
58          ENDDO          ENDDO
59         ENDDO         ENDDO
60        ENDDO         DO bj = myByLo(myThid), myByHi(myThid)
61        DO bj = myByLo(myThid), myByHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
62         DO bi = myBxLo(myThid), myBxHi(myThid)           DO j=1,sNy
63          DO j=1,sNy            DO i=1,sNx
64           DO i=1,sNx             distY = (yC(i,j,bi,bj)-yC0)/lY
65            distY = (yC(i,j,bi,bj)-yC0)/lY             tauX  = -tauMax*cos(2. _d 0*PI*distY)
66            tauX  = -tauMax*cos(2. _d 0*PI*distY)             tauX  = tauMax*sin(PI*distY)
67            tauX  = tauMax*sin(PI*distY)             fu(i,j,bi,bj) = tauX/(delZ(1)*rhonil)
68            fu(i,j,bi,bj) = tauX/(delZ(1)*rhonil)  C          fu(i,j,bi,bj) = -tauMax
69  C         fu(i,j,bi,bj) = -tauMax            ENDDO
70             ENDDO
71            ENDDO
72           ENDDO
73           fu(4,4,1,1) = fu(4,4,1,1)*0.917d0
74          ELSE
75           _BEGIN_MASTER(myThid)
76           CALL READ_FLD_XY_RS( zonalWindFile, ' ', fu, 0, myThid )
77           _END_MASTER(myThid)
78          ENDIF
79    C
80          _BARRIER
81    
82          IF ( meridWindFile .EQ. ' ' ) THEN
83           DO bj = myByLo(myThid), myByHi(myThid)
84            DO bi = myBxLo(myThid), myBxHi(myThid)
85             DO j=1-OLy,sNy+OLy
86              DO i=1-OLx,sNx+OLx
87               fv(i,j,bi,bj) = 0. _d 0
88              ENDDO
89           ENDDO           ENDDO
90          ENDDO          ENDDO
91         ENDDO         ENDDO
92        ENDDO        ELSE
93        fu(4,4,1,1) = fu(4,4,1,1)*0.917d0         _BEGIN_MASTER(myThid)
94           CALL READ_FLD_XY_RS( meridWindFile, ' ', fv, 0, myThid )
95           _END_MASTER(myThid)
96          ENDIF
97  C  C
98        _EXCH_XY_R4(fu  , myThid )        _EXCH_XY_R4(fu  , myThid )
99        _EXCH_XY_R4(fv  , myThid )        _EXCH_XY_R4(fv  , myThid )
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22