vertical interpolation

Dear Sir/Madam:   I have some model-level output data from a terrain-following model. I want to interpolate the model level data to a level surface, say: 5000 m. I did as:  do i = 0, 1350-1  do j = 0, 1200-1  w(i,j) = linint1(sht(:,i,j),ww(:,i,j),False,5000.,0)  p(i,j) = linint1(sht(:,i,j),th(:,i,j),False,5000.,0)  end do  end do Where the "sht(:,I,j)" is the height of model levels at a grid point.   This way takes a lot of CPU time. I am wonder if there is a better way to do this. Thanks very much!   Ed

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Apr 22 16:25:26 2014