Ignore:
Timestamp:
2020-07-01T15:42:06+02:00 (17 months ago)
Author:
smasson
Message:

better e3: update with trunk@13136 see #2385

Location:
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3

    • Property svn:externals
      •  

        old new  
        88 
        99# SETTE 
        10 ^/utils/CI/sette@HEAD         sette 
         10^/utils/CI/sette@12931        sette 
  • NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/DIA/diaar5.F90

    r13156 r13193  
    7676      REAL(wp) ::   zaw, zbw, zrw 
    7777      ! 
    78       REAL(wp), ALLOCATABLE, DIMENSION(:,:)     :: zarea_ssh , zbotpres       ! 2D workspace  
     78      REAL(wp), ALLOCATABLE, DIMENSION(:,:)     :: zarea_ssh, zbotpres        ! 2D workspace  
    7979      REAL(wp), ALLOCATABLE, DIMENSION(:,:)     :: z2d, zpe                   ! 2D workspace  
    80       REAL(wp), ALLOCATABLE, DIMENSION(:,:,:)   :: z3d, zrhd , zrhop, ztpot, zgdept   ! 3D workspace (zgdept: needed to use the substitute) 
     80      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:)   :: z3d, zrhd, ztpot, zgdept   ! 3D workspace (zgdept: needed to use the substitute) 
    8181      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:,:) :: ztsn                       ! 4D workspace 
    8282 
     
    8888      IF( l_ar5 ) THEN  
    8989         ALLOCATE( zarea_ssh(jpi,jpj), zbotpres(jpi,jpj), z2d(jpi,jpj) ) 
    90          ALLOCATE( zrhd(jpi,jpj,jpk) , zrhop(jpi,jpj,jpk) ) 
     90         ALLOCATE( zrhd(jpi,jpj,jpk) ) 
    9191         ALLOCATE( ztsn(jpi,jpj,jpk,jpts) ) 
    9292         zarea_ssh(:,:) = e1e2t(:,:) * ssh(:,:,Kmm) 
     
    163163       
    164164         !                                         ! steric sea surface height 
    165          CALL eos( ts(:,:,:,:,Kmm), zrhd, zrhop, zgdept )                 ! now in situ and potential density 
    166          zrhop(:,:,jpk) = 0._wp 
    167          CALL iom_put( 'rhop', zrhop ) 
    168          ! 
    169165         zbotpres(:,:) = 0._wp                        ! no atmospheric surface pressure, levitating sea-ice 
    170166         DO jk = 1, jpkm1 
    171             zbotpres(:,:) = zbotpres(:,:) + e3t(:,:,jk,Kmm) * zrhd(:,:,jk) 
     167            zbotpres(:,:) = zbotpres(:,:) + e3t(:,:,jk,Kmm) * rhd(:,:,jk) 
    172168         END DO 
    173169         IF( ln_linssh ) THEN 
     
    176172                  DO jj = 1,jpj 
    177173                     iks = mikt(ji,jj) 
    178                      zbotpres(ji,jj) = zbotpres(ji,jj) + ssh(ji,jj,Kmm) * zrhd(ji,jj,iks) + riceload(ji,jj) 
     174                     zbotpres(ji,jj) = zbotpres(ji,jj) + ssh(ji,jj,Kmm) * rhd(ji,jj,iks) + riceload(ji,jj) 
    179175                  END DO 
    180176               END DO 
    181177            ELSE 
    182                zbotpres(:,:) = zbotpres(:,:) + ssh(:,:,Kmm) * zrhd(:,:,1) 
     178               zbotpres(:,:) = zbotpres(:,:) + ssh(:,:,Kmm) * rhd(:,:,1) 
    183179            END IF 
    184180         END IF 
     
    303299      IF( l_ar5 ) THEN 
    304300        DEALLOCATE( zarea_ssh , zbotpres, z2d ) 
    305         DEALLOCATE( zrhd      , zrhop    ) 
    306301        DEALLOCATE( ztsn                 ) 
    307302      ENDIF 
     
    377372      IF(   iom_use( 'voltot'  ) .OR. iom_use( 'sshtot'    )  .OR. iom_use( 'sshdyn' )  .OR.  &  
    378373         &  iom_use( 'masstot' ) .OR. iom_use( 'temptot'   )  .OR. iom_use( 'saltot' ) .OR.  &     
    379          &  iom_use( 'botpres' ) .OR. iom_use( 'sshthster' )  .OR. iom_use( 'sshsteric' )  ) L_ar5 = .TRUE. 
     374         &  iom_use( 'botpres' ) .OR. iom_use( 'sshthster' )  .OR. iom_use( 'sshsteric' ) .OR. & 
     375         &  iom_use( 'rhop' )  ) L_ar5 = .TRUE. 
    380376   
    381377      IF( l_ar5 ) THEN 
Note: See TracChangeset for help on using the changeset viewer.