New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 5208 for branches/2014/dev_r4650_UKMO11_restart_functionality/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_eiv.F90 – NEMO

Ignore:
Timestamp:
2015-04-13T15:08:59+02:00 (9 years ago)
Author:
davestorkey
Message:

Merge in changes from trunk up to 5021.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2014/dev_r4650_UKMO11_restart_functionality/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_eiv.F90

    r3787 r5208  
    2525   USE phycst          ! physical constants 
    2626   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
    27    USE diaar5, ONLY:   lk_diaar5 
    2827# endif   
    2928   USE wrk_nemo        ! Memory Allocation 
     
    161160         CALL iom_put( "voce_eiv", v_eiv )    ! j-eiv current 
    162161         CALL iom_put( "woce_eiv", w_eiv )    ! vert. eiv current 
    163          IF( lk_diaar5 ) THEN 
     162         IF( iom_use('ueiv_heattr') ) THEN 
    164163            zztmp = 0.5 * rau0 * rcp  
    165164            z2d(:,:) = 0.e0  
     
    167166               DO jj = 2, jpjm1 
    168167                  DO ji = fs_2, fs_jpim1   ! vector opt. 
    169                      z2d(ji,jj) = z2d(ji,jj) + zztmp * u_eiv(ji,jj,jk) & 
     168                     z2d(ji,jj) = z2d(ji,jj) + u_eiv(ji,jj,jk) & 
    170169                       &         * (tsn(ji,jj,jk,jp_tem)+tsn(ji+1,jj,jk,jp_tem)) * e2u(ji,jj) * fse3u(ji,jj,jk)  
    171170                  END DO 
     
    173172            END DO 
    174173            CALL lbc_lnk( z2d, 'U', -1. ) 
    175             CALL iom_put( "ueiv_heattr", z2d )                  ! heat transport in i-direction 
     174            CALL iom_put( "ueiv_heattr", zztmp * z2d )                  ! heat transport in i-direction 
     175         ENDIF 
     176             
     177         IF( iom_use('veiv_heattr') ) THEN 
     178            zztmp = 0.5 * rau0 * rcp  
    176179            z2d(:,:) = 0.e0  
    177180            DO jk = 1, jpkm1 
    178181               DO jj = 2, jpjm1 
    179182                  DO ji = fs_2, fs_jpim1   ! vector opt. 
    180                      z2d(ji,jj) = z2d(ji,jj) + zztmp * v_eiv(ji,jj,jk) & 
     183                     z2d(ji,jj) = z2d(ji,jj) + v_eiv(ji,jj,jk) & 
    181184                     &           * (tsn(ji,jj,jk,jp_tem)+tsn(ji,jj+1,jk,jp_tem)) * e1v(ji,jj) * fse3v(ji,jj,jk)  
    182185                  END DO 
     
    184187            END DO 
    185188            CALL lbc_lnk( z2d, 'V', -1. ) 
    186             CALL iom_put( "veiv_heattr", z2d )                  !  heat transport in i-direction 
     189            CALL iom_put( "veiv_heattr", zztmp * z2d )                  !  heat transport in i-direction 
    187190         ENDIF 
    188191    END IF 
Note: See TracChangeset for help on using the changeset viewer.