<div class="gmail_quote">On Thu, Mar 27, 2008 at 12:27 AM, Hamza Varikoden <<a href="mailto:hamzavarikoden@gmail.com">hamzavarikoden@gmail.com</a>> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div>Dear Arlindo<br> </div>
<div>I gone through the discussion on stream function and velocity potential in the grads user forum. Based on that I tried it and got the results..</div>
<div> </div>
<div>But when I am trying for a limited domain (MM5 out put, it latitutes are not linear) it shows some problem. can you suggest a suitable solution to rectify this ? .. I here with enclose the ctl file and desciption of my grads out put. The output which i got is error from fish:lat/lon coords must be LINEAR.</div>
<div> </div></blockquote><div><br></div><div>The current poisson solver in fish() is for constant lat/lon on a global domain only. We need someone to help develop a regional version. Any volunteers? I just need someone to contribute a *self-contained* C or Fortran routine that takes vorticity/divergence and the coordinates and returns the streamfunction/vel potential on a regional domain. I'll then wrap it up as an extension.</div>
<div><br></div><div> Arlindo</div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div></div>
<div>CTL file</div>
<div>======</div><font size="2">
<p>dset ^20070115t12_12km.dat</p>
<p>byteswapped</p>
<p>title MM5 data</p>
<p>undef -9999.</p>
<p>xdef 219 linear 97.6112 0.1079</p>
<p>ydef 129 levels</p>
<p>-2.094353 </p>
<p>-1.986481 </p>
<p>-1.878609 </p>
<p>-1.770721 </p>
<p>-1.662834 </p>
<p>-1.554947 </p>
<p>-1.447044 </p>
<p>-1.339142 </p>
<p>-1.231232 </p>
<p>-1.123322 </p>
<p>-1.015404 </p>
<p>-0.9074860 </p>
<p>-0.7995605 </p>
<p>-0.6916351 </p>
<p>-0.5837021 </p>
<p>-0.4757767 </p>
<p>-0.3678436 </p>
<p>-0.2599030 </p>
<p>-0.1519775 </p>
<p>-4.4036865E-02</p>
<p>6.3896179E-02</p>
<p>0.1718369 </p>
<p>0.2797699 </p>
<p>0.3877029 </p>
<p>0.4956360 </p>
<p>0.6035690 </p>
<p>0.7114868 </p>
<p>0.8194199 </p>
<p>0.9273453 </p>
<p>1.035263 </p>
<p>1.143173 </p>
<p>1.251091 </p>
<p>1.359001 </p>
<p>1.466904 </p>
<p>1.574799 </p>
<p>1.682686 </p>
<p>1.790573 </p>
<p>1.898453 </p>
<p>2.006332 </p>
<p>2.114189 </p>
<p>2.222046 </p>
<p>2.329903 </p>
<p>2.437744 </p>
<p>2.545570 </p>
<p>2.653404 </p>
<p>2.761215 </p>
<p>2.869019 </p>
<p>2.976814 </p>
<p>3.084602 </p>
<p>3.192383 </p>
<p>3.300140 </p>
<p>3.407883 </p>
<p>3.515625 </p>
<p>3.623352 </p>
<p>3.731071 </p>
<p>3.838768 </p>
<p>3.946449 </p>
<p>4.054131 </p>
<p>4.161781 </p>
<p>4.269424 </p>
<p>4.377052 </p>
<p>4.484665 </p>
<p>4.592270 </p>
<p>4.699844 </p>
<p>4.807411 </p>
<p>4.914955 </p>
<p>5.022484 </p>
<p>5.130005 </p>
<p>5.237495 </p>
<p>5.344971 </p>
<p>5.452431 </p>
<p>5.559860 </p>
<p>5.667274 </p>
<p>5.774673 </p>
<p>5.882050 </p>
<p>5.989410 </p>
<p>6.096748 </p>
<p>6.204063 </p>
<p>6.311356 </p>
<p>6.418625 </p>
<p>6.525871 </p>
<p>6.633102 </p>
<p>6.740303 </p>
<p>6.847481 </p>
<p>6.954628 </p>
<p>7.061760 </p>
<p>7.168869 </p>
<p>7.275948 </p>
<p>7.382996 </p>
<p>7.490021 </p>
<p>7.597023 </p>
<p>7.704002 </p>
<p>7.810951 </p>
<p>7.917870 </p>
<p>8.024757 </p>
<p>8.131622 </p>
<p>8.238464 </p>
<p>8.345268 </p>
<p>8.452049 </p>
<p>8.558800 </p>
<p>8.665512 </p>
<p>8.772202 </p>
<p>8.878860 </p>
<p>8.985489 </p>
<p>9.092087 </p>
<p>9.198647 </p>
<p>9.305183 </p>
<p>9.411674 </p>
<p>9.518143 </p>
<p>9.624573 </p>
<p>9.730972 </p>
<p>9.837341 </p>
<p>9.943665 </p>
<p>10.04996 </p>
<p>10.15623 </p>
<p>10.26246 </p>
<p>10.36864 </p>
<p>10.47480 </p>
<p>10.58092 </p>
<p>10.68700 </p>
<p>10.79305 </p>
<p>10.89906 </p>
<p>11.00502 </p>
<p>11.11095 </p>
<p>11.21685 </p>
<p>11.32270 </p>
<p>11.42852 </p>
<p>11.53429 </p>
<p>11.64002 </p>
<p>zdef 20 levels </p>
<p>1000.00</p>
<p>950.00</p>
<p>925.00</p>
<p>900.00</p>
<p>850.00</p>
<p>800.00</p>
<p>750.00</p>
<p>700.00</p>
<p>650.00</p>
<p>600.00</p>
<p>550.00</p>
<p>500.00</p>
<p>450.00</p>
<p>400.00</p>
<p>350.00</p>
<p>300.00</p>
<p>250.00</p>
<p>200.00</p>
<p>150.00</p>
<p>100.00</p>
<p>tdef 73 linear 12:00Z15JAN2007 60MN</p>
<p>vars 62</p>
<p>ps 0 99 pstar (Pa) </p>
<p>tg 0 99 ground temp (K) </p>
<p>rc 0 99 accum conv pcn (cm) </p>
<p>rn 0 99 accum non-c pcn (cm)</p>
<p>ter 0 99 ter elevation (m) </p>
<p>xmf 0 99 cross map factor </p>
<p>dmf 0 99 dot map factor </p>
<p>cor 0 99 coriolis (s-1) </p>
<p>tr 0 99 reservoir temp (K) </p>
<p>xlat 0 99 cross lat (degree) </p>
<p>xlon 0 99 cross lon (degree) </p>
<p>lu 0 99 land use </p>
<p>pblh 0 99 pbl height (m) </p>
<p>pblr 0 99 pbl regime </p>
<p>shf 0 99 sen heat flux (W/m2)</p>
<p>lhf 0 99 lat heat flux (W/m2)</p>
<p>ust 0 99 friction vel (m/s) </p>
<p>swd 0 99 down sw rad (W/m2) </p>
<p>lwd 0 99 down lw rad (W/m2) </p>
<p>swo 0 99 out sw rad (W/m2) </p>
<p>lwo 0 99 out lw rad (W/m2) </p>
<p>sst 0 99 sea sfc temp </p>
<p>st1 0 99 soil temp 1 (K) </p>
<p>st2 0 99 soil temp 2 (K) </p>
<p>st3 0 99 soil temp 3 (K) </p>
<p>st4 0 99 soil temp 4 (K) </p>
<p>st5 0 99 soil temp 5 (K) </p>
<p>st6 0 99 soil temp 6 (K) </p>
<p>t2m 0 99 2 m temperature (K) </p>
<p>q2m 0 99 2m mix ratio (kg/kg)</p>
<p>u10 0 99 10 m u wind (m/sec) </p>
<p>v10 0 99 10 m v wind (m/sec) </p>
<p>pslv 0 99 sea level prs (mb) </p>
<p>cref 0 99 composite refl (dbz)</p>
<p>iclw 0 99 integrat clw (cm) </p>
<p>irnw 0 99 integrat rnw (cm) </p>
<p>pwat 0 99 precipit water (cm) </p>
<p>clfrlo 0 99 low cloud fraction </p>
<p>clfrmi 0 99 mid cloud fraction </p>
<p>clfrhi 0 99 high cloud fraction </p>
<p>u 20 99 u wind (m/s) </p>
<p>v 20 99 v wind (m/s) </p>
<p>w 20 99 vertical vel (m/s) </p>
<p>pp 20 99 prs pert (Pa) </p>
<p>t 20 99 temperature (C) </p>
<p>q 20 99 mixing ratio (kg/kg)</p>
<p>clw 20 99 cloud water (kg/kg) </p>
<p>rnw 20 99 rain water (kg/kg) </p>
<p>rtnd 20 99 rad tend (K/day) </p>
<p>z 20 99 height AGL (m) </p>
<p>h 20 99 geopot height (m) </p>
<p>td 20 99 dewpoint temp (C) </p>
<p>rh 20 99 rel humidity (%) </p>
<p>th 20 99 potential temp (K) </p>
<p>the 20 99 theta-e (K) </p>
<p>prs 20 99 pressure (Pa) </p>
<p>vor 20 99 vorticity (s-1) </p>
<p>pv 20 99 potential vort (pvu)</p>
<p>dbz 20 99 reflectivity (dbz) </p>
<p>div 20 99 divergence (s-1) </p>
<p>dir 20 99 dir (degrees) </p>
<p>tadv 20 99 temp advection (K/s)</p>
<p>endvars</p></font>
</blockquote></div><br><br clear="all"><br>-- <br>Arlindo da Silva<br><a href="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu</a>