<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>
</head>
<body ocsi="0" fpstyle="1">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Hi Users<br>
<br>
I would like to apply a multi-model ensemble drying mask to various CMIP5 model fields.<br>
This is how I am currently going about it in the code below. I work out the MME and then mask out wetting to only have regions of drying and then I multiply the mask to the original CMIP5 model files I would like to apply it to.<br>
Just would like to check if this is how to go about it or another better way. I feel my output is not 100% correct.<br>
<br>
Many thanks for your help in advance.<br>
<br>
Kindest Regards<br>
Melissa<br>
<br>
if (gg.eq.0) then<br>
<br>
dP_ENS = new ((/dimsizes(model),dimsizes(lat),dimsizes(lon)/), typeof(dP))<br>
dP_ENS!0="model"<br>
dP_ENS&model=model <br>
dP_ENS!1="lat"<br>
dP_ENS&lat=lat<br>
dP_ENS!2="lon"<br>
dP_ENS&lon=lon<br>
end if<br>
<br>
dP_ENS(gg,:,:) = dP(0,:,:)<br>
<br>
dP_ENSMEAN = dim_avg_n_Wrap(dP_ENS,0) ; dP_ENS(lat, lon)<br>
<br>
printVarSummary(dP_ENSMEAN)<br>
<br>
dP_ENS_dry = mask(dP_ENSMEAN,dP_ENSMEAN.gt.0, False) ; mask out all the wetting regions i.e. any values above 0 for dP<br>
<br>
dP_dry = new ((/dimsizes(lat),dimsizes(lon)/), typeof(dP_ENSMEAN))<br>
dP_dry!0="lat"<br>
dP_dry&lat=lat<br>
dP_dry!1="lon"<br>
dP_dry&lon=lon<br>
<br>
dP_dry = dP_ENS_dry(:,:) * dP(0,:,:) ; apply mask to original files dP<br>
<br>
</div>
</body>
</html>