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

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

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

revision 1.9 by mlosch, Wed Jun 24 08:55:43 2009 UTC revision 1.10 by jmc, Wed Mar 24 03:05:11 2010 UTC
# Line 82  c     == external functions == Line 82  c     == external functions ==
82  c     == end of interface ==  c     == end of interface ==
83    
84        if ( myTime .GT. (endTime - lastinterval) ) then        if ( myTime .GT. (endTime - lastinterval) ) then
85           tempVar = 1. /           tempVar = 1. _d 0/
86       &             ( ( 1. + min(endTime-startTime,lastinterval) )       &             ( ( 1. _d 0 + min(endTime-startTime,lastinterval) )
87       &             / deltaTClock )       &             / deltaTClock )
88    
89           kSrf = 1           kSrf = 1
# Line 146  c     heat content of top level plus lat Line 146  c     heat content of top level plus lat
146                 do i =  1,snx                 do i =  1,snx
147                  objf_ice(bi,bj) = objf_ice(bi,bj) +                  objf_ice(bi,bj) = objf_ice(bi,bj) +
148       &                 tempVar * rA(i,j,bi,bj) * (       &                 tempVar * rA(i,j,bi,bj) * (
149       &                 (THETA(i,j,kSrf,bi,bj) + 1.96 ) *       &                 (THETA(i,j,kSrf,bi,bj) + 1.96 _d 0 ) *
150       &                 drF(1) * 3996 * 1026 -       &                 drF(1) * 3996. _d 0 * 1026. _d 0 -
151       &                 HEFF(i,j,bi,bj) * 334000 * 910 )       &                 HEFF(i,j,bi,bj) * 334000. _d 0 * 910. _d 0 )
152                 enddo                 enddo
153                enddo                enddo
154               enddo               enddo
# Line 162  c     heat content of top level Line 162  c     heat content of top level
162                 do i =  1,snx                 do i =  1,snx
163                  objf_ice(bi,bj) = objf_ice(bi,bj) +                  objf_ice(bi,bj) = objf_ice(bi,bj) +
164       &                 tempVar * rA(i,j,bi,bj) * (       &                 tempVar * rA(i,j,bi,bj) * (
165       &                 (THETA(i,j,kSrf,bi,bj) + 1.96 ) *       &                 (THETA(i,j,kSrf,bi,bj) + 1.96 _d 0 ) *
166       &                 drF(1) * 3996 * 1026 )       &                 drF(1) * 3996. _d 0 * 1026. _d 0 )
167                 enddo                 enddo
168                enddo                enddo
169               enddo               enddo
# Line 177  c     heat content of top level plus sea Line 177  c     heat content of top level plus sea
177                 do i =  1,snx                 do i =  1,snx
178                  objf_ice(bi,bj) = objf_ice(bi,bj) +                  objf_ice(bi,bj) = objf_ice(bi,bj) +
179       &                 tempVar * rA(i,j,bi,bj) * (       &                 tempVar * rA(i,j,bi,bj) * (
180       &                 (THETA(i,j,kSrf,bi,bj) + 1.96 ) *       &                 (THETA(i,j,kSrf,bi,bj) + 1.96 _d 0 ) *
181       &                 drF(1) * 3996 * 1026 +       &                 drF(1) * 3996. _d 0 * 1026. _d 0 +
182       &                 (TICE(i,j,bi,bj) - 273.15 + 1.96 ) *       &                 (TICE(i,j,bi,bj) - 273.15 _d 0 + 1.96 _d 0 ) *
183       &                 HEFF(i,j,bi,bj) * 2090 * 910 -       &                 HEFF(i,j,bi,bj) * 2090. _d 0 * 910. _d 0 -
184       &                 HEFF(i,j,bi,bj) * 334000 * 910 -       &                 HEFF(i,j,bi,bj) * 334000. _d 0 * 910. _d 0 -
185       &                 HSNOW(i,j,bi,bj) * 334000 * 330 )       &                 HSNOW(i,j,bi,bj) * 334000. _d 0 * 330. _d 0 )
186                 enddo                 enddo
187                enddo                enddo
188               enddo               enddo
# Line 198  c     assumed that measurements are AREA Line 198  c     assumed that measurements are AREA
198                    do j = 1,sny                    do j = 1,sny
199                       do i =  1,snx                       do i =  1,snx
200                          objf_ice(bi,bj) = objf_ice(bi,bj) +                          objf_ice(bi,bj) = objf_ice(bi,bj) +
201       &                       ( AREA(i,j,bi,bj) - 0.5 ) *       &                       ( AREA(i,j,bi,bj) - 0.5 _d 0 ) *
202       &                       ( AREA(i,j,bi,bj) - 0.5 )       &                       ( AREA(i,j,bi,bj) - 0.5 _d 0 )
203                       enddo                       enddo
204                    enddo                    enddo
205                 enddo                 enddo
# Line 220  cph) Line 220  cph)
220    
221  #endif /* ALLOW_COST_ICE */  #endif /* ALLOW_COST_ICE */
222    
223          return
224        end        end
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22