/[MITgcm]/MITgcm/model/src/ini_theta.F
ViewVC logotype

Diff of /MITgcm/model/src/ini_theta.F

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

revision 1.6 by cnh, Mon Jun 15 05:13:56 1998 UTC revision 1.7 by cnh, Sat Aug 22 17:51:08 1998 UTC
# Line 55  C--    Initialise temperature field to a Line 55  C--    Initialise temperature field to a
55         rad = MIN(Ny/8,Nx/8)         rad = MIN(Ny/8,Nx/8)
56         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
57          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
58           DO K=1,Nz           DO K=1,Nr
59            DO J=1,sNy            DO J=1,sNy
60             DO I=1,sNx             DO I=1,sNx
61              theta(I,J,K,bi,bj) = 0. _d 0              theta(I,J,K,bi,bj) = 0. _d 0
# Line 74  C--    Example 2 Line 74  C--    Example 2
74  C--    Initialise temperature field to the vertical reference profile  C--    Initialise temperature field to the vertical reference profile
75         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
76          DO bi = myBxLo(myThid), myBxHi(myThid)          DO bi = myBxLo(myThid), myBxHi(myThid)
77           DO K=1,Nz           DO K=1,Nr
78            DO J=1,sNy            DO J=1,sNy
79             DO I=1,sNx             DO I=1,sNx
80              theta(I,J,K,bi,bj) = tRef(K)              theta(I,J,K,bi,bj) = tRef(K)
# Line 91  C--    Initialise temperature field to t Line 91  C--    Initialise temperature field to t
91  C     Set initial tendency terms  C     Set initial tendency terms
92        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
93         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
94          DO K=1,Nz          DO K=1,Nr
95           DO J=1,sNy           DO J=1,sNy
96            DO I=1,sNx            DO I=1,sNx
97             gt   (I,J,K,bi,bj) = 0. _d 0             gt   (I,J,K,bi,bj) = 0. _d 0
# Line 106  C Line 106  C
106        _EXCH_XYZ_R8(gt , myThid )        _EXCH_XYZ_R8(gt , myThid )
107        _EXCH_XYZ_R8(gtNM1 , myThid )        _EXCH_XYZ_R8(gtNM1 , myThid )
108    
109        CALL PLOT_FIELD_XYZRL( theta, 'Initial Temperature' , Nz, 1, myThid )        CALL PLOT_FIELD_XYZRL( theta, 'Initial Temperature' , Nr, 1, myThid )
110    
111        RETURN        RETURN
112        END        END
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22