/[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.12 by cnh, Fri Nov 6 22:44:47 1998 UTC revision 1.13 by adcroft, Tue Dec 8 18:03:21 1998 UTC
# Line 53  C--   Initialise surface bc arrays Line 53  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.
55  C      In spherical polar yc, delY and ly are degrees  C      In spherical polar yc, delY and ly are degrees
56         tauMax = 0.1 _d 0         tauMax = 0. _d 0
57         lY = 0. _d 0         lY = 0. _d 0
58         DO j=1,nY-1         DO j=2,nY-1
59          lY = lY + delY(j)          lY = lY + delY(j)
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)
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+delY(1)))/lY
66  C          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)
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22