Changeset 11424


Ignore:
Timestamp:
2019-08-08T18:08:53+02:00 (15 months ago)
Author:
mattmartin
Message:

Remove the MAX statementsin zdftke which were added as part of the STOPACK changes in order to reproduce the results from the standard GO6 branch.

Location:
branches/UKMO/dev_r5518_GO6_package_FOAMv14_STOPACK/NEMOGCM/NEMO/OPA_SRC
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package_FOAMv14_STOPACK/NEMOGCM/NEMO/OPA_SRC/TRA/trabbl.F90

    r11419 r11424  
    566566         DO ji = 1, jpim1 
    567567            mgrhu(ji,jj) = INT(  SIGN( 1.e0, gdept_0(ji+1,jj,mbkt(ji+1,jj)) - gdept_0(ji,jj,mbkt(ji,jj)) )  ) 
    568             ! 
    569568            mgrhv(ji,jj) = INT(  SIGN( 1.e0, gdept_0(ji,jj+1,mbkt(ji,jj+1)) - gdept_0(ji,jj,mbkt(ji,jj)) )  ) 
    570569         END DO 
  • branches/UKMO/dev_r5518_GO6_package_FOAMv14_STOPACK/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90

    r11417 r11424  
    358358                  zwlc = zind * rn_lc0(ji,jj) * zus * SIN( rpi * fsdepw(ji,jj,jk) / zhlc(ji,jj) ) 
    359359                  !                                           ! TKE Langmuir circulation source term 
    360                   en(ji,jj,jk) = en(ji,jj,jk) + rdt * MAX(0.,1._wp - fr_i(ji,jj) ) * ( zwlc * zwlc * zwlc ) /   & 
     360                  en(ji,jj,jk) = en(ji,jj,jk) + rdt * ( 1._wp - fr_i(ji,jj) )* ( zwlc * zwlc * zwlc ) /   & 
    361361                     &   zhlc(ji,jj) * wmask(ji,jj,jk) * tmask(ji,jj,1) 
    362362 
     
    487487               DO ji = fs_2, fs_jpim1   ! vector opt. 
    488488                  en(ji,jj,jk) = en(ji,jj,jk) + rn_efr0(ji,jj) * en(ji,jj,1) * EXP( -fsdepw(ji,jj,jk) / htau(ji,jj) )   & 
    489                      &                                 * MAX(0.,1._wp - fr_i(ji,jj) )  * wmask(ji,jj,jk) * tmask(ji,jj,1) 
     489                     &                                 * ( 1._wp - fr_i(ji,jj) )  * wmask(ji,jj,jk) * tmask(ji,jj,1) 
    490490               END DO 
    491491            END DO 
     
    496496               jk = nmln(ji,jj) 
    497497               en(ji,jj,jk) = en(ji,jj,jk) + rn_efr0(ji,jj) * en(ji,jj,1) * EXP( -fsdepw(ji,jj,jk) / htau(ji,jj) )   & 
    498                   &                                 * MAX(0.,1._wp - fr_i(ji,jj) )  * wmask(ji,jj,jk) * tmask(ji,jj,1) 
     498                  &                                 * ( 1._wp - fr_i(ji,jj) )  * wmask(ji,jj,jk) * tmask(ji,jj,1) 
    499499            END DO 
    500500         END DO 
     
    512512                  zdif = rhftau_scl * MAX( 0._wp, zdif + rhftau_add )  ! apply some modifications... 
    513513                  en(ji,jj,jk) = en(ji,jj,jk) + zbbrau(ji,jj) * zdif * EXP( -fsdepw(ji,jj,jk) / htau(ji,jj) )   & 
    514                      &                        * MAX(0.,1._wp - fr_i(ji,jj) ) * wmask(ji,jj,jk) * tmask(ji,jj,1) 
     514                     &                        * ( 1._wp - fr_i(ji,jj) ) * wmask(ji,jj,jk) * tmask(ji,jj,1) 
    515515               END DO 
    516516            END DO 
     
    885885      ALLOCATE( rn_efr0  (jpi,jpj) ) ; rn_efr0   = rn_efr 
    886886       
     887      IF( nn_etau == 2  ) THEN 
     888          ierr = zdf_mxl_alloc() 
     889          nmln(:,:) = nlb10           ! Initialization of nmln 
     890      ENDIF 
    887891 
    888892      IF( nn_etau /= 0 .and. nn_htau == 2 ) THEN 
Note: See TracChangeset for help on using the changeset viewer.