/[MITgcm]/MITgcm/pkg/ctrl/ctrl_map_ini.F
ViewVC logotype

Diff of /MITgcm/pkg/ctrl/ctrl_map_ini.F

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

revision 1.44 by gforget, Tue Oct 21 13:13:40 2014 UTC revision 1.45 by heimbach, Wed Feb 18 12:31:10 2015 UTC
# Line 124  c--   Temperature field. Line 124  c--   Temperature field.
124       $          tmpfld3d(i,j,k,bi,bj)=       $          tmpfld3d(i,j,k,bi,bj)=
125       $          sign(2.0/sqrt(wtheta(k,bi,bj)),tmpfld3d(i,j,k,bi,bj))       $          sign(2.0/sqrt(wtheta(k,bi,bj)),tmpfld3d(i,j,k,bi,bj))
126  #endif  #endif
 #ifdef ALLOW_OPENAD  
                 theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) +  
      &                               fac*xx_theta(i,j,k,bi,bj) +  
      &                               fac*tmpfld3d(i,j,k,bi,bj)  
 #else  
127                  theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) +                  theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) +
128       &                               fac*tmpfld3d(i,j,k,bi,bj)       &                               fac*tmpfld3d(i,j,k,bi,bj)
 #endif  
129  #ifndef DISABLE_CTRL_THETA_LIMIT  #ifndef DISABLE_CTRL_THETA_LIMIT
130                  if(theta(i,j,k,bi,bj).lt.-2.0)                  if(theta(i,j,k,bi,bj).lt.-2.0)
131       &               theta(i,j,k,bi,bj)= -2.0       &               theta(i,j,k,bi,bj)= -2.0
# Line 164  c--   Temperature field. Line 158  c--   Temperature field.
158       $          tmpfld3d(i,j,k,bi,bj)=       $          tmpfld3d(i,j,k,bi,bj)=
159       $          sign(2.0/sqrt(wsalt(k,bi,bj)),tmpfld3d(i,j,k,bi,bj))       $          sign(2.0/sqrt(wsalt(k,bi,bj)),tmpfld3d(i,j,k,bi,bj))
160  #endif  #endif
 #ifdef ALLOW_OPENAD  
161                  salt(i,j,k,bi,bj) = salt(i,j,k,bi,bj) +                  salt(i,j,k,bi,bj) = salt(i,j,k,bi,bj) +
      &                               fac*xx_salt(i,j,k,bi,bj) +  
162       &                               fac*tmpfld3d(i,j,k,bi,bj)       &                               fac*tmpfld3d(i,j,k,bi,bj)
 #else  
                 salt(i,j,k,bi,bj) = salt(i,j,k,bi,bj) +  
      &                               fac*tmpfld3d(i,j,k,bi,bj)  
 #endif  
163    
164                enddo                enddo
165              enddo              enddo
# Line 267  c--   diffkr. Line 255  c--   diffkr.
255            do k = 1,nr            do k = 1,nr
256              do j = jmin,jmax              do j = jmin,jmax
257                do i = imin,imax                do i = imin,imax
 #ifdef ALLOW_OPENAD  
258                  diffkr(i,j,k,bi,bj) = diffkr(i,j,k,bi,bj) +                  diffkr(i,j,k,bi,bj) = diffkr(i,j,k,bi,bj) +
      &                                xx_diffkr(i,j,k,bi,bj) +  
259       &                                tmpfld3d(i,j,k,bi,bj)       &                                tmpfld3d(i,j,k,bi,bj)
 #else  
                 diffkr(i,j,k,bi,bj) = diffkr(i,j,k,bi,bj) +  
      &                                tmpfld3d(i,j,k,bi,bj)  
 #endif  
260                enddo                enddo
261              enddo              enddo
262            enddo            enddo
Legend:
Removed from v.1.44  
changed lines
  Added in v.1.45
  ViewVC Help
Powered by ViewVC 1.1.22