<div dir="ltr">Dear WRF-experts,<div><br></div><div>I am trying to calculate cloud top temperature from a raw wrfout file and save it as a separate netcdf file that can be processed by GRADS and CDO.</div><div><br></div><div><br></div><div>Here is the script that I am using:</div><div><br></div><div><div>a = addfile("<a href="http://wrfout_d01_2013-08-16_00_00_00.nc">wrfout_d01_2013-08-16_00_00_00.nc</a>","r")</div><div><br></div><div>ctt   = wrf_user_getvar(a,"ctt",-1)    ; Time X Lat X Lon</div><div><br></div><div>Lat = wrf_user_getvar(a,"XLAT",0)</div><div>Lon = wrf_user_getvar(a,"XLONG",0)</div><div><br></div><div>printVarSummary(ctt)</div><div><br></div><div>system("/bin/rm -f <a href="http://test.nc">test.nc</a>")   ; remove any pre-existing file</div><div>   ncdf = addfile("<a href="http://test.nc">test.nc</a>" ,"c")  ; open output netCDF file</div><div>   fAtt               = True            ; assign file attributes</div><div>   fAtt@title         = "Cloud Top Temperature degC"</div><div>   fAtt@source_file   =  "wrfoutfile"</div><div>   fAtt@Conventions   = "None"</div><div>   fAtt@creation_date = systemfunc ("date")</div><div>   fileattdef( ncdf, fAtt )            ; copy file attributes</div><div>   filedimdef(ncdf,"time",-1,True)</div><div><br></div><div>ncdf->ctt = ctt</div><div>ncdf->Lat = Lat</div><div>ncdf->Lon = Lon</div></div><div><br></div><div><br></div><div><br></div><div><b>Problem</b></div><div><b>[1] </b>The output file cannot be read by GRADS/CDO. Am I missing something important when writing out the variable?</div><div><br></div><div><b>[2]</b> Any suggestion on how to do this correctly?</div><div><br></div><div><br></div><div>I'll appreciate any help.</div><div><br></div><div>Sincerely,</div><div><br></div><div><br></div><div><br clear="all"><div><div class="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"><font face="comic sans ms, sans-serif"><b><u>Lyndon Mark P. Olaguera</u></b></font><div>PhD Student </div><div>Laboratory of Climatology</div><div>Department of Geography</div><div>Faculty of Urban Environmental Science</div><div><span style="font-size:12.8px">Minami-Osawa Campus</span><br></div><div>Tokyo Metropolitan University</div><div>Hachioji Tokyo, Japan</div><div><br></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
</div></div>