- Timestamp:
- 2020-04-09T16:57:47+02:00 (4 years ago)
- Location:
- NEMO/branches/2020/r4.0-HEAD_ticket2425/src/OCE/DIA
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/r4.0-HEAD_ticket2425/src/OCE/DIA/diaar5.F90
r12712 r12727 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 … … 185 181 CALL iom_put( 'botpres', zbotpres ) 186 182 ! 187 ELSE IF( iom_use('rhop') ) THEN ! we want just the density field, not steric SSH or botpres188 CALL eos( tsn, zrhd, zrhop, gdept_n(:,:,:) ) ! now in situ and potential density189 zrhop(:,:,jpk) = 0._wp190 CALL iom_put( 'rhop', zrhop )191 183 ENDIF 192 184 -
NEMO/branches/2020/r4.0-HEAD_ticket2425/src/OCE/DIA/diawri.F90
r12494 r12727 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.