- Timestamp:
- 2017-12-19T09:47:17+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/TRP/trdmxl_trc.F90
r9050 r9125 31 31 USE prtctl ! print control 32 32 USE sms_pisces ! PISCES bio-model 33 USE wrk_nemo ! Memory allocation34 33 35 34 IMPLICIT NONE … … 97 96 ! 98 97 INTEGER :: ji, jj, jk, isum 99 REAL(wp), POINTER, DIMENSION(:,:) :: zvlmsk 100 !!---------------------------------------------------------------------- 101 102 CALL wrk_alloc( jpi, jpj, zvlmsk ) 98 REAL(wp), DIMENSION(jpi,jpj) :: zvlmsk 99 !!---------------------------------------------------------------------- 103 100 104 101 ! I. Definition of control surface and integration weights … … 182 179 tmltrd_trc(:,:,ktrd,kjn) = tmltrd_trc(:,:,ktrd,kjn) + ptrc_trdmxl(:,:,1) * wkx_trc(:,:,1) ! non penetrative 183 180 END SELECT 184 !185 CALL wrk_dealloc( jpi, jpj, zvlmsk )186 181 ! 187 182 END SUBROUTINE trd_mxl_trc_zint … … 241 236 REAL(wp) :: zavt, zfn, zfn2 242 237 ! 243 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztmltot ! d(trc)/dt over the anlysis window (incl. Asselin)244 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztmlres ! residual = dh/dt entrainment term245 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztmlatf ! for storage only246 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztmlrad ! for storage only (for trb<0 corr in trcrad)247 ! 248 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztmltot2 ! -+249 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztmlres2 ! | working arrays to diagnose the trends250 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztmltrdm2 ! | associated with the time meaned ML251 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztmlatf2 ! | passive tracers252 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztmlrad2 ! | (-> for trb<0 corr in trcrad)238 REAL(wp), DIMENSION(jpi,jpj,jptra) :: ztmltot ! d(trc)/dt over the anlysis window (incl. Asselin) 239 REAL(wp), DIMENSION(jpi,jpj,jptra) :: ztmlres ! residual = dh/dt entrainment term 240 REAL(wp), DIMENSION(jpi,jpj,jptra) :: ztmlatf ! for storage only 241 REAL(wp), DIMENSION(jpi,jpj,jptra) :: ztmlrad ! for storage only (for trb<0 corr in trcrad) 242 ! 243 REAL(wp), DIMENSION(jpi,jpj,jptra) :: ztmltot2 ! -+ 244 REAL(wp), DIMENSION(jpi,jpj,jptra) :: ztmlres2 ! | working arrays to diagnose the trends 245 REAL(wp), DIMENSION(jpi,jpj,jptra) :: ztmltrdm2 ! | associated with the time meaned ML 246 REAL(wp), DIMENSION(jpi,jpj,jptra) :: ztmlatf2 ! | passive tracers 247 REAL(wp), DIMENSION(jpi,jpj,jptra) :: ztmlrad2 ! | (-> for trb<0 corr in trcrad) 253 248 ! 254 249 CHARACTER (LEN=10) :: clvar 255 250 !!---------------------------------------------------------------------- 256 251 257 ! Set-up pointers into sub-arrays of workspaces258 CALL wrk_alloc( jpi, jpj, jptra, ztmltot , ztmlres , ztmlatf , ztmlrad )259 CALL wrk_alloc( jpi, jpj, jptra, ztmltot2, ztmlres2, ztmlatf2, ztmlrad2, ztmltrdm2 )260 252 261 253 IF( nn_dttrc /= 1 ) CALL ctl_stop( " Be careful, trends diags never validated " ) … … 734 726 735 727 IF( lrst_trc ) CALL trd_mxl_trc_rst_write( kt ) ! this must be after the array swap above (III.3) 736 737 CALL wrk_dealloc( jpi, jpj, jptra, ztmltot , ztmlres , ztmlatf , ztmlrad )738 CALL wrk_dealloc( jpi, jpj, jptra, ztmltot2, ztmlres2, ztmlatf2, ztmlrad2, ztmltrdm2 )739 728 ! 740 729 END SUBROUTINE trd_mxl_trc
Note: See TracChangeset
for help on using the changeset viewer.