Changeset 1193


Ignore:
Timestamp:
2008-09-18T14:45:43+02:00 (12 years ago)
Author:
cetlod
Message:

Correction of transport module to ensure reproductibility for TOP configurations, see ticket:253

Location:
trunk/NEMO/TOP_SRC/TRP
Files:
11 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/TOP_SRC/TRP/trcadv_cen2.F90

    r1175 r1193  
    284284         ! 3. Save the horizontal advective trends for diagnostics 
    285285         ! ------------------------------------------------------- 
    286 !CDIR BEGIN COLLAPSE 
    287286         TRDTRC_XY : IF( l_trdtrc )THEN 
    288287 
     
    342341 
    343342         ENDIF TRDTRC_XY 
    344 !CDIR END 
    345343         !                                                    ! =========== 
    346344      END DO                                                  ! tracer loop 
     
    415413         ! ---------------------------------------------------- 
    416414 
    417 !CDIR BEGIN COLLAPSE 
    418415         TRDTRC_Z : IF( l_trdtrc )THEN 
    419416            ztrtrd(:,:,:) = 0.e0 
     
    435432 
    436433         ENDIF TRDTRC_Z 
    437 !CDIR END 
    438434         !                                                    ! =========== 
    439435      END DO                                                  ! tracer loop 
  • trunk/NEMO/TOP_SRC/TRP/trcadv_muscl.F90

    r1175 r1193  
    238238         ! 3. Save the horizontal advective trends for diagnostics 
    239239         ! ------------------------------------------------------- 
    240 !CDIR BEGIN COLLAPSE 
    241240         TRDTRC_XY : IF( l_trdtrc ) THEN 
    242241 
     
    282281 
    283282         ENDIF TRDTRC_XY 
    284 !CDIR END 
    285283 
    286284      ENDDO 
     
    375373         ! 3. Save the vertical advective trends for diagnostic 
    376374         ! ---------------------------------------------------- 
    377 !CDIR BEGIN COLLAPSE 
    378375         TRDTRC_Z : IF( l_trdtrc )THEN 
    379376 
     
    394391 
    395392         END IF TRDTRC_Z 
    396 !CDIR END 
    397393      END DO 
    398394 
  • trunk/NEMO/TOP_SRC/TRP/trcadv_muscl2.F90

    r1175 r1193  
    293293         TRDTRC_XY : IF( l_trdtrc ) THEN 
    294294 
    295 !CDIRR COLLAPSE 
    296295            ! 3.1) Passive tracer ZONAL advection trends 
    297296            DO jk = 1, jpkm1 
  • trunk/NEMO/TOP_SRC/TRP/trcadv_smolar.F90

    r1189 r1193  
    117117          
    118118         !--1.1 Horizontal advection 
    119 !CDIRR COLLAPSE 
    120119         IF( l_trdtrc ) ztrtrd(:,:,:,:) = 0.e0             ! trends 
    121120 
     
    127126            zz  (:,:,jk) = 0.e0 
    128127             
    129 !CDIRR COLLAPSE 
    130128            IF( l_trdtrc ) ztrtrd(:,:,:,:) = 0.e0             ! trends 
    131129 
  • trunk/NEMO/TOP_SRC/TRP/trcadv_tvd.F90

    r1189 r1193  
    190190         ! ------------------------------------------------------------------- 
    191191 
    192 !CDIR BEGIN COLLAPSE 
    193192         IF( l_trdtrc ) THEN 
    194193 
     
    252251 
    253252         ENDIF 
    254 !CDIR END 
    255253 
    256254         ! 4. Update and guess with monotonic sheme 
     
    369367         ! terms in the computations below ??? 
    370368 
    371 !CDIR BEGIN COLLAPSE 
    372369         IF( l_trdtrc ) THEN 
    373370             
     
    412409             
    413410         ENDIF 
    414 !CDIR END 
    415411 
    416412 
  • trunk/NEMO/TOP_SRC/TRP/trcbbc.F90

    r1175 r1193  
    103103         DO jn = 1, jptra                                        ! tracer loop 
    104104            !                                                    ! =========== 
    105 !CDIR COLLAPSE 
    106105            IF( l_trdtrc )  ztrtrd(:,:,:) = tra(:,:,:,jn)  ! save trends 
    107106 
     
    120119             
    121120            IF( l_trdtrc ) THEN 
    122 !CDIR COLLAPSE 
    123121               ztrtrd(:,:,:) = tra(:,:,:,jn) - ztrtrd(:,:,:) 
    124122               IF (luttrd(jn)) CALL trd_mod_trc(ztrtrd, jn, jptrc_trd_bbc, kt) 
  • trunk/NEMO/TOP_SRC/TRP/trcldf_iso.F90

    r1175 r1193  
    119119         !                                                       ! =========== 
    120120 
    121 !CDIR COLLAPSE 
    122121         IF( l_trdtrc ) ztrtrd(:,:,:) = tra(:,:,:,jn) 
    123122 
     
    239238               END DO 
    240239                
    241 !CDIRR COLLAPSE 
    242240               DO jj = 2 , jpjm1 
    243241                  DO ji = fs_2, fs_jpim1   ! vector opt. 
  • trunk/NEMO/TOP_SRC/TRP/trcldf_iso_zps.F90

    r1175 r1193  
    126126      DO jn = 1, jptra                                           ! tracer loop 
    127127         !                                                       ! =========== 
    128 !CDIR COLLAPSE 
    129128         IF( l_trdtrc ) ztrtrd(:,:,:) = tra(:,:,:,jn)   ! save trends 
    130129 
     
    253252         ! 4. Save the horizontal diffusive and advective (eiv) trends for diagnostics 
    254253         ! --------------------------------------------------------------------------- 
    255 !CDIR BEGIN COLLAPSE 
    256254         IF( l_trdtrc ) THEN 
    257255 
  • trunk/NEMO/TOP_SRC/TRP/trcrad.F90

    r1175 r1193  
    125125         DO jn = jp_sms0, jp_sms1 
    126126         !                                                        ! =========== 
    127          ztrcorb = 0.e0   ;   ztrmasb = 0.e0 
    128          ztrcorn = 0.e0   ;   ztrmasn = 0.e0 
    129  
    130 !CDIR COLLAPSE 
    131          IF( l_trdtrc ) THEN 
    132             ztrtrdb(:,:,:) = ptrb(:,:,:,jn)                        ! save input trb for trend computation 
    133             ztrtrdn(:,:,:) = ptrn(:,:,:,jn)                        ! save input trn for trend computation 
    134          ENDIF 
     127            ztrcorb = 0.e0   ;   ztrmasb = 0.e0 
     128            ztrcorn = 0.e0   ;   ztrmasn = 0.e0 
     129 
     130           IF( l_trdtrc ) THEN 
     131              ztrtrdb(:,:,:) = ptrb(:,:,:,jn)                        ! save input trb for trend computation 
     132              ztrtrdn(:,:,:) = ptrn(:,:,:,jn)                        ! save input trn for trend computation 
     133           ENDIF 
    135134 
    136135 
  • trunk/NEMO/TOP_SRC/TRP/trcsbc.F90

    r1175 r1193  
    8484      DO jn = 1, jptra 
    8585         ! 1. Concentration dillution effect on tra 
    86 !CDIR COLLAPSE 
    8786         IF( l_trdtrc ) ztrtrd(:,:,:) = tra(:,:,:,jn)  ! save trends 
    8887 
     
    9998          
    10099         IF( l_trdtrc ) THEN 
    101 !CDIR COLLAPSE 
    102100            ztrtrd(:,:,:) = tra(:,:,:,jn) - ztrtrd(:,:,:) 
    103101            IF (luttrd(jn)) CALL trd_mod_trc(ztrtrd, jn, jptrc_trd_sbc, kt) 
  • trunk/NEMO/TOP_SRC/TRP/trczdf_iso_vopt.F90

    r1175 r1193  
    7979      IF( l_trdtrc ) THEN 
    8080         ALLOCATE( ztrcavg(jpi,jpj,jpk,jptra) ) 
    81 !CDIR COLLAPSE 
    8281         ztrcavg(:,:,:,:) = 0.e0          ! initialisation step 
    8382      ENDIF 
     
    187186         !                                                       ! =========== 
    188187          
    189 !CDIR COLLAPSE 
    190188         IF( l_trdtrc ) ztrtrd(:,:,:) = tra(:,:,:,jn)          ! save trends 
    191189          
     
    345343            !   of the rotated tensor (from trc_zdf_iso). 
    346344            IF( ln_trcldf_iso ) THEN 
    347 !CDIR COLLAPSE 
    348345               DO jk = 1, jpkm1 
    349346                  ztrtrd(:,:,jk) = ( (tra(:,:,jk,jn) - trb(:,:,jk,jn))/rdttrc(jk) ) - ztrtrd(:,:,jk)  & 
     
    351348               END DO 
    352349            ELSE 
    353 !CDIR COLLAPSE 
    354350               DO jk = 1, jpkm1 
    355351                  ztrtrd(:,:,jk) = ( (tra(:,:,jk,jn) - trb(:,:,jk,jn))/rdttrc(jk) ) - ztrtrd(:,:,jk) 
     
    449445         ztfw(1,:,:) = 0.e0    ;     ztfw(jpi,:,:) = 0.e0 
    450446 
    451 !CDIRR COLLAPSE 
    452447         IF( l_trdtrc ) ztrtrd(:,:,:) = tra(:,:,:,jn)          ! save trends 
    453448 
     
    589584 
    590585            ! II.1) Compute the eiv VERTICAL trend 
    591 !CDIRR COLLAPSE 
    592586            DO jk = 1, jpkm1 
    593587               DO jj = 2, jpjm1 
Note: See TracChangeset for help on using the changeset viewer.