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 8191 for branches/UKMO/dev_r5518_GO6_package_fix_rnf_MOCI_TEST_SUITE/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_tvd.F90 – NEMO

Ignore:
Timestamp:
2017-06-19T18:16:30+02:00 (7 years ago)
Author:
andmirek
Message:

merge with XIOS restart read branch

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package_fix_rnf_MOCI_TEST_SUITE/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_tvd.F90

    r7771 r8191  
    8686      REAL(wp) ::   zfp_ui, zfp_vj, zfp_wk   !   -      - 
    8787      REAL(wp) ::   zfm_ui, zfm_vj, zfm_wk   !   -      - 
     88<<<<<<< .working 
    8889      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: zwi, zwz 
    8990      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: ztrdx, ztrdy, ztrdz, zptry 
    9091      REAL(wp), ALLOCATABLE, DIMENSION(:,:)   :: z2d 
     92======= 
     93      REAL(wp), POINTER, DIMENSION(:,:,:) :: zwi, zwz 
     94      REAL(wp), POINTER, DIMENSION(:,:,:) :: ztrdx, ztrdy, ztrdz, zptry 
     95      REAL(wp), POINTER, DIMENSION(:,:)   :: z2d 
     96>>>>>>> .merge-right.r7923 
    9197      !!---------------------------------------------------------------------- 
    9298      ! 
     
    108114      IF( cdtype == 'TRA' .AND. (iom_use("uadv_heattr") .OR. iom_use("vadv_heattr") ) ) l_trans = .TRUE. 
    109115      ! 
     116<<<<<<< .working 
    110117      IF( l_trd .OR. l_trans )  THEN 
    111118         ALLOCATE(ztrdx(1:jpi, 1:jpj, 1:jpk)) 
    112119         ALLOCATE(ztrdy(1:jpi, 1:jpj, 1:jpk)) 
    113120         ALLOCATE(ztrdz(1:jpi, 1:jpj, 1:jpk)) 
     121======= 
     122      IF( l_trd .OR. l_trans )  THEN 
     123         CALL wrk_alloc( jpi, jpj, jpk, ztrdx, ztrdy, ztrdz ) 
     124>>>>>>> .merge-right.r7923 
    114125         ztrdx(:,:,:) = 0.e0   ;    ztrdy(:,:,:) = 0.e0   ;   ztrdz(:,:,:) = 0.e0 
     126<<<<<<< .working 
    115127         ALLOCATE(z2d(1:jpi, 1:jpj)) 
    116       ENDIF 
    117       ! 
     128======= 
     129         CALL wrk_alloc( jpi, jpj, z2d ) 
     130>>>>>>> .merge-right.r7923 
     131      ENDIF 
     132      ! 
     133<<<<<<< .working 
    118134      IF( cdtype == 'TRA' .AND. ln_diaptr ) THEN   
    119135         ALLOCATE(zptry(1:jpi, 1:jpj, 1:jpk)) 
     
    121137      ENDIF 
    122138      ! 
     139======= 
     140      IF( cdtype == 'TRA' .AND. ln_diaptr ) THEN   
     141         CALL wrk_alloc( jpi, jpj, jpk, zptry ) 
     142         zptry(:,:,:) = 0._wp 
     143      ENDIF 
     144      ! 
     145>>>>>>> .merge-right.r7923 
    123146      zwi(:,:,:) = 0.e0 ;  
    124147      ! 
     
    308331      END DO 
    309332      ! 
     333<<<<<<< .working 
    310334      DEALLOCATE( zwi ) 
    311335      DEALLOCATE( zwz ) 
     
    317341      ENDIF 
    318342      IF( cdtype == 'TRA' .AND. ln_diaptr ) DEALLOCATE( zptry ) 
     343======= 
     344      CALL wrk_dealloc( jpi, jpj, jpk, zwi, zwz ) 
     345      IF( l_trd .OR. l_trans )  THEN  
     346         CALL wrk_dealloc( jpi, jpj, jpk, ztrdx, ztrdy, ztrdz ) 
     347         CALL wrk_dealloc( jpi, jpj, z2d ) 
     348      ENDIF 
     349      IF( cdtype == 'TRA' .AND. ln_diaptr ) CALL wrk_dealloc( jpi, jpj, jpk, zptry ) 
     350>>>>>>> .merge-right.r7923 
    319351      ! 
    320352      IF( nn_timing == 1 )  CALL timing_stop('tra_adv_tvd') 
     
    360392      REAL(wp) ::   zfp_ui, zfp_vj, zfp_wk   !   -      - 
    361393      REAL(wp) ::   zfm_ui, zfm_vj, zfm_wk   !   -      - 
     394<<<<<<< .working 
    362395      REAL(wp), ALLOCATABLE, DIMENSION(:,:  ) :: zwx_sav , zwy_sav 
    363396      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: zwi, zwz, zhdiv, zwz_sav, zwzts 
     
    365398      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: zptry 
    366399      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:,:) :: ztrs 
     400======= 
     401      REAL(wp), POINTER, DIMENSION(:,:  ) :: zwx_sav , zwy_sav 
     402      REAL(wp), POINTER, DIMENSION(:,:,:) :: zwi, zwz, zhdiv, zwz_sav, zwzts 
     403      REAL(wp), POINTER, DIMENSION(:,:,:) :: ztrdx, ztrdy, ztrdz 
     404      REAL(wp), POINTER, DIMENSION(:,:,:) :: zptry 
     405      REAL(wp), POINTER, DIMENSION(:,:,:,:) :: ztrs 
     406>>>>>>> .merge-right.r7923 
    367407      !!---------------------------------------------------------------------- 
    368408      ! 
     
    394434      ENDIF 
    395435      ! 
     436<<<<<<< .working 
    396437      IF( cdtype == 'TRA' .AND. ln_diaptr ) THEN   
    397438         ALLOCATE(zptry(1:jpi, 1:jpj, 1:jpk))        
     
    399440      ENDIF 
    400441      ! 
     442======= 
     443      IF( cdtype == 'TRA' .AND. ln_diaptr ) THEN   
     444         CALL wrk_alloc( jpi, jpj,jpk, zptry ) 
     445         zptry(:,:,:) = 0._wp 
     446      ENDIF 
     447      ! 
     448>>>>>>> .merge-right.r7923 
    401449      zwi(:,:,:) = 0._wp 
    402450      z_rzts = 1._wp / REAL( jnzts, wp ) 
     
    633681 
    634682      IF( cdtype == 'TRA' .AND. ln_diaptr ) DEALLOCATE(zptry ) 
     683      IF( cdtype == 'TRA' .AND. ln_diaptr ) CALL wrk_dealloc( jpi, jpj, jpk, zptry ) 
    635684      ! 
    636685      IF( nn_timing == 1 )  CALL timing_stop('tra_adv_tvd_zts') 
Note: See TracChangeset for help on using the changeset viewer.