<div dir="ltr">Hi all,<div><br></div><div>I just wanted to make sure that I used this function correctly</div><div><br></div><div><br></div><div>here is my script for CAM models:</div><div><br></div><div><div><br></div><div>ana=addfile("<a href="http://temp.nc">temp.nc</a>","r")</div><div>T  = ana->T</div><div>TS  = ana->TS</div><div>ps=ana->ps</div><div>delete(ana)</div><div><br></div><div><br></div><div><br></div><div>ana=addfile("<a href="http://F_2000_CAM5_PM_fine_no_ds.cam.h0.0008-09.nc">F_2000_CAM5_PM_fine_no_ds.cam.h0.0008-09.nc</a>","r")</div><div>P0mb=ana->P0/100</div><div>;ps=ana->PS</div><div>hyam=ana->hyam</div><div>hybm=ana->hybm</div><div>;lev=ana->lev</div><div>lat=ana->lat</div><div>delete(ana)</div><div><br></div><div>  lev           = (/  10, 20, 30, 50, 70,100,150,200,250 \</div><div>                     , 300,400,500,600,700,850,925,1000 /) </div><div><br></div><div>  lev!0         = "lev"                  ; variable/dim name </div><div>  lev&lev     =  lev                   ; create coordinate variable</div><div>  lev@long_name = "pressure"               ; attach some attributes</div><div>  lev@units     = "hPa"</div><div>  lev@positive  = "down"</div><div><br></div><div> t_pres   = vinth2p(T, hyam, hybm, lev, ps, 1, P0mb, 1, False )</div></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div>and for GFDL I changed it to this:</div><div><br></div><div><div><br></div><div>P0mb=1000</div><div>hyam=pk/100000.</div><div>hybm=bk</div><div><br></div><div>  lev           = (/  10, 20, 30, 50, 70,100,150,200,250 \</div><div>                     , 300,400,500,600,700,850,925,1000 /) </div><div><br></div><div>  lev!0         = "lev"                  ; variable/dim name </div><div>  lev&lev     =  lev                   ; create coordinate variable</div><div>  lev@long_name = "pressure"               ; attach some attributes</div><div>  lev@units     = "hPa"</div><div>  lev@positive  = "down"</div><div><br></div><div> t_pres   = vinth2p(T, hyam, hybm, lev, ps, 1, P0mb, 1, False )</div></div><div><br></div><div><br></div></div>