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 10036 for NEMO/trunk/src/OCE/TRD – NEMO

Ignore:
Timestamp:
2018-08-07T09:44:56+02:00 (6 years ago)
Author:
davestorkey
Message:

Remove erroneous duplicated code in trdtra.F90. Ticket #2125.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/TRD/trdtra.F90

    r9598 r10036  
    358358         SELECT CASE( ktrd ) 
    359359         CASE( jptra_atf  )   ;   CALL iom_put( "ttrd_atf" , ptrdx )        ! asselin time Filter 
    360             CALL iom_put( "strd_atf" , ptrdy ) 
     360                                  CALL iom_put( "strd_atf" , ptrdy ) 
    361361         END SELECT 
    362362      END IF 
    363       ! 
    364       ! These trends are done every second time step. When 1ts output is selected must go different (2ts) file from standard T-file 
    365       IF( MOD( kt, 2 ) == 0 ) THEN 
    366          SELECT CASE( ktrd ) 
    367          CASE( jptra_xad  )   ;   CALL iom_put( "ttrd_xad"  , ptrdx )         ! x- horizontal advection 
    368                                   CALL iom_put( "strd_xad"  , ptrdy ) 
    369          CASE( jptra_yad  )   ;   CALL iom_put( "ttrd_yad"  , ptrdx )         ! y- horizontal advection 
    370                                   CALL iom_put( "strd_yad"  , ptrdy ) 
    371          CASE( jptra_zad  )   ;   CALL iom_put( "ttrd_zad"  , ptrdx )         ! z- vertical   advection 
    372                                   CALL iom_put( "strd_zad"  , ptrdy ) 
    373                                   IF( ln_linssh ) THEN                   ! cst volume : adv flux through z=0 surface 
    374                                      ALLOCATE( z2dx(jpi,jpj) , z2dy(jpi,jpj) ) 
    375                                      z2dx(:,:) = wn(:,:,1) * tsn(:,:,1,jp_tem) / e3t_n(:,:,1) 
    376                                      z2dy(:,:) = wn(:,:,1) * tsn(:,:,1,jp_sal) / e3t_n(:,:,1) 
    377                                      CALL iom_put( "ttrd_sad", z2dx ) 
    378                                      CALL iom_put( "strd_sad", z2dy ) 
    379                                      DEALLOCATE( z2dx, z2dy ) 
    380                                   ENDIF 
    381          CASE( jptra_totad  ) ;   CALL iom_put( "ttrd_totad", ptrdx )         ! total   advection 
    382                                   CALL iom_put( "strd_totad", ptrdy ) 
    383          CASE( jptra_ldf  )   ;   CALL iom_put( "ttrd_ldf"  , ptrdx )         ! lateral diffusion 
    384                                   CALL iom_put( "strd_ldf"  , ptrdy ) 
    385          CASE( jptra_zdf  )   ;   CALL iom_put( "ttrd_zdf"  , ptrdx )         ! vertical diffusion (including Kz contribution) 
    386                                   CALL iom_put( "strd_zdf"  , ptrdy ) 
    387          CASE( jptra_zdfp )   ;   CALL iom_put( "ttrd_zdfp" , ptrdx )         ! PURE vertical diffusion (no isoneutral contribution) 
    388                                   CALL iom_put( "strd_zdfp" , ptrdy ) 
    389          CASE( jptra_evd )    ;   CALL iom_put( "ttrd_evd"  , ptrdx )         ! EVD trend (convection) 
    390                                   CALL iom_put( "strd_evd"  , ptrdy ) 
    391          CASE( jptra_dmp  )   ;   CALL iom_put( "ttrd_dmp"  , ptrdx )         ! internal restoring (damping) 
    392                                   CALL iom_put( "strd_dmp"  , ptrdy ) 
    393          CASE( jptra_bbl  )   ;   CALL iom_put( "ttrd_bbl"  , ptrdx )         ! bottom boundary layer 
    394                                   CALL iom_put( "strd_bbl"  , ptrdy ) 
    395          CASE( jptra_npc  )   ;   CALL iom_put( "ttrd_npc"  , ptrdx )         ! static instability mixing 
    396                                   CALL iom_put( "strd_npc"  , ptrdy ) 
    397          CASE( jptra_bbc  )   ;   CALL iom_put( "ttrd_bbc"  , ptrdx )         ! geothermal heating   (only on temperature) 
    398          CASE( jptra_nsr  )   ;   CALL iom_put( "ttrd_qns"  , ptrdx(:,:,1) )  ! surface forcing + runoff (ln_rnf=T) 
    399                                   CALL iom_put( "strd_cdt"  , ptrdy(:,:,1) )        ! output as 2D surface fields 
    400          CASE( jptra_qsr  )   ;   CALL iom_put( "ttrd_qsr"  , ptrdx )         ! penetrative solar radiat. (only on temperature) 
    401          END SELECT 
    402          ! the Asselin filter trend  is also every other time step but needs to be lagged one time step 
    403          ! Even when 1ts output is selected can go to the same (2ts) file as the trends plotted every even time step. 
    404       ELSEIF( MOD( kt, 2 ) == 1 ) THEN 
    405          SELECT CASE( ktrd ) 
    406          CASE( jptra_atf  )   ;   CALL iom_put( "ttrd_atf" , ptrdx )        ! asselin time Filter 
    407                                   CALL iom_put( "strd_atf" , ptrdy ) 
    408          END SELECT 
    409       ENDIF 
    410363      ! 
    411364   END SUBROUTINE trd_tra_iom 
Note: See TracChangeset for help on using the changeset viewer.