Changeset 8077
- Timestamp:
- 2017-05-26T11:13:07+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/OPA_SRC/DIA/diaar5.F90
r7515 r8077 206 206 ! Exclude points where rn2 is negative as convection kicks in here and 207 207 ! work is not being done against stratification 208 pe(:,:) = 0._wp 209 IF( lk_zdfddm ) THEN 210 DO ji=1,jpi 211 DO jj=1,jpj 212 DO jk=1,jpk 213 zrw = ( fsdepw(ji,jj,jk ) - fsdept(ji,jj,jk) ) & 214 & / ( fsdept(ji,jj,jk-1) - fsdept(ji,jj,jk) ) 215 ! 216 zaw = rab_n(ji,jj,jk,jp_tem) * (1. - zrw) + rab_n(ji,jj,jk-1,jp_tem)* zrw 217 zbw = rab_n(ji,jj,jk,jp_sal) * (1. - zrw) + rab_n(ji,jj,jk-1,jp_sal)* zrw 218 ! 219 pe(ji, jj) = pe(ji, jj) - MIN(0._wp, rn2(ji,jj,jk)) * & 220 & grav * (avt(ji,jj,jk) * zaw * (tsn(ji,jj,jk-1,jp_tem) - tsn(ji,jj,jk,jp_tem) ) & 221 & - fsavs(ji,jj,jk) * zbw * (tsn(ji,jj,jk-1,jp_sal) - tsn(ji,jj,jk,jp_sal) ) ) 222 223 ENDDO 224 ENDDO 225 ENDDO 208 pe(:,:) = 0._wp 209 IF( ln_zdfddm ) THEN 210 DO jk = 2, jpk 211 DO jj = 1, jpj 212 DO ji = 1, jpi 213 IF( rn2(ji,jj,jk) > 0._wp ) THEN 214 zrw = ( gdepw_n(ji,jj,jk ) - gdept_n(ji,jj,jk) ) & 215 & / ( gdept_n(ji,jj,jk-1) - gdept_n(ji,jj,jk) ) 216 ! 217 zaw = rab_n(ji,jj,jk,jp_tem) * (1. - zrw) + rab_n(ji,jj,jk-1,jp_tem)* zrw 218 zbw = rab_n(ji,jj,jk,jp_sal) * (1. - zrw) + rab_n(ji,jj,jk-1,jp_sal)* zrw 219 ! 220 zpe(ji, jj) = zpe(ji, jj) & 221 & - grav * ( avt(ji,jj,jk) * zaw * (tsn(ji,jj,jk-1,jp_tem) - tsn(ji,jj,jk,jp_tem) ) & 222 & - fsavs(ji,jj,jk) * zbw * (tsn(ji,jj,jk-1,jp_sal) - tsn(ji,jj,jk,jp_sal) ) ) 223 ENDIF 224 END DO 225 END DO 226 END DO 226 227 ELSE 227 DO ji=1,jpi228 DO jj=1,jpj229 DO jk=1,jpk230 pe(ji,jj) = pe(ji,jj) + avt(ji, jj, jk) * MIN(0._wp,rn2(ji, jj, jk)) * rau0 * fse3w(ji, jj, jk)231 ENDDO232 ENDDO233 ENDDO234 235 236 228 DO jk = 1, jpk 229 DO ji = 1, jpi 230 DO jj = 1, jpj 231 zpe(ji,jj) = zpe(ji,jj) + avt(ji, jj, jk) * MIN(0._wp,rn2(ji, jj, jk)) * rau0 * e3w_n(ji, jj, jk) 232 END DO 233 END DO 234 END DO 235 ENDIF 236 CALL lbc_lnk(pe, 'T', 1._wp) 237 CALL iom_put( 'tnpeo', pe ) 237 238 ENDIF 238 239 !
Note: See TracChangeset
for help on using the changeset viewer.