Ignore:
Timestamp:
2019-09-30T11:07:57+02:00 (14 months ago)
Author:
davestorkey
Message:

UKMO/NEMO_4.0_momentum_trends branch : first set of code changes.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0_momentum_trends/src/OCE/DYN/dynspg.F90

    r10888 r11613  
    8181      IF( ln_timing )   CALL timing_start('dyn_spg') 
    8282      ! 
    83       IF( l_trddyn )   THEN                      ! temporary save of ta and sa trends 
     83      IF( l_trddyn .AND. nspg == np_EXP )   THEN                      ! temporary save of ta and sa trends 
    8484         ALLOCATE( ztrdu(jpi,jpj,jpk) , ztrdv(jpi,jpj,jpk) )  
    8585         ztrdu(:,:,:) = ua(:,:,:) 
     
    165165      END SELECT 
    166166      !                     
    167       IF( l_trddyn )   THEN                  ! save the surface pressure gradient trends for further diagnostics 
     167      IF( l_trddyn .AND. nspg == np_EXP )   THEN     ! save the surface pressure gradient trends for further diagnostics 
    168168         ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:) 
    169169         ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:) 
    170170         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 ) 
    171175         DEALLOCATE( ztrdu , ztrdv )  
    172176      ENDIF 
Note: See TracChangeset for help on using the changeset viewer.