Changeset 2528 for trunk/NEMOGCM/NEMO/OPA_SRC/DOM/closea.F90
- Timestamp:
- 2010-12-27T18:33:53+01:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/DOM/closea.F90
- Property svn:eol-style deleted
r1601 r2528 44 44 # include "vectopt_loop_substitute.h90" 45 45 !!---------------------------------------------------------------------- 46 !! OPA 9.0 , LOCEAN-IPSL (2006)46 !! NEMO/OPA 3.3 , NEMO Consortium (2010) 47 47 !! $Id$ 48 !! Software governed by the CeCILL licence ( modipsl/doc/NEMO_CeCILL.txt)48 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 49 49 !!---------------------------------------------------------------------- 50 50 … … 180 180 INTEGER :: ji, jj, jc, jn ! dummy loop indices 181 181 REAL(wp) :: zze2 182 REAL(wp), DIMENSION (jpncs) :: zemp 182 REAL(wp), DIMENSION (jpncs) :: zfwf 183 183 184 !!---------------------------------------------------------------------- 184 185 ! … … 216 217 ! !--------------------! 217 218 ! ! update emp, emps ! 218 z emp= 0.e0 !--------------------!219 zfwf = 0.e0 !--------------------! 219 220 DO jc = 1, jpncs 220 221 DO jj = ncsj1(jc), ncsj2(jc) 221 222 DO ji = ncsi1(jc), ncsi2(jc) 222 z emp(jc) = zemp(jc) + e1t(ji,jj) * e2t(ji,jj) * emp(ji,jj) * tmask_i(ji,jj)223 zfwf(jc) = zfwf(jc) + e1t(ji,jj) * e2t(ji,jj) * ( emp(ji,jj)-rnf(ji,jj) ) * tmask_i(ji,jj) 223 224 END DO 224 225 END DO 225 226 END DO 226 IF( lk_mpp ) CALL mpp_sum ( z emp(:) , jpncs ) ! mpp: sum over all the global domain227 IF( lk_mpp ) CALL mpp_sum ( zfwf(:) , jpncs ) ! mpp: sum over all the global domain 227 228 228 229 IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN ! Black Sea case for ORCA_R2 configuration 229 zze2 = ( z emp(3) + zemp(4) ) / 2.230 z emp(3) = zze2231 z emp(4) = zze2230 zze2 = ( zfwf(3) + zfwf(4) ) / 2. 231 zfwf(3) = zze2 232 zfwf(4) = zze2 232 233 ENDIF 233 234 … … 236 237 IF( ncstt(jc) == 0 ) THEN 237 238 ! water/evap excess is shared by all open ocean 238 emp (:,:) = emp (:,:) + z emp(jc) / surf(jpncs+1)239 emps(:,:) = emps(:,:) + z emp(jc) / surf(jpncs+1)239 emp (:,:) = emp (:,:) + zfwf(jc) / surf(jpncs+1) 240 emps(:,:) = emps(:,:) + zfwf(jc) / surf(jpncs+1) 240 241 ELSEIF( ncstt(jc) == 1 ) THEN 241 242 ! Excess water in open sea, at outflow location, excess evap shared 242 IF ( z emp(jc) <= 0.e0 ) THEN243 IF ( zfwf(jc) <= 0.e0 ) THEN 243 244 DO jn = 1, ncsnr(jc) 244 245 ji = mi0(ncsir(jc,jn)) … … 246 247 IF ( ji > 1 .AND. ji < jpi & 247 248 .AND. jj > 1 .AND. jj < jpj ) THEN 248 emp (ji,jj) = emp (ji,jj) + z emp(jc) / &249 emp (ji,jj) = emp (ji,jj) + zfwf(jc) / & 249 250 (FLOAT(ncsnr(jc)) * e1t(ji,jj) * e2t(ji,jj)) 250 emps(ji,jj) = emps(ji,jj) + z emp(jc) / &251 emps(ji,jj) = emps(ji,jj) + zfwf(jc) / & 251 252 (FLOAT(ncsnr(jc)) * e1t(ji,jj) * e2t(ji,jj)) 252 253 END IF 253 254 END DO 254 255 ELSE 255 emp (:,:) = emp (:,:) + z emp(jc) / surf(jpncs+1)256 emps(:,:) = emps(:,:) + z emp(jc) / surf(jpncs+1)256 emp (:,:) = emp (:,:) + zfwf(jc) / surf(jpncs+1) 257 emps(:,:) = emps(:,:) + zfwf(jc) / surf(jpncs+1) 257 258 ENDIF 258 259 ELSEIF( ncstt(jc) == 2 ) THEN … … 263 264 ji = mi0(ncsir(jc,jn)) 264 265 jj = mj0(ncsjr(jc,jn)) ! Location of outflow in open ocean 265 emp (ji,jj) = emp (ji,jj) + z emp(jc) &266 emp (ji,jj) = emp (ji,jj) + zfwf(jc) & 266 267 / (FLOAT(ncsnr(jc)) * e1t(ji,jj) * e2t(ji,jj) ) 267 emps(ji,jj) = emps(ji,jj) + z emp(jc) &268 emps(ji,jj) = emps(ji,jj) + zfwf(jc) & 268 269 / (FLOAT(ncsnr(jc)) * e1t(ji,jj) * e2t(ji,jj) ) 269 270 END DO … … 273 274 DO jj = ncsj1(jc), ncsj2(jc) 274 275 DO ji = ncsi1(jc), ncsi2(jc) 275 emp (ji,jj) = emp (ji,jj) - z emp(jc) / surf(jc)276 emps(ji,jj) = emps(ji,jj) - z emp(jc) / surf(jc)276 emp (ji,jj) = emp (ji,jj) - zfwf(jc) / surf(jc) 277 emps(ji,jj) = emps(ji,jj) - zfwf(jc) / surf(jc) 277 278 END DO 278 279 END DO
Note: See TracChangeset
for help on using the changeset viewer.