/[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.14 by adcroft, Wed Dec 9 16:11:52 1998 UTC revision 1.15 by adcroft, Tue Dec 15 00:20:34 1998 UTC
# Line 26  C     == Local variables == Line 26  C     == Local variables ==
26  C     iG, jG - Global coordinate index  C     iG, jG - Global coordinate index
27  C     bi,bj  - Loop counters  C     bi,bj  - Loop counters
28  C     I,J,K  C     I,J,K
 C     distY  - Distance accross domain of point j  
 C     lY     - Width of the basin ( last point is land )  
 C     tauMax - Peak zonal wind stress ( N/m^2 )  
 C     tauX   - Local zonal wind stress  
29        INTEGER iG, jG        INTEGER iG, jG
30        INTEGER bi, bj        INTEGER bi, bj
31        INTEGER  I,  J, K        INTEGER  I,  J, K
       _RL distY  
       _RL lY  
       _RL tauX, tauMax  
32    
33        _BARRIER        _BARRIER
34    
# Line 43  C     tauX   - Local zonal wind stress Line 36  C     tauX   - Local zonal wind stress
36         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
37          DO j=1-OLy,sNy+OLy          DO j=1-OLy,sNy+OLy
38           DO i=1-OLx,sNx+OLx           DO i=1-OLx,sNx+OLx
39            fu(i,j,bi,bj) = 0. _d 0            fu   (i,j,bi,bj) = 0. _d 0
40            fv(i,j,bi,bj) = 0. _d 0            fv   (i,j,bi,bj) = 0. _d 0
41              Qnet (i,j,bi,bj) = 0. _d 0
42              EmPmR(i,j,bi,bj) = 0. _d 0
43              SST  (i,j,bi,bj) = 0. _d 0
44              SSS  (i,j,bi,bj) = 0. _d 0
45           ENDDO           ENDDO
46          ENDDO          ENDDO
47         ENDDO         ENDDO
48        ENDDO        ENDDO
   
 C--   Initialise surface bc arrays  
       IF ( zonalWindFile .EQ. ' ' ) THEN  
 C      In cartesian yc, delY and ly are meters.  
 C      In spherical polar yc, delY and ly are degrees  
        tauMax = 0. _d 0  
        lY = 0. _d 0  
        DO j=2,nY-1  
         lY = lY + delY(j)  
        ENDDO  
        DO bj = myByLo(myThid), myByHi(myThid)  
         DO bi = myBxLo(myThid), myBxHi(myThid)  
          DO j=1,sNy  
           DO i=1,sNx  
            distY = (yC(i,j,bi,bj)-(yC0+delY(1)))/lY  
 C          tauX  = -tauMax*cos(2. _d 0*PI*distY)  
            tauX  = tauMax*sin(PI*distY)  
            fu(i,j,bi,bj) = tauX/(delZ(1)*rhonil)  
           ENDDO  
          ENDDO  
         ENDDO  
        ENDDO  
 C      fu(4,4,1,1) = fu(4,4,1,1)*0.917d0  
       ENDIF  
49  C  C
50        _BARRIER        _EXCH_XY_R4(fu   , myThid )
51          _EXCH_XY_R4(fv   , myThid )
52          _EXCH_XY_R4(Qnet , myThid )
53          _EXCH_XY_R4(EmPmR, myThid )
54          _EXCH_XY_R4(SST  , myThid )
55          _EXCH_XY_R4(SSS  , myThid )
56    
57    C     CALL PLOT_FIELD_XYRS( fu, 'S/R INI_FORCING FU',1,myThid)
58    C     CALL PLOT_FIELD_XYRS( fv, 'S/R INI_FORCING FV',1,myThid)
59    
60        IF ( meridWindFile .EQ. ' ' ) THEN        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-OLy,sNy+OLy
63           DO j=1-OLy,sNy+OLy           DO i=1-OLx,sNx+OLx
64            DO i=1-OLx,sNx+OLx            fudot   (i,j,bi,bj) = 0. _d 0
65             fv(i,j,bi,bj) = 0./(delZ(1)*rhonil)            fvdot   (i,j,bi,bj) = 0. _d 0
66            ENDDO            Qnetdot (i,j,bi,bj) = 0. _d 0
67              EmPmRdot(i,j,bi,bj) = 0. _d 0
68              SSTdot  (i,j,bi,bj) = 0. _d 0
69              SSSdot  (i,j,bi,bj) = 0. _d 0
70           ENDDO           ENDDO
71          ENDDO          ENDDO
72         ENDDO         ENDDO
73        ENDIF        ENDDO
74  C  C
75        _EXCH_XY_R4(fu  , myThid )        _EXCH_XY_R4(fudot   , myThid )
76        _EXCH_XY_R4(fv  , myThid )        _EXCH_XY_R4(fvdot   , myThid )
77          _EXCH_XY_R4(Qnetdot , myThid )
78        CALL PLOT_FIELD_XYRS( fu, 'S/R INI_FORCING FU',1,myThid)        _EXCH_XY_R4(EmPmRdot, myThid )
79        CALL PLOT_FIELD_XYRS( fv, 'S/R INI_FORCING FV',1,myThid)        _EXCH_XY_R4(SSTdot  , myThid )
80          _EXCH_XY_R4(SSSdot  , myThid )
81    
82        RETURN        RETURN
83        END        END
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.22