[gradsusr] staggered Arakawa C-grid

Charles Seman (NOAA Federal) charles.seman at noaa.gov
Tue Nov 6 13:44:59 EST 2012


Gabriel,

In the Arakawa C-grid, because the u and v components are on grids that 
are staggered from one another and also from the mass grid, one ctl file 
can't represent them all.  One thing to try is to make 3 ctl files; one 
for "u", one for "v", and one for "mass" grid variables...  This will 
allow you to view each wind component and the mass variables 
separately... But to do any calculations in GrADS with the "mass" 
(center) grid point, I don't know of any way to do this except if all 
variables are on a common grid.  For this, you can average both "u" and 
"v" to the "mass" grid (assuming salt and temperature are defined on the 
"mass" grid).  One way: if you're using GrADS v2, you can define "um", 
"vm" components averaged to the mass grid in GrADS, and use sdfwrite to 
write out component "um" and "vm" files.  This can also be done in a 
Fortran (or C) program if you like that better.

Hope this helps,
Chuck

On 11/06/2012 12:53 PM, Gabriel Bonow Munchow wrote:
> Hello,
>
> Can someone give some helṕ on how to open a nc that is in staggered
> Arakawa C-grid ? I have tried this ctl:
>
>
> DSET ^tsa1_avg_old600.nc <http://tsa1_avg_old600.nc>
> DTYPE netcdf
> TITLE 4-D Ocean Variables: Tracers Components
> UNDEF 1.e+37f _FillValue
> OPTIONS zrev
>
> XDEF 262 LEVELS -65.50  -65.3333333333333  -65.1666666666667  -65.00
>   -64.8333333333333      -64.66666666
> 66667  -64.50  -64.3333333333333  -64.1666666666667  -64.00
>   -63.8333333333333   -63.6666666666667  -63.5
> 0  -63.3333333333333  -63.1666666666667  -63.00  -62.8333333333333
> -62.6666666666667  -62.50  -62.33333
> 33333333  -62.1666666666667  -62.00  -61.8333333333333
>   -61.6666666666667  -61.50  -61.3333333333333  -61
> .1666666666667  -61.00  -60.8333333333333   -60.6666666666667  -60.50
>   -60.3333333333333  -60.16666666666
> 67  -60.00  -59.8333333333333   -59.6666666666667  -59.50
>   -59.3333333333333  -59.1666666666667  -59.00
> -58.8333333333333  -58.6666666666667  -58.50  -58.3333333333333
>   -58.1666666666667  -58.00  -57.833333333
> 3333   -57.6666666666667  -57.50  -57.3333333333333  -57.1666666666667
>   -57.00  -56.8333333333333   -56.6
> 666666666667  -56.50  -56.3333333333333  -56.1666666666667  -56.00
>   -55.8333333333333   -55.6666666666667
>    -55.50  -55.3333333333333  -55.1666666666667  -55.00
>   -54.8333333333333   -54.6666666666667  -54.50  -5
> 4.3333333333333  -54.1666666666667  -54.00  -53.8333333333333
> -53.6666666666667  -53.50  -53.3333333333
> 333  -53.1666666666667  -53.00  -52.8333333333333   -52.6666666666667
>   -52.50  -52.3333333333333  -52.166
> 6666666667  -52.00  -51.8333333333333   -51.6666666666667  -51.50
>   -51.3333333333333  -51.1666666666667
> -51.00  -50.8333333333333   -50.6666666666667  -50.50  -50.3333333333333
>   -50.1666666666667  -50.00  -49.
> 8333333333333   -49.6666666666667  -49.50  -49.3333333333333
>   -49.1666666666667  -49.00  -48.833333333333
> 3   -48.6666666666667  -48.50  -48.3333333333333  -48.1666666666667
>   -48.00  -47.8333333333333   -47.6666
> 666666667  -47.50  -47.3333333333333  -47.1666666666667  -47.00
>   -46.8333333333333   -46.6666666666667  -
> 46.50  -46.3333333333333  -46.1666666666667  -46.00  -45.8333333333333
> -45.6666666666667  -45.50  -45.3
> 333333333333  -45.1666666666667  -45.00  -44.8333333333333
> -44.6666666666667  -44.50  -44.3333333333333
>    -44.1666666666667  -44.00  -43.8333333333333   -43.6666666666667
>   -43.50  -43.3333333333333  -43.166666
> 6666667  -43.00  -42.8333333333333   -42.6666666666667  -42.50
>   -42.3333333333333  -42.1666666666667  -42
> .00  -41.8333333333333   -41.6666666666667  -41.50  -41.3333333333333
>   -41.1666666666667  -41.00  -40.833
> 3333333333   -40.6666666666667  -40.50  -40.3333333333333
>   -40.1666666666667  -40.00  -39.8333333333333
>   -39.6666666666667  -39.50  -39.3333333333333  -39.1666666666667
>   -39.00  -38.8333333333333   -38.6666666
> 666667  -38.50  -38.3333333333333  -38.1666666666667  -38.00
>   -37.8333333333333   -37.6666666666667  -37.
> 50  -37.3333333333333  -37.1666666666667  -37.00  -36.8333333333333
> -36.6666666666667  -36.50  -36.3333
> 333333333  -36.1666666666667  -36.00  -35.8333333333333
> -35.6666666666667  -35.50  -35.3333333333333  -
> 35.1666666666667  -35.00  -34.8333333333333   -34.6666666666667  -34.50
>   -34.3333333333333  -34.166666666
> 6667  -34.00  -33.8333333333333   -33.6666666666667  -33.50
>   -33.3333333333333  -33.1666666666667  -33.00
>    -32.8333333333333   -32.6666666666667  -32.50  -32.3333333333333
>   -32.1666666666667  -32.00  -31.833333
> 3333333   -31.6666666666667  -31.50  -31.3333333333333
>   -31.1666666666667  -31.00  -30.8333333333333   -3
> 0.6666666666667  -30.50  -30.3333333333333  -30.1666666666667  -30.00
>   -29.8333333333333   -29.6666666666
> 667  -29.50  -29.3333333333333  -29.1666666666667  -29.00
>   -28.8333333333333   -28.6666666666667  -28.50
>   -28.3333333333333  -28.1666666666667  -28.00  -27.8333333333333
> -27.6666666666667  -27.50  -27.3333333
> 333333  -27.1666666666667  -27.00  -26.8333333333333   -26.6666666666667
>   -26.50  -26.3333333333333  -26.
> 1666666666667  -26.00  -25.8333333333333   -25.6666666666667  -25.50
>   -25.3333333333333  -25.166666666666
> 7  -25.00  -24.8333333333333   -24.6666666666667  -24.50
>   -24.3333333333333  -24.1666666666667  -24.00  -
> 23.8333333333333   -23.6666666666667  -23.50  -23.3333333333333
>   -23.1666666666667  -23.00  -22.833333333
> 3333   -22.6666666666667  -22.50  -22.3333333333333  -22.1666666666667
>   -22.00
>
> YDEF 128 LEVELS -42.8  -42.6777116892495  -42.5551819648098
>   -42.4324109099214  -42.3093986124075  -42.18
> 61451646983  -42.0626506638555  -41.9389152115963  -41.8149389143181
>   -41.6907218831225  -41.566264233839
> 1  -41.4415660870501  -41.3166275681137  -41.1914488071884
>   -41.0660299392564  -40.9403711041476  -40.814
> 4724465629  -40.688334116098  -40.5619562672665  -40.4353390595236
>   -40.3084826572887  -40.1813872299691
>   -40.0540529519827  -39.9264800027809  -39.7986685668712
>   -39.6706188338402  -39.5423309983758  -39.41380
> 52602898  -39.28504182454  -39.1560409012524  -39.0268027057435
>   -38.8973274585417  -38.7676153854092  -3
> 8.6376667173637  -38.5074816907001  -38.377060547011  -38.2464035332088
>   -38.1155109015461  -37.984382909
> 6366  -37.8530198204761  -37.7214219024627  -37.5895894294174
>   -37.4575226806042  -37.3252219407498  -37.
> 192687500064
> -37.0599196542591  -36.9269187045694  -36.7936849577703
>   -36.6602187261978  -36.526520327767  -36.3925900
> 859914  -36.2584283300007  -36.1240353945596  -35.9894116200861
>   -35.8545573526688  -35.7194729440851  -3
> 5.5841587518186  -35.4486151390761  -35.3128424748053  -35.1768411337108
>   -35.0406114962714  -34.90415394
> 8756  -34.7674688832401  -34.6305566976216  -34.4934177956363
>   -34.3560525868735  -34.2184614867912  -34.
> 0806449167309  -33.9426033039323  -33.8043370815481  -33.6658466886573
>   -33.5271325702802  -33.3881951773
> 911  -33.2490349669321  -33.1096524018264 -32.9700479509905
>   -32.8302220893474  -32.6901752978385  -32.54
> 99080634354  -32.4094208791522  -32.2687142440563  -32.1277886632798
>   -31.9866446480301  -31.845282715600
> 8  -31.7037033893814  -31.5619071988675  -31.4198946796703
>   -31.2776663735259  -31.1352228283041  -30.992
> 5645980173  -30.8496922428287  -30.70660632906  -30.5633074291995
>   -30.4197961219092  -30.2760729920316
> -30.1321386305969  -29.9879936348286  -29.8436386081501
>   -29.6990741601898  -29.5543009067868  -29.409319
> 4699954  -29.2641304780902  -29.1187345655695  -28.97313237316
>   -28.8273245478197  -28.6813117427412  -28
> .5350946173545  -28.3886738373291  -28.2420500745764  -28.0952240072509
>   -27.9481963197516  -27.800967702
> 7231  -27.6535388530555  -27.5059104738849  -27.3580832745931
>   -27.2100579708065  -27.0618352843955  -26.
> 9134159434725  -26.7648006823905  -26.6159902417403  -26.466985368348
>   -26.317786815272  -26.168395341799
> 4  -26.0188117134416  -25.8690367019309  -25.7190710852147
>   -25.5689156474509  -25.4185711790023
>
> ZDEF  30 LEVELS  -0.9875  -0.9625  -0.9375  -0.9125  -0.8875  -0.8625
>   -0.8375  -0.8125  -0.7875   -0.762
> 5  -0.7375  -0.7125  -0.6875  -0.6625  -0.6375  -0.6125  -0.5875
>   -0.5625  -0.5375  -0.5125  -0.4875  -0.
> 4625  -0.4375  -0.4125  -0.3875  -0.3625  -0.3375  -0.3125  -0.2875
>   -0.2625  -0.2375  -0.2125  -0.1875
> -0.1625  -0.1375  -0.1125  -0.0875  -0.0625  -0.0375  -0.0125
>
> TDEF   60 LINEAR 06Z13dec1979 12hr
> VARS 4
> temp=>temp 30  t z y x  Potential Temperature (C)
> salt=>salt 30  t z y x  Salinity
> u=>u       30  t z y x  time-averaged u-momentum component (m/s)
> v=>v       30  t z y x  time-averaged v-momentum component  (m/s)
> ENDVARS
>
>
> Thanks,
>
> Gabriel
>
> --
> \/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\
>   Gabriel Bonow Münchow
>   LMQA/CEPSRM - UFRGS
>   Av. Bento Gonçalves, 9500 - Campus do Vale
>   91501-970 - Porto Alegre - RS - Brasil
>   Fone: (051) 3308-6966/7478      Fax: (051) 3308-7477
>   E-mails: gabriel.munchow at ufrgs.br <mailto:gabriel.munchow at ufrgs.br> /
> gabrielbmw at gmail.com <mailto:gabrielbmw at gmail.com>
> \/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\/\
>
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>

-- 

Please note that Charles.Seman at noaa.gov should be considered my NOAA
email address, not cjs at gfdl.noaa.gov.

********************************************************************
  Charles Seman                                Charles.Seman at noaa.gov
  U.S. Department of Commerce / NOAA / OAR
  Geophysical Fluid Dynamics Laboratory         voice: (609) 452-6547
  201 Forrestal Road                              fax: (609) 987-5063
  Princeton, NJ  08540-6649            http://www.gfdl.noaa.gov/~cjs/
********************************************************************

"The contents of this message are mine personally and do not reflect any
official or unofficial position of the United States Federal Government,
the United States Department of Commerce, or NOAA."



More information about the gradsusr mailing list