Changeset 11417


Ignore:
Timestamp:
2019-08-07T13:04:25+02:00 (15 months ago)
Author:
mattmartin
Message:

Fixed bug in zdftke code for stopack.

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/ZDF/zdftke.F90

    r11404 r11417  
    222222      IF( .NOT.Agrif_Root() )   CALL Agrif_Update_Tke( kt )      ! children only 
    223223#endif       
    224       IF (  kt .eq. nitend ) THEN 
     224      IF (  kt == nitend ) THEN 
    225225        DEALLOCATE ( rn_lc0, rn_ediff0, rn_ediss0, rn_ebb0, rn_efr0 ) 
    226226      ENDIF 
     
    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 - 2.*fr_i(ji,jj) ) * ( zwlc * zwlc * zwlc ) /   & 
     360                  en(ji,jj,jk) = en(ji,jj,jk) + rdt * MAX(0.,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 - 2.*fr_i(ji,jj) )  * wmask(ji,jj,jk) * tmask(ji,jj,1) 
     489                     &                                 * MAX(0.,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 - 2.*fr_i(ji,jj) )  * wmask(ji,jj,jk) * tmask(ji,jj,1) 
     498                  &                                 * MAX(0.,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 - 2.*fr_i(ji,jj) ) * wmask(ji,jj,jk) * tmask(ji,jj,1) 
     514                     &                        * MAX(0.,1._wp - fr_i(ji,jj) ) * wmask(ji,jj,jk) * tmask(ji,jj,1) 
    515515               END DO 
    516516            END DO 
  • branches/UKMO/dev_r5518_GO6_package_FOAMv14_STOPACK/NEMOGCM/NEMO/OPA_SRC/step.F90

    r11404 r11417  
    176176      ! 
    177177      IF( lk_ldfslp ) THEN                            ! slope of lateral mixing 
    178          IF(ln_sto_eos ) CALL sto_pts( tsn )          ! Random T/S fluctuations 
    179178                         CALL eos( tsb, rhd, gdept_0(:,:,:) )               ! before in situ density 
    180179         IF( ln_zps .AND. .NOT. ln_isfcav)                               & 
     
    421420         IF( lwm.AND.numoni /= -1 ) CALL FLUSH    ( numoni )     ! flush output namelist ice 
    422421      ENDIF 
    423       IF( ln_sto_eos       )   CALL sto_rst_write( kstp )   ! write restart file for stochastic parameters 
    424  
    425422 
    426423      IF( lrst_bias )          CALL bias_wrt     ( kstp ) 
Note: See TracChangeset for help on using the changeset viewer.