/[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.7 by cnh, Fri Nov 6 22:44:48 1998 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    
3  #include "CPP_EEOPTIONS.h"  #include "CPP_OPTIONS.h"
4    
5  CStartOfInterface  CStartOfInterface
6        SUBROUTINE INI_VERTICAL_GRID( myThid )        SUBROUTINE INI_VERTICAL_GRID( myThid )
# 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)     = delR(K)
33         zFace(K+1) = zFace(K)-delZ(K)         rF(K+1) = rF(K)-rkFac*delR(K)
34    C      drF(K)     = delZ(K)
35    C      rF(K+1)    = rF(K)-delZ(K)
36        ENDDO        ENDDO
37        dzC(1)      = delZ(1) * 0.5 _d 0        drC(1)      = delR(1) * 0.5 _d 0
38        zC(1)       = -delZ(1) * 0.5 _d 0        rC(1)       = rf(1)-rkFac*delR(1) * 0.5 _d 0
39        DO K=2,Nz  C     drC(1)      = delZ(1) * 0.5 _d 0
40         dzC(K)     = 0.5 _d 0 *(delZ(K-1)+delZ(K))  C     rC(1)       = -delZ(1) * 0.5 _d 0
41         zC(K)      = zC(K-1) - dzC(K)        DO K=2,Nr
42           drC(K)     = 0.5 _d 0 *(delR(K-1)+delR(K))
43           rC(K)      = rC(K-1) - rkFac*drC(K)
44    C      drC(K)     = 0.5 _d 0 *(delZ(K-1)+delZ(K))
45    C      rC(K)      = rC(K-1) - drC(K)
46        ENDDO        ENDDO
47        DO K=1,Nz        DO K=1,Nr
48         saFac(K)  = 1. D0         saFac(K)  = 1. _d 0
49         rDzC(K)   = 1. _d 0/dzC(K)         recip_drC(K)   = 1. _d 0/drC(K)
50         rDzF(K)   = 1. _d 0/dzF(K)         recip_drF(K)   = 1. _d 0/drF(K)
51        ENDDO        ENDDO
52  C  C
53        RETURN        RETURN
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22