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

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

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

revision 1.39 by dimitri, Wed Dec 5 15:37:12 2007 UTC revision 1.40 by dimitri, Thu Dec 6 05:23:41 2007 UTC
# Line 24  C     === Global variables === Line 24  C     === Global variables ===
24  #ifdef ALLOW_EXF  #ifdef ALLOW_EXF
25  # include "EXF_OPTIONS.h"  # include "EXF_OPTIONS.h"
26  # include "EXF_FIELDS.h"  # include "EXF_FIELDS.h"
27    # include "EXF_PARAM.h"
28  #endif  #endif
29  #ifdef ALLOW_SALT_PLUME  #ifdef ALLOW_SALT_PLUME
30  # include "SALT_PLUME.h"  # include "SALT_PLUME.h"
# Line 253  C FIRST DO OPEN WATER Line 254  C FIRST DO OPEN WATER
254       I       bi, bj, myThid )       I       bi, bj, myThid )
255    
256  C NOW DO ICE  C NOW DO ICE
257            IF (useRelativeWind) THEN
258    C     Compute relative wind speed over sea ice.
259             DO J=1,sNy
260              DO I=1,sNx
261               SPEED_SQ =
262         &          (uWind(I,J,bi,bj)
263         &          +0.5 _d 0*(uVel(i,j,1,bi,bj)+uVel(i+1,j,1,bi,bj))
264         &          -0.5 _d 0*(uice(i,j,1,bi,bj)+uice(i+1,j,1,bi,bj)))**2
265         &          +(vWind(I,J,bi,bj)
266         &          +0.5 _d 0*(vVel(i,j,1,bi,bj)+vVel(i,j+1,1,bi,bj))
267         &          -0.5 _d 0*(vice(i,j,1,bi,bj)+vice(i,j+1,1,bi,bj)))**2
268               IF ( SPEED_SQ .LE. SEAICE_EPS_SQ ) THEN
269                 UG(I,J)=SEAICE_EPS
270               ELSE
271                 UG(I,J)=SQRT(SPEED_SQ)
272               ENDIF
273              ENDDO
274             ENDDO
275            ENDIF
276  #ifdef SEAICE_MULTICATEGORY  #ifdef SEAICE_MULTICATEGORY
277  C--  Start loop over muli-categories  C--  Start loop over muli-categories
278          DO IT=1,MULTDIM          DO IT=1,MULTDIM
Legend:
Removed from v.1.39  
changed lines
  Added in v.1.40
  ViewVC Help
Powered by ViewVC 1.1.22