<div dir="ltr"><div><div>Dear Sir <br></div> Presently I am using NCAR Command Language Version 6.1.2. I have made the anomaly from this version and from pyferret and after compression we have got a lot of differences.<br></div>I am mentioning the script use for NCAR to make anomaly<br><br><br><div>;<<<<<<<<<<<<<<<<<<<<<<<<<<program for making the anomaly <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<!<br>load "/usr/share/ncarg/nclscripts/csm/gsn_code.ncl"<br>load "/usr/share/ncarg/nclscripts/csm/gsn_csm.ncl"<br>load "/usr/share/ncarg/nclscripts/csm/contributed.ncl"<br>load "/usr/share/ncarg/nclscripts/csm/shea_util.ncl"<br><br>;<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><br>begin<br>;<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><br> fili = "/home/abhishek/Desktop/work/phd/Ridge_SCTR/<a href="http://soda_d20_1958_2010_ind.nc">soda_d20_1958_2010_ind.nc</a>" ; define filename<br> f = addfile (fili, "r")<br> d20=f->D20<br>; print(d20)<br>; printVarSummary(d20)<br> anom = d20<br> anom = dim_rmvmean(d20)<br>; printVarSummary(anom )<br><br>;<<<<<<<<<<<<<<<<<<<<<<remove any preexicting file <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<!<br> system("/bin/rm -f <a href="http://soda_d20_anm.nc">soda_d20_anm.nc</a>")<br><br>;<<<<<<<<<<<<<<<<<<<< open outputnetcdf file <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<br>ncdf = addfile("<a href="http://soda_d20_anm.nc">soda_d20_anm.nc</a>" ,"c") ; open output netCDF file<br><br>;<<<<<<<<<<<<<<<<<<<<<<<<<for giveing the information of the file <<<<<<<<<<<<<<<<<<<<<<<<<<<!<br>fAtt = True ; assign file attributes<br> fAtt@title = "NCL Simple Approach to netCDF Creation"<br> fAtt@source_file = "<a href="http://original-file.nc">original-file.nc</a>"<br> fAtt@Conventions = "None"<br> fAtt@creation_date = systemfunc ("date")<br> fileattdef( ncdf, fAtt ) ; copy file attributes<br><br>;<<<<<<<<<<<<<<<<<<<<<<< make time an UNLIMITED dimension; recommended for most applications<<<<<<<<<<<<<<<<<<<<<<<<<br> filedimdef(ncdf,"time",-1,True)<br>;<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><br> ; output variables directly; NCL will call appropriate function to write the meta data associated with each variable<br> ncdf->anom = anom<br>end<br><br></div><div>Eagerly waiting for your reply<br></div><div>regards<br></div><div>Abhishek<br></div><div><br><div><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><span style="color:rgb(0,0,127);font-family:"times new roman","new york",times,serif;font-size:16px;font-weight:bold">Regards</span><br style="color:rgb(0,0,127);font-family:"times new roman","new york",times,serif;font-size:16px;font-weight:bold"><br style="color:rgb(0,0,127);font-family:"times new roman","new york",times,serif;font-size:16px;font-weight:bold"><font face="times new roman, new york, times, serif" size="3" color="#00007f"><b>Abhishek Savita</b></font></div><div dir="ltr"><font face="times new roman, new york, times, serif" size="3" color="#00007f"><b><br></b></font></div><div dir="ltr"><span style="color:rgb(0,0,127);font-family:"times new roman","new york",times,serif;font-size:16px;font-weight:bold">Research Scholar (Earth </span><span style="color:rgb(0,0,127);font-family:"times new roman","new york",times,serif;font-size:16px;font-weight:bold">System Science Technology)</span></div><div dir="ltr"><font face="times new roman, new york, times, serif" size="3" color="#00007f"><b>Center For Oceans, Revers, Atmoshphere & Land Science Technology</b></font><br style="color:rgb(0,0,127);font-family:"times new roman","new york",times,serif;font-size:16px;font-weight:bold"><span style="color:rgb(0,0,127);font-family:"times new roman","new york",times,serif;font-size:16px;font-weight:bold">Indian Institute of Technology, Kharagpur</span><br style="color:rgb(0,0,127);font-family:"times new roman","new york",times,serif;font-size:16px;font-weight:bold"><span style="color:rgb(0,0,127);font-family:"times new roman","new york",times,serif;font-size:16px;font-weight:bold">+91-8609704619</span><div><font face="times new roman, new york, times, serif" size="3" color="#00007f"><b><br></b></font><div><span style="color:rgb(0,0,127);font-family:"times new roman","new york",times,serif;font-size:16px;font-weight:bold"> </span></div></div></div></div></div></div>
</div></div></div></div>