- Timestamp:
- 2020-06-10T10:59:29+02:00 (4 years ago)
- Location:
- NEMO/releases/r4.0/r4.0-HEAD/src/OCE/DIA
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/releases/r4.0/r4.0-HEAD/src/OCE/DIA/diaar5.F90
r12631 r13085 156 156 157 157 ! ! steric sea surface height 158 CALL eos( tsn, zrhd, zrhop, gdept_n(:,:,:) ) ! now in situ and potential density159 zrhop(:,:,jpk) = 0._wp160 CALL iom_put( 'rhop', zrhop )161 !162 158 zbotpres(:,:) = 0._wp ! no atmospheric surface pressure, levitating sea-ice 163 159 DO jk = 1, jpkm1 164 zbotpres(:,:) = zbotpres(:,:) + e3t_n(:,:,jk) * zrhd(:,:,jk)160 zbotpres(:,:) = zbotpres(:,:) + e3t_n(:,:,jk) * rhd(:,:,jk) 165 161 END DO 166 162 IF( ln_linssh ) THEN … … 169 165 DO jj = 1,jpj 170 166 iks = mikt(ji,jj) 171 zbotpres(ji,jj) = zbotpres(ji,jj) + sshn(ji,jj) * zrhd(ji,jj,iks) + riceload(ji,jj)167 zbotpres(ji,jj) = zbotpres(ji,jj) + sshn(ji,jj) * rhd(ji,jj,iks) + riceload(ji,jj) 172 168 END DO 173 169 END DO 174 170 ELSE 175 zbotpres(:,:) = zbotpres(:,:) + sshn(:,:) * zrhd(:,:,1)171 zbotpres(:,:) = zbotpres(:,:) + sshn(:,:) * rhd(:,:,1) 176 172 END IF 177 173 END IF … … 392 388 IF( iom_use( 'voltot' ) .OR. iom_use( 'sshtot' ) .OR. iom_use( 'sshdyn' ) .OR. & 393 389 & iom_use( 'masstot' ) .OR. iom_use( 'temptot' ) .OR. iom_use( 'saltot' ) .OR. & 394 & iom_use( 'botpres' ) .OR. iom_use( 'sshthster' ) .OR. iom_use( 'sshsteric' ) ) L_ar5 = .TRUE. 390 & iom_use( 'botpres' ) .OR. iom_use( 'sshthster' ) .OR. iom_use( 'sshsteric' ) .OR. & 391 & iom_use( 'rhop' ) ) L_ar5 = .TRUE. 395 392 396 393 IF( l_ar5 ) THEN -
NEMO/releases/r4.0/r4.0-HEAD/src/OCE/DIA/diawri.F90
r12494 r13085 167 167 CALL iom_put( "sbs", z2d ) ! bottom salinity 168 168 ENDIF 169 170 CALL iom_put( "rhop", rhop(:,:,:) ) ! 3D potential density (sigma0) 169 171 170 172 IF ( iom_use("taubot") ) THEN ! bottom stress
Note: See TracChangeset
for help on using the changeset viewer.