/[MITgcm]/MITgcm/model/src/ini_depths.F
ViewVC logotype

Diff of /MITgcm/model/src/ini_depths.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.5 by cnh, Wed Jun 10 01:44:03 1998 UTC revision 1.6 by cnh, Fri Jun 12 19:33:33 1998 UTC
# Line 105  C         Inverse of depth Line 105  C         Inverse of depth
105  C--   Now calculate "lopping" factors hFac.  C--   Now calculate "lopping" factors hFac.
106        zG = delZ(1)*0.5 D0        zG = delZ(1)*0.5 D0
107        zFace(1) = 0        zFace(1) = 0
108          dzC(1)   = delz(1)/2. _d 0
109        rDzC(1)  = 2. _d 0/delZ(1)        rDzC(1)  = 2. _d 0/delZ(1)
110        DO K=1,Nz-1        DO K=1,Nz-1
111         saFac(K) = 1. D0         saFac(K) = 1. D0
# Line 112  C--   Now calculate "lopping" factors hF Line 113  C--   Now calculate "lopping" factors hF
113         zG         = zG + (delZ(K)+delZ(K+1))*0.5 D0         zG         = zG + (delZ(K)+delZ(K+1))*0.5 D0
114         zFace(K+1) = zFace(K)+delZ(K)         zFace(K+1) = zFace(K)+delZ(K)
115         rDzC(K+1)  = 2. _d 0/(delZ(K)+delZ(K+1))         rDzC(K+1)  = 2. _d 0/(delZ(K)+delZ(K+1))
116           dzC(K+1)   = (delZ(K)+delZ(K+1))/2. _d 0
117        ENDDO        ENDDO
118        zC(Nz)      = zG        zC(Nz)      = zG
119        zFace(Nz+1) = zFace(Nz)+delZ(Nz)        zFace(Nz+1) = zFace(Nz)+delZ(Nz)
# Line 188  C--   Calculate recipricols of hFacC, hF Line 190  C--   Calculate recipricols of hFacC, hF
190        _EXCH_XYZ_R4(rhFacS , myThid )        _EXCH_XYZ_R4(rhFacS , myThid )
191    
192  C  C
193        CALL PLOT_FIELD_XYR4( H, 'Model depths' , 1, myThid )        CALL PLOT_FIELD_XYR8( H, 'Model depths' , 1, myThid )
194  C  C
195        RETURN        RETURN
196        END        END
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22