/[MITgcm]/MITgcm/verification/tidal_basin_2d/code/external_forcing.F
ViewVC logotype

Diff of /MITgcm/verification/tidal_basin_2d/code/external_forcing.F

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

revision 1.1 by adcroft, Thu Jan 30 18:33:57 2003 UTC revision 1.2 by adcroft, Tue Mar 25 15:24:25 2003 UTC
# Line 74  c     tidal_freq=2.*pi/(43200.+25.*60.) Line 74  c     tidal_freq=2.*pi/(43200.+25.*60.)
74  C But for convenience we are using 12 hour period  C But for convenience we are using 12 hour period
75        tidal_freq=2.*pi/(43200.)        tidal_freq=2.*pi/(43200.)
76  C Make the tide relatively strong (about 1 m)  C Make the tide relatively strong (about 1 m)
77        tidal_Hscale=100.        tidal_Hscale=10.
78        IF ( usingCartesianGrid ) THEN        IF ( usingCartesianGrid ) THEN
79         Coord2longitude=1./rSphere         Coord2longitude=1./rSphere
80        ELSEIF ( usingSphericalPolarGrid ) THEN        ELSEIF ( usingSphericalPolarGrid ) THEN
# Line 88  C Make the tide relatively strong (about Line 88  C Make the tide relatively strong (about
88           longitud2=XC(i,j,bi,bj)*Coord2longitude           longitud2=XC(i,j,bi,bj)*Coord2longitude
89           gU(i,j,kLev,bi,bj) = gU(i,j,kLev,bi,bj)           gU(i,j,kLev,bi,bj) = gU(i,j,kLev,bi,bj)
90       &   +gravity*tidal_Hscale*       &   +gravity*tidal_Hscale*
91       &    ( SIN( tidal_freq*myCurrentTime +longitud2 )       &    ( SIN( tidal_freq*myCurrentTime + 2.*longitud2 )
92       &     -SIN( tidal_freq*myCurrentTime +longitud1 )       &     -SIN( tidal_freq*myCurrentTime + 2.*longitud1 )
93       &    )*recip_DXC(i,j,bi,bj)       &    )*recip_DXC(i,j,bi,bj)
94       &   *_maskW(i,j,kLev,bi,bj)       &   *_maskW(i,j,kLev,bi,bj)
95         &   *min( myCurrentTime/86400. , 1.)
96         ENDDO         ENDDO
97        ENDDO        ENDDO
98    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22