On Fri, Apr 4, 2008 at 1:10 PM, Semyon Grodsky <<a href="mailto:senya@atmos.umd.edu">senya@atmos.umd.edu</a>> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Arlindo,<br>
thank you very much. It works.<div class="Ih2E3d"><br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
ga-> open ../sst/sst_month<br>
ga-> set x 1 360<br>
ga-> sst1 = sst - 273<br>
ga-> lats4d -o qq -func sst1 -time nov1981 dec2007 -vars sst<br>
<br>
</blockquote></div>
Hopefuly, the last question. Suppose I have two defined variables, sst1<br>
and sst2. What modifications should be done to the above lats4d command in<br>
order to output two variables in a single netcdf.<br>
</blockquote><div><br>Create a fake ctl file:<br><br>dset ^fake.data<br>options template<br>xdef ...<br>ydef ...<br>zdef ...<br>tdef ...<br>vars 2<br>sst1 0 0 whatever<br>sst2 0 0 whatever<br>endvars<br><br>If you open this file you always get missing values for whatever you try to plot; this is OK. Define the *def records above to match your original file; the easiest way of doing this is to open your original file and do a<br>
<br>
ga-> q ctlinfo<br><br>Once you create a fake file with the structure and variable names of what you want on output, then<br><br>ga-> open fake.ctl<br>
ga-> open ../sst/sst_month<br>
ga-> set x 1 360<br>
ga-> dsst1 = sst.2 - 273<br>ga-> dsst2 = sst.2 / 3<br>ga-> set dfile 1 (to make sure)<br>ga-> lats4d -o qq -func d@<br><br>So it loops over the time/variables on the fake file but saves the defined variables dsst1 and dsst2; the output file will have variable names ss1, sst2. If you want more than 1 time level on output try this variant:<br>
<br>ga-> open fake.ctl<br>
ga-> open ../sst/sst_month<br>
ga-> set x 1 360<br>ga-> set t 1 12<br>ga-> dsst1 = sst.2 - 273 (this takes memory)<br>
ga-> dsst2 = sst.2 / 3<br>ga-> set t 1<br>
ga-> set dfile 1 (to make sure)<br>
ga-> lats4d -o qq -func d@<br><br>and have the TDEF match the number of times above.<br><br> Arlindo<br><br><br><br><br><br><br><br><br><br><br> </div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
--Senya<br>
</blockquote></div><br><br clear="all"><br>-- <br>Arlindo da Silva<br><a href="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu</a>