Changeset 13085


Ignore:
Timestamp:
2020-06-10T10:59:29+02:00 (4 months ago)
Author:
davestorkey
Message:

r4.0-HEAD : fix for ticket #2425.

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  
    156156       
    157157         !                                         ! steric sea surface height 
    158          CALL eos( tsn, zrhd, zrhop, gdept_n(:,:,:) )                 ! now in situ and potential density 
    159          zrhop(:,:,jpk) = 0._wp 
    160          CALL iom_put( 'rhop', zrhop ) 
    161          ! 
    162158         zbotpres(:,:) = 0._wp                        ! no atmospheric surface pressure, levitating sea-ice 
    163159         DO jk = 1, jpkm1 
    164             zbotpres(:,:) = zbotpres(:,:) + e3t_n(:,:,jk) * zrhd(:,:,jk) 
     160            zbotpres(:,:) = zbotpres(:,:) + e3t_n(:,:,jk) * rhd(:,:,jk) 
    165161         END DO 
    166162         IF( ln_linssh ) THEN 
     
    169165                  DO jj = 1,jpj 
    170166                     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) 
    172168                  END DO 
    173169               END DO 
    174170            ELSE 
    175                zbotpres(:,:) = zbotpres(:,:) + sshn(:,:) * zrhd(:,:,1) 
     171               zbotpres(:,:) = zbotpres(:,:) + sshn(:,:) * rhd(:,:,1) 
    176172            END IF 
    177173         END IF 
     
    392388      IF(   iom_use( 'voltot'  ) .OR. iom_use( 'sshtot'    )  .OR. iom_use( 'sshdyn' )  .OR.  &  
    393389         &  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. 
    395392   
    396393      IF( l_ar5 ) THEN 
  • NEMO/releases/r4.0/r4.0-HEAD/src/OCE/DIA/diawri.F90

    r12494 r13085  
    167167         CALL iom_put( "sbs", z2d )                ! bottom salinity 
    168168      ENDIF 
     169 
     170      CALL iom_put( "rhop", rhop(:,:,:) )          ! 3D potential density (sigma0) 
    169171 
    170172      IF ( iom_use("taubot") ) THEN                ! bottom stress 
Note: See TracChangeset for help on using the changeset viewer.