/[MITgcm]/MITgcm/model/src/taueddy_init_varia.F
ViewVC logotype

Diff of /MITgcm/model/src/taueddy_init_varia.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by heimbach, Wed Feb 15 03:52:53 2006 UTC revision 1.2 by gforget, Fri May 30 02:45:43 2008 UTC
# Line 29  C     == Routine arguments == Line 29  C     == Routine arguments ==
29  C     myThid -  Number of this instance of INI_UVEL  C     myThid -  Number of this instance of INI_UVEL
30        INTEGER myThid        INTEGER myThid
31    
32  #ifdef ALLOW_TAU_EDDY  #ifdef ALLOW_EDDYPSI
33    
34  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
35  C     == Local variables ==  C     == Local variables ==
# Line 45  C--   Over all tiles Line 45  C--   Over all tiles
45          DO K=1,Nr          DO K=1,Nr
46           DO J=1-Oly,sNy+Oly           DO J=1-Oly,sNy+Oly
47            DO I=1-Olx,sNx+Olx            DO I=1-Olx,sNx+Olx
48             eddyTauX(i,j,k,bi,bj)=0.             eddyPsiX(i,j,k,bi,bj)=0.
49             eddyTauY(i,j,k,bi,bj)=0.             eddyPsiY(i,j,k,bi,bj)=0.
50            ENDDO            ENDDO
51           ENDDO           ENDDO
52          ENDDO          ENDDO
53         ENDDO         ENDDO
54        ENDDO        ENDDO
55    
56        IF ( eddyTauxFile .NE. ' ' ) THEN        IF ( eddyPsixFile .NE. ' ' ) THEN
57         CALL MDSREADFIELD(eddyTauxFile,precFloat32,         CALL MDSREADFIELD(eddyPsiXFile,precFloat32,
58       &                  'RS',Nr,eddyTauX,1,myThid)       &                  'RS',Nr,eddyPsiX,1,myThid)
59        ENDIF        ENDIF
60        IF ( eddyTauyFile .NE. ' ' ) THEN        IF ( eddyPsiyFile .NE. ' ' ) THEN
61         CALL MDSREADFIELD(eddyTauyFile,precFloat32,         CALL MDSREADFIELD(eddyPsiYFile,precFloat32,
62       &                  'RS',Nr,eddyTauY,1,myThid)       &                  'RS',Nr,eddyPsiY,1,myThid)
63        ENDIF        ENDIF
64    
65  cph(        CALL EXCH_UV_XYZ_RS(eddyPsiX,eddyPsiY,.TRUE.,myThid)
 cph divide by rhoConst to suppress it in momentum eqn.  
 cph see taueddy_external_forcing  
       DO bj = myByLo(myThid), myByHi(myThid)  
        DO bi = myBxLo(myThid), myBxHi(myThid)  
         DO K=1,Nr  
          DO J=1-Oly,sNy+Oly  
           DO I=1-Olx,sNx+Olx  
            eddyTauX(i,j,k,bi,bj)=eddyTauX(i,j,k,bi,bj)  
      &      *recip_rhoConst  
            eddyTauY(i,j,k,bi,bj)=eddyTauY(i,j,k,bi,bj)  
      &      *recip_rhoConst  
           ENDDO  
          ENDDO  
         ENDDO  
        ENDDO  
       ENDDO  
 cph)  
   
       CALL EXCH_UV_XYZ_RS(eddyTauX,eddyTauY,.TRUE.,myThid)  
66    
67  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
68    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22