Changeset 12917


Ignore:
Timestamp:
2020-05-13T09:56:41+02:00 (7 months ago)
Author:
davestorkey
Message:

UKMO/NEMO_4.0.1_biharmonic_GM : remove redundant code (which breaks reproducibility).

Location:
NEMO/branches/UKMO/NEMO_4.0.1_biharmonic_GM/src/OCE
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0.1_biharmonic_GM/src/OCE/DYN/dynldf_lap_blp.F90

    r12788 r12917  
    174174 
    175175!CW: set boundary conditions: d/dz(del^2 u) = 0 at top and bottom, so that eddy-induced velocity, w*=0 
     176!DS: note that w*=0 next to topography is already set because of the use of wumask and wvmask above. 
    176177!     Surface 
    177178      zulapdz(:,:,1)   = 0._wp   ;   zvlapdz(:,:,1)   = 0._wp  
    178179!     Flat bottom case 
    179180      zulapdz(:,:,jpk) = 0._wp   ;   zvlapdz(:,:,jpk) = 0._wp 
    180 !     Variable bathymetry case, including z-partial steps, as in dynhpg.F90, subroutine hpg_zps 
    181         DO jj = 2, jpjm1 
    182             DO ji = 2, jpim1 
    183                 iku = mbku(ji,jj)+1 
    184                 ikv = mbkv(ji,jj)+1 
    185                 zulapdz(:,:,iku) = 0._wp 
    186                 zvlapdz(:,:,ikv) = 0._wp 
    187             ENDDO 
    188         ENDDO 
    189181 
    190182!!    calculate d/dz(-bhm * d/dz(del^2 u)) 
     
    199191         ENDDO 
    200192      ENDDO 
    201  
    202       ! ----- 
    203  
    204193    
    205194   END SUBROUTINE dyn_ldf_bgm 
  • NEMO/branches/UKMO/NEMO_4.0.1_biharmonic_GM/src/OCE/LDF/ldfdyn.F90

    r12788 r12917  
    380380            ENDDO 
    381381 
    382             !     Surface BC : set to zero 
    383             bhm(:,:,1)   = 0._wp 
    384             !     Flat bottom BC : set to zero 
    385             bhm(:,:,jpk) = 0._wp  
    386             !     Variable bathymetry case: diffusive fluxes masked in dyn_ldf_bgm 
    387  
    388382            !CW Test criterion for two grid-length mode for BGM scheme with leapfrog timestepping 
    389383            ! \frac{{\kappa \Delta t}{(\Delta z)^2 (\Delta x)^2}< 1/32  
     
    460454            ENDDO 
    461455 
    462             !     Surface BC : set to zero 
    463             bhm(:,:,1)   = 0._wp  
    464             !     Flat bottom BC : set to zero 
    465             bhm(:,:,jpk) = 0._wp  
    466             !     Variable bathymetry case: diffusive fluxes masked in dyn_ldf_bgm 
    467  
    468456            !-------------------------------- 
    469457         CASE(  13  )      !==  fixed profile: steady profile of the form bhm0 X dx^2 X dz^2 in interior, zero at bottom and top  ==! 
     
    488476               ENDDO 
    489477            ENDDO 
    490  
    491             !     Surface BC : set to zero 
    492             bhm(:,:,1)   = 0._wp  
    493             !     Flat bottom BC : set to zero 
    494             bhm(:,:,jpk) = 0._wp  
    495             !     Variable bathymetry case: diffusive fluxes masked in dyn_ldf_bgm 
    496478 
    497479            !-------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.