/[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.18 by heimbach, Tue Jan 3 17:10:35 2006 UTC revision 1.19 by heimbach, Thu Jan 5 17:48:01 2006 UTC
# Line 404  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) +                  uVel(i,j,k,bi,bj) = uVel(i,j,k,bi,bj) +
409       &                              tmpfld3d(i,j,k,bi,bj)       &                                  fac*xx_uvel(i,j,k,bi,bj)
410    #else
411                    uVel(i,j,k,bi,bj) = uVel(i,j,k,bi,bj) +
412         &                                  fac*tmpfld3d(i,j,k,bi,bj)
413    #endif
414                enddo                enddo
415              enddo              enddo
416            enddo            enddo
# Line 426  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 447  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.18  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22