/[MITgcm]/MITgcm/pkg/seaice/seaice_ocean_stress.F
ViewVC logotype

Diff of /MITgcm/pkg/seaice/seaice_ocean_stress.F

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

revision 1.6 by mlosch, Thu Mar 16 14:25:40 2006 UTC revision 1.28 by mlosch, Wed Nov 24 15:36:31 2010 UTC
# Line 4  C $Name$ Line 4  C $Name$
4  #include "SEAICE_OPTIONS.h"  #include "SEAICE_OPTIONS.h"
5    
6  CStartOfInterface  CStartOfInterface
7        SUBROUTINE SEAICE_OCEAN_STRESS(        SUBROUTINE SEAICE_OCEAN_STRESS(
8       I     myTime, myIter, myThid )       I     myTime, myIter, myThid )
9  C     /==========================================================\  C     /==========================================================\
10  C     | SUBROUTINE SEAICE_OCEAN_STRESS                           |  C     | SUBROUTINE SEAICE_OCEAN_STRESS                           |
# Line 18  C     === Global variables === Line 18  C     === Global variables ===
18  #include "SIZE.h"  #include "SIZE.h"
19  #include "EEPARAMS.h"  #include "EEPARAMS.h"
20  #include "PARAMS.h"  #include "PARAMS.h"
21    #include "DYNVARS.h"
22  #include "GRID.h"  #include "GRID.h"
23  #include "FFIELDS.h"  #include "FFIELDS.h"
24  #include "SEAICE.h"  #include "SEAICE.h"
25  #include "SEAICE_PARAMS.h"  #include "SEAICE_PARAMS.h"
 #include "SEAICE_FFIELDS.h"  
26    
27  C     === Routine arguments ===  C     === Routine arguments ===
28  C     myTime - Simulation time  C     myTime - Simulation time
# Line 33  C     myThid - Thread no. that called th Line 33  C     myThid - Thread no. that called th
33        INTEGER myThid        INTEGER myThid
34  CEndOfInterface  CEndOfInterface
35    
36  #ifdef SEAICE_CGRID  #ifdef SEAICE_CGRID
37  C     === Local variables ===  C     === Local variables ===
38  C     i,j,bi,bj - Loop counters  C     i,j,bi,bj - Loop counters
39    C     kSrf      - vertical index of surface layer
40        INTEGER i, j, bi, bj        INTEGER i, j, bi, bj
41        _RL  SINWAT, COSWAT, SINWIN, COSWIN        INTEGER kSrf
42        _RL  fuIce, fvIce, FX, FY        _RL  COSWAT
43          _RS  SINWAT
44          _RL  fuIceLoc, fvIceLoc
45        _RL  areaW, areaS        _RL  areaW, areaS
46    
47        _RL press       (1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy)  C     surrface level
48        _RL etaPlusZeta (1-Olx:sNx+Olx,1-Oly:sNy+Oly)        kSrf = 1
49        _RL zetaMinusEta(1-Olx:sNx+Olx,1-Oly:sNy+Oly)  C     introduce turning angle (default is zero)
       _RL etaMeanZ    (1-Olx:sNx+Olx,1-Oly:sNy+Oly)  
       _RL etaMeanU    (1-Olx:sNx+Olx,1-Oly:sNy+Oly)  
       _RL etaMeanV    (1-Olx:sNx+Olx,1-Oly:sNy+Oly)  
       _RL dVdx        (1-Olx:sNx+Olx,1-Oly:sNy+Oly)  
       _RL dVdy        (1-Olx:sNx+Olx,1-Oly:sNy+Oly)  
       _RL dUdx        (1-Olx:sNx+Olx,1-Oly:sNy+Oly)  
       _RL dUdy        (1-Olx:sNx+Olx,1-Oly:sNy+Oly)  
   
 c     introduce turning angle (default is zero)  
50        SINWAT=SIN(SEAICE_waterTurnAngle*deg2rad)        SINWAT=SIN(SEAICE_waterTurnAngle*deg2rad)
51        COSWAT=COS(SEAICE_waterTurnAngle*deg2rad)        COSWAT=COS(SEAICE_waterTurnAngle*deg2rad)
       SINWIN=SIN(SEAICE_airTurnAngle*deg2rad)  
       COSWIN=COS(SEAICE_airTurnAngle*deg2rad)  
   
 C--   Update overlap regions  
       CALL EXCH_UV_XY_RL(WINDX, WINDY, .TRUE., myThid)  
   
 #ifndef SEAICE_EXTERNAL_FLUXES  
 C--   Interpolate wind stress (N/m^2) from C-points of C-grid  
 C     to U and V points of C-grid for forcing the ocean model.  
       DO bj=myByLo(myThid),myByHi(myThid)  
        DO bi=myBxLo(myThid),myBxHi(myThid)  
         DO j=1,sNy  
          DO i=1,sNx  
           fu(I,J,bi,bj)=0.5*(WINDX(I,J,bi,bj) + WINDX(I-1,J,bi,bj))  
           fv(I,J,bi,bj)=0.5*(WINDY(I,J,bi,bj) + WINDY(I,J-1,bi,bj))  
          ENDDO  
         ENDDO  
        ENDDO  
       ENDDO  
 #endif /* ifndef SEAICE_EXTERNAL_FLUXES */  
52    
53        IF ( useHB87StressCoupling ) THEN        IF ( useHB87StressCoupling ) THEN
54  C  C
55  C     use an intergral over ice and ocean surface layer to define  C     use an intergral over ice and ocean surface layer to define
56  C     surface stresses on ocean following Hibler and Bryan (1987, JPO)  C     surface stresses on ocean following Hibler and Bryan (1987, JPO)
57  C      C
 C     recompute viscosities from updated ice velocities  
        CALL SEAICE_CALC_VISCOSITIES(  
      I      uIce(1-Olx,1-Oly,1,1,1), vIce(1-Olx,1-Oly,1,1,1),  
      I      zMin, zMax, hEffM, press0,  
      O      eta, zeta, press,  
      I      myThid )  
 C     re-compute internal stresses with updated ice velocities  
58         DO bj=myByLo(myThid),myByHi(myThid)         DO bj=myByLo(myThid),myByHi(myThid)
59          DO bi=myBxLo(myThid),myBxHi(myThid)          DO bi=myBxLo(myThid),myBxHi(myThid)
60           DO j=1-Oly+1,sNy+Oly-1           DO J=1,sNy
61            DO i=1-Olx+1,sNx+Olx-1            DO I=1,sNx
62             etaPlusZeta (I,J) =  eta(I,J,bi,bj) + zeta(I,J,bi,bj)  C     average wind stress over ice and ocean and apply averaged wind
63             zetaMinusEta(I,J) = zeta(I,J,bi,bj) -  eta(I,J,bi,bj)  C     stress and internal ice stresses to surface layer of ocean
64             etaMeanU (I,J) =             areaW = 0.5 * (AREA(I,J,bi,bj) + AREA(I-1,J,bi,bj))
65       &          HALF*(ETA (I,J,bi,bj) + ETA (I-1,J  ,bi,bj))       &          * SEAICEstressFactor
66             etaMeanV (I,J) =             fu(I,J,bi,bj)=(ONE-areaW)*fu(I,J,bi,bj)
67       &          HALF*(ETA (I,J,bi,bj) + ETA (I  ,J-1,bi,bj))       &          + areaW*taux(I,J,bi,bj)
68             etaMeanZ (I,J) = QUART *         &          + stressDivergenceX(I,J,bi,bj) * SEAICEstressFactor
69       &          ( eta(I  ,J,bi,bj) + eta(I  ,J-1,bi,bj)            ENDDO
70       &          + eta(I-1,J,bi,bj) + eta(I-1,J-1,bi,bj) )           ENDDO
71             dUdx(I,J) = ( uIce(I+1,J,1,bi,bj) - uIce(I,J,1,bi,bj) )           DO J=1,sNy
72       &          * _recip_dxF(I,J,bi,bj)            DO I=1,sNx
73             dUdy(I,J) = ( uIce(I,J+1,1,bi,bj) - uIce(I,J,1,bi,bj) )             areaS = 0.5 * (AREA(I,J,bi,bj) + AREA(I,J-1,bi,bj))
74       &          * _recip_dyU(I,J+1,bi,bj)       &          * SEAICEstressFactor
75             dVdx(I,J) = ( vIce(I+1,J,1,bi,bj) - vIce(I,J,1,bi,bj) )             fv(I,J,bi,bj)=(ONE-areaS)*fv(I,J,bi,bj)
76       &          * _recip_dxV(I+1,J,bi,bj)       &          + areaS*tauy(I,J,bi,bj)
77             dVdy(I,J) = ( vIce(I,J+1,1,bi,bj) - vIce(I,J,1,bi,bj) )       &          + stressDivergenceY(I,J,bi,bj) * SEAICEstressFactor
      &          * _recip_dyF(I,J,bi,bj)  
78            ENDDO            ENDDO
79           ENDDO           ENDDO
          DO J = 1,sNy  
           DO I = 1,sNx  
 C     First FX = (d/dx)*sigam  
 C     + d/dx[ eta+zeta d/dx ] U  
            FX = _recip_dxC(I,J,bi,bj) *  
      &            ( etaPlusZeta(I  ,J) * dUdx(I  ,J)  
      &            - etaPlusZeta(I-1,J) * dUdx(I-1,J) )  
 C     + (d/dy)[eta*(d/dy + tanphi/a)] U (also on UVRT1/2)  
            FX = FX + _recip_dyG(I,J,bi,bj) * (  
      &          ( etaMeanZ(I,J+1) * dUdy(I,J+1)  
      &          - etaMeanZ(I,J  ) * dUdy(I,J  )  
      &          )  
      &          - ( etaMeanZ(I,J+1)  
      &            * ( uIce(I,J+1,1,bi,bj)+uIce(I,J,1,bi,bj) )  
      &            - etaMeanZ(I,J  )  
      &            * ( uIce(I,J-1,1,bi,bj)+uIce(I,J,1,bi,bj) ) )  
      &          * 0.5 _d 0 * _tanPhiAtU(I,J,bi,bj)  
      &          * recip_rSphere )  
 C     - 2*eta*(tanphi/a) * ( tanphi/a ) U  
            FX = FX - TWO * uIce(I,J,1,bi,bj)  
      &          * etaMeanU(I,J)*recip_rSphere*recip_rSphere  
      &          * _tanPhiAtU(I,J,bi,bj)  * _tanPhiAtU(I,J,bi,bj)  
 C     + d/dx[ (zeta-eta) dV/dy]  
            FX = FX +  
      &          ( zetaMinusEta(I  ,J  ) * dVdy(I  ,J  )  
      &          - zetaMinusEta(I-1,J  ) * dVdy(I-1,J  )  
      &          ) * _recip_dxC(I,J,bi,bj)  
 C     + d/dy[ eta dV/x ]  
            FX = FX + (  
      &            etaMeanZ(I,J+1)  
      &          * ( vIce(I  ,J+1,1,bi,bj) - vIce(I-1,J+1,1,bi,bj) )  
      &          * _recip_dxV(I,J+1,bi,bj)  
      &          - etaMeanZ(I,J  )  
      &          * ( vIce(I  ,J,1,bi,bj) - vIce(I-1,J,1,bi,bj) )  
      &          * _recip_dxV(I,J,bi,bj)  
      &          ) * _recip_dyG(I,J,bi,bj)  
 C     - d/dx[ (eta+zeta) * v * (tanphi/a) ]  
            FX = FX - (  
      &            etaPlusZeta(I  ,J)  
      &          * 0.5 * (vIce(I  ,J,1,bi,bj)+vIce(I  ,J+1,1,bi,bj))  
      &          * 0.5 * ( _tanPhiAtU(I  ,J,bi,bj)  
      &          + _tanPhiAtU(I+1,J,bi,bj) )  
      &          - etaPlusZeta(I-1,J) *  
      &          * 0.5 * (vIce(I-1,J,1,bi,bj)+vIce(I-1,J+1,1,bi,bj))  
      &          * 0.5 * ( _tanPhiAtU(I-1,J,bi,bj)  
      &          + _tanPhiAtU(I  ,J,bi,bj) )  
      &          )* _recip_dxC(I,J,bi,bj)*recip_rSphere  
 C     - 2*eta*(tanphi/a) * dV/dx  
            FX = FX  
      &          -TWO * etaMeanU(I,J) * _tanPhiAtV(I,J,bi,bj)  
      &          *recip_rSphere  
      &          *(vIce(I  ,J,1,bi,bj) + vIce(I  ,J+1,1,bi,bj)  
      &          - vIce(I-1,J,1,bi,bj) - vIce(I-1,J+1,1,bi,bj))  
      &          * _recip_dxC(I,J,bi,bj)  
 C     - (d/dx) P/2  
            FX = _maskW(I,J,1,bi,bj) * ( FX - _recip_dxC(I,J,bi,bj)  
      &          * ( press(I,J,bi,bj) - press(I-1,J,bi,bj) ) )  
 C  
 C     then FY = (d/dy)*sigam  
 C     + d/dy [(eta+zeta) d/dy] V  
            FY = _recip_dyC(I,J,bi,bj) *  
      &          ( dVdy(I,J  ) * etaPlusZeta(I,J  )  
      &          - dVdy(I,J-1) * etaPlusZeta(I,J-1) )  
 C     + d/dx [eta d/dx] V  
            FY = FY +  _recip_dxC(I,J,bi,bj) *  
      &          ( eta(I  ,J,bi,bj) * dVdx(I  ,J)  
      &          - eta(I-1,J,bi,bj) * dVdx(I-1,J) )  
 C     - d/dy [(zeta-eta) tanphi/a] V  
            FY = FY - _recip_dyC(I,J,bi,bj) * recip_rSphere * (  
      &            zetaMinusEta(I,J  ) * tanPhiAtU(I,J  ,bi,bj)  
      &          * 0.5 * ( vIce(I,J,1,bi,bj) + vIce(I,J+1,1,bi,bj))  
      &          - zetaMinusEta(I,J-1) * tanPhiAtU(I,J-1,bi,bj)  
      &          * 0.5 * ( vIce(I,J,1,bi,bj) + vIce(I,J-1,1,bi,bj)) )  
 C     2*eta tanphi/a ( - tanphi/a - d/dy) V  
            FY = FY - TWO*etaMeanV(I,J) * recip_rSphere  
      &          * _tanPhiAtV(I,J,bi,bj) * (  
      &            _tanPhiAtV(I,J,bi,bj) * recip_rSphere  
      &          + _recip_dyC(I,J,bi,bj) *  
      &          ( 0.5 * ( vIce(I,J,1,bi,bj) + vIce(I,J+1,1,bi,bj))  
      &          - 0.5 * ( vIce(I,J,1,bi,bj) + vIce(I,J-1,1,bi,bj)) ) )  
 C     + d/dy[ (zeta-eta) dU/dx ]  
            FY = FY +  
      &          ( zetaMinusEta(I,J  )*dUdx(I,J  )  
      &          - zetaMinusEta(I,J-1)*dUdx(I,J-1) )  
      &          * _recip_dyC(I,J,bi,bj)  
 C     + d/dx[ eta dU/dy ]  
            FY = FY + _recip_dxG(I,J,bi,bj) *  
      &          ( etaMeanZ(I+1,J) * dUdy(I+1,J)  
      &          - etaMeanZ(I  ,J) * dUdy(I  ,J) )  
 C     + d/dx[ eta * (tanphi/a) * U ]  
            FY = FY + (  
      &            etaMeanZ(I+1,J) * 0.5 *  
      &          ( uIce(I+1,J  ,1,bi,bj) * _tanPhiAtU(I+1,J  ,bi,bj)  
      &          + uIce(I+1,J-1,1,bi,bj) * _tanPhiAtU(I+1,J-1,bi,bj) )  
      &          - etaMeanZ(I  ,J) * 0.5 *  
      &          ( uIce(I  ,J  ,1,bi,bj) * _tanPhiAtU(I  ,J  ,bi,bj)  
      &          + uIce(I  ,J-1,1,bi,bj) * _tanPhiAtU(I  ,J  ,bi,bj) )  
      &          ) *  _recip_dxG(I,J,bi,bj)*recip_rSphere  
 C     + 2*eta*(tanphi/a) dU/dx  
            FY = FY +  
      &          TWO * etaMeanV(I,J)*TWO  * _tanPhiAtV(I,J,bi,bj)  
      &          * ( uIce(I+1,J,1,bi,bj)+uIce(I+1,J-1,1,bi,bj)  
      &            - uIce(I  ,J,1,bi,bj)-uIce(I  ,J-1,1,bi,bj) )  
      &          * _recip_dxG(I,J,bi,bj) * recip_rSphere  
 C     - (d/dy) P/2  
            FY = _maskS(I,J,1,bi,bj) * ( FY - _recip_dyC(I,J,bi,bj)  
      &          * ( press(I,J,bi,bj) - press(I,J-1,bi,bj) ) )  
 C      
 C     recompute wind stress over ice (done already in seaice_dynsolver,  
 C     but not saved)  
            fuIce = 0.5 _d 0 *  
      &          ( DAIRN(I  ,J,bi,bj)*(  
      &          COSWIN*uWind(I  ,J,bi,bj)  
      &          -SIGN(SINWIN, _fCori(I  ,J,bi,bj))*vWind(I  ,J,bi,bj) )  
      &          + DAIRN(I-1,J,bi,bj)*(  
      &          COSWIN*uWind(I-1,J,bi,bj)  
      &          -SIGN(SINWIN, _fCori(I-1,J,bi,bj))*vWind(I-1,J,bi,bj) )  
      &          )  
            fvIce = 0.5 _d 0 *  
      &          ( DAIRN(I,J  ,bi,bj)*(  
      &          SIGN(SINWIN, _fCori(I  ,J,bi,bj))*uWind(I,J  ,bi,bj)  
      &          +COSWIN*vWind(I,J  ,bi,bj) )  
      &          + DAIRN(I,J-1,bi,bj)*(  
      &          SIGN(SINWIN, _fCori(I,J-1,bi,bj))*uWind(I,J-1,bi,bj)  
      &          +COSWIN*vWind(I,J-1,bi,bj) )  
      &          )  
 C     average wind stress over ice and ocean and apply averaged wind  
 C     stress and internal ice stresses to surface layer of ocean  
            areaW = 0.5 * (AREA(I,J,1,bi,bj) + AREA(I-1,J,1,bi,bj))  
            areaS = 0.5 * (AREA(I,J,1,bi,bj) + AREA(I,J-1,1,bi,bj))  
            fu(I,J,bi,bj)=(ONE-areaW)*fu(I,J,bi,bj)+areaW*(fuIce + FX)  
            fv(I,J,bi,bj)=(ONE-areaS)*fv(I,J,bi,bj)+areaS*(fvIce + FY)  
           END DO  
          END DO  
80          ENDDO          ENDDO
81         ENDDO         ENDDO
82    
83        ELSE        ELSE
84    C     else: useHB87StressCoupling=F
85    
86  C--   Compute ice-affected wind stress (interpolate to U/V-points)  C--   Compute ice-affected wind stress (interpolate to U/V-points)
87  C     by averaging wind stress and ice-ocean stress according to  C     by averaging wind stress and ice-ocean stress according to
88  C     ice cover  C     ice cover
89        DO bj=myByLo(myThid),myByHi(myThid)        DO bj=myByLo(myThid),myByHi(myThid)
90         DO bi=myBxLo(myThid),myBxHi(myThid)         DO bi=myBxLo(myThid),myBxHi(myThid)
91          DO j=1,sNy          DO j=1,sNy
92           DO i=1,sNx           DO i=1,sNx
93            fuIce=HALF*( DWATN(I,J,bi,bj)+DWATN(I,J+1,bi,bj) )*            fuIceLoc=HALF*( DWATN(I,J,bi,bj)+DWATN(I-1,J,bi,bj) )*
94       &         COSWAT *       &         COSWAT *
95       &         ( UICE(I,J,1,bi,bj)-GWATX(I,J,bi,bj) )       &         ( uIce(I,J,bi,bj)-uVel(I,J,kSrf,bi,bj) )
96       &         - SIGN(SINWAT, _fCori(I,J,bi,bj)) * 0.5 _d 0 *       &         - SIGN(SINWAT, _fCori(I,J,bi,bj)) * 0.5 _d 0 *
97       &         ( DWATN(I  ,J,bi,bj) *       &         ( DWATN(I  ,J,bi,bj) *
98       &         0.5 _d 0*(vIce(I  ,J  ,1,bi,bj)-GWATY(I  ,J  ,bi,bj)       &         0.5 _d 0*(vIce(I  ,J  ,bi,bj)-vVel(I  ,J  ,kSrf,bi,bj)
99       &                  +vIce(I  ,J+1,1,bi,bj)-GWATY(I  ,J+1,bi,bj))       &                  +vIce(I  ,J+1,bi,bj)-vVel(I  ,J+1,kSrf,bi,bj))
100       &         + DWATN(I-1,J,bi,bj) *       &         + DWATN(I-1,J,bi,bj) *
101       &         0.5 _d 0*(vIce(I-1,J  ,1,bi,bj)-GWATY(I-1,J  ,bi,bj)       &         0.5 _d 0*(vIce(I-1,J  ,bi,bj)-vVel(I-1,J  ,kSrf,bi,bj)
102       &                  +vIce(I-1,J+1,1,bi,bj)-GWATY(I-1,J+1,bi,bj))       &                  +vIce(I-1,J+1,bi,bj)-vVel(I-1,J+1,kSrf,bi,bj))
103       &         )       &         )
104            fvIce=HALF*( DWATN(I,J,bi,bj)+DWATN(I+1,J,bi,bj) )*            fvIceLoc=HALF*( DWATN(I,J,bi,bj)+DWATN(I,J-1,bi,bj) )*
105       &         COSWAT *       &         COSWAT *
106       &         ( VICE(I,J,1,bi,bj)-GWATY(I,J,bi,bj) )       &         ( vIce(I,J,bi,bj)-vVel(I,J,kSrf,bi,bj) )
107       &         + SIGN(SINWAT,  _fCori(I,J,bi,bj)) * 0.5 _d 0 *       &         + SIGN(SINWAT,  _fCori(I,J,bi,bj)) * 0.5 _d 0 *
108       &         ( DWATN(I,J  ,bi,bj) *       &         ( DWATN(I,J  ,bi,bj) *
109       &         0.5 _d 0*(uIce(I  ,J  ,1,bi,bj)-GWATX(I  ,J  ,bi,bj)       &         0.5 _d 0*(uIce(I  ,J  ,bi,bj)-uVel(I  ,J  ,kSrf,bi,bj)
110       &                  +uIce(I+1,J  ,1,bi,bj)-GWATX(I+1,J  ,bi,bj))       &                  +uIce(I+1,J  ,bi,bj)-uVel(I+1,J  ,kSrf,bi,bj))
111       &         + DWATN(I,J-1,bi,bj) *       &         + DWATN(I,J-1,bi,bj) *
112       &         0.5 _d 0*(uIce(I  ,J-1,1,bi,bj)-GWATX(I  ,J-1,bi,bj)       &         0.5 _d 0*(uIce(I  ,J-1,bi,bj)-uVel(I  ,J-1,kSrf,bi,bj)
113       &                  +uIce(I+1,J-1,1,bi,bj)-GWATX(I+1,J-1,bi,bj))       &                  +uIce(I+1,J-1,bi,bj)-uVel(I+1,J-1,kSrf,bi,bj))
114       &         )       &         )
115            areaW = 0.5 _d 0 * (AREA(I,J,1,bi,bj) + AREA(I-1,J,1,bi,bj))            areaW = 0.5 _d 0 * (AREA(I,J,bi,bj) + AREA(I-1,J,bi,bj))
116            areaS = 0.5 _d 0 * (AREA(I,J,1,bi,bj) + AREA(I,J-1,1,bi,bj))       &         * SEAICEstressFactor
117            fu(I,J,bi,bj)=(ONE-areaW)*fu(I,J,bi,bj)+areaW*fuIce            areaS = 0.5 _d 0 * (AREA(I,J,bi,bj) + AREA(I,J-1,bi,bj))
118            fv(I,J,bi,bj)=(ONE-areaS)*fv(I,J,bi,bj)+areaS*fvIce       &         * SEAICEstressFactor
119              fu(I,J,bi,bj)=(ONE-areaW)*fu(I,J,bi,bj)+areaW*fuIceLoc
120              fv(I,J,bi,bj)=(ONE-areaS)*fv(I,J,bi,bj)+areaS*fvIceLoc
121           ENDDO           ENDDO
122          ENDDO          ENDDO
123         ENDDO         ENDDO
# Line 290  C     ice cover Line 125  C     ice cover
125        ENDIF        ENDIF
126        CALL EXCH_UV_XY_RS(fu, fv, .TRUE., myThid)        CALL EXCH_UV_XY_RS(fu, fv, .TRUE., myThid)
127    
128  #endif /* not SEAICE_CGRID */  #endif /* SEAICE_CGRID */
129    
130        RETURN        RETURN
131        END        END
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.28
  ViewVC Help
Powered by ViewVC 1.1.22