/[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.3 by adcroft, Wed Jul 29 21:37:21 1998 UTC revision 1.4 by cnh, Sat Aug 22 17:51:09 1998 UTC
# Line 27  C     K Line 27  C     K
27        INTEGER K        INTEGER K
28    
29  C     Calculate depths of centers and interfaces  C     Calculate depths of centers and interfaces
30        zFace(1) = 0. _d 0        rF(1) = 0. _d 0
31        DO K=1,Nz        DO K=1,Nr
32         dzF(K)     = delZ(K)         drF(K)     = delZ(K)
33         zFace(K+1) = zFace(K)-delZ(K)         rF(K+1) = rF(K)-delZ(K)
34        ENDDO        ENDDO
35        dzC(1)      = delZ(1) * 0.5 _d 0        drC(1)      = delZ(1) * 0.5 _d 0
36        zC(1)       = -delZ(1) * 0.5 _d 0        rC(1)       = -delZ(1) * 0.5 _d 0
37        DO K=2,Nz        DO K=2,Nr
38         dzC(K)     = 0.5 _d 0 *(delZ(K-1)+delZ(K))         drC(K)     = 0.5 _d 0 *(delZ(K-1)+delZ(K))
39         zC(K)      = zC(K-1) - dzC(K)         rC(K)      = rC(K-1) - drC(K)
40        ENDDO        ENDDO
41        DO K=1,Nz        DO K=1,Nr
42         saFac(K)  = 1. D0         saFac(K)  = 1. D0
43         rDzC(K)   = 1. _d 0/dzC(K)         recip_drC(K)   = 1. _d 0/drC(K)
44         rDzF(K)   = 1. _d 0/dzF(K)         recip_drF(K)   = 1. _d 0/drF(K)
45        ENDDO        ENDDO
46  C  C
47        RETURN        RETURN
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22