/[MITgcm]/MITgcm_contrib/jscott/pkg_atm2d/fixed_flux_add.F
ViewVC logotype

Diff of /MITgcm_contrib/jscott/pkg_atm2d/fixed_flux_add.F

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

revision 1.1 by jscott, Fri Aug 11 18:55:50 2006 UTC revision 1.2 by jscott, Tue Aug 22 20:21:38 2006 UTC
# Line 119  C--   Interpolate and add to anomaly Line 119  C--   Interpolate and add to anomaly
119       &                (wght0*qnet0(i,j,1,1) + wght1*qnet1(i,j,1,1))       &                (wght0*qnet0(i,j,1,1) + wght1*qnet1(i,j,1,1))
120    
121  c note below is different from older code...  c note below is different from older code...
122            IF (iceMask(i,j,1,1).NE.0.D0)            IF (iceMask(i,j,1,1) .NE. 0. _d 0)
123       &      qneti_2D(i,j)= qneti_2D(i,j) +       &      qneti_2D(i,j)= qneti_2D(i,j) +
124       &                     (wght0*qnet0(i,j,1,1) + wght1*qnet1(i,j,1,1))       &                     (wght0*qnet0(i,j,1,1) + wght1*qnet1(i,j,1,1))
125    
126            IF (useObsEmP) THEN            IF (useObsEmP) THEN
127              evapo_2D(i,j)= (wght0*evap0(i,j,1,1) + wght1*evap1(i,j,1,1))              evapo_2D(i,j)= (wght0*evap0(i,j,1,1) + wght1*evap1(i,j,1,1))
128              precipo_2D(i,j)= (wght0*precip0(i,j,1,1) + wght1*precip1(i,j,1,1))              precipo_2D(i,j)= (wght0*precip0(i,j,1,1) + wght1*precip1(i,j,1,1))
129              IF (iceMask(i,j,1,1).NE.0.D0) THEN              IF (iceMask(i,j,1,1) .NE. 0. _d 0) THEN
130                evapi_2D(i,j)= (wght0*evap0(i,j,1,1) + wght1*evap1(i,j,1,1))                evapi_2D(i,j)= (wght0*evap0(i,j,1,1) + wght1*evap1(i,j,1,1))
131                precipi_2D(i,j)= (wght0*precip0(i,j,1,1) +                precipi_2D(i,j)= (wght0*precip0(i,j,1,1) +
132       &                          wght1*precip1(i,j,1,1))       &                          wght1*precip1(i,j,1,1))
# Line 136  c note below is different from older cod Line 136  c note below is different from older cod
136       &                (wght0*evap0(i,j,1,1) + wght1*evap1(i,j,1,1))       &                (wght0*evap0(i,j,1,1) + wght1*evap1(i,j,1,1))
137              precipo_2D(i,j)= precipo_2D(i,j) +              precipo_2D(i,j)= precipo_2D(i,j) +
138       &                (wght0*precip0(i,j,1,1) + wght1*precip1(i,j,1,1))       &                (wght0*precip0(i,j,1,1) + wght1*precip1(i,j,1,1))
139              IF (iceMask(i,j,1,1).NE.0.D0) THEN              IF (iceMask(i,j,1,1) .NE. 0. _d 0) THEN
140                evapi_2D(i,j)= evapi_2D(i,j) +                evapi_2D(i,j)= evapi_2D(i,j) +
141       &                (wght0*evap0(i,j,1,1) + wght1*evap1(i,j,1,1))       &                (wght0*evap0(i,j,1,1) + wght1*evap1(i,j,1,1))
142                precipi_2D(i,j)= precipi_2D(i,j) +                precipi_2D(i,j)= precipi_2D(i,j) +
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22