problem using fish function

Thu Mar 27 09:08:10 EDT 2008

On Thu, Mar 27, 2008 at 12:27 AM, Hamza Varikoden <hamzavarikoden at gmail.com>
wrote:

> Dear  Arlindo
>
> 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..
>
> 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.
>
>

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.

        Arlindo

> CTL file
> ======
>
> dset ^20070115t12_12km.dat
>
> byteswapped
>
> title MM5 data
>
> undef -9999.
>
> xdef 219 linear 97.6112 0.1079
>
> ydef 129 levels
>
> -2.094353
>
> -1.986481
>
> -1.878609
>
> -1.770721
>
> -1.662834
>
> -1.554947
>
> -1.447044
>
> -1.339142
>
> -1.231232
>
> -1.123322
>
> -1.015404
>
> -0.9074860
>
> -0.7995605
>
> -0.6916351
>
> -0.5837021
>
> -0.4757767
>
> -0.3678436
>
> -0.2599030
>
> -0.1519775
>
> -4.4036865E-02
>
> 6.3896179E-02
>
> 0.1718369
>
> 0.2797699
>
> 0.3877029
>
> 0.4956360
>
> 0.6035690
>
> 0.7114868
>
> 0.8194199
>
> 0.9273453
>
> 1.035263
>
> 1.143173
>
> 1.251091
>
> 1.359001
>
> 1.466904
>
> 1.574799
>
> 1.682686
>
> 1.790573
>
> 1.898453
>
> 2.006332
>
> 2.114189
>
> 2.222046
>
> 2.329903
>
> 2.437744
>
> 2.545570
>
> 2.653404
>
> 2.761215
>
> 2.869019
>
> 2.976814
>
> 3.084602
>
> 3.192383
>
> 3.300140
>
> 3.407883
>
> 3.515625
>
> 3.623352
>
> 3.731071
>
> 3.838768
>
> 3.946449
>
> 4.054131
>
> 4.161781
>
> 4.269424
>
> 4.377052
>
> 4.484665
>
> 4.592270
>
> 4.699844
>
> 4.807411
>
> 4.914955
>
> 5.022484
>
> 5.130005
>
> 5.237495
>
> 5.344971
>
> 5.452431
>
> 5.559860
>
> 5.667274
>
> 5.774673
>
> 5.882050
>
> 5.989410
>
> 6.096748
>
> 6.204063
>
> 6.311356
>
> 6.418625
>
> 6.525871
>
> 6.633102
>
> 6.740303
>
> 6.847481
>
> 6.954628
>
> 7.061760
>
> 7.168869
>
> 7.275948
>
> 7.382996
>
> 7.490021
>
> 7.597023
>
> 7.704002
>
> 7.810951
>
> 7.917870
>
> 8.024757
>
> 8.131622
>
> 8.238464
>
> 8.345268
>
> 8.452049
>
> 8.558800
>
> 8.665512
>
> 8.772202
>
> 8.878860
>
> 8.985489
>
> 9.092087
>
> 9.198647
>
> 9.305183
>
> 9.411674
>
> 9.518143
>
> 9.624573
>
> 9.730972
>
> 9.837341
>
> 9.943665
>
> 10.04996
>
> 10.15623
>
> 10.26246
>
> 10.36864
>
> 10.47480
>
> 10.58092
>
> 10.68700
>
> 10.79305
>
> 10.89906
>
> 11.00502
>
> 11.11095
>
> 11.21685
>
> 11.32270
>
> 11.42852
>
> 11.53429
>
> 11.64002
>
> zdef 20 levels
>
> 1000.00
>
> 950.00
>
> 925.00
>
> 900.00
>
> 850.00
>
> 800.00
>
> 750.00
>
> 700.00
>
> 650.00
>
> 600.00
>
> 550.00
>
> 500.00
>
> 450.00
>
> 400.00
>
> 350.00
>
> 300.00
>
> 250.00
>
> 200.00
>
> 150.00
>
> 100.00
>
> tdef 73 linear 12:00Z15JAN2007 60MN
>
> vars 62
>
> ps 0 99 pstar (Pa)
>
> tg 0 99 ground temp (K)
>
> rc 0 99 accum conv pcn (cm)
>
> rn 0 99 accum non-c pcn (cm)
>
> ter 0 99 ter elevation (m)
>
> xmf 0 99 cross map factor
>
> dmf 0 99 dot map factor
>
> cor 0 99 coriolis (s-1)
>
> tr 0 99 reservoir temp (K)
>
> xlat 0 99 cross lat (degree)
>
> xlon 0 99 cross lon (degree)
>
> lu 0 99 land use
>
> pblh 0 99 pbl height (m)
>
> pblr 0 99 pbl regime
>
> shf 0 99 sen heat flux (W/m2)
>
> lhf 0 99 lat heat flux (W/m2)
>
> ust 0 99 friction vel (m/s)
>
> swd 0 99 down sw rad (W/m2)
>
> lwd 0 99 down lw rad (W/m2)
>
> swo 0 99 out sw rad (W/m2)
>
> lwo 0 99 out lw rad (W/m2)
>
> sst 0 99 sea sfc temp
>
> st1 0 99 soil temp 1 (K)
>
> st2 0 99 soil temp 2 (K)
>
> st3 0 99 soil temp 3 (K)
>
> st4 0 99 soil temp 4 (K)
>
> st5 0 99 soil temp 5 (K)
>
> st6 0 99 soil temp 6 (K)
>
> t2m 0 99 2 m temperature (K)
>
> q2m 0 99 2m mix ratio (kg/kg)
>
> u10 0 99 10 m u wind (m/sec)
>
> v10 0 99 10 m v wind (m/sec)
>
> pslv 0 99 sea level prs (mb)
>
> cref 0 99 composite refl (dbz)
>
> iclw 0 99 integrat clw (cm)
>
> irnw 0 99 integrat rnw (cm)
>
> pwat 0 99 precipit water (cm)
>
> clfrlo 0 99 low cloud fraction
>
> clfrmi 0 99 mid cloud fraction
>
> clfrhi 0 99 high cloud fraction
>
> u 20 99 u wind (m/s)
>
> v 20 99 v wind (m/s)
>
> w 20 99 vertical vel (m/s)
>
> pp 20 99 prs pert (Pa)
>
> t 20 99 temperature (C)
>
> q 20 99 mixing ratio (kg/kg)
>
> clw 20 99 cloud water (kg/kg)
>
> rnw 20 99 rain water (kg/kg)
>
> rtnd 20 99 rad tend (K/day)
>
> z 20 99 height AGL (m)
>
> h 20 99 geopot height (m)
>
> td 20 99 dewpoint temp (C)
>
> rh 20 99 rel humidity (%)
>
> th 20 99 potential temp (K)
>
> the 20 99 theta-e (K)
>
> prs 20 99 pressure (Pa)
>
> vor 20 99 vorticity (s-1)
>
> pv 20 99 potential vort (pvu)
>
> dbz 20 99 reflectivity (dbz)
>
> div 20 99 divergence (s-1)
>
> dir 20 99 dir (degrees)
>
> tadv 20 99 temp advection (K/s)
>
> endvars
>

--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20080327/fd57a7e5/attachment.html 