/[MITgcm]/MITgcm/pkg/aim/aim_do_atmos_physics.F
ViewVC logotype

Diff of /MITgcm/pkg/aim/aim_do_atmos_physics.F

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

revision 1.6 by jmc, Wed Aug 15 13:50:43 2001 UTC revision 1.7 by jmc, Mon Aug 27 18:45:47 2001 UTC
# Line 211  c---+----1----+----2----+----3----+----4 Line 211  c---+----1----+----2----+----3----+----4
211  c_jmc: add square of surface wind speed (center of C grid) = 2 * KE_surf  c_jmc: add square of surface wind speed (center of C grid) = 2 * KE_surf
212        DO J = 1, sNy        DO J = 1, sNy
213         DO I = 1, sNx         DO I = 1, sNx
214           I2 = I+(J-1)*sNx          I2 = I+(J-1)*sNx
215           K = k_surf(i,j,bi,bj)          K = ksurfC(i,j,bi,bj)
216            IF (K.LE.Nr) THEN
217           Vsurfsq(I2,myThid) = 0.5 * (           Vsurfsq(I2,myThid) = 0.5 * (
218       &                 uVel(I,J,K,bi,bj)*uVel(I,J,K,bi,bj)       &                 uVel(I,J,K,bi,bj)*uVel(I,J,K,bi,bj)
219       &               + uVel(I+1,J,K,bi,bj)*uVel(I+1,J,K,bi,bj)       &               + uVel(I+1,J,K,bi,bj)*uVel(I+1,J,K,bi,bj)
# Line 226  c - to reproduce old results : Line 227  c - to reproduce old results :
227       &     UG1(I2,Katm,myThid)*UG1(I2,Katm,myThid)       &     UG1(I2,Katm,myThid)*UG1(I2,Katm,myThid)
228       &   + VG1(I2,Katm,myThid)*VG1(I2,Katm,myThid)       &   + VG1(I2,Katm,myThid)*VG1(I2,Katm,myThid)
229  #endif /* OLD_AIM_GRIG_MAPPING */  #endif /* OLD_AIM_GRIG_MAPPING */
230            ELSE
231             Vsurfsq(I2,myThid) = 0.
232            ENDIF
233         ENDDO         ENDDO
234        ENDDO        ENDDO
235  c---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  c---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22