/[MITgcm]/MITgcm_contrib/ksnow/press_release/code/update_cg2d.F
ViewVC logotype

Diff of /MITgcm_contrib/ksnow/press_release/code/update_cg2d.F

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

revision 1.3 by dgoldberg, Sat Feb 4 18:55:11 2017 UTC revision 1.4 by dgoldberg, Sat Mar 4 11:57:16 2017 UTC
# Line 101  C     aS2d: integral in Z Ay/dY Line 101  C     aS2d: integral in Z Ay/dY
101    
102                            
103               drag_fac = _recip_hFacW(i,j,k,bi,bj)**2*recip_drF(k)**2*               drag_fac = _recip_hFacW(i,j,k,bi,bj)**2*recip_drF(k)**2*
104       &        12000. * viscArNr(k) *       &        pReleaseDamp * viscArNr(k) *
105       &        1./(1+exp(-10./cg2dminColumnEps*       &        1./(1+exp(-10./cg2dminColumnEps*
106       &         (-1.)*(depthColW(i,j,bi,bj)-cg2dminColumnEps/2)))       &         (-1.)*(depthColW(i,j,bi,bj)-cg2dminColumnEps/2)))
107    
108               aW2d(i,j,bi,bj) = aW2d(i,j,bi,bj)               aW2d(i,j,bi,bj) = aW2d(i,j,bi,bj)
109       &               + faceArea*dU_psFacX(i,j,k,bi,bj)       &               + faceArea*dU_psFacX(i,j,k,bi,bj)
110       &               / (1 + drag_fac*deltaTmom)       &               / (1 + drag_fac*deltaTmom)
111         &                         *recip_dxC(i,j,bi,bj)
112              ELSE              ELSE
113  #endif  #endif
114               aW2d(i,j,bi,bj) = aW2d(i,j,bi,bj)               aW2d(i,j,bi,bj) = aW2d(i,j,bi,bj)
# Line 123  C     aS2d: integral in Z Ay/dY Line 124  C     aS2d: integral in Z Ay/dY
124              IF (depthColS(i,j,bi,bj).lt.cg2dminColumnEps) THEN                      IF (depthColS(i,j,bi,bj).lt.cg2dminColumnEps) THEN        
125    
126               drag_fac = _recip_hFacS(i,j,k,bi,bj)**2*recip_drF(k)**2*               drag_fac = _recip_hFacS(i,j,k,bi,bj)**2*recip_drF(k)**2*
127       &        12000. * viscArNr(k)*       &        pReleaseDamp * viscArNr(k)*
128       &        1./(1+exp(-10./cg2dminColumnEps*       &        1./(1+exp(-10./cg2dminColumnEps*
129       &         (-1.)*(depthColS(i,j,bi,bj)-cg2dminColumnEps/2)))       &         (-1.)*(depthColS(i,j,bi,bj)-cg2dminColumnEps/2)))
130    
# Line 131  C     aS2d: integral in Z Ay/dY Line 132  C     aS2d: integral in Z Ay/dY
132               aS2d(i,j,bi,bj) = aS2d(i,j,bi,bj)               aS2d(i,j,bi,bj) = aS2d(i,j,bi,bj)
133       &               + faceArea*dV_psFacY(i,j,k,bi,bj)       &               + faceArea*dV_psFacY(i,j,k,bi,bj)
134       &               / (1 + drag_fac*deltaTmom)       &               / (1 + drag_fac*deltaTmom)
135         &                         *recip_dyC(i,j,bi,bj)
136              ELSE              ELSE
137  #endif  #endif
138               aS2d(i,j,bi,bj) = aS2d(i,j,bi,bj)               aS2d(i,j,bi,bj) = aS2d(i,j,bi,bj)
# Line 152  C  deep-model: *deepFacC (faceArea), /de Line 154  C  deep-model: *deepFacC (faceArea), /de
154  #ifdef ALLOW_PRESSURE_RELEASE_CODE  #ifdef ALLOW_PRESSURE_RELEASE_CODE
155              IF (depthColW(i,j,bi,bj).lt.cg2dminColumnEps) THEN                      IF (depthColW(i,j,bi,bj).lt.cg2dminColumnEps) THEN        
156               drag_fac = _recip_hFacW(i,j,k,bi,bj)**2*recip_drF(k)**2*               drag_fac = _recip_hFacW(i,j,k,bi,bj)**2*recip_drF(k)**2*
157       &        12000. * viscArNr(k)*       &        pReleaseDamp * viscArNr(k)*
158       &        1./(1+exp(-10./cg2dminColumnEps*       &        1./(1+exp(-10./cg2dminColumnEps*
159       &         (-1.)*(depthColW(i,j,bi,bj)-cg2dminColumnEps/2)))       &         (-1.)*(depthColW(i,j,bi,bj)-cg2dminColumnEps/2)))
160              ELSE              ELSE
# Line 171  C  deep-model: *deepFacC (faceArea), /de Line 173  C  deep-model: *deepFacC (faceArea), /de
173  #ifdef ALLOW_PRESSURE_RELEASE_CODE  #ifdef ALLOW_PRESSURE_RELEASE_CODE
174              IF (depthColS(i,j,bi,bj).lt.cg2dminColumnEps) THEN                      IF (depthColS(i,j,bi,bj).lt.cg2dminColumnEps) THEN        
175               drag_fac = _recip_hFacS(i,j,k,bi,bj)**2*recip_drF(k)**2*               drag_fac = _recip_hFacS(i,j,k,bi,bj)**2*recip_drF(k)**2*
176       &        12000. * viscArNr(k) *       &        pReleaseDamp * viscArNr(k) *
177       &        1./(1+exp(-10./cg2dminColumnEps*       &        1./(1+exp(-10./cg2dminColumnEps*
178       &         (-1.)*(depthColS(i,j,bi,bj)-cg2dminColumnEps/2)))       &         (-1.)*(depthColS(i,j,bi,bj)-cg2dminColumnEps/2)))
179              ELSE              ELSE
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22