Opened 6 years 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.
Note: See
TracTickets for help on using
tickets.