Missing upper-air pressure levels in output, and missing all geopotential height!

Hi all

I've got a setup of OpenIFS cy40r1v2 where I'm writing both surface, pressure levels and model levels to output at the same time. This works fine for all variables except one: geopotential height. 

My namelist has the following in it:

&NAMFPC
CFPFMT="MODEL",
NFP3DFS=4,
NFP3DFT=0,
NFP3DFV=0,
MFP3DFS(:)=131,132,130,133 !! u,v,t,q on model levels
!MFP3DFT=60,
!MFP3DFV=133,
NFP2DF=3,
MFP2DF=129,134,152, !z, sp, lnsp
NFPPHY=39,
NFP3DFP=7,
MFP3DFP=131,132,130,133,203,135,156 !! u,v,t,q,o3,w,zh on pressure levels
MFPPHY(:)=167,235,51,52,151, 134,165,166,207,168, 129,228,144,143,182, 44,180,181,146,147, 175,176,169,177,210, 211,164,212,178,179, 20\
8,209,157,135,78, 79,8,9,206
NRFP3S(:)=1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45\
,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,
RFP3P(:)=100000, 92500, 85000, 77500, 70000, 60000, 50000, 40000, 30000, 25000, 20000, 15000, 10000, 7000, 5000, 3000, 2000, 1000, 700,\
500, 300, 200, 100, 50, 20, 10, 5, 2, 1,
NFPCLI=0,
LFPQ=false,
LTRACEFP=false,
RFPCORR=60000.,
/

Does anyone know why the model does not write geopotential height on pressure levels? 

Also, I've noticed that some upper-atmosphere pressure levels are NaN in the output (see attached plot from ncview of temperature on a latitude-pressure plane). 
In particular, some pressure levels above 20 hPa are NaNs in the tropics. Does OpenIFS not interpolate/extrapolate between model levels to get data on a pressure level? Or is there a setting in the namelist I'm missing? 

Any help would be much appreciated! 

Best regards

Joakim