/[MITgcm]/MITgcm/pkg/obcs/obcs_init_variables.F
ViewVC logotype

Diff of /MITgcm/pkg/obcs/obcs_init_variables.F

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

revision 1.4 by adcroft, Tue May 29 14:01:39 2001 UTC revision 1.9 by adcroft, Thu Jul 8 16:54:43 2004 UTC
# Line 28  C     myThid -  Number of this instance Line 28  C     myThid -  Number of this instance
28  C     == Local variables ==  C     == Local variables ==
29        INTEGER bi, bj        INTEGER bi, bj
30        INTEGER I, J, K        INTEGER I, J, K
31          CHARACTER*(10) suff
32          INTEGER prec
33    
34        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
35         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
# Line 63  C     == Local variables == Line 65  C     == Local variables ==
65           ENDDO           ENDDO
66          ENDDO          ENDDO
67    
68    #ifdef NONLIN_FRSURF
69            DO I=1-Olx,sNx+Olx
70              OBNeta(I,bi,bj)=0.
71              OBSeta(I,bi,bj)=0.
72            ENDDO
73            DO J=1-Oly,sNy+Oly
74              OBEeta(J,bi,bj)=0.
75              OBWeta(J,bi,bj)=0.
76            ENDDO
77    #endif /* NONLIN_FRSURF */
78    
79  #ifdef ALLOW_ORLANSKI  #ifdef ALLOW_ORLANSKI
80          IF (useOrlanskiNorth.OR.useOrlanskiSouth.OR.          IF (useOrlanskiNorth.OR.useOrlanskiSouth.OR.
81       &      useOrlanskiEast.OR.useOrlanskiWest) THEN       &      useOrlanskiEast.OR.useOrlanskiWest) THEN
# Line 73  C     == Local variables == Line 86  C     == Local variables ==
86         ENDDO         ENDDO
87        ENDDO            ENDDO    
88    
89    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
90    C jmc: here is the logical place to read OBCS-pickup files
91    C      but a) without Orlanski: pass the test 1+1=2 without reading pickup.
92    C          b) with Orlanski: 1+1=2 fail even with this bit of code
93          IF ( nIter0.NE.0 ) THEN
94            prec = precFloat64
95            IF (pickupSuff.EQ.' ') THEN
96             WRITE(suff,'(I10.10)') nIter0
97            ELSE
98             WRITE(suff,'(A10)') pickupSuff
99            ENDIF
100    c       CALL OBCS_READ_CHECKPOINT(prec, nIter0, suff, myThid)
101          ENDIF
102    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
103    
104  C--   Apply OBCS values to initial conditions for consistancy  C--   Apply OBCS values to initial conditions for consistancy
105        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
106         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
107          CALL OBCS_CALC( bi, bj, startTime,          CALL OBCS_CALC( bi, bj, startTime, nIter0,
108       &        uVel, vVel, wVel, theta, salt, myThid )       &        uVel, vVel, wVel, theta, salt, myThid )
109          DO K=1,Nr          DO K=1,Nr
110           CALL OBCS_APPLY_UV( bi, bj, k, uVel, vVel, myThid )           CALL OBCS_APPLY_UV( bi, bj, k, uVel, vVel, myThid )
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22