/[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.22 by heimbach, Mon May 14 22:02:34 2007 UTC revision 1.24 by heimbach, Thu Jun 21 04:06:21 2007 UTC
# Line 108  c--   Temperature field. Line 108  c--   Temperature field.
108       $          tmpfld3d(i,j,k,bi,bj)=       $          tmpfld3d(i,j,k,bi,bj)=
109       $          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))
110  #endif  #endif
111  #ifdef ALLOW_OPENAD  #ifdef ALLOW_AUTODIFF_OPENAD
112                  theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) +                  theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) +
113       &                               fac*xx_theta(i,j,k,bi,bj) +       &                               fac*xx_theta(i,j,k,bi,bj) +
114       &                               fac*tmpfld3d(i,j,k,bi,bj)       &                               fac*tmpfld3d(i,j,k,bi,bj)
# Line 146  c--   Temperature field. Line 146  c--   Temperature field.
146       $          tmpfld3d(i,j,k,bi,bj)=       $          tmpfld3d(i,j,k,bi,bj)=
147       $          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))
148  #endif  #endif
149  #ifdef ALLOW_OPENAD  #ifdef ALLOW_AUTODIFF_OPENAD
150                  salt(i,j,k,bi,bj) = salt(i,j,k,bi,bj) +                  salt(i,j,k,bi,bj) = salt(i,j,k,bi,bj) +
151       &                               fac*xx_salt(i,j,k,bi,bj) +       &                               fac*xx_salt(i,j,k,bi,bj) +
152       &                               fac*tmpfld3d(i,j,k,bi,bj)       &                               fac*tmpfld3d(i,j,k,bi,bj)
# Line 406  c-- initial zonal velocity Line 406  c-- initial zonal velocity
406            do k = 1,nr            do k = 1,nr
407              do j = jmin,jmax              do j = jmin,jmax
408                do i = imin,imax                do i = imin,imax
409  #ifdef ALLOW_OPENAD  #ifdef ALLOW_AUTODIFF_OPENAD
410                  uVel(i,j,k,bi,bj) = uVel(i,j,k,bi,bj) +                  uVel(i,j,k,bi,bj) = uVel(i,j,k,bi,bj) +
411       &                                  fac*xx_uvel(i,j,k,bi,bj)       &                                  fac*xx_uvel(i,j,k,bi,bj)
412  #else  #else
# Line 433  c-- initial merid. velocity Line 433  c-- initial merid. velocity
433            do k = 1,nr            do k = 1,nr
434              do j = jmin,jmax              do j = jmin,jmax
435                do i = imin,imax                do i = imin,imax
436  #ifdef ALLOW_OPENAD  #ifdef ALLOW_AUTODIFF_OPENAD
437                  vVel(i,j,k,bi,bj) = vVel(i,j,k,bi,bj) +                  vVel(i,j,k,bi,bj) = vVel(i,j,k,bi,bj) +
438       &                                  fac*xx_vvel(i,j,k,bi,bj)       &                                  fac*xx_vvel(i,j,k,bi,bj)
439  #else  #else
# Line 459  c--   initial Eta. Line 459  c--   initial Eta.
459          do bi = itlo,ithi          do bi = itlo,ithi
460            do j = jmin,jmax            do j = jmin,jmax
461              do i = imin,imax              do i = imin,imax
462  #ifdef ALLOW_OPENAD  #ifdef ALLOW_AUTODIFF_OPENAD
463                etaN(i,j,bi,bj) = etaN(i,j,bi,bj) +                etaN(i,j,bi,bj) = etaN(i,j,bi,bj) +
464       &                              fac*xx_etan(i,j,bi,bj)       &                              fac*xx_etan(i,j,bi,bj)
465  #else  #else
# Line 514  c--   SSS relaxation coefficient. Line 514  c--   SSS relaxation coefficient.
514        enddo        enddo
515  #endif  #endif
516    
517    #ifdef ALLOW_SEAICE
518          call seaice_ctrl_map_ini( mythid )
519    #endif
520    
521  c--   Update the tile edges.  c--   Update the tile edges.
522    
523  #if (defined (ALLOW_THETA0_CONTROL) || defined (ALLOW_SST0_CONTROL))  #if (defined (ALLOW_THETA0_CONTROL) || defined (ALLOW_SST0_CONTROL))
Legend:
Removed from v.1.22  
changed lines
  Added in v.1.24
  ViewVC Help
Powered by ViewVC 1.1.22