/[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.21 by jmc, Thu Aug 18 20:29:49 2011 UTC revision 1.22 by jmc, Fri Jul 26 14:22:22 2013 UTC
# Line 105  C        when trying to optimise this lo Line 105  C        when trying to optimise this lo
105    
106        IF (setCenterDr) THEN        IF (setCenterDr) THEN
107  C--   Cell Center r-distances are defined:  C--   Cell Center r-distances are defined:
108         DO k=1,Nr         DO k=1,Nr+1
109           drC(k) = delRc(k)           drC(k) = delRc(k)
110         ENDDO         ENDDO
111  C-    Check that all thickness are > 0 :  C-    Check that all thickness are > 0 :
# Line 127  C     assume Center at middle between 2 Line 127  C     assume Center at middle between 2
127         DO k=2,Nr         DO k=2,Nr
128           drC(k) = 0.5 _d 0 *(delR(k-1)+delR(k))           drC(k) = 0.5 _d 0 *(delR(k-1)+delR(k))
129         ENDDO         ENDDO
130           drC(Nr+1) = 0.5 _d 0 *delR(Nr)
131        ENDIF        ENDIF
132    
133  C---  Set r-position of  interFace (rF) and cell-Center (rC):  C---  Set r-position of  interFace (rF) and cell-Center (rC):
# Line 164  c       write(0,*) 'rC=',rC Line 165  c       write(0,*) 'rC=',rC
165        ENDDO        ENDDO
166    
167  C-    Calculate reciprol vertical grid spacing :  C-    Calculate reciprol vertical grid spacing :
168        DO k=1,Nr        DO k=1,Nr+1
169         recip_drC(k)   = 1. _d 0/drC(k)         recip_drC(k)   = 1. _d 0/drC(k)
170          ENDDO
171          DO k=1,Nr
172         recip_drF(k)   = 1. _d 0/drF(k)         recip_drF(k)   = 1. _d 0/drF(k)
173        ENDDO        ENDDO
174    
Legend:
Removed from v.1.21  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.22