/[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.1.1.1 by cnh, Wed Apr 22 19:15:30 1998 UTC revision 1.9.2.1 by cnh, Fri Jun 19 04:55:13 1998 UTC
# Line 1  Line 1 
1  C $Id$  C $Header$
2    
3  #include "CPP_EEOPTIONS.h"  #include "CPP_EEOPTIONS.h"
4    
# 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
58            ENDDO
59           ENDDO
60           DO bj = myByLo(myThid), myByHi(myThid)
61            DO bi = myBxLo(myThid), myBxHi(myThid)
62             DO j=1,sNy
63              DO i=1,sNx
64               distY = (yC(i,j,bi,bj)-yC0)/lY
65               tauX  = -tauMax*cos(2. _d 0*PI*distY)
66               tauX  = tauMax*sin(PI*distY)
67               fu(i,j,bi,bj) = tauX/(delZ(1)*rhonil)
68    C          fu(i,j,bi,bj) = -tauMax
69              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           DO bj = myByLo(myThid), myByHi(myThid)
78            DO bi = myBxLo(myThid), myBxHi(myThid)
79             DO j=1,sNy
80              DO i=1,sNx
81               fu(i,j,bi,bj) = fu(i,j,bi,bj)/(delZ(1)*rhonil)
82    C          fu(i,j,bi,bj) = 0.1/(delZ(1)*rhonil)
83              ENDDO
84             ENDDO
85            ENDDO
86           ENDDO
87           _END_MASTER(myThid)
88          ENDIF
89    C
90          _BARRIER
91    
92          IF ( meridWindFile .EQ. ' ' ) THEN
93           DO bj = myByLo(myThid), myByHi(myThid)
94            DO bi = myBxLo(myThid), myBxHi(myThid)
95             DO j=1-OLy,sNy+OLy
96              DO i=1-OLx,sNx+OLx
97               fv(i,j,bi,bj) = 0. _d 0
98              ENDDO
99           ENDDO           ENDDO
100          ENDDO          ENDDO
101         ENDDO         ENDDO
102        ENDDO        ELSE
103        DO bj = myByLo(myThid), myByHi(myThid)         _BEGIN_MASTER(myThid)
104         DO bi = myBxLo(myThid), myBxHi(myThid)         CALL READ_FLD_XY_RS( meridWindFile, ' ', fv, 0, myThid )
105          DO j=1,sNy         DO bj = myByLo(myThid), myByHi(myThid)
106           DO i=1,sNx          DO bi = myBxLo(myThid), myBxHi(myThid)
107            distY = yC(i,j,bi,bj)/lY           DO j=1,sNy
108  C         tauX  = -tauMax*cos(2. _d 0*PI*distY)            DO i=1,sNx
109            tauX  = tauMax*sin(PI*distY)             fv(i,j,bi,bj) = fv(i,j,bi,bj)/(delZ(1)*rhonil)
110            fu(i,j,bi,bj) = tauX/(delZ(1)*rhonil)            ENDDO
 C         fu(i,j,bi,bj) = -tauMax  
111           ENDDO           ENDDO
112          ENDDO          ENDDO
113         ENDDO         ENDDO
114        ENDDO         _END_MASTER(myThid)
115          ENDIF
116  C  C
117        _EXCH_XY_R4(fu  , myThid )        _EXCH_XY_R4(fu  , myThid )
118        _EXCH_XY_R4(fv  , myThid )        _EXCH_XY_R4(fv  , myThid )
119    
120  CcnhDebugStarts        CALL PLOT_FIELD_XYRS( fu, 'S/R INI_FORCING FU',1,myThid)
121  Cdbg  WRITE(0,*) ' distY = ', distY        CALL PLOT_FIELD_XYRS( fv, 'S/R INI_FORCING FV',1,myThid)
 Cdbg  WRITE(0,*) '    ly = ', lY  
 Cdbg  WRITE(0,*) ' tauMax= ', tauMax  
 Cdbg  CALL PLOT_FIELD_XYR8( fu, 'INI_FORCING FU',1,myThid)  
 Cdbg  STOP 'INI_FORCING'  
 CcnhDebugEnds  
122    
123        RETURN        RETURN
124        END        END
Legend:
Removed from v.1.1.1.1  
changed lines
  Added in v.1.9.2.1
  ViewVC Help
Powered by ViewVC 1.1.22