/[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.2 by ksnow, Mon Jan 30 16:32:17 2017 UTC revision 1.3 by dgoldberg, Sat Feb 4 18:55:11 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       &        30000. * viscArNr(k) *       &        12000. * 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    
# Line 123  C     aS2d: integral in Z Ay/dY Line 123  C     aS2d: integral in Z Ay/dY
123              IF (depthColS(i,j,bi,bj).lt.cg2dminColumnEps) THEN                      IF (depthColS(i,j,bi,bj).lt.cg2dminColumnEps) THEN        
124    
125               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*
126       &        30000. * viscArNr(k)*       &        12000. * viscArNr(k)*
127       &        1./(1+exp(-10./cg2dminColumnEps*       &        1./(1+exp(-10./cg2dminColumnEps*
128       &         (-1.)*(depthColS(i,j,bi,bj)-cg2dminColumnEps/2)))       &         (-1.)*(depthColS(i,j,bi,bj)-cg2dminColumnEps/2)))
129    
# Line 152  C  deep-model: *deepFacC (faceArea), /de Line 152  C  deep-model: *deepFacC (faceArea), /de
152  #ifdef ALLOW_PRESSURE_RELEASE_CODE  #ifdef ALLOW_PRESSURE_RELEASE_CODE
153              IF (depthColW(i,j,bi,bj).lt.cg2dminColumnEps) THEN                      IF (depthColW(i,j,bi,bj).lt.cg2dminColumnEps) THEN        
154               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*
155       &        30000. * viscArNr(k)*       &        12000. * viscArNr(k)*
156       &        1./(1+exp(-10./cg2dminColumnEps*       &        1./(1+exp(-10./cg2dminColumnEps*
157       &         (-1.)*(depthColW(i,j,bi,bj)-cg2dminColumnEps/2)))       &         (-1.)*(depthColW(i,j,bi,bj)-cg2dminColumnEps/2)))
158              ELSE              ELSE
# Line 171  C  deep-model: *deepFacC (faceArea), /de Line 171  C  deep-model: *deepFacC (faceArea), /de
171  #ifdef ALLOW_PRESSURE_RELEASE_CODE  #ifdef ALLOW_PRESSURE_RELEASE_CODE
172              IF (depthColS(i,j,bi,bj).lt.cg2dminColumnEps) THEN                      IF (depthColS(i,j,bi,bj).lt.cg2dminColumnEps) THEN        
173               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*
174       &        30000. * viscArNr(k) *       &        12000. * viscArNr(k) *
175       &        1./(1+exp(-10./cg2dminColumnEps*       &        1./(1+exp(-10./cg2dminColumnEps*
176       &         (-1.)*(depthColS(i,j,bi,bj)-cg2dminColumnEps/2)))       &         (-1.)*(depthColS(i,j,bi,bj)-cg2dminColumnEps/2)))
177              ELSE              ELSE
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22