/[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.20 by heimbach, Wed Feb 15 03:52:54 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 350  c-- zonal eddy stress : edtaux Line 360  c-- zonal eddy stress : edtaux
360            do k = 1,nr            do k = 1,nr
361              do j = jmin,jmax              do j = jmin,jmax
362                do i = imin,imax                do i = imin,imax
363                  Eddytaux(i,j,k,bi,bj) = Eddytaux(i,j,k,bi,bj) +                  eddyTauX(i,j,k,bi,bj) = eddyTauX(i,j,k,bi,bj) +
364       &                                tmpfld3d(i,j,k,bi,bj)       &            fCori(i,j,bi,bj)*tmpfld3d(i,j,k,bi,bj)
365                enddo                enddo
366              enddo              enddo
367            enddo            enddo
# Line 372  c-- meridional eddy stress : edtauy Line 382  c-- meridional eddy stress : edtauy
382            do k = 1,nr            do k = 1,nr
383              do j = jmin,jmax              do j = jmin,jmax
384                do i = imin,imax                do i = imin,imax
385                  Eddytauy(i,j,k,bi,bj) = Eddytauy(i,j,k,bi,bj) +                  eddyTauY(i,j,k,bi,bj) = eddyTauY(i,j,k,bi,bj) +
386       &                                tmpfld3d(i,j,k,bi,bj)       &            fCoriG(i,j,bi,bj)*tmpfld3d(i,j,k,bi,bj)
387                enddo                enddo
388              enddo              enddo
389            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) +                  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 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
# Line 516  c--   Update the tile edges. Line 542  c--   Update the tile edges.
542  #endif  #endif
543    
544  #if (defined (ALLOW_EDTAUX_CONTROL) && defined (ALLOW_EDTAUY_CONTROL))  #if (defined (ALLOW_EDTAUX_CONTROL) && defined (ALLOW_EDTAUY_CONTROL))
545         CALL EXCH_UV_XYZ_RS(Eddytaux,Eddytauy,.TRUE.,myThid)         CALL EXCH_UV_XYZ_RS(eddyTauX,eddyTauY,.TRUE.,myThid)
546  #elif (defined (ALLOW_EDTAUX_CONTROL) || defined (ALLOW_EDTAUY_CONTROL))  #elif (defined (ALLOW_EDTAUX_CONTROL) || defined (ALLOW_EDTAUY_CONTROL))
547         STOP 'ctrl_map_forcing: need BOTH ALLOW_EDTAU[X,Y]_CONTROL'         STOP 'ctrl_map_forcing: need BOTH ALLOW_EDTAU[X,Y]_CONTROL'
548  #endif  #endif
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.20
  ViewVC Help
Powered by ViewVC 1.1.22