[gradsusr] Problem with RCM Projections
Roberto Mera
rjmera at coas.oregonstate.edu
Fri Oct 22 18:35:58 EDT 2010
Hi,
Yes, I'm having a similar problem so I'd be interested if anyone out
there can help out with this dataset. Adam, can you attach a ctl file so
I can look at it?
Here is the ncdump -c of the file:
netcdf spdmax_RCM3_cgcm3_1968010103 {
dimensions:
yc = 104 ;
xc = 134 ;
time = UNLIMITED ; // (1095 currently)
variables:
double lat(yc, xc) ;
lat:units = "degrees_north" ;
lat:long_name = "latitude" ;
lat:standard_name = "latitude" ;
lat:axis = "Y" ;
double lon(yc, xc) ;
lon:units = "degrees_east" ;
lon:long_name = "longitude" ;
lon:standard_name = "longitude" ;
lon:axis = "X" ;
float spdmax(time, yc, xc) ;
spdmax:grid_mapping = "Transverse_Mercator" ;
spdmax:standard_name = "wind_speed_of_gust" ;
spdmax:long_name = "Maximum Daily 10-Meter Wind Speed" ;
spdmax:units = "m s-1" ;
spdmax:coordinates = "lon lat level" ;
spdmax:cell_methods = "time: maximum (interval: 1 day)" ;
spdmax:missing_value = 1.e+20f ;
spdmax:original_name = "W10MX" ;
spdmax:original_units = "m s-1" ;
spdmax:_FillValue = 1.e+20f ;
double time(time) ;
time:standard_name = "time" ;
time:long_name = "time" ;
time:axis = "T" ;
time:units = "days since 1968-01-01 00:00:0.0" ;
time:calendar = "noleap" ;
double xc(xc) ;
xc:long_name = "x-coordinate in Cartesian system" ;
xc:standard_name = "projection_x_coordinate" ;
xc:units = "m" ;
xc:axis = "X" ;
double yc(yc) ;
yc:long_name = "y-coordinate in Cartesian system" ;
yc:standard_name = "projection_y_coordinate" ;
yc:units = "m" ;
yc:axis = "Y" ;
char Transverse_Mercator ;
Transverse_Mercator:grid_mapping_name = "transverse_mercator" ;
Transverse_Mercator:longitude_of_central_meridian = -97. ;
Transverse_Mercator:latitude_of_projection_origin = 47.5 ;
Transverse_Mercator:scale_factor_at_central_meridian = 1. ;
Transverse_Mercator:false_easting = 3925000. ;
Transverse_Mercator:false_northing = 3175000. ;
double level ;
level:long_name = "height" ;
level:units = "m" ;
level:postive = "up" ;
level:axis = "Z" ;
// global attributes:
:title = "UC Santa Cruz RegCM3 model output prepared for
NARCCAP historical climate using CGCM" ;
:Conventions = "CF-1.0" ;
:experiment_id = "historical climate using CGCM" ;
:realization = "1" ;
:table_id = "Table 1" ;
:project_id = "NARCCAP" ;
:source = "RegCM3 (2006) atmosphere: RegCM3v3.1 hydrostatic,
split-explicit, horizontal equal area grid resolution of 50 km by 50 km,
vertical sigma-coordinates levels ={0.025, 0.075, 0.13, 0.195, 0.27,
0.35, 0.43, 0.51, 0.59, 0.67, 0.745, 0.8099999, 0.865, 0.91, 0.945,
0.97, 0.985, 0.995} land: BATSv1e Biosphere-Atmosphere Transfer Scheme" ;
:institution = "CCIL, UC Santa Cruz (Climate Change and Impacts
Laboratory, University of California, Santa Cruz, CA, USA)" ;
:NCO = "3.9.9" ;
data:
time = 0.375, 1.375, 2.375, 3.375, 4.375, 5.375, 6.375, 7.375, 8.375,
9.375,
10.375, 11.375, 12.375, 13.375, 14.375, 15.375, 16.375, 17.375,
18.375,
19.375, 20.375, 21.375, 22.375, 23.375, 24.375, 25.375, 26.375,
27.375,
28.375, 29.375, 30.375, 31.375, 32.375, 33.375, 34.375, 35.375,
36.375,
37.375, 38.375, 39.375, 40.375, 41.375, 42.375, 43.375, 44.375,
45.375,
46.375, 47.375, 48.375, 49.375, 50.375, 51.375, 52.375, 53.375,
54.375,
55.375, 56.375, 57.375, 58.375, 59.375, 60.375, 61.375, 62.375,
63.375,
64.375, 65.375, 66.375, 67.375, 68.375, 69.375, 70.375, 71.375,
72.375,
73.375, 74.375, 75.375, 76.375, 77.375, 78.375, 79.375, 80.375,
81.375,
82.375, 83.375, 84.375, 85.375, 86.375, 87.375, 88.375, 89.375,
90.375,
91.375, 92.375, 93.375, 94.375, 95.375, 96.375, 97.375, 98.375,
99.375,
100.375, 101.375, 102.375, 103.375, 104.375, 105.375, 106.375,
107.375,
108.375, 109.375, 110.375, 111.375, 112.375, 113.375, 114.375,
115.375,
116.375, 117.375, 118.375, 119.375, 120.375, 121.375, 122.375,
123.375,
124.375, 125.375, 126.375, 127.375, 128.375, 129.375, 130.375,
131.375,
132.375, 133.375, 134.375, 135.375, 136.375, 137.375, 138.375,
139.375,
140.375, 141.375, 142.375, 143.375, 144.375, 145.375, 146.375,
147.375,
148.375, 149.375, 150.375, 151.375, 152.375, 153.375, 154.375,
155.375,
156.375, 157.375, 158.375, 159.375, 160.375, 161.375, 162.375,
163.375,
164.375, 165.375, 166.375, 167.375, 168.375, 169.375, 170.375,
171.375,
172.375, 173.375, 174.375, 175.375, 176.375, 177.375, 178.375,
179.375,
180.375, 181.375, 182.375, 183.375, 184.375, 185.375, 186.375,
187.375,
188.375, 189.375, 190.375, 191.375, 192.375, 193.375, 194.375,
195.375,
196.375, 197.375, 198.375, 199.375, 200.375, 201.375, 202.375,
203.375,
204.375, 205.375, 206.375, 207.375, 208.375, 209.375, 210.375,
211.375,
212.375, 213.375, 214.375, 215.375, 216.375, 217.375, 218.375,
219.375,
220.375, 221.375, 222.375, 223.375, 224.375, 225.375, 226.375,
227.375,
228.375, 229.375, 230.375, 231.375, 232.375, 233.375, 234.375,
235.375,
236.375, 237.375, 238.375, 239.375, 240.375, 241.375, 242.375,
243.375,
244.375, 245.375, 246.375, 247.375, 248.375, 249.375, 250.375,
251.375,
252.375, 253.375, 254.375, 255.375, 256.375, 257.375, 258.375,
259.375,
260.375, 261.375, 262.375, 263.375, 264.375, 265.375, 266.375,
267.375,
268.375, 269.375, 270.375, 271.375, 272.375, 273.375, 274.375,
275.375,
276.375, 277.375, 278.375, 279.375, 280.375, 281.375, 282.375,
283.375,
284.375, 285.375, 286.375, 287.375, 288.375, 289.375, 290.375,
291.375,
292.375, 293.375, 294.375, 295.375, 296.375, 297.375, 298.375,
299.375,
300.375, 301.375, 302.375, 303.375, 304.375, 305.375, 306.375,
307.375,
308.375, 309.375, 310.375, 311.375, 312.375, 313.375, 314.375,
315.375,
316.375, 317.375, 318.375, 319.375, 320.375, 321.375, 322.375,
323.375,
324.375, 325.375, 326.375, 327.375, 328.375, 329.375, 330.375,
331.375,
332.375, 333.375, 334.375, 335.375, 336.375, 337.375, 338.375,
339.375,
340.375, 341.375, 342.375, 343.375, 344.375, 345.375, 346.375,
347.375,
348.375, 349.375, 350.375, 351.375, 352.375, 353.375, 354.375,
355.375,
356.375, 357.375, 358.375, 359.375, 360.375, 361.375, 362.375,
363.375,
364.375, 365.375, 366.375, 367.375, 368.375, 369.375, 370.375,
371.375,
372.375, 373.375, 374.375, 375.375, 376.375, 377.375, 378.375,
379.375,
380.375, 381.375, 382.375, 383.375, 384.375, 385.375, 386.375,
387.375,
388.375, 389.375, 390.375, 391.375, 392.375, 393.375, 394.375,
395.375,
396.375, 397.375, 398.375, 399.375, 400.375, 401.375, 402.375,
403.375,
404.375, 405.375, 406.375, 407.375, 408.375, 409.375, 410.375,
411.375,
412.375, 413.375, 414.375, 415.375, 416.375, 417.375, 418.375,
419.375,
420.375, 421.375, 422.375, 423.375, 424.375, 425.375, 426.375,
427.375,
428.375, 429.375, 430.375, 431.375, 432.375, 433.375, 434.375,
435.375,
436.375, 437.375, 438.375, 439.375, 440.375, 441.375, 442.375,
443.375,
444.375, 445.375, 446.375, 447.375, 448.375, 449.375, 450.375,
451.375,
452.375, 453.375, 454.375, 455.375, 456.375, 457.375, 458.375,
459.375,
460.375, 461.375, 462.375, 463.375, 464.375, 465.375, 466.375,
467.375,
468.375, 469.375, 470.375, 471.375, 472.375, 473.375, 474.375,
475.375,
476.375, 477.375, 478.375, 479.375, 480.375, 481.375, 482.375,
483.375,
484.375, 485.375, 486.375, 487.375, 488.375, 489.375, 490.375,
491.375,
492.375, 493.375, 494.375, 495.375, 496.375, 497.375, 498.375,
499.375,
500.375, 501.375, 502.375, 503.375, 504.375, 505.375, 506.375,
507.375,
508.375, 509.375, 510.375, 511.375, 512.375, 513.375, 514.375,
515.375,
516.375, 517.375, 518.375, 519.375, 520.375, 521.375, 522.375,
523.375,
524.375, 525.375, 526.375, 527.375, 528.375, 529.375, 530.375,
531.375,
532.375, 533.375, 534.375, 535.375, 536.375, 537.375, 538.375,
539.375,
540.375, 541.375, 542.375, 543.375, 544.375, 545.375, 546.375,
547.375,
548.375, 549.375, 550.375, 551.375, 552.375, 553.375, 554.375,
555.375,
556.375, 557.375, 558.375, 559.375, 560.375, 561.375, 562.375,
563.375,
564.375, 565.375, 566.375, 567.375, 568.375, 569.375, 570.375,
571.375,
572.375, 573.375, 574.375, 575.375, 576.375, 577.375, 578.375,
579.375,
580.375, 581.375, 582.375, 583.375, 584.375, 585.375, 586.375,
587.375,
588.375, 589.375, 590.375, 591.375, 592.375, 593.375, 594.375,
595.375,
596.375, 597.375, 598.375, 599.375, 600.375, 601.375, 602.375,
603.375,
604.375, 605.375, 606.375, 607.375, 608.375, 609.375, 610.375,
611.375,
612.375, 613.375, 614.375, 615.375, 616.375, 617.375, 618.375,
619.375,
620.375, 621.375, 622.375, 623.375, 624.375, 625.375, 626.375,
627.375,
628.375, 629.375, 630.375, 631.375, 632.375, 633.375, 634.375,
635.375,
636.375, 637.375, 638.375, 639.375, 640.375, 641.375, 642.375,
643.375,
644.375, 645.375, 646.375, 647.375, 648.375, 649.375, 650.375,
651.375,
652.375, 653.375, 654.375, 655.375, 656.375, 657.375, 658.375,
659.375,
660.375, 661.375, 662.375, 663.375, 664.375, 665.375, 666.375,
667.375,
668.375, 669.375, 670.375, 671.375, 672.375, 673.375, 674.375,
675.375,
676.375, 677.375, 678.375, 679.375, 680.375, 681.375, 682.375,
683.375,
684.375, 685.375, 686.375, 687.375, 688.375, 689.375, 690.375,
691.375,
692.375, 693.375, 694.375, 695.375, 696.375, 697.375, 698.375,
699.375,
700.375, 701.375, 702.375, 703.375, 704.375, 705.375, 706.375,
707.375,
708.375, 709.375, 710.375, 711.375, 712.375, 713.375, 714.375,
715.375,
716.375, 717.375, 718.375, 719.375, 720.375, 721.375, 722.375,
723.375,
724.375, 725.375, 726.375, 727.375, 728.375, 729.375, 730.375,
731.375,
732.375, 733.375, 734.375, 735.375, 736.375, 737.375, 738.375,
739.375,
740.375, 741.375, 742.375, 743.375, 744.375, 745.375, 746.375,
747.375,
748.375, 749.375, 750.375, 751.375, 752.375, 753.375, 754.375,
755.375,
756.375, 757.375, 758.375, 759.375, 760.375, 761.375, 762.375,
763.375,
764.375, 765.375, 766.375, 767.375, 768.375, 769.375, 770.375,
771.375,
772.375, 773.375, 774.375, 775.375, 776.375, 777.375, 778.375,
779.375,
780.375, 781.375, 782.375, 783.375, 784.375, 785.375, 786.375,
787.375,
788.375, 789.375, 790.375, 791.375, 792.375, 793.375, 794.375,
795.375,
796.375, 797.375, 798.375, 799.375, 800.375, 801.375, 802.375,
803.375,
804.375, 805.375, 806.375, 807.375, 808.375, 809.375, 810.375,
811.375,
812.375, 813.375, 814.375, 815.375, 816.375, 817.375, 818.375,
819.375,
820.375, 821.375, 822.375, 823.375, 824.375, 825.375, 826.375,
827.375,
828.375, 829.375, 830.375, 831.375, 832.375, 833.375, 834.375,
835.375,
836.375, 837.375, 838.375, 839.375, 840.375, 841.375, 842.375,
843.375,
844.375, 845.375, 846.375, 847.375, 848.375, 849.375, 850.375,
851.375,
852.375, 853.375, 854.375, 855.375, 856.375, 857.375, 858.375,
859.375,
860.375, 861.375, 862.375, 863.375, 864.375, 865.375, 866.375,
867.375,
868.375, 869.375, 870.375, 871.375, 872.375, 873.375, 874.375,
875.375,
876.375, 877.375, 878.375, 879.375, 880.375, 881.375, 882.375,
883.375,
884.375, 885.375, 886.375, 887.375, 888.375, 889.375, 890.375,
891.375,
892.375, 893.375, 894.375, 895.375, 896.375, 897.375, 898.375,
899.375,
900.375, 901.375, 902.375, 903.375, 904.375, 905.375, 906.375,
907.375,
908.375, 909.375, 910.375, 911.375, 912.375, 913.375, 914.375,
915.375,
916.375, 917.375, 918.375, 919.375, 920.375, 921.375, 922.375,
923.375,
924.375, 925.375, 926.375, 927.375, 928.375, 929.375, 930.375,
931.375,
932.375, 933.375, 934.375, 935.375, 936.375, 937.375, 938.375,
939.375,
940.375, 941.375, 942.375, 943.375, 944.375, 945.375, 946.375,
947.375,
948.375, 949.375, 950.375, 951.375, 952.375, 953.375, 954.375,
955.375,
956.375, 957.375, 958.375, 959.375, 960.375, 961.375, 962.375,
963.375,
964.375, 965.375, 966.375, 967.375, 968.375, 969.375, 970.375,
971.375,
972.375, 973.375, 974.375, 975.375, 976.375, 977.375, 978.375,
979.375,
980.375, 981.375, 982.375, 983.375, 984.375, 985.375, 986.375,
987.375,
988.375, 989.375, 990.375, 991.375, 992.375, 993.375, 994.375,
995.375,
996.375, 997.375, 998.375, 999.375, 1000.375, 1001.375, 1002.375,
1003.375, 1004.375, 1005.375, 1006.375, 1007.375, 1008.375, 1009.375,
1010.375, 1011.375, 1012.375, 1013.375, 1014.375, 1015.375, 1016.375,
1017.375, 1018.375, 1019.375, 1020.375, 1021.375, 1022.375, 1023.375,
1024.375, 1025.375, 1026.375, 1027.375, 1028.375, 1029.375, 1030.375,
1031.375, 1032.375, 1033.375, 1034.375, 1035.375, 1036.375, 1037.375,
1038.375, 1039.375, 1040.375, 1041.375, 1042.375, 1043.375, 1044.375,
1045.375, 1046.375, 1047.375, 1048.375, 1049.375, 1050.375, 1051.375,
1052.375, 1053.375, 1054.375, 1055.375, 1056.375, 1057.375, 1058.375,
1059.375, 1060.375, 1061.375, 1062.375, 1063.375, 1064.375, 1065.375,
1066.375, 1067.375, 1068.375, 1069.375, 1070.375, 1071.375, 1072.375,
1073.375, 1074.375, 1075.375, 1076.375, 1077.375, 1078.375, 1079.375,
1080.375, 1081.375, 1082.375, 1083.375, 1084.375, 1085.375, 1086.375,
1087.375, 1088.375, 1089.375, 1090.375, 1091.375, 1092.375, 1093.375,
1094.375 ;
xc = 600000, 650000, 700000, 750000, 800000, 850000, 900000, 950000,
1000000, 1050000, 1100000, 1150000, 1200000, 1250000, 1300000,
1350000,
1400000, 1450000, 1500000, 1550000, 1600000, 1650000, 1700000,
1750000,
1800000, 1850000, 1900000, 1950000, 2000000, 2050000, 2100000,
2150000,
2200000, 2250000, 2300000, 2350000, 2400000, 2450000, 2500000,
2550000,
2600000, 2650000, 2700000, 2750000, 2800000, 2850000, 2900000,
2950000,
3000000, 3050000, 3100000, 3150000, 3200000, 3250000, 3300000,
3350000,
3400000, 3450000, 3500000, 3550000, 3600000, 3650000, 3700000,
3750000,
3800000, 3850000, 3900000, 3950000, 4000000, 4050000, 4100000,
4150000,
4200000, 4250000, 4300000, 4350000, 4400000, 4450000, 4500000,
4550000,
4600000, 4650000, 4700000, 4750000, 4800000, 4850000, 4900000,
4950000,
5000000, 5050000, 5100000, 5150000, 5200000, 5250000, 5300000,
5350000,
5400000, 5450000, 5500000, 5550000, 5600000, 5650000, 5700000,
5750000,
5800000, 5850000, 5900000, 5950000, 6000000, 6050000, 6100000,
6150000,
6200000, 6250000, 6300000, 6350000, 6400000, 6450000, 6500000,
6550000,
6600000, 6650000, 6700000, 6750000, 6800000, 6850000, 6900000,
6950000,
7000000, 7050000, 7100000, 7150000, 7200000, 7250000 ;
yc = 600000, 650000, 700000, 750000, 800000, 850000, 900000, 950000,
1000000, 1050000, 1100000, 1150000, 1200000, 1250000, 1300000,
1350000,
1400000, 1450000, 1500000, 1550000, 1600000, 1650000, 1700000,
1750000,
1800000, 1850000, 1900000, 1950000, 2000000, 2050000, 2100000,
2150000,
2200000, 2250000, 2300000, 2350000, 2400000, 2450000, 2500000,
2550000,
2600000, 2650000, 2700000, 2750000, 2800000, 2850000, 2900000,
2950000,
3000000, 3050000, 3100000, 3150000, 3200000, 3250000, 3300000,
3350000,
3400000, 3450000, 3500000, 3550000, 3600000, 3650000, 3700000,
3750000,
3800000, 3850000, 3900000, 3950000, 4000000, 4050000, 4100000,
4150000,
4200000, 4250000, 4300000, 4350000, 4400000, 4450000, 4500000,
4550000,
4600000, 4650000, 4700000, 4750000, 4800000, 4850000, 4900000,
4950000,
5000000, 5050000, 5100000, 5150000, 5200000, 5250000, 5300000,
5350000,
5400000, 5450000, 5500000, 5550000, 5600000, 5650000, 5700000,
5750000 ;
}
On 10/22/10 6:48 AM, Adam Cinderich wrote:
>
> Hello,
>
> I am trying to plot model output from the various RCMs of NARCCAP in
> GrADS. Because the default projection in GrADS differs from those of
> the RCMs, when I try to plot temperature (for example), the plot does
> not overlay the map correctly. The RCM output is netcdf format and I
> am using a .ctl file in GrADS. In the .ctl file, I have tried
> altering the pdef line several times but still can't get it to work
> correctly. Can anyone offer any suggestions for what to enter in the
> pdef line for polar stereographic and/or lambert conformal? Or, can
> anybody offer any codes that will take care of this reprojecting?
>
> Thanks,
>
> *Adam B. Cinderich*
>
> Graduate Student
>
> Michigan State University
>
> Department of Geography
>
> 1E Geography Building
>
> East Lansing, MI 48824
>
> cinderic at msu.edu|(517)432-6218
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20101022/df322aa5/attachment-0003.html
More information about the gradsusr
mailing list