Ignore:
Timestamp:
2017-09-22T17:23:46+02:00 (3 years ago)
Author:
jpalmier
Message:

JPALM — trdtrc — Silly! SBC trend is 2D but trends routines are written for 3D variables only — So make SBC trend 3D again : all 0 but surface layer

File:
1 edited

Legend:

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

    r8557 r8560  
    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,     ztrtrd ) 
     78      IF( l_trdtrc )  CALL wrk_alloc( jpi, jpj, jpk, ztrtrd ) 
    7979      ! 
    8080      zrtrn = 1.e-15_wp 
     
    140140      DO jn = 1, jptra 
    141141         ! 
    142          IF( l_trdtrc ) ztrtrd(:,:) = tra(:,:,1,jn)  ! save trends 
    143          !                                             ! add the trend to the general tracer trend 
     142         IF( l_trdtrc ) THEN 
     143            ztrtrd(:,:,:) = 0.0 
     144            ztrtrd(:,:,1) = tra(:,:,1,jn)  ! save surface trends 
     145         !                                 ! add the trend to the general tracer trend 
     146         ENDIF 
    144147 
    145148         IF ( nn_ice_tr == -1 ) THEN  ! No tracers in sea ice (null concentration in sea ice) 
     
    184187         ! 
    185188         IF( l_trdtrc ) THEN 
    186             ztrtrd(:,:) = tra(:,:,1,jn) - ztrtrd(:,:) 
     189            ztrtrd(:,:,1) = tra(:,:,1,jn) - ztrtrd(:,:,1) 
    187190            CALL trd_tra( kt, 'TRC', jn, jptra_nsr, ztrtrd ) 
    188191         END IF 
     
    207210                                           CALL prt_ctl_trc( tab4d=tra, mask=tmask, clinfo=ctrcnm, clinfo2='trd' ) 
    208211      ENDIF 
    209                       CALL wrk_dealloc( jpi, jpj, zsfx   ) 
    210       IF( l_trdtrc )  CALL wrk_dealloc( jpi, jpj, ztrtrd ) 
     212                      CALL wrk_dealloc( jpi, jpj,      zsfx   ) 
     213      IF( l_trdtrc )  CALL wrk_dealloc( jpi, jpj, jpk, ztrtrd ) 
    211214      ! 
    212215      IF( nn_timing == 1 )  CALL timing_stop('trc_sbc') 
Note: See TracChangeset for help on using the changeset viewer.