Based on cy38t1_bf03 Additional fullpos parameters for convection forecasting - MLCAPE -> another option in computation of CAPE (NFPCAPE=6): mix the lowest GCAPEPSD Pa of atmosphere before computing CAPE - storm motion diagnostics -> average wind between 0 and 6km - vertical wind shear diagonostic -> simple vector difference for wind at 6km - 2m (Deep Layer Shear) and at 1km -2m (Low Level Shear) - added for convenience reasons - storm relative helicity -> vertical integral of dot product between storm relative wind and vorticity vector - lightning diagnostics -> acording to Price and Rind (1992) and McCaul et al (2009) - temperature lapse rate (vertical temperature gradient) Except for the last one, all fields are 2D. -- authors: Jure Cedilnik (ARSO) and Christoph Wittmann (ZAMG) -- modified routines: arp/fullpos/fpcica.F90 arp/fullpos/sufpc.F90 arp/fullpos/vpos.F90 arp/fullpos/phymfpos.F90 arp/fullpos/endpos.F90 arp/module/yomcape.F90 arp/module/yomfpc.F90 arp/module/yomafn.F90 arp/module/parfpos.F90 arp/namelist/namafn.h arp/namelist/namfpc.h arp/pp_obs/pos.F90 arp/setup/sucape.F90 arp/setup/suafn3.F90 arp/setup/suafn1.F90 arp/setup/suafn2.F90 -- new routines arp/fullpos/fpstrmm.F90 <- computation of storm motion arp/fullpos/fpdiagflash.F90 <- computation of diagnostics of flashes arp/fullpos/fpsrh.F90 <- computation of storm relative helicity arp/fullpos/fpshear.F90 <- computation of vertical wind shear -- doc: SUFPC.F90 <- - Additional option for computation of CAPE (Mixed Layer CAPE) SUFPC <- allow for NFPCAPE=6 - defaults for some parameters related to new variables SUCAPE.F90 < - check value of GCAPEPSD in case of NCAPE=6 SUAFN1.F90, SUAFN2.F90, SUAFN3.F90 < - addition of new fullpos_type variables ENDPOS.F90 < - changed the order of storing of data in buffer for variables LCL,FLC and EL (zoomed fullpos didn't work), probably because all 1D should be before 2D - additional fullpos_type variables - from YOMAFN + their corresoponding local arrays + local logical switches - modified call to FPCICA (to include KAPPA) PHYMFPOS.F90 < - additional fullpos_type variables - from YOMAFN + their corresoponding local arrays + local logical switches - modified call to FPCICA (to include KAPPA) - added global array for geometric vertical velocity - coming as input - to pass to ligthning diagnostics routine VPOS.F90 < - addition of ZWW (vertical velocity) - in the calls to POS and PHYMFPOS - to make it available FPDIAGFLASH POS.F90 < - added vertical velocity in output (PWW) - added VTG (vertical temperature gradient) and FLASH variables (from module yomafn and their Z...,ZPP.... and L.. counterparts)A - compute VTG (3D field) FPCICA.F90 < - added another call to FPCINCAPE, after mixing the lowest GCAPEPSD Pa thick layer YOMFPC.F90 <- added variables - NMTFLASH - choice of lighntning diagnostics (0,1,2,3) but only 0 and 2 are fully coded - RFLASHCAL - calibration factor for diagnostics of flashes - RFLASHQIL - some parameters for lightning diagnostics - RFLASHTMXPL - -ii- - RFLASHTCRIT - -ii- - RSTRMMH - height of the top for computation of storm motion vector (default 6 km) - RSRHH - height of the top for computation of storm relative helicity (default 3 km) - RDLSH - height (above orography) for computation of deep layer shear (default 6 km) - RLLSH - height (above orography) for computation of low layer shear (default 1 km) YOMAFN.F90 <- added fullpos_type variables YOMCAPE.F90 <- only additional comments PARFPOS.F90 <- increased JPOS3DF and JPOS2DF, according to additonal number of variables -- namelist example &NAMFPC CFPDOM(1)='LALON', CFPFMT='LALON', CFP3DF(1)='GEOPOTENTIEL', CFP3DF(2)='RAIN', CFP3DF(3)='VERTTEMPGRAD', CFPPHY(1)='SURFGEOPOTENTIEL', CFPPHY(2)='SURFIND.TERREMER' CFPXFU(1)='CLSVENT.ZONAL' CFPXFU(2)='CLSVENT.MERIDIEN' CFP2DF(1)='MSLPRESSURE', CFP2DF(2)='SURFCAPE.POS.F00', CFP2DF(3)='SURFFLASHDIAG' CFP2DF(4)='ICAOTROP.PRESSUR', CFP2DF(5)='ICAOTROP.TEMPERA', CFP2DF(6)='SURFTOT.WAT.VAPO', CFP2DF(7)='SURFPRESSION', CFP2DF(8)='SURFLIFTCONDLEV ' CFP2DF(9)='SURFUSTORMMOTION' CFP2DF(10)='SURFVSTORMMOTION' CFP2DF(11)='SURFDEEPLAYERSHR', CFP2DF(12)='SURFLOWLEVELSHR ', CFP2DF(13)='SURFHELICITY ', CFP2DF(14)='SURFLIFTCONDLEV', CFP2DF(15)='SURFFREECONVLEV', CFP2DF(16)='SURFEQUILIBRLEV', NRFP3S=1,10,20,30,40,50,60,70,80,81,82,83,84,85,86,87 RFP3P=30000.,50000.,55000.,60000.,65000.,70000.,75000.,80000.,85000.,90000.,92500.,100000., NFPCAPE=6 NMTFLASH=2 LFPCAPEX =.TRUE. LFPMOIS=.TRUE., NFPCLI=0, / &NAMFPD RLONC(1)=13.2643, RLATC(1)=46.3881, RDELX(1)=0.028, RDELY(1)=0.018, NLAT=151, NLON=171, NFPLUX=171, NFPGUX=151, /