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

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

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

revision 1.7 by cnh, Fri Nov 6 22:44:48 1998 UTC revision 1.9 by adcroft, Tue May 18 17:42:35 1999 UTC
# Line 10  C     | o Initialise vertical gridding a Line 10  C     | o Initialise vertical gridding a
10  C     |==========================================================|  C     |==========================================================|
11  C     |                                                          |  C     |                                                          |
12  C     \==========================================================/  C     \==========================================================/
13          IMPLICIT NONE
14    
15  C     === Global variables ===  C     === Global variables ===
16  #include "SIZE.h"  #include "SIZE.h"
# Line 31  C     Calculate depths of centers and in Line 32  C     Calculate depths of centers and in
32        DO K=1,Nr        DO K=1,Nr
33         drF(K)     = delR(K)         drF(K)     = delR(K)
34         rF(K+1) = rF(K)-rkFac*delR(K)         rF(K+1) = rF(K)-rkFac*delR(K)
 C      drF(K)     = delZ(K)  
 C      rF(K+1)    = rF(K)-delZ(K)  
35        ENDDO        ENDDO
36        drC(1)      = delR(1) * 0.5 _d 0        drC(1)      = delR(1) * 0.5 _d 0
37        rC(1)       = rf(1)-rkFac*delR(1) * 0.5 _d 0        rC(1)       = rf(1)-rkFac*delR(1) * 0.5 _d 0
 C     drC(1)      = delZ(1) * 0.5 _d 0  
 C     rC(1)       = -delZ(1) * 0.5 _d 0  
38        DO K=2,Nr        DO K=2,Nr
39         drC(K)     = 0.5 _d 0 *(delR(K-1)+delR(K))         drC(K)     = 0.5 _d 0 *(delR(K-1)+delR(K))
40         rC(K)      = rC(K-1) - rkFac*drC(K)         rC(K)      = rC(K-1) - rkFac*drC(K)
 C      drC(K)     = 0.5 _d 0 *(delZ(K-1)+delZ(K))  
 C      rC(K)      = rC(K-1) - drC(K)  
41        ENDDO        ENDDO
42        DO K=1,Nr        DO K=1,Nr
43         saFac(K)  = 1. _d 0         saFac(K)  = 1. _d 0
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22