- Timestamp:
- 2019-11-20T10:47:16+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0.1_momentum_trends/src/OCE/DYN/dynspg.F90
r11715 r11934 81 81 IF( ln_timing ) CALL timing_start('dyn_spg') 82 82 ! 83 IF( l_trddyn ) THEN ! temporary save of ta and sa trends83 IF( l_trddyn .AND. nspg == np_EXP ) THEN ! temporary save of ta and sa trends 84 84 ALLOCATE( ztrdu(jpi,jpj,jpk) , ztrdv(jpi,jpj,jpk) ) 85 85 ztrdu(:,:,:) = ua(:,:,:) … … 165 165 END SELECT 166 166 ! 167 IF( l_trddyn ) THEN! save the surface pressure gradient trends for further diagnostics167 IF( l_trddyn .AND. nspg == np_EXP ) THEN ! save the surface pressure gradient trends for further diagnostics 168 168 ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:) 169 169 ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:) 170 170 CALL trd_dyn( ztrdu, ztrdv, jpdyn_spg, kt ) 171 ! In this case we also need to finalise the write-out of the PVO 3D diagnostic with zero correction 172 ztrdu(:,:,1) = 0._wp 173 ztrdv(:,:,1) = 0._wp 174 CALL trd_dyn( ztrdu(:,:,1), ztrdv(:,:,1), jpdyn_pvo, kt ) 171 175 DEALLOCATE( ztrdu , ztrdv ) 172 176 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.