/[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.17 by heimbach, Thu Apr 7 23:38:43 2005 UTC revision 1.19 by heimbach, Thu Jan 5 17:48:01 2006 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
112                    theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) +
113         &                               fac*xx_theta(i,j,k,bi,bj)
114    #else
115                  theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) +                  theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) +
116       &                               fac*tmpfld3d(i,j,k,bi,bj)       &                               fac*tmpfld3d(i,j,k,bi,bj)
117    #endif
118                  if(theta(i,j,k,bi,bj).lt.-2.0)                  if(theta(i,j,k,bi,bj).lt.-2.0)
119       &               theta(i,j,k,bi,bj)= -2.0       &               theta(i,j,k,bi,bj)= -2.0
120                enddo                enddo
# Line 140  c--   Temperature field. Line 145  c--   Temperature field.
145       $          tmpfld3d(i,j,k,bi,bj)=       $          tmpfld3d(i,j,k,bi,bj)=
146       $          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))
147  #endif  #endif
148    #ifdef ALLOW_OPENAD
149                    salt(i,j,k,bi,bj) = salt(i,j,k,bi,bj) +
150         &                               fac*xx_salt(i,j,k,bi,bj)
151    #else
152                  salt(i,j,k,bi,bj) = salt(i,j,k,bi,bj) +                  salt(i,j,k,bi,bj) = salt(i,j,k,bi,bj) +
153       &                               fac*tmpfld3d(i,j,k,bi,bj)       &                               fac*tmpfld3d(i,j,k,bi,bj)
154    #endif
155    
156                enddo                enddo
157              enddo              enddo
# Line 394  c-- initial zonal velocity Line 404  c-- initial zonal velocity
404            do k = 1,nr            do k = 1,nr
405              do j = jmin,jmax              do j = jmin,jmax
406                do i = imin,imax                do i = imin,imax
407    #ifdef ALLOW_OPENAD
408                    uVel(i,j,k,bi,bj) = uVel(i,j,k,bi,bj) +
409         &                                  fac*xx_uvel(i,j,k,bi,bj)
410    #else
411                  uVel(i,j,k,bi,bj) = uVel(i,j,k,bi,bj) +                  uVel(i,j,k,bi,bj) = uVel(i,j,k,bi,bj) +
412       &                              tmpfld3d(i,j,k,bi,bj)       &                                  fac*tmpfld3d(i,j,k,bi,bj)
413    #endif
414                enddo                enddo
415              enddo              enddo
416            enddo            enddo
# Line 416  c-- initial merid. velocity Line 431  c-- initial merid. velocity
431            do k = 1,nr            do k = 1,nr
432              do j = jmin,jmax              do j = jmin,jmax
433                do i = imin,imax                do i = imin,imax
434    #ifdef ALLOW_OPENAD
435                    vVel(i,j,k,bi,bj) = vVel(i,j,k,bi,bj) +
436         &                                  fac*xx_vvel(i,j,k,bi,bj)
437    #else
438                  vVel(i,j,k,bi,bj) = vVel(i,j,k,bi,bj) +                  vVel(i,j,k,bi,bj) = vVel(i,j,k,bi,bj) +
439       &                              tmpfld3d(i,j,k,bi,bj)       &                                  fac*tmpfld3d(i,j,k,bi,bj)
440    #endif
441                enddo                enddo
442              enddo              enddo
443            enddo            enddo
# Line 437  c--   initial Eta. Line 457  c--   initial Eta.
457          do bi = itlo,ithi          do bi = itlo,ithi
458            do j = jmin,jmax            do j = jmin,jmax
459              do i = imin,imax              do i = imin,imax
460                etaN(i,j,bi,bj) = etaN(i,j,bi,bj) + tmpfld2d(i,j,bi,bj)  #ifdef ALLOW_OPENAD
461                  etaN(i,j,bi,bj) = etaN(i,j,bi,bj) +
462         &                              fac*xx_etan(i,j,bi,bj)
463    #else
464                  etaN(i,j,bi,bj) = etaN(i,j,bi,bj) +
465         &                              fac*tmpfld2d(i,j,bi,bj)
466    #endif
467              enddo              enddo
468            enddo            enddo
469          enddo          enddo
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22