- Timestamp:
- 2016-11-08T16:30:34+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/v3_6_extra_CMIP6_diagnostics/NEMOGCM/NEMO/OPA_SRC/DIA/diaar5.F90
r7063 r7209 196 196 CALL iom_put( 'temptot', ztemp ) 197 197 CALL iom_put( 'saltot' , zsal ) 198 ! 199 CALL wrk_dealloc( jpi , jpj , zarea_ssh , zbotpres ) 198 199 IF( iom_use( 'tnpeo' )) THEN 200 ! Work done against stratification by vertical mixing 201 ! Exclude points where rn2 is negative as convection kicks in here and 202 ! work is not being done against stratification 203 pe(:,:) = 0._wp 204 IF( lk_zdfddm ) THEN 205 DO ji=1,jpi 206 DO jj=1,jpj 207 DO jk=1,jpk 208 zrw = ( fsdepw(ji,jj,jk ) - fsdept(ji,jj,jk) ) & 209 & / ( fsdept(ji,jj,jk-1) - fsdept(ji,jj,jk) ) 210 ! 211 zaw = rab_n(ji,jj,jk,jp_tem) * (1. - zrw) + rab_n(ji,jj,jk-1,jp_tem)* zrw 212 zbw = rab_n(ji,jj,jk,jp_sal) * (1. - zrw) + rab_n(ji,jj,jk-1,jp_sal)* zrw 213 ! 214 pe(ji, jj) = pe(ji, jj) - MIN(0._wp, rn2(ji,jj,jk)) * & 215 & grav * (avt(ji,jj,jk) * zaw * (tsn(ji,jj,jk-1,jp_tem) - tsn(ji,jj,jk,jp_tem) ) & 216 & - fsavs(ji,jj,jk) * zbw * (tsn(ji,jj,jk-1,jp_sal) - tsn(ji,jj,jk,jp_sal) ) ) 217 218 ENDDO 219 ENDDO 220 ENDDO 221 ELSE 222 DO ji=1,jpi 223 DO jj=1,jpj 224 DO jk=1,jpk 225 pe(ji,jj) = pe(ji,jj) + avt(ji, jj, jk) * MIN(0._wp,rn2(ji, jj, jk)) * rau0 * fse3w(ji, jj, jk) 226 ENDDO 227 ENDDO 228 ENDDO 229 ENDIF 230 CALL lbc_lnk(pe, 'T', 1._wp) 231 CALL iom_put( 'tnpeo', pe ) 232 ENDIF 233 ! 234 CALL wrk_dealloc( jpi , jpj , zarea_ssh , zbotpres, pe ) 200 235 CALL wrk_dealloc( jpi , jpj , jpk , zrhd , zrhop ) 201 236 CALL wrk_dealloc( jpi , jpj , jpk , jpts , ztsn )
Note: See TracChangeset
for help on using the changeset viewer.