Ignore:
Timestamp:
2020-03-24T15:08:27+01:00 (8 months ago)
Author:
gsamson
Message:

correct mxld_abl & mxlm_abl init to avoid division by zero (ticket #2419)

Location:
NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/src/ABL
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/src/ABL/ablmod.F90

    r12588 r12592  
    860860            ! 
    861861            DO ji = 1, jpi 
    862                mxld_abl ( ji, jj,    1 ) = 0._wp 
     862               mxld_abl ( ji, jj,    1 ) = mxl_min 
    863863               mxld_abl ( ji, jj, jpka ) = mxl_min 
    864                mxlm_abl ( ji, jj,    1 ) = 0._wp 
     864               mxlm_abl ( ji, jj,    1 ) = mxl_min 
    865865               mxlm_abl ( ji, jj, jpka ) = mxl_min 
    866866               zldw     ( ji,        1 ) = zrough(ji,jj) * rn_Lsfc 
  • NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/src/ABL/sbcabl.F90

    r12588 r12592  
    289289         tq_abl (:,:,:,nt_a,:   ) = 0._wp 
    290290         tke_abl(:,:,:,nt_a     ) = 0._wp 
     291 
     292         mxlm_abl(:,:,:         ) = mxl_min 
     293         mxld_abl(:,:,:         ) = mxl_min 
    291294      ENDIF 
    292295 
Note: See TracChangeset for help on using the changeset viewer.