<div dir="ltr">Dear Sir Dennis,<div><br></div><div>Thank you for your help.</div><div><br></div><div>I got this working already but I am still not sure whether the output is correct or not.</div><div>Some areas have no CAPE values. Maybe I am just forgetting something basic here.</div><div><br></div><div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Attached files are my script and output file.</div><div><br></div><div>I'll appreciate any comments about this.</div></div></div></div></div></div></div></div></div></div></div></div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jul 26, 2019 at 12:34 AM Dennis Shea <<a href="mailto:shea@ucar.edu">shea@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><b><font color="#ff0000">printVarSummary(qf) ; [time | 365] x [level | 8] x [lat | 73] x [lon | 144]</font></b><br><br> p = (/100000,92500,85000,70000,60000,50000,40000,30000/)<br> z = (/111,762,1457,3011,4205,5572,7182,9160/)<br><br><div> zconform = <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/conform.shtml" target="_blank"><b>conform</b></a>(qf,z,1) ; <======<br></div><div> pconform = <b>conform</b>(qf,p,1) ; <======<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Jul 23, 2019 at 11:15 PM Lyndz via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div dir="ltr" class="gmail-m_2798422837249268762gmail-m_-1012073796851807859m_2922388677776988710gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Dear NCL-experts,</div><div><br></div><div>I am calculating CAPE based on NCEP reanalysis.</div><div>I encountered the following errors when conforming the data.</div><div><br></div><div>fatal:conform: The array to be conformed must have the same number of dimensions as indicated by the length of the last argument<br></div><div><br></div><div>Here's the part of the script:</div><div><br></div><div> begin<br> year = "1979"<br> rh_f = addfile("../rh/rhum."+year+".nc", "r")<br> temp_f = addfile("../temp/air."+year+".nc", "r")<br> q_f = addfile("../shum/shum_"+year+".nc", "r")<br> psfc_f= addfile("../pres/pres.sfc.gauss."+year+".nc","r")<br> zsfc_f= addfile("<a href="http://hgt.sfc.nc" target="_blank">hgt.sfc.nc</a>","r")<br><br> rhf = short2flt(rh_f->rhum(:,{1000:300},:,:))<br> tempf = short2flt(temp_f->air(:,{1000:300},:,:))<br> presf = short2flt(psfc_f->pres(:,:,:))<br> qf = q_f->shum(:,{1000:300},:,:)<br> zsfcf = short2flt(zsfc_f->hgt(:,:,:))<br><br><b><font color="#ff0000"> printVarSummary(qf) ; [time | 365] x [level | 8] x [lat | 73] x [lon | 144]</font></b><br><br> p = (/100000,92500,85000,70000,60000,50000,40000,30000/)<br> z = (/111,762,1457,3011,4205,5572,7182,9160/)<br><br> zconform = conform(qf,z,(/0,1,2,3/))<br> pconform = conform(qf,p,(/0,1,2,3/))<br></div><div><br></div><div><br></div><div>;;**************************************************<br>;;Calculates CAPE(J/kg)<br>;;**************************************************<br> cinfo = wrf_cape_3d(pconform,tempf,qf,zconform,zsfcf,presf,False)<br></div><div><br></div><div><br></div><div>Any ideas on how I should do this correctly in NCL?</div><div><br></div><div>I'll appreciate any help.</div><div><br></div><div>Sincerely,</div><div>Lyndz</div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>
</blockquote></div>