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

Ignore:
Timestamp:
2017-06-20T10:46:27+02:00 (7 years ago)
Author:
andmirek
Message:

ticket #1914 back to 8188

File:
1 edited

Legend:

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

    r8191 r8193  
    8282      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: zun, zvn, zwn 
    8383      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: ztrdt, ztrds   ! 3D workspace 
    84       REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztrdt, ztrds   ! 3D workspace 
    8584      !!---------------------------------------------------------------------- 
    8685      ! 
     
    126125      IF( ln_diaptr )   CALL dia_ptr( zvn )                                     ! diagnose the effective MSF  
    127126      ! 
    128 <<<<<<< .working 
    129127      IF( l_trdtra )   THEN                    !* Save ta and sa trends 
    130128         ALLOCATE(ztrdt( 1:jpi, 1:jpj, 1:jpk) ) 
     
    134132      ENDIF 
    135133      ! 
    136 ======= 
    137       IF( l_trdtra )   THEN                    !* Save ta and sa trends 
    138          CALL wrk_alloc( jpi, jpj, jpk, ztrdt, ztrds ) 
    139          ztrdt(:,:,:) = tsa(:,:,:,jp_tem) 
    140          ztrds(:,:,:) = tsa(:,:,:,jp_sal) 
    141       ENDIF 
    142       ! 
    143 >>>>>>> .merge-right.r7923 
    144134      SELECT CASE ( nadv )                            !==  compute advection trend and add it to general trend  ==! 
    145135      CASE ( 1 )   ;    CALL tra_adv_cen2   ( kt, nit000, 'TRA',         zun, zvn, zwn, tsb, tsn, tsa, jpts )   !  2nd order centered 
     
    172162      END SELECT 
    173163      ! 
    174 <<<<<<< .working 
    175164      IF( l_trdtra )   THEN                      ! save the advective trends for further diagnostics 
    176165         DO jk = 1, jpkm1 
     
    183172         DEALLOCATE (ztrds) 
    184173      ENDIF 
    185 ======= 
    186       IF( l_trdtra )   THEN                      ! save the advective trends for further diagnostics 
    187          DO jk = 1, jpkm1 
    188             ztrdt(:,:,jk) = tsa(:,:,jk,jp_tem) - ztrdt(:,:,jk) 
    189             ztrds(:,:,jk) = tsa(:,:,jk,jp_sal) - ztrds(:,:,jk) 
    190          END DO 
    191          CALL trd_tra( kt, 'TRA', jp_tem, jptra_totad, ztrdt ) 
    192          CALL trd_tra( kt, 'TRA', jp_sal, jptra_totad, ztrds ) 
    193          CALL wrk_dealloc( jpi, jpj, jpk, ztrdt, ztrds ) 
    194       ENDIF 
    195 >>>>>>> .merge-right.r7923 
    196174      !                                              ! print mean trends (used for debugging) 
    197175      IF(ln_ctl)   CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' adv  - Ta: ', mask1=tmask,               & 
Note: See TracChangeset for help on using the changeset viewer.