Opened 4 weeks ago

#213 new Defect

wrong unit for max convection level in gridcheck

Reported by: pesei Owned by:
Priority: major Milestone: FLEXPART 10
Component: FP coding/compilation Version: FLEXPART 9.0.3
Keywords: Cc:

Description

gridcheck.f90 (in Fp10, gridcheck_ecmwf.f90 and gridchec_gfs.f90 has the following code

  ! Determine the uppermost level for which the convection scheme shall be applied
  ! by assuming that there is no convection above 50 hPa (for standard SLP)
  !*****************************************************************************

  do i=1,nuvz-2
    pint=akz(i)+bkz(i)*101325.
    if (pint.lt.5000.) goto 96
  end do
96   nconvlev=i
  if (nconvlev.gt.nconvlevmax-1) then
    nconvlev=nconvlevmax-1
    write(*,*) 'Attention, convection only calculated up to ', &
         akz(nconvlev)+bkz(nconvlev)*1013.25,' hPa'
  endif

The output message will not give the pressure at nconvlev in hPa, because akz is in Pa.

furthermore, the purpose of all this is not clear as the convection subroutines only use nconvlevmax.

Note that these lines were added going from v4.4 to v5.1 and are still present in the latest v10.3 draft.

Change History (0)

Note: See TracTickets for help on using tickets.
hosted by ZAMG