Ignore:
Timestamp:
2019-12-07T12:40:06+01:00 (8 months ago)
Author:
cetlod
Message:

check out & merge dev_r11643_ENHANCE-11_CEthe_Shaconemo_diags branch onto dev_r12072_TOP-01_ENHANCE-11_CEthe

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r12072_TOP-01_ENHANCE-11_CEthe/src/OCE/LDF/ldftra.F90

    r11993 r12109  
    851851      CALL iom_put( "woce_eiv", zw3d ) 
    852852      ! 
     853      IF( iom_use('weiv_masstr') ) THEN   ! vertical mass transport & its square value 
     854         zw2d(:,:) = rau0 * e1e2t(:,:) 
     855         DO jk = 1, jpk 
     856            zw3d(:,:,jk) = zw3d(:,:,jk) * zw2d(:,:) 
     857         END DO 
     858         CALL iom_put( "weiv_masstr" , zw3d )   
     859      ENDIF 
     860      ! 
     861      IF( iom_use('ueiv_masstr') ) THEN 
     862         zw3d(:,:,:) = 0.e0 
     863         DO jk = 1, jpkm1 
     864            zw3d(:,:,jk) = rau0 * ( psi_uw(:,:,jk+1) - psi_uw(:,:,jk) )  
     865         END DO 
     866         CALL iom_put( "ueiv_masstr", zw3d )                  ! mass transport in i-direction 
     867      ENDIF 
    853868      ! 
    854869      zztmp = 0.5_wp * rau0 * rcp  
     
    870885        CALL iom_put( "ueiv_heattr3d", zztmp * zw3d )                  ! heat transport in i-direction 
    871886      ENDIF 
     887      ! 
     888      IF( iom_use('veiv_masstr') ) THEN 
     889         zw3d(:,:,:) = 0.e0 
     890         DO jk = 1, jpkm1 
     891            zw3d(:,:,jk) = rau0 * ( psi_vw(:,:,jk+1) - psi_vw(:,:,jk) )  
     892         END DO 
     893         CALL iom_put( "veiv_masstr", zw3d )                  ! mass transport in i-direction 
     894      ENDIF 
     895      ! 
    872896      zw2d(:,:)   = 0._wp  
    873897      zw3d(:,:,:) = 0._wp  
     
    885909      CALL iom_put( "veiv_heattr", zztmp * zw3d )                  !  heat transport in j-direction 
    886910      ! 
    887       IF( ln_diaptr )  CALL dia_ptr_hst( jp_tem, 'eiv', 0.5 * zw3d ) 
     911      IF( iom_use( 'sophteiv' ) )   CALL dia_ptr_hst( jp_tem, 'eiv', 0.5 * zw3d ) 
    888912      ! 
    889913      zztmp = 0.5_wp * 0.5 
     
    920944      CALL iom_put( "veiv_salttr", zztmp * zw3d )                  !  salt transport in j-direction 
    921945      ! 
    922       IF( ln_diaptr ) CALL dia_ptr_hst( jp_sal, 'eiv', 0.5 * zw3d ) 
     946      IF( iom_use( 'sopsteiv' ) ) CALL dia_ptr_hst( jp_sal, 'eiv', 0.5 * zw3d ) 
    923947      ! 
    924948      ! 
Note: See TracChangeset for help on using the changeset viewer.