/[MITgcm]/MITgcm/pkg/atm2d/sum_ocn_fluxes.F
ViewVC logotype

Diff of /MITgcm/pkg/atm2d/sum_ocn_fluxes.F

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

revision 1.1 by jscott, Wed Sep 6 15:32:40 2006 UTC revision 1.2 by jscott, Tue Sep 25 19:58:04 2007 UTC
# Line 35  C     LOCAL VARIABLES: Line 35  C     LOCAL VARIABLES:
35        DO j=1, sNy        DO j=1, sNy
36          DO i=1,sNx          DO i=1,sNx
37                    
38            IF (maskC(i,j,1,1,1).EQ.1.) THEN            IF (maskC(i,j,1,1,1) .EQ. 1.) THEN
39    
40  C         Ad hoc phase out wind stress if sufficient ice coverage  C         Ad hoc phase out wind stress if sufficient ice coverage
41            IF (iceMask(i,j,1,1).GT.windice_thres) THEN  C         similar idea as stressReduction used in thsice_main
42              IF (iceMask(i,j,1,1) .GT. windice_thres) THEN
43    
44             fu_2D(i,j)= fu_2D(i,j)*(1. _d 0-iceMask(i,j,1,1))             fu_2D(i,j)= fu_2D(i,j) * (1. _d 0 -
45       &                 /(1.d0-windice_thres)       &          (iceMask(i,j,1,1) + iceMask(i-1,j,1,1)) * 0.5 _d 0 )
46             fv_2D(i,j)= fv_2D(i,j)*(1. _d 0-iceMask(i,j,1,1))       &                 / (1. _d 0 - windice_thres)
47       &                 /(1.d0-windice_thres)             fv_2D(i,j)= fv_2D(i,j) * (1. _d 0 -
48         &          (iceMask(i,j,1,1) + iceMask(i,j-1,1,1)) * 0.5 _d 0 )
49         &                 / (1. _d 0 - windice_thres)
50               wspeed_2D(i,j)= wspeed_2D(i,j)
51         &                    * (1. _d 0 - iceMask(i,j,1,1))
52         &                    / (1. _d 0 - windice_thres)
53            ENDIF            ENDIF
54    
55            sum_runoff(i,j)= sum_runoff(i,j) + runoff_2D(i,j)            sum_runoff(i,j)= sum_runoff(i,j) + runoff_2D(i,j)
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22