Ignore:
Timestamp:
2012-01-28T17:44:18+01:00 (9 years ago)
Author:
rblod
Message:

Merge of 3.4beta into the trunk

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/TRD/trdicp.F90

    r2781 r3294  
    2626   USE eosbn2          ! equation of state 
    2727   USE phycst          ! physical constants 
     28   USE wrk_nemo        ! Memory allocation 
     29 
    2830 
    2931   IMPLICIT NONE 
     
    106108         ! 
    107109      CASE( 'TRA' )              ! Tracers 
    108          t2(ktrd) = SUM( ptrd2dx(:,:) * e1e2t(:,:) * fse3t(:,:,1) * tn(:,:,1) ) 
    109          s2(ktrd) = SUM( ptrd2dy(:,:) * e1e2t(:,:) * fse3t(:,:,1) * sn(:,:,1) ) 
     110         t2(ktrd) = SUM( ptrd2dx(:,:) * e1e2t(:,:) * fse3t(:,:,1) * tsn(:,:,1,jp_tem) ) 
     111         s2(ktrd) = SUM( ptrd2dy(:,:) * e1e2t(:,:) * fse3t(:,:,1) * tsn(:,:,1,jp_sal) ) 
    110112         !       
    111113      END SELECT 
     
    184186         s2(ktrd) = 0._wp 
    185187         DO jk = 1, jpkm1 
    186             t2(ktrd) = t2(ktrd) + SUM( ptrd3dx(:,:,jk) * tn(:,:,jk) * e1e2t(:,:) * fse3t(:,:,jk) ) 
    187             s2(ktrd) = s2(ktrd) + SUM( ptrd3dy(:,:,jk) * sn(:,:,jk) * e1e2t(:,:) * fse3t(:,:,jk) ) 
     188            t2(ktrd) = t2(ktrd) + SUM( ptrd3dx(:,:,jk) * tsn(:,:,jk,jp_tem) * e1e2t(:,:) * fse3t(:,:,jk) ) 
     189            s2(ktrd) = s2(ktrd) + SUM( ptrd3dy(:,:,jk) * tsn(:,:,jk,jp_sal) * e1e2t(:,:) * fse3t(:,:,jk) ) 
    188190         END DO 
    189191         ! 
     
    251253      !! ** Purpose :  write dynamic trends in ocean.output  
    252254      !!---------------------------------------------------------------------- 
    253       USE wrk_nemo, ONLY:   wrk_in_use, wrk_not_released 
    254       USE wrk_nemo, ONLY:   zkepe => wrk_3d_1 , zkx => wrk_3d_2   ! 3D workspace 
    255       USE wrk_nemo, ONLY:   zky   => wrk_3d_3 , zkz => wrk_3d_4   !  -      - 
    256255      ! 
    257256      INTEGER, INTENT(in) ::   kt   ! ocean time-step index 
     
    259258      INTEGER  ::   ji, jj, jk   ! dummy loop indices 
    260259      REAL(wp) ::   zcof         ! local scalar 
    261       !!---------------------------------------------------------------------- 
    262  
    263       IF( wrk_in_use(3, 1,2,3,4) ) THEN 
    264          CALL ctl_stop('trd_dwr: requested workspace arrays unavailable')   ;   RETURN 
    265       ENDIF 
     260      REAL(wp), POINTER, DIMENSION(:,:,:)  ::  zkx, zky, zkz, zkepe   
     261      !!---------------------------------------------------------------------- 
     262 
     263      CALL wrk_alloc( jpi, jpj, jpk, zkx, zky, zkz, zkepe ) 
    266264 
    267265      ! I. Momentum trends 
     
    439437      ENDIF 
    440438      ! 
    441       IF( wrk_not_released(3, 1,2,3,4) )   CALL ctl_stop('trd_dwr: failed to release workspace arrays') 
     439      CALL wrk_dealloc( jpi, jpj, jpk, zkx, zky, zkz, zkepe ) 
    442440      ! 
    443441   END SUBROUTINE trd_dwr 
Note: See TracChangeset for help on using the changeset viewer.