Dear all,<div> I'd like to use CCSM POP SST output for driving the wrf model. When I wrote the fortran program to read the SST data, some details of the data were poster below.I found there are some parameters I was not make sure, I hope I can get some help.</div>
<div>CCSM POP SST file (<a href="http://b30.030e.pop.h.SST.1870-01_cat_1999-12.nc">b30.030e.pop.h.SST.1870-01_cat_1999-12.nc</a>)</div><div>My little fortran program like that:</div><div><div><span class="Apple-style-span" style="text-decoration: underline; font-family: Times; font-size: 14px; font-style: italic; -webkit-text-decorations-in-effect: underline; "><br>
</span></div><div><p style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 12px/normal Times; "><br></p></div></div><div><div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; ">integer :: iproj ! Code for projection of data in array:</span></span></div>
<div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; "> ! 0 = cylindrical equidistant</span></span></div>
<div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; "> ! 1 = Mercator</span></span></div>
<div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; "> ! 3 = Lambert conformal conic</span></span></div>
<div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; "> ! 4 = Gaussian</span></span></div>
<div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; "> ! 5 = Polar stereographic</span></span></div>
<div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; ">real :: nlats ! Number of latitudes north of equator</span></span></div>
<div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; "> ! (whole number, for Gaussian grids)</span></span></div>
<div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; ">real :: xfcst ! Forecast hour of data</span></span></div>
<div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; ">real :: xlvl ! Vertical level of data in 2-d array</span></span></div>
<div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; ">real :: startlat, startlon ! Lat/lon of point in array indicated by</span></span></div>
<div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; "> ! startloc string</span></span></div>
<div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; ">real :: dx, dy ! Grid spacing, km</span></span></div>
<div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; ">real :: xlonc ! Standard longitude of projection</span></span></div>
<div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; ">real :: earth_radius ! Earth radius, km<br>
</span></span></div><div><div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; ">real :: deltalat, deltalon ! Grid spacing, degrees</span></span></div>
<div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; "><br></span></span></div><div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; ">....</span></span></div>
<div><div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; "> write (unit= iunit) version</span></span></div>
<div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; "> write (unit= iunit) hdate_output,xfcst,map_source, field,&</span></span></div>
<div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; "> units, desc,xlvl, nx, ny, iproj</span></span></div>
<div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; "> write (unit=iunit) startloc,startlat,startlon,nlats,deltalon,earth_radius</span></span></div>
<div><span class="Apple-style-span" style="font-style: italic; "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; "> close (iunit)</span></span><br></div><div>
<br></div></div></div><div><span class="Apple-style-span" style="background-color: rgb(255, 153, 102);">I'd like to know how can I get the referenced parameters like deltalat,deltalon since the grid projection is not so regular; and how should I choose the startlat, startlon, that mean the first tlat, tlong in my ducument? what is the earth_radius? Thanks a lot.</span></div>
<div><br></div><div><br></div></div><div>Information about the data</div><div><div> nlat = 384 ;</div><div> nlon = 320 ;</div><div> time = UNLIMITED ; // (1560 currently)</div><div> z_t = 1 ;</div>
<div> z_w = 40 ;</div><div><div> float SST(time, z_t, nlat, nlon) ;</div><div> SST:long_name = "Potential Temperature" ;</div><div> SST:units = "degC" ;</div>
<div> SST:coordinates = "TLONG TLAT z_t time" ;</div><div> SST:cell_methods = "time: mean" ;</div><div> SST:_FillValue = 9.96921e+36f ;</div><div> SST:missing_value = 9.96921e+36f ;</div>
<div> float TLAT(nlat, nlon) ;</div><div> TLAT:long_name = "array of t-grid latitudes" ;</div><div> TLAT:units = "degrees_north" ;</div><div> float TLONG(nlat, nlon) ;</div>
<div> TLONG:long_name = "array of t-grid longitudes" ;</div><div> TLONG:units = "degrees_east" </div><div><div> float z_t(z_t) ;</div><div> z_t:long_name = "depth from surface to midpoint of layer" ;</div>
<div> z_t:units = "centimeters" ;</div><div> z_t:positive = "down" ;</div><div> z_t:valid_min = 500.622f ;</div><div> z_t:valid_max = 537499.9f ;</div>
<div><br></div></div><div>in South hemisphere, I have TLAT like this:</div><div><div>TLAT =</div><div> -79.22052, -79.22052, -79.22052, -79.22052, -79.22052, -79.22052...</div><div> -78.68631, -78.68631,...</div><div> -78.15209, -78.15209,...</div>
<div> -77.61787, -77.61787,...</div><div>...</div><div> -0.9349651, -0.9349651, -0.9349651, -0.9349651, -0.9349651, -0.9349651,...</div><div>-0.1335664, -0.1335664, -0.1335664, -0.1335664, -0.1335664, -0.1335664,...<br>
</div><div> while in north hemisphere, the TLAT like this:</div><div>0.1334779, 0.1334779, 0.133478, 0.1334781, 0.1334782, 0.1334784, 0.1334786,...<br></div><div>0.4003012, 0.4003014, 0.4003017, 0.4003021, 0.4003028, 0.4003035,...</div>
<div>0.4003044, 0.4003055, 0.4003067, 0.4003081, 0.4003096, 0.4003113,...</div><div>...</div><div><div> The last line like this :</div><div>72.18597, 72.18933, 72.19604, 72.20609, 72.21945, 72.23611, 72.25603, 72.27917, 72.30551, 72.33498, 72.36753, 72.40311, 72.44166, 72.48309,</div>
<div> 72.52736, 72.57437, 72.62406, 72.67633, 72.73111, 72.78831, 72.84784,</div><div> 72.90961, 72.97351, 73.03947, 73.10739, 73.17717, 73.24872, 73.32194,</div><div> 73.39672, 73.47299, 73.55065, 73.6296, 73.70975, 73.791, 73.87328,</div>
<div> 73.95649, 74.04054, 74.12536, 74.21085, 74.29695, 74.38358, 74.47064,</div><div> 74.55808, 74.64583, 74.7338, 74.82195, 74.91021, 74.9985, 75.08678,</div><div> 75.175, 75.26308, 75.351, 75.43868, 75.52609, 75.61319, 75.69992,</div>
<div> 75.78625, 75.87215, 75.95757, 76.04248, 76.12685, 76.21063, 76.29382,</div><div> 76.37639, 76.45829, 76.53952, 76.62004, 76.69984, 76.77891, 76.85721,</div><div> 76.93474, 77.01147, 77.08739, 77.16249, 77.23676, 77.3102, 77.38277,</div>
<div> 77.45448, 77.52531, 77.59527, 77.66434, 77.73251, 77.79979, 77.86617,</div><div> 77.93163, 77.99619, 78.05982, 78.12255, 78.18435, 78.24523, 78.3052,</div><div> 78.36423, 78.42235, 78.47954, 78.53581, 78.59116, 78.64558, 78.69908,</div>
<div> 78.75166, 78.80332, 78.85407, 78.9039, 78.95282, 79.00082, 79.04792,</div><div> 79.09412, 79.1394, 79.18378, 79.22727, 79.26986, 79.31155, 79.35236,</div><div> 79.39228, 79.43131, 79.46946, 79.50673, 79.54312, 79.57864, 79.61328,</div>
<div> 79.64706, 79.67998, 79.71203, 79.74322, 79.77355, 79.80303, 79.83166,</div><div> 79.85944, 79.88637, 79.91245, 79.9377, 79.9621, 79.98568, 80.00842,</div><div> 80.03032, 80.05139, 80.07163, 80.09105, 80.10965, 80.12742, 80.14437,</div>
<div> 80.16051, 80.17582, 80.19032, 80.20401, 80.21688, 80.22895, 80.2402,</div><div> 80.25065, 80.26028, 80.26911, 80.27714, 80.28436, 80.29077, 80.29639,</div><div> 80.30119, 80.3052, 80.3084, 80.31081, 80.31241, 80.31321, 80.31321,</div>
<div> 80.31241, 80.31081, 80.3084, 80.3052, 80.30119, 80.29639, 80.29077,</div><div> 80.28436, 80.27714, 80.26911, 80.26028, 80.25065, 80.2402, 80.22895,</div><div> 80.21688, 80.20401, 80.19032, 80.17582, 80.16051, 80.14437, 80.12742,</div>
<div> 80.10965, 80.09105, 80.07163, 80.05139, 80.03032, 80.00842, 79.98568,</div><div> 79.9621, 79.9377, 79.91245, 79.88637, 79.85944, 79.83166, 79.80303,</div><div> 79.77355, 79.74322, 79.71203, 79.67998, 79.64706, 79.61328, 79.57864,</div>
<div> 79.54312, 79.50673, 79.46946, 79.43131, 79.39228, 79.35236, 79.31155,</div><div> 79.26986, 79.22727, 79.18378, 79.1394, 79.09412, 79.04792, 79.00082,</div><div> 78.95282, 78.9039, 78.85407, 78.80332, 78.75166, 78.69908, 78.64558,</div>
<div> 78.59116, 78.53581, 78.47954, 78.42235, 78.36423, 78.3052, 78.24523,</div><div> 78.18435, 78.12255, 78.05982, 77.99619, 77.93163, 77.86617, 77.79979,</div><div> 77.73251, 77.66434, 77.59527, 77.52531, 77.45448, 77.38277, 77.3102,</div>
<div><div> 77.23676, 77.16249, 77.08739, 77.01147, 76.93474, 76.85721, 76.77891,</div><div> 76.69984, 76.62004, 76.53952, 76.45829, 76.37639, 76.29382, 76.21063,</div><div> 76.12685, 76.04248, 75.95757, 75.87215, 75.78625, 75.69992, 75.61319,</div>
<div> 75.52609, 75.43868, 75.351, 75.26308, 75.175, 75.08678, 74.9985,</div><div> 74.91021, 74.82195, 74.7338, 74.64583, 74.55808, 74.47064, 74.38358,</div><div> 74.29695, 74.21085, 74.12536, 74.04054, 73.95649, 73.87328, 73.791,</div>
<div> 73.70975, 73.6296, 73.55065, 73.47299, 73.39672, 73.32194, 73.24872,</div><div> 73.17717, 73.10739, 73.03947, 72.97351, 72.90961, 72.84784, 72.78831,</div><div> 72.73111, 72.67633, 72.62406, 72.57437, 72.52736, 72.48309, 72.44166,</div>
<div> 72.40311, 72.36753, 72.33498, 72.30551, 72.27917, 72.25603, 72.23611,</div><div> 72.21945, 72.20609, 72.19604, 72.18933, 72.18597 ;</div><div><br></div><div>TLONG =</div><div> 320.5625, 321.6875, 322.8125, 323.9375, 325.0625, 326.1875, 327.3125,</div>
<div> 328.4375, 329.5625, 330.6875, 331.8125, 332.9375, 334.0625, 335.1875,</div><div> 336.3125, 337.4375, 338.5625, 339.6875, 340.8125, 341.9375, 343.0625,</div><div> 344.1875, 345.3125, 346.4375, 347.5625, 348.6875, 349.8125, 350.9375,</div>
</div></div></div></div></div>