/[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.4 by cnh, Sat Aug 22 17:51:09 1998 UTC revision 1.11 by cnh, Sun Feb 4 14:38:48 2001 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #include "CPP_EEOPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6  CStartOfInterface  CStartOfInterface
7        SUBROUTINE INI_VERTICAL_GRID( myThid )        SUBROUTINE INI_VERTICAL_GRID( myThid )
# Line 10  C     | o Initialise vertical gridding a Line 11  C     | o Initialise vertical gridding a
11  C     |==========================================================|  C     |==========================================================|
12  C     |                                                          |  C     |                                                          |
13  C     \==========================================================/  C     \==========================================================/
14          IMPLICIT NONE
15    
16  C     === Global variables ===  C     === Global variables ===
17  #include "SIZE.h"  #include "SIZE.h"
# Line 27  C     K Line 29  C     K
29        INTEGER K        INTEGER K
30    
31  C     Calculate depths of centers and interfaces  C     Calculate depths of centers and interfaces
32        rF(1) = 0. _d 0        rF(1) = Ro_SeaLevel
33        DO K=1,Nr        DO K=1,Nr
34         drF(K)     = delZ(K)         drF(K)     = delR(K)
35         rF(K+1) = rF(K)-delZ(K)         rF(K+1) = rF(K)-rkFac*delR(K)
36        ENDDO        ENDDO
37        drC(1)      = delZ(1) * 0.5 _d 0        drC(1)      = delR(1) * 0.5 _d 0
38        rC(1)       = -delZ(1) * 0.5 _d 0        rC(1)       = rf(1)-rkFac*delR(1) * 0.5 _d 0
39        DO K=2,Nr        DO K=2,Nr
40         drC(K)     = 0.5 _d 0 *(delZ(K-1)+delZ(K))         drC(K)     = 0.5 _d 0 *(delR(K-1)+delR(K))
41         rC(K)      = rC(K-1) - drC(K)         rC(K)      = rC(K-1) - rkFac*drC(K)
42        ENDDO        ENDDO
43        DO K=1,Nr        DO K=1,Nr
44         saFac(K)  = 1. D0         saFac(K)  = 1. _d 0
45         recip_drC(K)   = 1. _d 0/drC(K)         recip_drC(K)   = 1. _d 0/drC(K)
46         recip_drF(K)   = 1. _d 0/drF(K)         recip_drF(K)   = 1. _d 0/drF(K)
47        ENDDO        ENDDO
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22