Changeset 12712 for NEMO

Timestamp:

2020-04-08T11:42:02+02:00 (4 years ago)

Author:

davestorkey

Message:

branches/2020/r4.0-HEAD_ticket2425 : first guess solution

File:

• NEMO/branches/2020/r4.0-HEAD_ticket2425/src/OCE/DIA/diaar5.F90 (modified) (2 diffs)

NEMO/branches/2020/r4.0-HEAD_ticket2425/src/OCE/DIA/diaar5.F90

@@ -185 +185 @@
          CALL iom_put( 'botpres', zbotpres )
          !
+      ELSE IF( iom_use('rhop') ) THEN ! we want just the density field, not steric SSH or botpres
+         CALL eos( tsn, zrhd, zrhop, gdept_n(:,:,:) )                 ! now in situ and potential density
+         zrhop(:,:,jpk) = 0._wp
+         CALL iom_put( 'rhop', zrhop )
       ENDIF
 
@@ -392 +396 @@
       IF(   iom_use( 'voltot'  ) .OR. iom_use( 'sshtot'    )  .OR. iom_use( 'sshdyn' )  .OR.  & 
          &  iom_use( 'masstot' ) .OR. iom_use( 'temptot'   )  .OR. iom_use( 'saltot' ) .OR.  &    
-         &  iom_use( 'botpres' ) .OR. iom_use( 'sshthster' )  .OR. iom_use( 'sshsteric' )  ) L_ar5 = .TRUE.
+         &  iom_use( 'botpres' ) .OR. iom_use( 'sshthster' )  .OR. iom_use( 'sshsteric' ) .OR. &
+         &  iom_use( 'rhop' )  ) L_ar5 = .TRUE.
   
       IF( l_ar5 ) THEN