/[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.13 by jmc, Thu Mar 7 14:09:02 2002 UTC revision 1.14 by jmc, Wed Jun 22 00:25:32 2005 UTC
# Line 53  C-    Check that all thickness are > 0 : Line 53  C-    Check that all thickness are > 0 :
53    
54  C-    Calculate depths of centers and interfaces  C-    Calculate depths of centers and interfaces
55          rF(1)  = Ro_SeaLevel          rF(1)  = Ro_SeaLevel
56          rC(1)  = rF(1) - rkFac*delRc(1)          rC(1)  = rF(1) + rkSign*delRc(1)
57          drC(1) = delRc(1)          drC(1) = delRc(1)
58          drF(1) = delRc(1)          drF(1) = delRc(1)
59         DO K=2,Nr         DO K=2,Nr
60          drC(K)   = delRc(K)          drC(K)   = delRc(K)
61          drF(K-1) =  drF(K-1) + 0.5 _d 0*delRc(K)          drF(K-1) =  drF(K-1) + 0.5 _d 0*delRc(K)
62          drF(K)   = 0.5 _d 0*delRc(K)          drF(K)   = 0.5 _d 0*delRc(K)
63          rC(K)    = rC(K-1) - rkFac*drC(K)          rC(K)    = rC(K-1) + rkSign*drC(K)
64          rF(K)    = rF(K-1) - rkFac*drF(K-1)          rF(K)    = rF(K-1) + rkSign*drF(K-1)
65         ENDDO         ENDDO
66          drF(Nr)  = drF(Nr) + delRc(Nr+1)          drF(Nr)  = drF(Nr) + delRc(Nr+1)
67          rF(Nr+1) = rF(Nr) - rkFac*drF(Nr)          rF(Nr+1) = rF(Nr) + rkSign*drF(Nr)
68    
69        ELSE        ELSE
70  C-- Center at middle between 2 interfaces :  C-- Center at middle between 2 interfaces :
# Line 86  C-    Calculate depths of interfaces and Line 86  C-    Calculate depths of interfaces and
86        rF(1) = Ro_SeaLevel        rF(1) = Ro_SeaLevel
87        DO K=1,Nr        DO K=1,Nr
88         drF(K)     = delR(K)         drF(K)     = delR(K)
89         rF(K+1) = rF(K)-rkFac*delR(K)         rF(K+1) = rF(K) + rkSign*delR(K)
90        ENDDO        ENDDO
91        drC(1)      = delR(1) * 0.5 _d 0        drC(1)      = delR(1) * 0.5 _d 0
92        rC(1)       = rf(1)-rkFac*delR(1) * 0.5 _d 0        rC(1)       = rf(1) + rkSign*delR(1) * 0.5 _d 0
93        DO K=2,Nr        DO K=2,Nr
94         drC(K)     = 0.5 _d 0 *(delR(K-1)+delR(K))         drC(K)     = 0.5 _d 0 *(delR(K-1)+delR(K))
95         rC(K)      = rC(K-1) - rkFac*drC(K)         rC(K)      = rC(K-1) + rkSign*drC(K)
96        ENDDO        ENDDO
97    
98  C--  C--
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22