<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META content="text/html; charset=gb2312" http-equiv=Content-Type>
<STYLE>
BLOCKQUOTE {
MARGIN-BOTTOM: 0px; MARGIN-LEFT: 2em; MARGIN-TOP: 0px
}
OL {
MARGIN-BOTTOM: 0px; MARGIN-TOP: 0px
}
UL {
MARGIN-BOTTOM: 0px; MARGIN-TOP: 0px
}
P {
MARGIN-BOTTOM: 0px; MARGIN-TOP: 0px
}
BODY {
FONT-SIZE: 10.5pt; FONT-FAMILY: Microsoft YaHei UI; COLOR: #000000; LINE-HEIGHT: 1.5
}
</STYLE>
<META name=GENERATOR content="MSHTML 10.00.9200.17492"></HEAD>
<BODY style="MARGIN: 10px">
<DIV>
<DIV>Dear all,</DIV>
<DIV> </DIV>
<DIV>I used the latlon_to_SCRIP function to
generate lat-lon grid file with SCRIP format which will then be used
as input files of the ESMF application "ESMF_RegridWeightGen".</DIV>
<DIV>When I applied latlon_to_SCRIP to 1x1, 0.5x0.5 and 0.25x0.25 grids,
the output results were right. However, when I applied latlon_to_SCRIP to
0.1x0.1 grids, the output results </DIV>
<DIV>are wrong. Do the latlon_to_SCRIP function support to generate a
0.1x0.1 grid file? Can you help me to solve this problem? </DIV>
<DIV>Thank you very much.</DIV>
<DIV> </DIV>
<DIV style="COLOR: #ff0000">The code used to generate the 0.25x0.25 grids
are as follows:</DIV>
<DIV>_______________</DIV>
<DIV>
<DIV>Opt
= True</DIV>
<DIV><A
href="mailto:Opt@LLCorner">Opt@LLCorner</A>
= (/ -89.875d, -179.875d/)</DIV>
<DIV><A
href="mailto:Opt@URCorner">Opt@URCorner</A>
= (/ 89.875d, 179.875d/)</DIV>
<DIV><A href="mailto:Opt@ForceOverwrite">Opt@ForceOverwrite</A> = True</DIV>
<DIV><A href="mailto:Opt@PrintTimings">Opt@PrintTimings</A> =
True</DIV>
<DIV>latlon_to_SCRIP(srcGridName,"0.25x0.25",Opt)</DIV>
<DIV>delete(Opt)</DIV>
<DIV>________________</DIV>
<DIV> </DIV></DIV>
<DIV style="COLOR: #ff0000">The header of the output netcdf file is as follows:
</DIV>
<DIV>________________</DIV>
<DIV>
<DIV>dimensions:</DIV>
<DIV> grid_size = 1036800 ;</DIV>
<DIV> grid_corners = 4 ;</DIV>
<DIV> grid_rank = 2 ;</DIV>
<DIV>variables:</DIV>
<DIV> int grid_dims(grid_rank) ;</DIV>
<DIV> double
grid_center_lat(grid_size) ;</DIV>
<DIV>
grid_center_lat:units = "degrees" ;</DIV>
<DIV> double
grid_center_lon(grid_size) ;</DIV>
<DIV>
grid_center_lon:units = "degrees" ;</DIV>
<DIV> int grid_imask(grid_size)
;</DIV>
<DIV>
grid_imask:units = "unitless" ;</DIV>
<DIV> double
grid_corner_lat(grid_size, grid_corners) ;</DIV>
<DIV>
grid_corner_lat:units = "degrees" ;</DIV>
<DIV> double
grid_corner_lon(grid_size, grid_corners) ;</DIV>
<DIV>
grid_corner_lon:units = "degrees" ;</DIV>
<DIV> </DIV>
<DIV>// global attributes:</DIV>
<DIV>
:date_created = "Thu Oct 15 19:49:54 CST 2015" ;</DIV>
<DIV>
:Createdby = "ESMF_regridding.ncl" ;</DIV>
<DIV>
:Conventions = "SCRIP" ;</DIV>
<DIV>
:title = "0.25x0.25 grid" ;</DIV>
<DIV>data:</DIV>
<DIV> </DIV>
<DIV>grid_dims = 1440, 720 ;</DIV></DIV>
<DIV>________________</DIV>
<DIV> </DIV>
<DIV> </DIV>
<DIV>
<DIV style="COLOR: #ff0000">The code used to generate the 0.1x0.1 grids
are as follows:</DIV></DIV>
<DIV>
<DIV>_______________</DIV>
<DIV>
<DIV>Opt
= True</DIV>
<DIV><A
href="mailto:Opt@LLCorner">Opt@LLCorner</A>
= (/ -89.95d, -179.95d/)</DIV>
<DIV><A
href="mailto:Opt@URCorner">Opt@URCorner</A>
= (/ 89.95d, 179.95d/)</DIV>
<DIV><A href="mailto:Opt@ForceOverwrite">Opt@ForceOverwrite</A> = True</DIV>
<DIV><A href="mailto:Opt@PrintTimings">Opt@PrintTimings</A> =
True</DIV>
<DIV>latlon_to_SCRIP(srcGridName,"0.1x0.1",Opt)</DIV>
<DIV>delete(Opt)</DIV>
<DIV>________________</DIV>
<DIV> </DIV>
<DIV style="COLOR: #ff0000">The header of the output netcdf file is as follows:
</DIV></DIV></DIV>
<DIV>_____________</DIV>
<DIV>
<DIV>dimensions:</DIV>
<DIV> grid_size = 6474601 ;</DIV>
<DIV> grid_corners = 4 ;</DIV>
<DIV> grid_rank = 2 ;</DIV>
<DIV>variables:</DIV>
<DIV> int grid_dims(grid_rank) ;</DIV>
<DIV> double
grid_center_lat(grid_size) ;</DIV>
<DIV>
grid_center_lat:units = "degrees" ;</DIV>
<DIV> double
grid_center_lon(grid_size) ;</DIV>
<DIV>
grid_center_lon:units = "degrees" ;</DIV>
<DIV> int grid_imask(grid_size)
;</DIV>
<DIV>
grid_imask:units = "unitless" ;</DIV>
<DIV> double
grid_corner_lat(grid_size, grid_corners) ;</DIV>
<DIV>
grid_corner_lat:units = "degrees" ;</DIV>
<DIV> double
grid_corner_lon(grid_size, grid_corners) ;</DIV>
<DIV>
grid_corner_lon:units = "degrees" ;</DIV>
<DIV> </DIV>
<DIV>// global attributes:</DIV>
<DIV>
:date_created = "Thu Oct 15 20:27:42 CST 2015" ;</DIV>
<DIV>
:Createdby = "ESMF_regridding.ncl" ;</DIV>
<DIV>
:Conventions = "SCRIP" ;</DIV>
<DIV>
:title = "0.1deg grid" ;</DIV>
<DIV>data:</DIV>
<DIV> </DIV>
<DIV>grid_dims = 3599, 1799 ;</DIV></DIV>
<DIV>_____________</DIV>
<DIV style="FONT-WEIGHT: bold; COLOR: #ff0000">Why the grid_dims are
3599x1799? I think it should be 3600x1800.</DIV>
<DIV> </DIV>
<DIV> </DIV>
<DIV>Yours Sincerely,</DIV>
<DIV>Huansheng Chen</DIV></DIV>
<DIV> </DIV>
<HR style="HEIGHT: 1px; WIDTH: 210px" align=left color=#b5c4df SIZE=1>
<DIV><SPAN>
<DIV style="FONT-SIZE: 10pt; FONT-FAMILY: verdana; MARGIN: 10px">
<DIV>
<DIV style="FONT-SIZE: 10pt; FONT-FAMILY: verdana; MARGIN: 10px">
<DIV>陈焕盛（LAPC/IAP）</DIV>
<DIV>大气边界层物理和大气化学国家重点实验室</DIV>
<DIV>中国科学院大气物理研究所</DIV></DIV></DIV></DIV></SPAN></DIV></BODY></HTML>