Changeset 8193 for branches/UKMO/dev_r5518_GO6_package_fix_rnf_MOCI_TEST_SUITE/NEMOGCM/NEMO/OPA_SRC/TRA/traadv.F90
- Timestamp:
- 2017-06-20T10:46:27+02:00 (7 years ago)
- 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 82 82 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: zun, zvn, zwn 83 83 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: ztrdt, ztrds ! 3D workspace 84 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztrdt, ztrds ! 3D workspace85 84 !!---------------------------------------------------------------------- 86 85 ! … … 126 125 IF( ln_diaptr ) CALL dia_ptr( zvn ) ! diagnose the effective MSF 127 126 ! 128 <<<<<<< .working129 127 IF( l_trdtra ) THEN !* Save ta and sa trends 130 128 ALLOCATE(ztrdt( 1:jpi, 1:jpj, 1:jpk) ) … … 134 132 ENDIF 135 133 ! 136 =======137 IF( l_trdtra ) THEN !* Save ta and sa trends138 CALL wrk_alloc( jpi, jpj, jpk, ztrdt, ztrds )139 ztrdt(:,:,:) = tsa(:,:,:,jp_tem)140 ztrds(:,:,:) = tsa(:,:,:,jp_sal)141 ENDIF142 !143 >>>>>>> .merge-right.r7923144 134 SELECT CASE ( nadv ) !== compute advection trend and add it to general trend ==! 145 135 CASE ( 1 ) ; CALL tra_adv_cen2 ( kt, nit000, 'TRA', zun, zvn, zwn, tsb, tsn, tsa, jpts ) ! 2nd order centered … … 172 162 END SELECT 173 163 ! 174 <<<<<<< .working175 164 IF( l_trdtra ) THEN ! save the advective trends for further diagnostics 176 165 DO jk = 1, jpkm1 … … 183 172 DEALLOCATE (ztrds) 184 173 ENDIF 185 =======186 IF( l_trdtra ) THEN ! save the advective trends for further diagnostics187 DO jk = 1, jpkm1188 ztrdt(:,:,jk) = tsa(:,:,jk,jp_tem) - ztrdt(:,:,jk)189 ztrds(:,:,jk) = tsa(:,:,jk,jp_sal) - ztrds(:,:,jk)190 END DO191 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 ENDIF195 >>>>>>> .merge-right.r7923196 174 ! ! print mean trends (used for debugging) 197 175 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.