<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hello,<div class=""><br class=""></div><div class="">I created a netcdf output file. But I am not able to get proper dimensions. I think I need proper dimensions to perform nco operations on the .nc file. I saw an NCL page on creating the output file, but it assumes the output variables already have dimensions and no meta data. <a href="https://www.ncl.ucar.edu/Applications/o-netcdf.shtml" class="">https://www.ncl.ucar.edu/Applications/o-netcdf.shtml</a></div><div class=""><br class=""></div><div class=""><b class="">Here I need to create dimensions and meta data.</b></div><div class="">Could you please help me to get proper dimensions with the output variable, so that I can use NCO tool. I have attached ncl code, ncdump file and NCO script for suggestion. </div><div class=""><br class=""></div><div class="">Here is my code:</div><div class=""><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">;================================================;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">; ================================================;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">begin</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">;=================================================;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124); min-height: 13px;" class=""><br class=""></div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">datadir_00 = "/glade/scratch/sharmaa/BLUE_WATERS_PROJECT/nco_wrf2000/"</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">FILES_00 = systemfunc ("ls -1 " + datadir_00 + "2000d01_day00* " )</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">numFILES_00 = dimsizes(FILES_00)</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">print(FILES_00)</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124); min-height: 13px;" class=""><br class=""></div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">sn = 515</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">we = 599</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124); min-height: 13px;" class=""><br class=""></div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">dailyT2maxS_00 = new((/numFILES_00,sn,we/),float)</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">dailyT2minS_00 = new((/numFILES_00,sn,we/),float)</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">dailyT2avS_00 = new((/numFILES_00,sn,we/),float)</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124); min-height: 13px;" class=""><br class=""></div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">do ifil = 0,numFILES_00-1</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> a_00 = addfile(FILES_00(ifil)+".nc","r") ; Open the next file</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> temp_00 = a_00->T2 ; perturbation potential temperature (theta)</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> temp_00 = temp_00 - 273.15</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124); min-height: 13px;" class=""><br class=""></div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> dailyT2maxS_00(ifil,:,:) = dim_max_n_Wrap(temp_00,0)</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> dailyT2minS_00(ifil,:,:) = dim_min_n_Wrap(temp_00,0)</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> dailyT2avS_00(ifil,:,:) = (dailyT2maxS_00(ifil,:,:)+dailyT2minS_00(ifil,:,:))/2</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124); min-height: 13px;" class=""><br class=""></div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">end do</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124); min-height: 13px;" class=""><br class=""></div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">;=================================================;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">f1 = addfile("YEARd01_day001.nc","r")</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">LAT = f1->XLAT(0,:,:)</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">LON = f1->XLONG(0,:,:)</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">;=================================================;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">;=================================================;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">;OPEN NEW NETCDF FOR WRITING</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">;-------------------------------------------------</div></div><div class=""><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124); min-height: 13px;" class=""><br class=""></div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> system("/bin/rm -f WRF_T2_d01_2000_test.nc") ; remove any pre-existing file</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> fout = addfile("WRF_T2_d01_2000_test.nc" ,"c") ; open output netCDF file</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">;===================================================================</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> ; create global attributes of the file (optional)</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> ;===================================================================</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> fAtt = True ; assign file attributes</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> fAtt@title = "NCL Simple Approach to netCDF Creation"</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> fAtt@source_file = "original-file.nc"</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> fAtt@Conventions = "None"</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> fAtt@creation_date = systemfunc ("date")</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> fileattdef( fout, fAtt ) ; copy file attributes</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124); min-height: 13px;" class=""><br class=""></div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> ;===================================================================</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> ; make time an UNLIMITED dimension; recommended for most applications</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> ;===================================================================</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> filedimdef(fout,"time",-1,True)</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124); min-height: 13px;" class=""><br class=""></div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> ;===================================================================</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> ; output variables directly; NCL will call appropriate functions</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> ; to write the meta data associated with each variable</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> ;===================================================================</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124); min-height: 13px;" class=""><br class=""></div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> XLAT = f1->XLAT</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> XLONG = f1->XLONG</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124); min-height: 13px;" class=""><br class=""></div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> fout->XLAT = XLAT</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> fout->XLONG = XLONG</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""> fout->T2_AV = dailyT2avS_00</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124); min-height: 13px;" class=""><br class=""></div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">end</div></div><div class=""><br class=""></div><div class="">Here is the ncdump file:</div><div class=""><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">[sharmaa@yslogin1 T2_SIM]$ ncdump -h WRF_T2_d01_2000_test.nc </div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">netcdf WRF_T2_d01_2000_test {</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">dimensions:</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span>time = UNLIMITED ; // (0 currently)</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span>Time = 24 ;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span>south_north = 515 ;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span>west_east = 599 ;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span>ncl4 = 9 ;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span>ncl5 = 515 ;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span>ncl6 = 599 ;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">variables:</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span>float XLAT(Time, south_north, west_east) ;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span>XLAT:FieldType = 104 ;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span>XLAT:MemoryOrder = "XY " ;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span>XLAT:description = "LATITUDE, SOUTH IS NEGATIVE" ;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span>XLAT:units = "degree_north" ;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span>XLAT:stagger = "" ;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span>float XLONG(Time, south_north, west_east) ;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span>XLONG:FieldType = 104 ;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span>XLONG:MemoryOrder = "XY " ;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span>XLONG:description = "LONGITUDE, WEST IS NEGATIVE" ;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span>XLONG:units = "degree_east" ;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span>XLONG:stagger = "" ;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span><b class="">float T2_AV(ncl4, ncl5, ncl6) ;</b></div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><b class=""><span class="Apple-tab-span" style="white-space:pre"> </span>T2_AV:_FillValue = 9.96921e+36f ;</b></div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124); min-height: 13px;" class=""><br class=""></div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">// global attributes:</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span>:creation_date = "Mon Nov 2 12:12:21 MST 2015" ;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span>:Conventions = "None" ;</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class=""><span class="Apple-tab-span" style="white-space:pre"> </span>:source_file = "original-file.nc" ;</div></div><div class=""><br class=""></div><div class="">Here is the issue with NCO operation:</div><div class=""><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">[sharmaa@yslogin1 T2_SIM]$ ncrcat -d ncl4,1,5 -v T2_AV WRF_T2_d01_2000_test.nc test.nc</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">ncrcat: ERROR no variables fit criteria for processing</div><div style="margin: 0px; font-size: 11px; font-family: Menlo; background-color: rgb(163, 152, 124);" class="">ncrcat: HINT Extraction list must contain a record variable which to concatenate</div></div><div class=""><br class=""></div><div class=""><br class=""></div></body></html>