/[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.29 by heimbach, Fri Mar 28 18:48:05 2008 UTC revision 1.30 by gforget, Fri May 30 02:48:28 2008 UTC
# Line 381  c--   bottom drag Line 381  c--   bottom drag
381        enddo        enddo
382  #endif  #endif
383    
384  #ifdef ALLOW_EDTAUX_CONTROL  #ifdef ALLOW_EDDYPSI_CONTROL
385  c-- zonal eddy stress : edtaux  c-- zonal eddy streamfunction : eddyPsiX
386        il=ilnblnk( xx_edtaux_file )        il=ilnblnk( xx_edtaux_file )
387        write(fnamegeneric(1:80),'(2a,i10.10)')        write(fnamegeneric(1:80),'(2a,i10.10)')
388       &     xx_edtaux_file(1:il),'.',optimcycle       &     xx_edtaux_file(1:il),'.',optimcycle
# Line 394  c-- zonal eddy stress : edtaux Line 394  c-- zonal eddy stress : edtaux
394            do k = 1,nr            do k = 1,nr
395              do j = jmin,jmax              do j = jmin,jmax
396                do i = imin,imax                do i = imin,imax
397                  eddyTauX(i,j,k,bi,bj) = eddyTauX(i,j,k,bi,bj) +                  eddyPsiX(i,j,k,bi,bj) = eddyPsiX(i,j,k,bi,bj) +
398       &            fCori(i,j,bi,bj)*tmpfld3d(i,j,k,bi,bj)       &            tmpfld3d(i,j,k,bi,bj)
399                enddo                enddo
400              enddo              enddo
401            enddo            enddo
402         enddo         enddo
403        enddo        enddo
404  #endif  c-- meridional eddy streamfunction : eddyPsiY
   
 #ifdef ALLOW_EDTAUY_CONTROL  
 c-- meridional eddy stress : edtauy  
405        il=ilnblnk( xx_edtauy_file )        il=ilnblnk( xx_edtauy_file )
406        write(fnamegeneric(1:80),'(2a,i10.10)')        write(fnamegeneric(1:80),'(2a,i10.10)')
407       &     xx_edtauy_file(1:il),'.',optimcycle       &     xx_edtauy_file(1:il),'.',optimcycle
# Line 416  c-- meridional eddy stress : edtauy Line 413  c-- meridional eddy stress : edtauy
413            do k = 1,nr            do k = 1,nr
414              do j = jmin,jmax              do j = jmin,jmax
415                do i = imin,imax                do i = imin,imax
416                  eddyTauY(i,j,k,bi,bj) = eddyTauY(i,j,k,bi,bj) +                  eddyPsiY(i,j,k,bi,bj) = eddyPsiY(i,j,k,bi,bj) +
417       &            fCoriG(i,j,bi,bj)*tmpfld3d(i,j,k,bi,bj)       &            tmpfld3d(i,j,k,bi,bj)
418                enddo                enddo
419              enddo              enddo
420            enddo            enddo
# Line 586  c--   Update the tile edges. Line 583  c--   Update the tile edges.
583        _EXCH_XY_R8( bottomdragfld, mythid )        _EXCH_XY_R8( bottomdragfld, mythid )
584  #endif  #endif
585    
586  #if (defined (ALLOW_EDTAUX_CONTROL) && defined (ALLOW_EDTAUY_CONTROL))  #ifdef ALLOW_EDDYPSI_CONTROL
587         CALL EXCH_UV_XYZ_RS(eddyTauX,eddyTauY,.TRUE.,myThid)         CALL EXCH_UV_XYZ_RS(eddyPsiX,eddyPsiY,.TRUE.,myThid)
 #elif (defined (ALLOW_EDTAUX_CONTROL) || defined (ALLOW_EDTAUY_CONTROL))  
        STOP 'ctrl_map_forcing: need BOTH ALLOW_EDTAU[X,Y]_CONTROL'  
588  #endif  #endif
589    
590  #ifdef ALLOW_UVEL0_CONTROL  #ifdef ALLOW_UVEL0_CONTROL
Legend:
Removed from v.1.29  
changed lines
  Added in v.1.30
  ViewVC Help
Powered by ViewVC 1.1.22