New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 8470 – NEMO

Changeset 8470


Ignore:
Timestamp:
2017-08-29T16:18:37+02:00 (7 years ago)
Author:
jpalmier
Message:

JPALM -- 29-08-2017 -- trends bug fix TOT, ZDF, ZDP

Location:
branches/NERC/dev_r5518_GO6_package_trdtrc/NEMOGCM/NEMO/TOP_SRC/TRP
Files:
3 edited

Legend:

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

    r8453 r8470  
    130130         ztrdt(:,:,jpk,:) = 0._wp 
    131131         IF( ln_trcldf_iso ) THEN                       ! diagnose the "pure" Kz diffusive trend  
    132             CALL trd_tra( kt, 'TRC', jn, jptra_zdfp, ztrdt ) 
     132            DO jn = 1, jptra 
     133               CALL trd_tra( kt, 'TRC', jn, jptra_zdfp, ztrdt(:,:,:,jn) ) 
     134            ENDDO 
    133135         ENDIF 
    134136         ! total trend for the non-time-filtered variables. 
     
    151153            END DO 
    152154         END IF 
    153          CALL trd_tra( kt, 'TRC', jn, jptra_tot, ztrdt ) 
     155         DO jn = 1, jptra 
     156            CALL trd_tra( kt, 'TRC', jn, jptra_tot, ztrdt(:,:,:,jn) ) 
     157         ENDDO 
    154158         IF( .NOT.lk_vvl )  THEN 
    155159            ! Store now fields before applying the Asselin filter  
     
    172176                                            ! this time step 
    173177            ztrdt(:,:,:,:) = 0._wp 
    174             CALL trd_tra( kt, 'TRC', jn, jptra_atf, ztrdt ) 
     178            DO jn = 1, jptra 
     179               CALL trd_tra( kt, 'TRC', jn, jptra_atf, ztrdt(:,:,:,jn) ) 
     180            ENDDO 
    175181         END IF 
    176182 
     
    191197               ztrdt(:,:,jk,jn) = ( trb(:,:,jk,jn) - ztrdt(:,:,jk,jn) ) * zfact  
    192198            END DO 
    193             CALL trd_tra( kt, 'TRC', jn, jptra_atf, ztrdt ) 
     199            CALL trd_tra( kt, 'TRC', jn, jptra_atf, ztrdt(:,:,:,jn) ) 
    194200         END DO 
    195201      END IF 
  • branches/NERC/dev_r5518_GO6_package_trdtrc/NEMOGCM/NEMO/TOP_SRC/TRP/trczdf.F90

    r8453 r8470  
    103103            DO jn = 1, jptra 
    104104               DO jk = 1, jpkm1 
    105                   ztrtrd(:,:,jk) = ( ( tra(:,:,jk,jn)*fse3t_a(:,:,jk) - & 
    106                                         trn(:,:,jk,jn)*fse3t_b(:,:,jk) ) & 
    107                                       / (fse3t_n(:,:,jk)*r2dtra(jk)) ) - ztrtrd(:,:,jk) 
     105                  ztrtrd(:,:,jk,jn) = ( ( tra(:,:,jk,jn)*fse3t_a(:,:,jk) - & 
     106                                          trn(:,:,jk,jn)*fse3t_b(:,:,jk) ) & 
     107                                       / (fse3t_n(:,:,jk)*r2dtra(jk)) ) - ztrtrd(:,:,jk,jn) 
    108108               END DO 
    109109            END DO 
  • branches/NERC/dev_r5518_GO6_package_trdtrc/NEMOGCM/NEMO/TOP_SRC/TRP/trdtrc.F90

    r8467 r8470  
    187187 
    188188 
    189                            IF  (iom_use(cltra)) THEN 
     189      IF  (iom_use(cltra)) THEN 
    190190# if defined key_debug_medusa 
    191                               IF(lwp) WRITE(numout,*) ' TREND stats (min, max,sum) kt = ',kt ,' jn = ',kjn 
    192                               CALL trc_rst_dia_stat( ptrtrd(:,:,1), cltra) 
     191         IF(lwp) WRITE(numout,*) ' TREND stats (min, max,sum) kt = ',kt ,' jn = ',kjn 
     192         CALL trc_rst_dia_stat( ptrtrd(:,:,1), cltra) 
    193193# endif 
    194                               CALL iom_put( cltra,  ptrtrd(:,:,:) ) 
     194         CALL iom_put( cltra,  ptrtrd(:,:,:) ) 
    195195# if defined key_debug_medusa 
    196                            ELSE 
    197                               IF(lwp) WRITE(numout,*) & 
    198                               ' TREND -- No output asked for ',cltra,' kt = ',kt,' jn = ',kjn 
    199                               CALL trc_rst_dia_stat( ptrtrd(:,:,1), cltra) 
     196      ELSE 
     197         IF(lwp) WRITE(numout,*) & 
     198                      ' TREND -- No output asked for ',cltra,' kt = ',kt,' jn = ',kjn 
     199         CALL trc_rst_dia_stat( ptrtrd(:,:,1), cltra) 
    200200# endif 
    201                            ENDIF 
     201      ENDIF 
    202202 
    203203   END SUBROUTINE trd_trc_iomput 
Note: See TracChangeset for help on using the changeset viewer.