[ncl-talk] Error with conform

Lyndz olagueralyndonmark429 at gmail.com
Tue Jul 23 23:14:27 MDT 2019


Dear NCL-experts,

I am calculating CAPE based on NCEP reanalysis.
I encountered the following errors when conforming the data.

fatal:conform: The array to be conformed must have the same number of
dimensions as indicated by the length of the last argument

Here's the part of the script:

 begin
   year = "1979"
   rh_f = addfile("../rh/rhum."+year+".nc", "r")
   temp_f = addfile("../temp/air."+year+".nc", "r")
   q_f = addfile("../shum/shum_"+year+".nc", "r")
   psfc_f= addfile("../pres/pres.sfc.gauss."+year+".nc","r")
   zsfc_f= addfile("hgt.sfc.nc","r")

   rhf = short2flt(rh_f->rhum(:,{1000:300},:,:))
   tempf = short2flt(temp_f->air(:,{1000:300},:,:))
   presf = short2flt(psfc_f->pres(:,:,:))
   qf = q_f->shum(:,{1000:300},:,:)
   zsfcf = short2flt(zsfc_f->hgt(:,:,:))

*   printVarSummary(qf) ; [time | 365] x [level | 8] x [lat | 73] x [lon |
144]*

   p  = (/100000,92500,85000,70000,60000,50000,40000,30000/)
   z = (/111,762,1457,3011,4205,5572,7182,9160/)

   zconform = conform(qf,z,(/0,1,2,3/))
   pconform = conform(qf,p,(/0,1,2,3/))


;;**************************************************
;;Calculates CAPE(J/kg)
;;**************************************************
 cinfo = wrf_cape_3d(pconform,tempf,qf,zconform,zsfcf,presf,False)


Any ideas on how I should do this correctly in NCL?

I'll appreciate any help.

Sincerely,
Lyndz
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190724/4f973029/attachment.html>


More information about the ncl-talk mailing list