Ignore:
Timestamp:
2017-09-22T13:06:31+02:00 (4 years ago)
Author:
jpalmier
Message:

JPALM — trdtrc — fix SBC trend (2D not 3D) — remove DMP from XIOS group to avoid Na Ns?

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/NERC/dev_r5518_GO6_package_trdtrc/NEMOGCM/NEMO/TOP_SRC/TRP/trcsbc.F90

    r8356 r8557  
    6868      CHARACTER (len=22) :: charout 
    6969      REAL(wp), POINTER, DIMENSION(:,:  ) :: zsfx 
    70       REAL(wp), POINTER, DIMENSION(:,:,:) :: ztrtrd 
     70      REAL(wp), POINTER, DIMENSION(:,:  ) :: ztrtrd 
    7171 
    7272      !!--------------------------------------------------------------------- 
     
    7676      ! Allocate temporary workspace 
    7777                      CALL wrk_alloc( jpi, jpj,      zsfx   ) 
    78       IF( l_trdtrc )  CALL wrk_alloc( jpi, jpj, jpk, ztrtrd ) 
     78      IF( l_trdtrc )  CALL wrk_alloc( jpi, jpj,     ztrtrd ) 
    7979      ! 
    8080      zrtrn = 1.e-15_wp 
     
    140140      DO jn = 1, jptra 
    141141         ! 
    142          IF( l_trdtrc ) ztrtrd(:,:,:) = tra(:,:,:,jn)  ! save trends 
     142         IF( l_trdtrc ) ztrtrd(:,:) = tra(:,:,1,jn)  ! save trends 
    143143         !                                             ! add the trend to the general tracer trend 
    144144 
     
    184184         ! 
    185185         IF( l_trdtrc ) THEN 
    186             ztrtrd(:,:,:) = tra(:,:,:,jn) - ztrtrd(:,:,:) 
     186            ztrtrd(:,:) = tra(:,:,1,jn) - ztrtrd(:,:) 
    187187            CALL trd_tra( kt, 'TRC', jn, jptra_nsr, ztrtrd ) 
    188188         END IF 
     
    207207                                           CALL prt_ctl_trc( tab4d=tra, mask=tmask, clinfo=ctrcnm, clinfo2='trd' ) 
    208208      ENDIF 
    209                       CALL wrk_dealloc( jpi, jpj,      zsfx   ) 
    210       IF( l_trdtrc )  CALL wrk_dealloc( jpi, jpj, jpk, ztrtrd ) 
     209                      CALL wrk_dealloc( jpi, jpj, zsfx   ) 
     210      IF( l_trdtrc )  CALL wrk_dealloc( jpi, jpj, ztrtrd ) 
    211211      ! 
    212212      IF( nn_timing == 1 )  CALL timing_stop('trc_sbc') 
Note: See TracChangeset for help on using the changeset viewer.