Vinth2p_ecmwf errors

  
Hello,
 
I am trying to convert CCSM 3.0 output variables to different pressure levels using the vinth2p_ecmwf command. I have attached the program (CCSM_2_MM5_levs.ncl) that I am using to carry out the conversion. I get the following errors while trying to run the program.
 
When I set tbot = t(:,nlev-1,:,:), It gives me the following error.

 "fatal:Number of subscripts on rhs do not match number of dimensions of
 variable,(4) Subscripts used, (3) Subscripts expected
 fatal:Execute: Error occurred at or near line 59 in file
 CCSM_2_MM5_levs.ncl"

 If I remove the last semi colon ( I know I should not but I just tried)
 and use tbot = t(:,nlev-1,:) it does not give me the previous error but instead gives the following error.

 "fatal:vinth2p_ecmwf: Surface pressure must have same number of time, lat
 and lon elements as input, number of dimensions does not match.
 fatal:Execute: Error occurred at or near line 75 in file
 CCSM_2_MM5_levs.ncl"

 The line where the error is occurring is this

 lev_array(i,:,:,:,:) = vinth2p_ecmwf(var,hyam,hybm,pnew,PS,interp(i),P0,1,True,varflg(i),tbot,PHIS)
 
I am using a program that someone else has written and when I checked with them, they informed me that they did not get this error.
 
I printed the variable summary using the printVarSummary command and I obtained the following output:
 
Variable: t
Type: float
Total Size: 15859712 bytes
            3964928 values
Number of Dimensions: 3
Dimensions and sizes: [time | 121] x [lat | 128] x [lon | 256]
Coordinates:
            time: [54750..54780]
            lat: [-88.9277353522959..88.9277353522959]
            lon: [ 0..358.59375]
Number Of Attributes: 2
  units : K
  long_name : Surface temperature (radiative)

Variable: tbot
Type: float
Total Size: 15859712 bytes
            3964928 values
Number of Dimensions: 3
Dimensions and sizes: [time | 121] x [lat | 128] x [lon | 256]
Coordinates:
            time: [54750..54780]
            lat: [-88.9277353522959..88.9277353522959]
            lon: [ 0..358.59375]
Number Of Attributes: 2
  long_name : Surface temperature (radiative)
  units : K

Variable: PS
Type: float
Total Size: 15859712 bytes
            3964928 values
Number of Dimensions: 3
Dimensions and sizes: [time | 121] x [lat | 128] x [lon | 256]
Coordinates:
            time: [54750..54780]
            lat: [-88.9277353522959..88.9277353522959]
            lon: [ 0..358.59375]
Number Of Attributes: 2
  units : Pa
  long_name : Surface pressure

Variable: PHIS
Type: float
Total Size: 15859712 bytes
            3964928 values
Number of Dimensions: 3
Dimensions and sizes: [time | 121] x [lat | 128] x [lon | 256]
Coordinates:
            time: [54750..54780]
            lat: [-88.9277353522959..88.9277353522959]
            lon: [ 0..358.59375]
Number Of Attributes: 2
  units : m2/s2
  long_name : Surface geopotential

 
I am not very familiar with NCL programming language and hence am not able to figure out why I am getting the above mentioned errors.
 
Can anyone please help me figure out the solution?
 
Thanks in advance.
 
Spandana Tummuri

 
 
Research Assistant,
Texas Tech University,
2717, 3rd St. # 315
Lubbock, Texas-79415
806-741-0901 (Home)
806-470-7732 (Cell).

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk