/[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.9 by cnh, Mon Jun 15 05:13:56 1998 UTC revision 1.10 by adcroft, Wed Jun 17 21:07:02 1998 UTC
# Line 38  C     tauX   - Local zonal wind stress Line 38  C     tauX   - Local zonal wind stress
38    
39        _BARRIER        _BARRIER
40    
41          DO bj = myByLo(myThid), myByHi(myThid)
42           DO bi = myBxLo(myThid), myBxHi(myThid)
43            DO j=1-OLy,sNy+OLy
44             DO i=1-OLx,sNx+OLx
45              fu(i,j,bi,bj) = 0. _d 0
46              fv(i,j,bi,bj) = 0. _d 0
47             ENDDO
48            ENDDO
49           ENDDO
50          ENDDO
51    
52  C--   Initialise surface bc arrays  C--   Initialise surface bc arrays
53        IF ( zonalWindFile .EQ. ' ' ) THEN        IF ( zonalWindFile .EQ. ' ' ) THEN
54  C      In cartesian yc, delY and ly are meters.  C      In cartesian yc, delY and ly are meters.
# Line 49  C      In spherical polar yc, delY and l Line 60  C      In spherical polar yc, delY and l
60         ENDDO         ENDDO
61         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
62          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
          DO j=1-OLy,sNy+OLy  
           DO i=1-OLx,sNx+OLx  
            fu(i,j,bi,bj) = 0. _d 0  
            fv(i,j,bi,bj) = 0. _d 0  
           ENDDO  
          ENDDO  
         ENDDO  
        ENDDO  
        DO bj = myByLo(myThid), myByHi(myThid)  
         DO bi = myBxLo(myThid), myBxHi(myThid)  
63           DO j=1,sNy           DO j=1,sNy
64            DO i=1,sNx            DO i=1,sNx
65             distY = (yC(i,j,bi,bj)-yC0)/lY             distY = (yC(i,j,bi,bj)-yC0)/lY
66             tauX  = -tauMax*cos(2. _d 0*PI*distY)  C          tauX  = -tauMax*cos(2. _d 0*PI*distY)
67             tauX  = tauMax*sin(PI*distY)             tauX  = tauMax*sin(PI*distY)
68             fu(i,j,bi,bj) = tauX/(delZ(1)*rhonil)             fu(i,j,bi,bj) = tauX/(delZ(1)*rhonil)
 C          fu(i,j,bi,bj) = -tauMax  
69            ENDDO            ENDDO
70           ENDDO           ENDDO
71          ENDDO          ENDDO
72         ENDDO         ENDDO
73         fu(4,4,1,1) = fu(4,4,1,1)*0.917d0         fu(4,4,1,1) = fu(4,4,1,1)*0.917d0
       ELSE  
        _BEGIN_MASTER(myThid)  
        CALL READ_FLD_XY_RS( zonalWindFile, ' ', fu, 0, myThid )  
        DO bj = myByLo(myThid), myByHi(myThid)  
         DO bi = myBxLo(myThid), myBxHi(myThid)  
          DO j=1,sNy  
           DO i=1,sNx  
            fu(i,j,bi,bj) = fu(i,j,bi,bj)/(delZ(1)*rhonil)*0.1  
            fu(i,j,bi,bj) = 0.1/(delZ(1)*rhonil)  
           ENDDO  
          ENDDO  
         ENDDO  
        ENDDO  
        _END_MASTER(myThid)  
74        ENDIF        ENDIF
75  C  C
76        _BARRIER        _BARRIER
# Line 94  C Line 80  C
80          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
81           DO j=1-OLy,sNy+OLy           DO j=1-OLy,sNy+OLy
82            DO i=1-OLx,sNx+OLx            DO i=1-OLx,sNx+OLx
83             fv(i,j,bi,bj) = 0. _d 0             fv(i,j,bi,bj) = 0./(delZ(1)*rhonil)
           ENDDO  
          ENDDO  
         ENDDO  
        ENDDO  
       ELSE  
        _BEGIN_MASTER(myThid)  
        CALL READ_FLD_XY_RS( meridWindFile, ' ', fv, 0, myThid )  
        DO bj = myByLo(myThid), myByHi(myThid)  
         DO bi = myBxLo(myThid), myBxHi(myThid)  
          DO j=1,sNy  
           DO i=1,sNx  
            fv(i,j,bi,bj) = fv(i,j,bi,bj)/(delZ(1)*rhonil)*0.  
84            ENDDO            ENDDO
85           ENDDO           ENDDO
86          ENDDO          ENDDO
87         ENDDO         ENDDO
        _END_MASTER(myThid)  
88        ENDIF        ENDIF
89  C  C
90        _EXCH_XY_R4(fu  , myThid )        _EXCH_XY_R4(fu  , myThid )
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22