爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 8110|回复: 10

[脚本编辑] 不同版本GrADS结果不一样,数据和底图不吻合?

[复制链接]

新浪微博达人勋

发表于 2020-6-12 19:05:50 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
刚接触NC格式数据和GrADS,遇到一些问题,请各位指教。
BCC_AVIM模式转出来的结果是nc格式的数据,用较新版本的GrADS sdfopen命令打开都显示
SDF Error: 365 day calendars are no longer supported by sdfopen.
To open this file with GrADS, use a descriptor file with a complete TDEF entry and OPTIONS 365_day_calendar.
去官网看了说要根据ncdump写一个数据描述文件。ncdump的结果在我另一个求问帖里有。尝试了用GrADS1.9能直接sdfopen打开nc数据,并且display其中一个变量分配系数dis1能出来合理的图(图1),在1.9里q ctlinfo,然后根据结果加了OPTIONS 365_day_calendar写了一个ctl文件(图2),在气象家园整合版的opengrads里open这个ctl,能打开nc数据了,但同样 d dis1,出现的图却不合理(图3),图上东西经两部分的数据和地图为什么会不吻合?而且会显示Contouring: -3e+38 to 3e+38 ,而在Grads1.9里就不会,分配系数dis1的范围在0-0.9是合理的。为什么会这样呢?是我的ctl文件写的有问题吗?有什么解决方法吗,恳请各位指教。


图1grads1.9合理结果

图1grads1.9合理结果

图2编写的ctl文件

图2编写的ctl文件

图3新版本grads不合理结果

图3新版本grads不合理结果
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2020-6-12 22:04:32 | 显示全部楼层
这个十有八九是缺测值设置有问题。你使用 q attr 命令看下missing value 是多少
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2020-6-15 19:22:47 | 显示全部楼层
river 发表于 2020-6-12 22:04
这个十有八九是缺测值设置有问题。你使用 q attr 命令看下missing value 是多少

用整合版的opengrad,open先前写的那个ctl文件后,q attr 显示
No Descriptor Attributes for File 1 : old1990-06
No Native Attributes for File 1 : old1990-06
在1.9下q attr的结果是图1这样,缺测值是1e+36的意思吧,是要把1.9中q ctlinfo 得到的ctl文件undef 9.96921e+36后面的9.96921e+36改成1e+36吗,我这样试过了,把ctl文件中undef那行数值修改成1e+36后,d dis1 依旧是原来错位的样子,求指教

图1

图1
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2020-6-15 21:42:54 | 显示全部楼层
本帖最后由 老A 于 2020-6-15 21:44 编辑

再看看,看提示感觉是缺测值的问题
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2020-6-15 21:45:13 | 显示全部楼层
用2.2打开看看,用q ctlinfo看看缺测值是多少
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2020-6-15 23:12:40 | 显示全部楼层
1、试一下加上 dtype netcdf
2、错位问题 我觉得旋转 180°刚刚好对上,所以资料排列不是从0°开始排列的,所以试一下 xdef 800 linear -180 0.45
3、从出图数据的变化范围来看,缺测值应该也是有问题的
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2020-6-16 17:26:49 | 显示全部楼层
river 发表于 2020-6-15 23:12
1、试一下加上 dtype netcdf
2、错位问题 我觉得旋转 180°刚刚好对上,所以资料排列不是从0°开始排列的 ...

感谢回复!我按照您说的进行了一些尝试
使用1.9q ctlinfo 结果undef 9.96921e+36ctl中添上一句dtype netcdf 会显示Open Error:  Invalid variable record(对应ctl文件、脚本、出图如图1
删掉dtype netcdf,能打开,就是原来那种错位的情况。
使用1.9q ctlinfo 结果undef 9.96921e+36,经度改成xdef 800 linear -180 0.45后,位置有改善,但还是一点偏移,分配系数依旧Contouring: -3e+38 to 3e+38 interval 5e+37。(对应ctl文件、脚本、出图如图2


undef改成undef 1e+36,经度用xdef 800 linear -180 0.45,图和上一种情况一样,有偏移,分配系数范围不合理(对应ctl文件、脚本、出图如图3)。
恳请看到帖的各位多加指教


图1

图1

图2

图2

图3

图3
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2020-6-16 17:36:46 | 显示全部楼层
老A 发表于 2020-6-15 21:45
用2.2打开看看,用q ctlinfo看看缺测值是多少

那个模式转出来的结果nc格式数据,用2.0以后的版本都没有办法sdfopen,显示SDF Error: 365 day calendars are no longer supported by sdfopen.To open this file with GrADS, use a descriptor file with a complete TDEF entry and OPTIONS 365_day_calendar.为了在新版本里打开它才去写的ctl文件,数据在1.9里能打开,目前写的ctl文件也是根据1.9里q ctlinfo的结果写的,在较新版本GrADS里还没打开数据,貌似没办法q ctlinfo
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2020-6-16 21:26:16 | 显示全部楼层
用ncdump.exe取nc文件的信息发出来看看
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2020-6-19 15:17:04 | 显示全部楼层
老A 发表于 2020-6-16 21:26
用ncdump.exe取nc文件的信息发出来看看

/cygdrive/E/GrADSv2.2.1/sample$>ncdump -c old1990-08.nc
netcdf old1990-08 {
dimensions:
        gridcell = 116421 ;
        landunit = 201491 ;
        column = 414348 ;
        pft = 414348 ;
        levsoi = 10 ;
        levlak = 10 ;
        numrad = 2 ;
        lon = 800 ;
        lat = 400 ;
        time = UNLIMITED ; // (1 currently)
        hist_interval = 2 ;
        string_length = 8 ;
variables:
        float lon(lon) ;
                lon:long_name = "coordinate longitude" ;
                lon:units = "degrees_east" ;
        float lat(lat) ;
                lat:long_name = "coordinate latitude" ;
                lat:units = "degrees_north" ;
        float time(time) ;
                time:logname = "time" ;
                time:units = "days since 1990-01-01 00:00:00" ;
                time:calendar = "noleap" ;
        float levsoi(levsoi) ;
                levsoi:long_name = "coordinate soil levels" ;
                levsoi:units = "m" ;
        float levlak(levlak) ;
                levlak:long_name = "coordinate lake levels" ;
                levlak:units = "m" ;
        float longxy(lat, lon) ;
                longxy:long_name = "longitude" ;
                longxy:units = "degrees_east" ;
        float latixy(lat, lon) ;
                latixy:long_name = "latitude" ;
                latixy:units = "degrees_north" ;
        float area(lat, lon) ;
                area:long_name = "grid cell areas" ;
                area:units = "km^2" ;
        int numlon(lat) ;
                numlon:long_name = "number of longitudes at each latitude" ;
        float landfrac(lat, lon) ;
                landfrac:long_name = "land fraction" ;
        float landmask(lat, lon) ;
                landmask:long_name = "land/ocean mask (0.=ocean and 1.=land)" ;
        int mcdate(time) ;
                mcdate:logname = "current date (YYYYMMDD)" ;
        int mcsec(time) ;
                mcsec:logname = "current seconds of current date" ;
                mcsec:units = "s" ;
        int mdcur(time) ;
                mdcur:logname = "current day (from base day)" ;
        int mscur(time) ;
                mscur:logname = "current seconds of current day" ;
        int nstep(time) ;
                nstep:logname = "time step" ;
        double time_bounds(time, hist_interval) ;
                time_bounds:logname = "history time interval endpoints" ;
        char date_written(time, string_length) ;
        char time_written(time, string_length) ;
        float ZSOI(levsoi, lat, lon) ;
                ZSOI:long_name = "soil depth" ;
                ZSOI:units = "m" ;
                ZSOI:_FillValue = 1.e+36f ;
                ZSOI:missing_value = 1.e+36f ;
        float DZSOI(levsoi, lat, lon) ;
                DZSOI:long_name = "soil thickness" ;
                DZSOI:units = "m" ;
                DZSOI:_FillValue = 1.e+36f ;
                DZSOI:missing_value = 1.e+36f ;
        float WATSAT(levsoi, lat, lon) ;
                WATSAT:long_name = "saturated soil water content (porosity)" ;
                WATSAT:units = "mm3/mm3" ;
                WATSAT:_FillValue = 1.e+36f ;
                WATSAT:missing_value = 1.e+36f ;
        float SUCSAT(levsoi, lat, lon) ;
                SUCSAT:long_name = "saturated soil matric potential" ;
                SUCSAT:units = "mm" ;
                SUCSAT:_FillValue = 1.e+36f ;
                SUCSAT:missing_value = 1.e+36f ;
        float BSW(levsoi, lat, lon) ;
                BSW:long_name = "slope of soil water retention curve" ;
                BSW:units = "unitless" ;
                BSW:_FillValue = 1.e+36f ;
                BSW:missing_value = 1.e+36f ;
        float SANDAM(levsoi, lat, lon) ;
                SANDAM:long_name = "soil sand content (percent)" ;
                SANDAM:units = "mm3/mm3" ;
                SANDAM:_FillValue = 1.e+36f ;
                SANDAM:missing_value = 1.e+36f ;
        float CLAYAM(levsoi, lat, lon) ;
                CLAYAM:long_name = "soil clay content (percent)" ;
                CLAYAM:units = "mm3/mm3" ;
                CLAYAM:_FillValue = 1.e+36f ;
                CLAYAM:missing_value = 1.e+36f ;
        float BM1(time, lat, lon) ;
                BM1:long_name = "biomass of leaf" ;
                BM1:units = "umol/m2s" ;
                BM1:cell_method = "time: mean" ;
                BM1:_FillValue = 1.e+36f ;
                BM1:missing_value = 1.e+36f ;
        float BM11(time, lat, lon) ;
                BM11:long_name = "biomass of leaf" ;
                BM11:units = "Kg C/m2" ;
                BM11:cell_method = "time: mean" ;
                BM11:_FillValue = 1.e+36f ;
                BM11:missing_value = 1.e+36f ;
        float BM2(time, lat, lon) ;
                BM2:long_name = "biomass of stem" ;
                BM2:units = "umol/m2s" ;
                BM2:cell_method = "time: mean" ;
                BM2:_FillValue = 1.e+36f ;
                BM2:missing_value = 1.e+36f ;
        float BM21(time, lat, lon) ;
                BM21:long_name = "biomass of stem" ;
                BM21:units = "KgC /m2" ;
                BM21:cell_method = "time: mean" ;
                BM21:_FillValue = 1.e+36f ;
                BM21:missing_value = 1.e+36f ;
        float BM3(time, lat, lon) ;
                BM3:long_name = "biomass of root" ;
                BM3:units = "umol/m2s" ;
                BM3:cell_method = "time: mean" ;
                BM3:_FillValue = 1.e+36f ;
                BM3:missing_value = 1.e+36f ;
        float BM31(time, lat, lon) ;
                BM31:long_name = "biomass of root" ;
                BM31:units = "Kg C/m2" ;
                BM31:cell_method = "time: mean" ;
                BM31:_FillValue = 1.e+36f ;
                BM31:missing_value = 1.e+36f ;
        float BTRAN(time, lat, lon) ;
                BTRAN:long_name = "transpiration beta factor" ;
                BTRAN:units = "unitless" ;
                BTRAN:cell_method = "time: mean" ;
                BTRAN:_FillValue = 1.e+36f ;
                BTRAN:missing_value = 1.e+36f ;
        float COLDDAYS(time, lat, lon) ;
                COLDDAYS:long_name = "days of too cold for leaffall" ;
                COLDDAYS:units = "days" ;
                COLDDAYS:cell_method = "time: instantaneous" ;
                COLDDAYS:_FillValue = 1.e+36f ;
                COLDDAYS:missing_value = 1.e+36f ;
        float DCO2(time, lat, lon) ;
                DCO2:long_name = "soil respiration" ;
                DCO2:units = "g/m2d" ;
                DCO2:cell_method = "time: mean" ;
                DCO2:_FillValue = 1.e+36f ;
                DCO2:missing_value = 1.e+36f ;
        float DIS1(time, lat, lon) ;
                DIS1:long_name = "distribution of npp to leaf" ;
                DIS1:units = "umol/m2s" ;
                DIS1:cell_method = "time: mean" ;
                DIS1:_FillValue = 1.e+36f ;
                DIS1:missing_value = 1.e+36f ;
        float DIS2(time, lat, lon) ;
                DIS2:long_name = "distribution of npp to stem" ;
                DIS2:units = "umol/m2s" ;
                DIS2:cell_method = "time: mean" ;
                DIS2:_FillValue = 1.e+36f ;
                DIS2:missing_value = 1.e+36f ;
        float DIS3(time, lat, lon) ;
                DIS3:long_name = "distribution of npp to root" ;
                DIS3:units = "umol/m2s" ;
                DIS3:cell_method = "time: mean" ;
                DIS3:_FillValue = 1.e+36f ;
                DIS3:missing_value = 1.e+36f ;
        float DSMC(time, lat, lon) ;
                DSMC:long_name = "total soil carbon" ;
                DSMC:units = "g/m2d" ;
                DSMC:cell_method = "time: mean" ;
                DSMC:_FillValue = 1.e+36f ;
                DSMC:missing_value = 1.e+36f ;
        float ELAI(time, lat, lon) ;
                ELAI:long_name = "exposed one-sided leaf area index" ;
                ELAI:units = "m^2/m^2" ;
                ELAI:cell_method = "time: mean" ;
                ELAI:_FillValue = 1.e+36f ;
                ELAI:missing_value = 1.e+36f ;
        float ERRH2O(time, lat, lon) ;
                ERRH2O:long_name = "total water conservation error" ;
                ERRH2O:units = "mm" ;
                ERRH2O:cell_method = "time: mean" ;
                ERRH2O:_FillValue = 1.e+36f ;
                ERRH2O:missing_value = 1.e+36f ;
        float ERRSEB(time, lat, lon) ;
                ERRSEB:long_name = "surface energy conservation error" ;
                ERRSEB:units = "watt/m^2" ;
                ERRSEB:cell_method = "time: mean" ;
                ERRSEB:_FillValue = 1.e+36f ;
                ERRSEB:missing_value = 1.e+36f ;
        float ERRSOI(time, lat, lon) ;
                ERRSOI:long_name = "soil/lake energy conservation error" ;
                ERRSOI:units = "watt/m^2" ;
                ERRSOI:cell_method = "time: mean" ;
                ERRSOI:_FillValue = 1.e+36f ;
                ERRSOI:missing_value = 1.e+36f ;
        float ERRSOL(time, lat, lon) ;
                ERRSOL:long_name = "solar radiation conservation error" ;
                ERRSOL:units = "watt/m^2" ;
                ERRSOL:cell_method = "time: mean" ;
                ERRSOL:_FillValue = 1.e+36f ;
                ERRSOL:missing_value = 1.e+36f ;
        float ESAI(time, lat, lon) ;
                ESAI:long_name = "exposed one-sided stem area index" ;
                ESAI:units = "m^2/m^2" ;
                ESAI:cell_method = "time: mean" ;
                ESAI:_FillValue = 1.e+36f ;
                ESAI:missing_value = 1.e+36f ;
        float FCEV(time, lat, lon) ;
                FCEV:long_name = "canopy evaporation" ;
                FCEV:units = "watt/m^2" ;
                FCEV:cell_method = "time: mean" ;
                FCEV:_FillValue = 1.e+36f ;
                FCEV:missing_value = 1.e+36f ;
        float FCTR(time, lat, lon) ;
                FCTR:long_name = "canopy transpiration" ;
                FCTR:units = "watt/m^2" ;
                FCTR:cell_method = "time: mean" ;
                FCTR:_FillValue = 1.e+36f ;
                FCTR:missing_value = 1.e+36f ;
        float FGEV(time, lat, lon) ;
                FGEV:long_name = "ground evaporation" ;
                FGEV:units = "watt/m^2" ;
                FGEV:cell_method = "time: mean" ;
                FGEV:_FillValue = 1.e+36f ;
                FGEV:missing_value = 1.e+36f ;
        float FGR(time, lat, lon) ;
                FGR:long_name = "heat flux into soil" ;
                FGR:units = "watt/m^2" ;
                FGR:cell_method = "time: mean" ;
                FGR:_FillValue = 1.e+36f ;
                FGR:missing_value = 1.e+36f ;
        float FIRA(time, lat, lon) ;
                FIRA:long_name = "net infrared (longwave) radiation" ;
                FIRA:units = "watt/m^2" ;
                FIRA:cell_method = "time: mean" ;
                FIRA:_FillValue = 1.e+36f ;
                FIRA:missing_value = 1.e+36f ;
        float FIRE(time, lat, lon) ;
                FIRE:long_name = "emitted infrared (longwave) radiation" ;
                FIRE:units = "watt/m^2" ;
                FIRE:cell_method = "time: mean" ;
                FIRE:_FillValue = 1.e+36f ;
                FIRE:missing_value = 1.e+36f ;
        float FLDS(time, lat, lon) ;
                FLDS:long_name = "atmospheric longwave radiation" ;
                FLDS:units = "watt/m^2" ;
                FLDS:cell_method = "time: mean" ;
                FLDS:_FillValue = 1.e+36f ;
                FLDS:missing_value = 1.e+36f ;
        float FORCECO2(time, lat, lon) ;
                FORCECO2:long_name = "atmospheric CO2 concentration" ;
                FORCECO2:units = "ppm" ;
                FORCECO2:cell_method = "time: mean" ;
                FORCECO2:_FillValue = 1.e+36f ;
                FORCECO2:missing_value = 1.e+36f ;
        float FPSN(time, lat, lon) ;
                FPSN:long_name = "photosynthesis" ;
                FPSN:units = "umol/m2s" ;
                FPSN:cell_method = "time: mean" ;
                FPSN:_FillValue = 1.e+36f ;
                FPSN:missing_value = 1.e+36f ;
        float FSA(time, lat, lon) ;
                FSA:long_name = "absorbed solar radiation" ;
                FSA:units = "watt/m^2" ;
                FSA:cell_method = "time: mean" ;
                FSA:_FillValue = 1.e+36f ;
                FSA:missing_value = 1.e+36f ;
        float FSDS(time, lat, lon) ;
                FSDS:long_name = "atmospheric incident solar radiation" ;
                FSDS:units = "watt/m^2" ;
                FSDS:cell_method = "time: mean" ;
                FSDS:_FillValue = 1.e+36f ;
                FSDS:missing_value = 1.e+36f ;
        float FSDSND(time, lat, lon) ;
                FSDSND:long_name = "direct nir incident solar radiation" ;
                FSDSND:units = "watt/m^2" ;
                FSDSND:cell_method = "time: mean" ;
                FSDSND:_FillValue = 1.e+36f ;
                FSDSND:missing_value = 1.e+36f ;
        float FSDSNDLN(time, lat, lon) ;
                FSDSNDLN:long_name = "direct nir incident solar radiation at local noon" ;
                FSDSNDLN:units = "watt/m^2" ;
                FSDSNDLN:cell_method = "time: mean" ;
                FSDSNDLN:_FillValue = 1.e+36f ;
                FSDSNDLN:missing_value = 1.e+36f ;
        float FSDSNI(time, lat, lon) ;
                FSDSNI:long_name = "diffuse nir incident solar radiation" ;
                FSDSNI:units = "watt/m^2" ;
                FSDSNI:cell_method = "time: mean" ;
                FSDSNI:_FillValue = 1.e+36f ;
                FSDSNI:missing_value = 1.e+36f ;
        float FSDSSNO(time, lat, lon) ;
                FSDSSNO:long_name = "snow diffuse nir reflected solar radiation" ;
                FSDSSNO:units = "watt/m^2" ;
                FSDSSNO:cell_method = "time: mean" ;
                FSDSSNO:_FillValue = 1.e+36f ;
                FSDSSNO:missing_value = 1.e+36f ;
        float FSDSVD(time, lat, lon) ;
                FSDSVD:long_name = "direct vis incident solar radiation" ;
                FSDSVD:units = "watt/m^2" ;
                FSDSVD:cell_method = "time: mean" ;
                FSDSVD:_FillValue = 1.e+36f ;
                FSDSVD:missing_value = 1.e+36f ;
        float FSDSVDLN(time, lat, lon) ;
                FSDSVDLN:long_name = "direct vis incident solar radiation at local noon" ;
                FSDSVDLN:units = "watt/m^2" ;
                FSDSVDLN:cell_method = "time: mean" ;
                FSDSVDLN:_FillValue = 1.e+36f ;
                FSDSVDLN:missing_value = 1.e+36f ;
        float FSDSVI(time, lat, lon) ;
                FSDSVI:long_name = "diffuse vis incident solar radiation" ;
                FSDSVI:units = "watt/m^2" ;
                FSDSVI:cell_method = "time: mean" ;
                FSDSVI:_FillValue = 1.e+36f ;
                FSDSVI:missing_value = 1.e+36f ;
        float FSH(time, lat, lon) ;
                FSH:long_name = "sensible heat" ;
                FSH:units = "watt/m^2" ;
                FSH:cell_method = "time: mean" ;
                FSH:_FillValue = 1.e+36f ;
                FSH:missing_value = 1.e+36f ;
        float FSH_G(time, lat, lon) ;
                FSH_G:long_name = "sensible heat from ground" ;
                FSH_G:units = "watt/m^2" ;
                FSH_G:cell_method = "time: mean" ;
                FSH_G:_FillValue = 1.e+36f ;
                FSH_G:missing_value = 1.e+36f ;
        float FSH_V(time, lat, lon) ;
                FSH_V:long_name = "sensible heat from veg" ;
                FSH_V:units = "watt/m^2" ;
                FSH_V:cell_method = "time: mean" ;
                FSH_V:_FillValue = 1.e+36f ;
                FSH_V:missing_value = 1.e+36f ;
        float FSM(time, lat, lon) ;
                FSM:long_name = "snow melt heat flux" ;
                FSM:units = "watt/m^2" ;
                FSM:cell_method = "time: mean" ;
                FSM:_FillValue = 1.e+36f ;
                FSM:missing_value = 1.e+36f ;
        float FSNO(time, lat, lon) ;
                FSNO:long_name = "fraction of ground covered by snow" ;
                FSNO:units = "unitless" ;
                FSNO:cell_method = "time: mean" ;
                FSNO:_FillValue = 1.e+36f ;
                FSNO:missing_value = 1.e+36f ;
        float FSR(time, lat, lon) ;
                FSR:long_name = "reflected solar radiation" ;
                FSR:units = "watt/m^2" ;
                FSR:cell_method = "time: mean" ;
                FSR:_FillValue = 1.e+36f ;
                FSR:missing_value = 1.e+36f ;
        float FSRND(time, lat, lon) ;
                FSRND:long_name = "direct nir reflected solar radiation" ;
                FSRND:units = "watt/m^2" ;
                FSRND:cell_method = "time: mean" ;
                FSRND:_FillValue = 1.e+36f ;
                FSRND:missing_value = 1.e+36f ;
        float FSRNDLN(time, lat, lon) ;
                FSRNDLN:long_name = "direct nir reflected solar radiation at local noon" ;
                FSRNDLN:units = "watt/m^2" ;
                FSRNDLN:cell_method = "time: mean" ;
                FSRNDLN:_FillValue = 1.e+36f ;
                FSRNDLN:missing_value = 1.e+36f ;
        float FSRNI(time, lat, lon) ;
                FSRNI:long_name = "diffuse nir reflected solar radiation" ;
                FSRNI:units = "watt/m^2" ;
                FSRNI:cell_method = "time: mean" ;
                FSRNI:_FillValue = 1.e+36f ;
                FSRNI:missing_value = 1.e+36f ;
        float FSRSNO(time, lat, lon) ;
                FSRSNO:long_name = "snow diffuse nir reflected solar radiation" ;
                FSRSNO:units = "watt/m^2" ;
                FSRSNO:cell_method = "time: mean" ;
                FSRSNO:_FillValue = 1.e+36f ;
                FSRSNO:missing_value = 1.e+36f ;
        float FSRVD(time, lat, lon) ;
                FSRVD:long_name = "direct vis reflected solar radiation" ;
                FSRVD:units = "watt/m^2" ;
                FSRVD:cell_method = "time: mean" ;
                FSRVD:_FillValue = 1.e+36f ;
                FSRVD:missing_value = 1.e+36f ;
        float FSRVDLN(time, lat, lon) ;
                FSRVDLN:long_name = "direct vis reflected solar radiation at local noon" ;
                FSRVDLN:units = "watt/m^2" ;
                FSRVDLN:cell_method = "time: mean" ;
                FSRVDLN:_FillValue = 1.e+36f ;
                FSRVDLN:missing_value = 1.e+36f ;
        float FSRVI(time, lat, lon) ;
                FSRVI:long_name = "diffuse vis reflected solar radiation" ;
                FSRVI:units = "watt/m^2" ;
                FSRVI:cell_method = "time: mean" ;
                FSRVI:_FillValue = 1.e+36f ;
                FSRVI:missing_value = 1.e+36f ;
        float FTOUT(time, lat, lon) ;
                FTOUT:long_name = "actual temperature adjustment" ;
                FTOUT:units = "g/m2d" ;
                FTOUT:cell_method = "time: mean" ;
                FTOUT:_FillValue = 1.e+36f ;
                FTOUT:missing_value = 1.e+36f ;
        float H2OCAN(time, lat, lon) ;
                H2OCAN:long_name = "intercepted water" ;
                H2OCAN:units = "mm" ;
                H2OCAN:cell_method = "time: mean" ;
                H2OCAN:_FillValue = 1.e+36f ;
                H2OCAN:missing_value = 1.e+36f ;
        float H2OSNO(time, lat, lon) ;
                H2OSNO:long_name = "snow depth (liquid water)" ;
                H2OSNO:units = "mm" ;
                H2OSNO:cell_method = "time: mean" ;
                H2OSNO:_FillValue = 1.e+36f ;
                H2OSNO:missing_value = 1.e+36f ;
        float H2OSOI(time, levsoi, lat, lon) ;
                H2OSOI:long_name = "volumetric soil water" ;
                H2OSOI:units = "mm3/mm3" ;
                H2OSOI:cell_method = "time: mean" ;
                H2OSOI:_FillValue = 1.e+36f ;
                H2OSOI:missing_value = 1.e+36f ;
        float LEAFSTATUS(time, lat, lon) ;
                LEAFSTATUS:long_name = "leaf status" ;
                LEAFSTATUS:units = "unitless" ;
                LEAFSTATUS:cell_method = "time: instantaneous" ;
                LEAFSTATUS:_FillValue = 1.e+36f ;
                LEAFSTATUS:missing_value = 1.e+36f ;
        float LOS1(time, lat, lon) ;
                LOS1:long_name = "loss of leaf" ;
                LOS1:units = "umol/m2s" ;
                LOS1:cell_method = "time: mean" ;
                LOS1:_FillValue = 1.e+36f ;
                LOS1:missing_value = 1.e+36f ;
        float LOS2(time, lat, lon) ;
                LOS2:long_name = "loss of stem" ;
                LOS2:units = "umol/m2s" ;
                LOS2:cell_method = "time: mean" ;
                LOS2:_FillValue = 1.e+36f ;
                LOS2:missing_value = 1.e+36f ;
        float LOS21(time, lat, lon) ;
                LOS21:long_name = "harvest loss of stem" ;
                LOS21:units = "umol/m2s" ;
                LOS21:cell_method = "time: mean" ;
                LOS21:_FillValue = 1.e+36f ;
                LOS21:missing_value = 1.e+36f ;
        float LOS3(time, lat, lon) ;
                LOS3:long_name = "loss of root" ;
                LOS3:units = "umol/m2s" ;
                LOS3:cell_method = "time: mean" ;
                LOS3:_FillValue = 1.e+36f ;
                LOS3:missing_value = 1.e+36f ;
        float LOS31(time, lat, lon) ;
                LOS31:long_name = "harvest loss of root" ;
                LOS31:units = "umol/m2s" ;
                LOS31:cell_method = "time: mean" ;
                LOS31:_FillValue = 1.e+36f ;
                LOS31:missing_value = 1.e+36f ;
        float MATUREDAYS(time, lat, lon) ;
                MATUREDAYS:long_name = "days of mature leaf" ;
                MATUREDAYS:units = "days" ;
                MATUREDAYS:cell_method = "time: instantaneous" ;
                MATUREDAYS:_FillValue = 1.e+36f ;
                MATUREDAYS:missing_value = 1.e+36f ;
        float MAXDAYS(time, lat, lon) ;
                MAXDAYS:long_name = "days of maximum growth" ;
                MAXDAYS:units = "days" ;
                MAXDAYS:cell_method = "time: instantaneous" ;
                MAXDAYS:_FillValue = 1.e+36f ;
                MAXDAYS:missing_value = 1.e+36f ;
        float NEP(time, lat, lon) ;
                NEP:long_name = "net ecosystem productivity" ;
                NEP:units = "g/m2d" ;
                NEP:cell_method = "time: mean" ;
                NEP:_FillValue = 1.e+36f ;
                NEP:missing_value = 1.e+36f ;
        float NORMALDAYS(time, lat, lon) ;
                NORMALDAYS:long_name = "days of normal growth" ;
                NORMALDAYS:units = "days" ;
                NORMALDAYS:cell_method = "time: instantaneous" ;
                NORMALDAYS:_FillValue = 1.e+36f ;
                NORMALDAYS:missing_value = 1.e+36f ;
        float NPP(time, lat, lon) ;
                NPP:long_name = "net primary productivity" ;
                NPP:units = "umol/m2s" ;
                NPP:cell_method = "time: mean" ;
                NPP:_FillValue = 1.e+36f ;
                NPP:missing_value = 1.e+36f ;
        float PANDAYS(time, lat, lon) ;
                PANDAYS:long_name = "days of mature leaf" ;
                PANDAYS:units = "days" ;
                PANDAYS:cell_method = "time: instantaneous" ;
                PANDAYS:_FillValue = 1.e+36f ;
                PANDAYS:missing_value = 1.e+36f ;
        float PHT(time, lat, lon) ;
                PHT:long_name = "photosynthesis" ;
                PHT:units = "umol/m2s" ;
                PHT:cell_method = "time: mean" ;
                PHT:_FillValue = 1.e+36f ;
                PHT:missing_value = 1.e+36f ;
        float POLC1(time, lat, lon) ;
                POLC1:long_name = "soil carbon pool 1" ;
                POLC1:units = "g/m2d" ;
                POLC1:cell_method = "time: mean" ;
                POLC1:_FillValue = 1.e+36f ;
                POLC1:missing_value = 1.e+36f ;
        float POLC2(time, lat, lon) ;
                POLC2:long_name = "soil carbon pool 2" ;
                POLC2:units = "g/m2d" ;
                POLC2:cell_method = "time: mean" ;
                POLC2:_FillValue = 1.e+36f ;
                POLC2:missing_value = 1.e+36f ;
        float POLC3(time, lat, lon) ;
                POLC3:long_name = "soil carbon pool 3" ;
                POLC3:units = "g/m2d" ;
                POLC3:cell_method = "time: mean" ;
                POLC3:_FillValue = 1.e+36f ;
                POLC3:missing_value = 1.e+36f ;
        float POLC4(time, lat, lon) ;
                POLC4:long_name = "soil carbon pool 4" ;
                POLC4:units = "g/m2d" ;
                POLC4:cell_method = "time: mean" ;
                POLC4:_FillValue = 1.e+36f ;
                POLC4:missing_value = 1.e+36f ;
        float POLC5(time, lat, lon) ;
                POLC5:long_name = "soil carbon pool 5" ;
                POLC5:units = "g/m2d" ;
                POLC5:cell_method = "time: mean" ;
                POLC5:_FillValue = 1.e+36f ;
                POLC5:missing_value = 1.e+36f ;
        float POLC6(time, lat, lon) ;
                POLC6:long_name = "soil carbon pool 6" ;
                POLC6:units = "g/m2d" ;
                POLC6:cell_method = "time: mean" ;
                POLC6:_FillValue = 1.e+36f ;
                POLC6:missing_value = 1.e+36f ;
        float POLC7(time, lat, lon) ;
                POLC7:long_name = "soil carbon pool 7" ;
                POLC7:units = "g/m2d" ;
                POLC7:cell_method = "time: mean" ;
                POLC7:_FillValue = 1.e+36f ;
                POLC7:missing_value = 1.e+36f ;
        float POLC8(time, lat, lon) ;
                POLC8:long_name = "soil carbon pool 8" ;
                POLC8:units = "g/m2d" ;
                POLC8:cell_method = "time: mean" ;
                POLC8:_FillValue = 1.e+36f ;
                POLC8:missing_value = 1.e+36f ;
        float Q2M(time, lat, lon) ;
                Q2M:long_name = "2m specific humidity" ;
                Q2M:units = "kg/kg" ;
                Q2M:cell_method = "time: mean" ;
                Q2M:_FillValue = 1.e+36f ;
                Q2M:missing_value = 1.e+36f ;
        float QBOT(time, lat, lon) ;
                QBOT:long_name = "atmospheric specific humidity" ;
                QBOT:units = "kg/kg" ;
                QBOT:cell_method = "time: mean" ;
                QBOT:_FillValue = 1.e+36f ;
                QBOT:missing_value = 1.e+36f ;
        float QDRAI(time, lat, lon) ;
                QDRAI:long_name = "sub-surface drainage" ;
                QDRAI:units = "mm/s" ;
                QDRAI:cell_method = "time: mean" ;
                QDRAI:_FillValue = 1.e+36f ;
                QDRAI:missing_value = 1.e+36f ;
        float QDRIP(time, lat, lon) ;
                QDRIP:long_name = "throughfall" ;
                QDRIP:units = "mm/s" ;
                QDRIP:cell_method = "time: mean" ;
                QDRIP:_FillValue = 1.e+36f ;
                QDRIP:missing_value = 1.e+36f ;
        float QINFL(time, lat, lon) ;
                QINFL:long_name = "infiltration" ;
                QINFL:units = "mm/s" ;
                QINFL:cell_method = "time: mean" ;
                QINFL:_FillValue = 1.e+36f ;
                QINFL:missing_value = 1.e+36f ;
        float QINTR(time, lat, lon) ;
                QINTR:long_name = "interception" ;
                QINTR:units = "mm/s" ;
                QINTR:cell_method = "time: mean" ;
                QINTR:_FillValue = 1.e+36f ;
                QINTR:missing_value = 1.e+36f ;
        float QMELT(time, lat, lon) ;
                QMELT:long_name = "snow melt" ;
                QMELT:units = "mm/s" ;
                QMELT:cell_method = "time: mean" ;
                QMELT:_FillValue = 1.e+36f ;
                QMELT:missing_value = 1.e+36f ;
        float QOVER(time, lat, lon) ;
                QOVER:long_name = "surface runoff" ;
                QOVER:units = "mm/s" ;
                QOVER:cell_method = "time: mean" ;
                QOVER:_FillValue = 1.e+36f ;
                QOVER:missing_value = 1.e+36f ;
        float QRGWL(time, lat, lon) ;
                QRGWL:long_name = "surface runoff at glaciers, wetlands, lakes" ;
                QRGWL:units = "mm/s" ;
                QRGWL:cell_method = "time: mean" ;
                QRGWL:_FillValue = 1.e+36f ;
                QRGWL:missing_value = 1.e+36f ;
        float QSOIL(time, lat, lon) ;
                QSOIL:long_name = "ground evaporation" ;
                QSOIL:units = "mm/s" ;
                QSOIL:cell_method = "time: mean" ;
                QSOIL:_FillValue = 1.e+36f ;
                QSOIL:missing_value = 1.e+36f ;
        float QSUBLIM(time, lat, lon) ;
                QSUBLIM:long_name = "sublimation rate from snow pack" ;
                QSUBLIM:units = "mmH2O/s" ;
                QSUBLIM:cell_method = "time: mean" ;
                QSUBLIM:_FillValue = 1.e+36f ;
                QSUBLIM:missing_value = 1.e+36f ;
        float QVEGE(time, lat, lon) ;
                QVEGE:long_name = "canopy evaporation" ;
                QVEGE:units = "mm/s" ;
                QVEGE:cell_method = "time: mean" ;
                QVEGE:_FillValue = 1.e+36f ;
                QVEGE:missing_value = 1.e+36f ;
        float QVEGT(time, lat, lon) ;
                QVEGT:long_name = "canopy transpiration" ;
                QVEGT:units = "mm/s" ;
                QVEGT:cell_method = "time: mean" ;
                QVEGT:_FillValue = 1.e+36f ;
                QVEGT:missing_value = 1.e+36f ;
        float RAIN(time, lat, lon) ;
                RAIN:long_name = "atmospheric rain" ;
                RAIN:units = "mm/s" ;
                RAIN:cell_method = "time: mean" ;
                RAIN:_FillValue = 1.e+36f ;
                RAIN:missing_value = 1.e+36f ;
        float RGP(time, lat, lon) ;
                RGP:long_name = "growth respiration" ;
                RGP:units = "umol/m2s" ;
                RGP:cell_method = "time: mean" ;
                RGP:_FillValue = 1.e+36f ;
                RGP:missing_value = 1.e+36f ;
        float RMP(time, lat, lon) ;
                RMP:long_name = "maintenance respiration" ;
                RMP:units = "umol/m2s" ;
                RMP:cell_method = "time: mean" ;
                RMP:_FillValue = 1.e+36f ;
                RMP:missing_value = 1.e+36f ;
        float ROTL(time, lat, lon) ;
                ROTL:long_name = "root carbon input to soil" ;
                ROTL:units = "umol/m2s" ;
                ROTL:cell_method = "time: mean" ;
                ROTL:_FillValue = 1.e+36f ;
                ROTL:missing_value = 1.e+36f ;
        float RRESIS(time, levsoi, lat, lon) ;
                RRESIS:long_name = "root resistance in each soil layer" ;
                RRESIS:units = "proportion" ;
                RRESIS:cell_method = "time: mean" ;
                RRESIS:_FillValue = 1.e+36f ;
                RRESIS:missing_value = 1.e+36f ;
        float RSSHA(time, lat, lon) ;
                RSSHA:long_name = "shaded leaf stomatal resistance" ;
                RSSHA:units = "s/m" ;
                RSSHA:cell_method = "time: minimum" ;
                RSSHA:_FillValue = 1.e+36f ;
                RSSHA:missing_value = 1.e+36f ;
        float RSSUN(time, lat, lon) ;
                RSSUN:long_name = "sunlit leaf stomatal resistance" ;
                RSSUN:units = "s/m" ;
                RSSUN:cell_method = "time: minimum" ;
                RSSUN:_FillValue = 1.e+36f ;
                RSSUN:missing_value = 1.e+36f ;
        float SABG(time, lat, lon) ;
                SABG:long_name = "solar rad absorbed by ground" ;
                SABG:units = "watt/m^2" ;
                SABG:cell_method = "time: mean" ;
                SABG:_FillValue = 1.e+36f ;
                SABG:missing_value = 1.e+36f ;
        float SABV(time, lat, lon) ;
                SABV:long_name = "solar rad absorbed by veg" ;
                SABV:units = "watt/m^2" ;
                SABV:cell_method = "time: mean" ;
                SABV:_FillValue = 1.e+36f ;
                SABV:missing_value = 1.e+36f ;
        float SNOW(time, lat, lon) ;
                SNOW:long_name = "atmospheric snow" ;
                SNOW:units = "mm/s" ;
                SNOW:cell_method = "time: mean" ;
                SNOW:_FillValue = 1.e+36f ;
                SNOW:missing_value = 1.e+36f ;
        float SNOWAGE(time, lat, lon) ;
                SNOWAGE:long_name = "snow age" ;
                SNOWAGE:units = "unitless" ;
                SNOWAGE:cell_method = "time: mean" ;
                SNOWAGE:_FillValue = 1.e+36f ;
                SNOWAGE:missing_value = 1.e+36f ;
        float SNOWDP(time, lat, lon) ;
                SNOWDP:long_name = "snow height" ;
                SNOWDP:units = "m" ;
                SNOWDP:cell_method = "time: mean" ;
                SNOWDP:_FillValue = 1.e+36f ;
                SNOWDP:missing_value = 1.e+36f ;
        float SNOWICE(time, lat, lon) ;
                SNOWICE:long_name = "snow ice" ;
                SNOWICE:units = "kg/m2" ;
                SNOWICE:cell_method = "time: mean" ;
                SNOWICE:_FillValue = 1.e+36f ;
                SNOWICE:missing_value = 1.e+36f ;
        float SNOWLIQ(time, lat, lon) ;
                SNOWLIQ:long_name = "snow liquid water" ;
                SNOWLIQ:units = "kg/m2" ;
                SNOWLIQ:cell_method = "time: mean" ;
                SNOWLIQ:_FillValue = 1.e+36f ;
                SNOWLIQ:missing_value = 1.e+36f ;
        float SOILICE(time, levsoi, lat, lon) ;
                SOILICE:long_name = "soil ice" ;
                SOILICE:units = "kg/m2" ;
                SOILICE:cell_method = "time: mean" ;
                SOILICE:_FillValue = 1.e+36f ;
                SOILICE:missing_value = 1.e+36f ;
        float SOILLIQ(time, levsoi, lat, lon) ;
                SOILLIQ:long_name = "soil liquid water" ;
                SOILLIQ:units = "kg/m2" ;
                SOILLIQ:cell_method = "time: mean" ;
                SOILLIQ:_FillValue = 1.e+36f ;
                SOILLIQ:missing_value = 1.e+36f ;
        float SURL(time, lat, lon) ;
                SURL:long_name = "surface carbon input to soil" ;
                SURL:units = "umol/m2s" ;
                SURL:cell_method = "time: mean" ;
                SURL:_FillValue = 1.e+36f ;
                SURL:missing_value = 1.e+36f ;
        float TAUX(time, lat, lon) ;
                TAUX:long_name = "zonal surface stress" ;
                TAUX:units = "kg/m/s^2" ;
                TAUX:cell_method = "time: mean" ;
                TAUX:_FillValue = 1.e+36f ;
                TAUX:missing_value = 1.e+36f ;
        float TAUY(time, lat, lon) ;
                TAUY:long_name = "meridional surface stress" ;
                TAUY:units = "kg/m/s^2" ;
                TAUY:cell_method = "time: mean" ;
                TAUY:_FillValue = 1.e+36f ;
                TAUY:missing_value = 1.e+36f ;
        float TBOT(time, lat, lon) ;
                TBOT:long_name = "atmospheric air temperature" ;
                TBOT:units = "K" ;
                TBOT:cell_method = "time: mean" ;
                TBOT:_FillValue = 1.e+36f ;
                TBOT:missing_value = 1.e+36f ;
        float TG(time, lat, lon) ;
                TG:long_name = "ground temperature" ;
                TG:units = "K" ;
                TG:cell_method = "time: mean" ;
                TG:_FillValue = 1.e+36f ;
                TG:missing_value = 1.e+36f ;
        float THBOT(time, lat, lon) ;
                THBOT:long_name = "atmospheric air potential temperature" ;
                THBOT:units = "K" ;
                THBOT:cell_method = "time: mean" ;
                THBOT:_FillValue = 1.e+36f ;
                THBOT:missing_value = 1.e+36f ;
        float TLAKE(time, levsoi, lat, lon) ;
                TLAKE:long_name = "lake temperature" ;
                TLAKE:units = "K" ;
                TLAKE:cell_method = "time: mean" ;
                TLAKE:_FillValue = 1.e+36f ;
                TLAKE:missing_value = 1.e+36f ;
        float TREFMNAV(time, lat, lon) ;
                TREFMNAV:long_name = "daily minimum of average 2-m temperature" ;
                TREFMNAV:units = "K" ;
                TREFMNAV:cell_method = "time: mean" ;
                TREFMNAV:_FillValue = 1.e+36f ;
                TREFMNAV:missing_value = 1.e+36f ;
        float TREFMXAV(time, lat, lon) ;
                TREFMXAV:long_name = "daily maximum of average 2-m temperature" ;
                TREFMXAV:units = "K" ;
                TREFMXAV:cell_method = "time: mean" ;
                TREFMXAV:_FillValue = 1.e+36f ;
                TREFMXAV:missing_value = 1.e+36f ;
        float TSA(time, lat, lon) ;
                TSA:long_name = "2m air temperature" ;
                TSA:units = "K" ;
                TSA:cell_method = "time: mean" ;
                TSA:_FillValue = 1.e+36f ;
                TSA:missing_value = 1.e+36f ;
        float TSNOW(time, lat, lon) ;
                TSNOW:long_name = "snow temperature" ;
                TSNOW:units = "K" ;
                TSNOW:cell_method = "time: mean" ;
                TSNOW:_FillValue = 1.e+36f ;
                TSNOW:missing_value = 1.e+36f ;
        float TSOI(time, levsoi, lat, lon) ;
                TSOI:long_name = "soil temperature" ;
                TSOI:units = "K" ;
                TSOI:cell_method = "time: mean" ;
                TSOI:_FillValue = 1.e+36f ;
                TSOI:missing_value = 1.e+36f ;
        float TV(time, lat, lon) ;
                TV:long_name = "vegetation temperature" ;
                TV:units = "K" ;
                TV:cell_method = "time: mean" ;
                TV:_FillValue = 1.e+36f ;
                TV:missing_value = 1.e+36f ;
        float VCMXOUT(time, lat, lon) ;
                VCMXOUT:long_name = "actual VCMX" ;
                VCMXOUT:units = "g/m2d" ;
                VCMXOUT:cell_method = "time: mean" ;
                VCMXOUT:_FillValue = 1.e+36f ;
                VCMXOUT:missing_value = 1.e+36f ;
        float WIND(time, lat, lon) ;
                WIND:long_name = "atmospheric wind velocity magnitude" ;
                WIND:units = "m/s" ;
                WIND:cell_method = "time: mean" ;
                WIND:_FillValue = 1.e+36f ;
                WIND:missing_value = 1.e+36f ;
        float ZBOT(time, lat, lon) ;
                ZBOT:long_name = "atmospheric reference height" ;
                ZBOT:units = "m" ;
                ZBOT:cell_method = "time: mean" ;
                ZBOT:_FillValue = 1.e+36f ;
                ZBOT:missing_value = 1.e+36f ;

// global attributes:
                :conventions = "CF1.0" ;
                :history = "created on 04/28/20 16:15:26" ;
                :logname = "zhouwy" ;
                :host = "cmb05n21" ;
                :source = "Community Land Model: CLM2" ;
                :version = "$Name: ccsm3_0_rel04 $" ;
                :revision_id = "$Id: histFileMod.F90,v 1.19.4.47 2004/06/10 16:41:39 mvertens Exp $" ;
                :case_title = "bccavim.T266crop" ;
                :case_id = "bccavim.T266crop" ;
                :Surface_dataset = "surface-data_T266_20180921.nc" ;
                :Initial_conditions_dataset = "arbitrary initialization" ;
                :PFT_physiological_constants_dataset = "pft-physiology.USGS.AVIM2" ;
                :RTM_input_datset = "rdirc.05" ;
data:

lon = 0, 0.45, 0.9, 1.35, 1.8, 2.25, 2.7, 3.15, 3.6, 4.05, 4.5, 4.95, 5.4,
    5.85, 6.3, 6.75, 7.2, 7.65, 8.1, 8.55, 9, 9.45, 9.9, 10.35, 10.8, 11.25,
    11.7, 12.15, 12.6, 13.05, 13.5, 13.95, 14.4, 14.85, 15.3, 15.75, 16.2,
    16.65, 17.1, 17.55, 18, 18.45, 18.9, 19.35, 19.8, 20.25, 20.7, 21.15,
    21.6, 22.05, 22.5, 22.95, 23.4, 23.85, 24.3, 24.75, 25.2, 25.65, 26.1,
    26.55, 27, 27.45, 27.9, 28.35, 28.8, 29.25, 29.7, 30.15, 30.6, 31.05,
    31.5, 31.95, 32.4, 32.85, 33.3, 33.75, 34.2, 34.65, 35.1, 35.55, 36,
    36.45, 36.9, 37.35, 37.8, 38.25, 38.7, 39.15, 39.6, 40.05, 40.5, 40.95,
    41.4, 41.85, 42.3, 42.75, 43.2, 43.65, 44.1, 44.55, 45, 45.45, 45.9,
    46.35, 46.8, 47.25, 47.7, 48.15, 48.6, 49.05, 49.5, 49.95, 50.4, 50.85,
    51.3, 51.75, 52.2, 52.65, 53.1, 53.55, 54, 54.45, 54.9, 55.35, 55.8,
    56.25, 56.7, 57.15, 57.6, 58.05, 58.5, 58.95, 59.4, 59.85, 60.3, 60.75,
    61.2, 61.65, 62.1, 62.55, 63, 63.45, 63.9, 64.35, 64.8, 65.25, 65.7,
    66.15, 66.6, 67.05, 67.5, 67.95, 68.4, 68.85, 69.3, 69.75, 70.2, 70.65,
    71.1, 71.55, 72, 72.45, 72.9, 73.35, 73.8, 74.25, 74.7, 75.15, 75.6,
    76.05, 76.5, 76.95, 77.4, 77.85, 78.3, 78.75, 79.2, 79.65, 80.1, 80.55,
    81, 81.45, 81.9, 82.35, 82.8, 83.25, 83.7, 84.15, 84.6, 85.05, 85.5,
    85.95, 86.4, 86.85, 87.3, 87.75, 88.2, 88.65, 89.1, 89.55, 90, 90.45,
    90.9, 91.35, 91.8, 92.25, 92.7, 93.15, 93.6, 94.05, 94.5, 94.95, 95.4,
    95.85, 96.3, 96.75, 97.2, 97.65, 98.1, 98.55, 99, 99.45, 99.9, 100.35,
    100.8, 101.25, 101.7, 102.15, 102.6, 103.05, 103.5, 103.95, 104.4,
    104.85, 105.3, 105.75, 106.2, 106.65, 107.1, 107.55, 108, 108.45, 108.9,
    109.35, 109.8, 110.25, 110.7, 111.15, 111.6, 112.05, 112.5, 112.95,
    113.4, 113.85, 114.3, 114.75, 115.2, 115.65, 116.1, 116.55, 117, 117.45,
    117.9, 118.35, 118.8, 119.25, 119.7, 120.15, 120.6, 121.05, 121.5,
    121.95, 122.4, 122.85, 123.3, 123.75, 124.2, 124.65, 125.1, 125.55, 126,
    126.45, 126.9, 127.35, 127.8, 128.25, 128.7, 129.15, 129.6, 130.05,
    130.5, 130.95, 131.4, 131.85, 132.3, 132.75, 133.2, 133.65, 134.1,
    134.55, 135, 135.45, 135.9, 136.35, 136.8, 137.25, 137.7, 138.15, 138.6,
    139.05, 139.5, 139.95, 140.4, 140.85, 141.3, 141.75, 142.2, 142.65,
    143.1, 143.55, 144, 144.45, 144.9, 145.35, 145.8, 146.25, 146.7, 147.15,
    147.6, 148.05, 148.5, 148.95, 149.4, 149.85, 150.3, 150.75, 151.2,
    151.65, 152.1, 152.55, 153, 153.45, 153.9, 154.35, 154.8, 155.25, 155.7,
    156.15, 156.6, 157.05, 157.5, 157.95, 158.4, 158.85, 159.3, 159.75,
    160.2, 160.65, 161.1, 161.55, 162, 162.45, 162.9, 163.35, 163.8, 164.25,
    164.7, 165.15, 165.6, 166.05, 166.5, 166.95, 167.4, 167.85, 168.3,
    168.75, 169.2, 169.65, 170.1, 170.55, 171, 171.45, 171.9, 172.35, 172.8,
    173.25, 173.7, 174.15, 174.6, 175.05, 175.5, 175.95, 176.4, 176.85,
    177.3, 177.75, 178.2, 178.65, 179.1, 179.55, 180, 180.45, 180.9, 181.35,
    181.8, 182.25, 182.7, 183.15, 183.6, 184.05, 184.5, 184.95, 185.4,
    185.85, 186.3, 186.75, 187.2, 187.65, 188.1, 188.55, 189, 189.45, 189.9,
    190.35, 190.8, 191.25, 191.7, 192.15, 192.6, 193.05, 193.5, 193.95,
    194.4, 194.85, 195.3, 195.75, 196.2, 196.65, 197.1, 197.55, 198, 198.45,
    198.9, 199.35, 199.8, 200.25, 200.7, 201.15, 201.6, 202.05, 202.5,
    202.95, 203.4, 203.85, 204.3, 204.75, 205.2, 205.65, 206.1, 206.55, 207,
    207.45, 207.9, 208.35, 208.8, 209.25, 209.7, 210.15, 210.6, 211.05,
    211.5, 211.95, 212.4, 212.85, 213.3, 213.75, 214.2, 214.65, 215.1,
    215.55, 216, 216.45, 216.9, 217.35, 217.8, 218.25, 218.7, 219.15, 219.6,
    220.05, 220.5, 220.95, 221.4, 221.85, 222.3, 222.75, 223.2, 223.65,
    224.1, 224.55, 225, 225.45, 225.9, 226.35, 226.8, 227.25, 227.7, 228.15,
    228.6, 229.05, 229.5, 229.95, 230.4, 230.85, 231.3, 231.75, 232.2,
    232.65, 233.1, 233.55, 234, 234.45, 234.9, 235.35, 235.8, 236.25, 236.7,
    237.15, 237.6, 238.05, 238.5, 238.95, 239.4, 239.85, 240.3, 240.75,
    241.2, 241.65, 242.1, 242.55, 243, 243.45, 243.9, 244.35, 244.8, 245.25,
    245.7, 246.15, 246.6, 247.05, 247.5, 247.95, 248.4, 248.85, 249.3,
    249.75, 250.2, 250.65, 251.1, 251.55, 252, 252.45, 252.9, 253.35, 253.8,
    254.25, 254.7, 255.15, 255.6, 256.05, 256.5, 256.95, 257.4, 257.85,
    258.3, 258.75, 259.2, 259.65, 260.1, 260.55, 261, 261.45, 261.9, 262.35,
    262.8, 263.25, 263.7, 264.15, 264.6, 265.05, 265.5, 265.95, 266.4,
    266.85, 267.3, 267.75, 268.2, 268.65, 269.1, 269.55, 270, 270.45, 270.9,
    271.35, 271.8, 272.25, 272.7, 273.15, 273.6, 274.05, 274.5, 274.95,
    275.4, 275.85, 276.3, 276.75, 277.2, 277.65, 278.1, 278.55, 279, 279.45,
    279.9, 280.35, 280.8, 281.25, 281.7, 282.15, 282.6, 283.05, 283.5,
    283.95, 284.4, 284.85, 285.3, 285.75, 286.2, 286.65, 287.1, 287.55, 288,
    288.45, 288.9, 289.35, 289.8, 290.25, 290.7, 291.15, 291.6, 292.05,
    292.5, 292.95, 293.4, 293.85, 294.3, 294.75, 295.2, 295.65, 296.1,
    296.55, 297, 297.45, 297.9, 298.35, 298.8, 299.25, 299.7, 300.15, 300.6,
    301.05, 301.5, 301.95, 302.4, 302.85, 303.3, 303.75, 304.2, 304.65,
    305.1, 305.55, 306, 306.45, 306.9, 307.35, 307.8, 308.25, 308.7, 309.15,
    309.6, 310.05, 310.5, 310.95, 311.4, 311.85, 312.3, 312.75, 313.2,
    313.65, 314.1, 314.55, 315, 315.45, 315.9, 316.35, 316.8, 317.25, 317.7,
    318.15, 318.6, 319.05, 319.5, 319.95, 320.4, 320.85, 321.3, 321.75,
    322.2, 322.65, 323.1, 323.55, 324, 324.45, 324.9, 325.35, 325.8, 326.25,
    326.7, 327.15, 327.6, 328.05, 328.5, 328.95, 329.4, 329.85, 330.3,
    330.75, 331.2, 331.65, 332.1, 332.55, 333, 333.45, 333.9, 334.35, 334.8,
    335.25, 335.7, 336.15, 336.6, 337.05, 337.5, 337.95, 338.4, 338.85,
    339.3, 339.75, 340.2, 340.65, 341.1, 341.55, 342, 342.45, 342.9, 343.35,
    343.8, 344.25, 344.7, 345.15, 345.6, 346.05, 346.5, 346.95, 347.4,
    347.85, 348.3, 348.75, 349.2, 349.65, 350.1, 350.55, 351, 351.45, 351.9,
    352.35, 352.8, 353.25, 353.7, 354.15, 354.6, 355.05, 355.5, 355.95,
    356.4, 356.85, 357.3, 357.75, 358.2, 358.65, 359.1, 359.55 ;

lat = -89.65597, -89.2103, -88.76199, -88.3131, -87.86398, -87.41474,
    -86.96545, -86.51612, -86.06676, -85.61739, -85.16801, -84.71862,
    -84.26922, -83.81981, -83.3704, -82.92099, -82.47157, -82.02216,
    -81.57273, -81.12331, -80.67389, -80.22446, -79.77503, -79.32561,
    -78.87617, -78.42675, -77.97732, -77.52789, -77.07846, -76.62903,
    -76.1796, -75.73016, -75.28073, -74.8313, -74.38187, -73.93243,
    -73.48299, -73.03356, -72.58413, -72.13469, -71.68526, -71.23582,
    -70.78639, -70.33695, -69.88752, -69.43808, -68.98865, -68.53922,
    -68.08978, -67.64034, -67.1909, -66.74147, -66.29203, -65.8426,
    -65.39316, -64.94373, -64.49429, -64.04485, -63.59542, -63.14598,
    -62.69654, -62.2471, -61.79767, -61.34823, -60.8988, -60.44936,
    -59.99992, -59.55049, -59.10105, -58.65161, -58.20218, -57.75274,
    -57.3033, -56.85386, -56.40443, -55.95499, -55.50555, -55.05611,
    -54.60668, -54.15724, -53.70781, -53.25837, -52.80893, -52.35949,
    -51.91006, -51.46062, -51.01118, -50.56174, -50.1123, -49.66287,
    -49.21343, -48.76399, -48.31456, -47.86512, -47.41568, -46.96624,
    -46.51681, -46.06737, -45.61793, -45.1685, -44.71906, -44.26962,
    -43.82018, -43.37074, -42.92131, -42.47187, -42.02243, -41.57299,
    -41.12356, -40.67412, -40.22468, -39.77524, -39.32581, -38.87637,
    -38.42693, -37.97749, -37.52806, -37.07862, -36.62918, -36.17974,
    -35.7303, -35.28087, -34.83143, -34.38199, -33.93255, -33.48312,
    -33.03368, -32.58424, -32.1348, -31.68537, -31.23593, -30.78649,
    -30.33705, -29.88762, -29.43818, -28.98874, -28.5393, -28.08986,
    -27.64043, -27.19099, -26.74155, -26.29211, -25.84267, -25.39324,
    -24.9438, -24.49436, -24.04492, -23.59549, -23.14605, -22.69661,
    -22.24717, -21.79774, -21.3483, -20.89886, -20.44942, -19.99998,
    -19.55055, -19.10111, -18.65167, -18.20223, -17.75279, -17.30336,
    -16.85392, -16.40448, -15.95504, -15.50561, -15.05617, -14.60673,
    -14.15729, -13.70785, -13.25842, -12.80898, -12.35954, -11.9101,
    -11.46066, -11.01123, -10.56179, -10.11235, -9.662913, -9.213476,
    -8.764038, -8.3146, -7.865162, -7.415724, -6.966287, -6.516849,
    -6.067411, -5.617973, -5.168535, -4.719098, -4.26966, -3.820222,
    -3.370784, -2.921346, -2.471908, -2.02247, -1.573032, -1.123595,
    -0.6741568, -0.2247189, 0.2247189, 0.6741568, 1.123595, 1.573032,
    2.02247, 2.471908, 2.921346, 3.370784, 3.820222, 4.26966, 4.719098,
    5.168535, 5.617973, 6.067411, 6.516849, 6.966287, 7.415724, 7.865162,
    8.3146, 8.764038, 9.213476, 9.662913, 10.11235, 10.56179, 11.01123,
    11.46066, 11.9101, 12.35954, 12.80898, 13.25842, 13.70785, 14.15729,
    14.60673, 15.05617, 15.50561, 15.95504, 16.40448, 16.85392, 17.30336,
    17.75279, 18.20223, 18.65167, 19.10111, 19.55055, 19.99998, 20.44942,
    20.89886, 21.3483, 21.79774, 22.24717, 22.69661, 23.14605, 23.59549,
    24.04492, 24.49436, 24.9438, 25.39324, 25.84267, 26.29211, 26.74155,
    27.19099, 27.64043, 28.08986, 28.5393, 28.98874, 29.43818, 29.88762,
    30.33705, 30.78649, 31.23593, 31.68537, 32.1348, 32.58424, 33.03368,
    33.48312, 33.93255, 34.38199, 34.83143, 35.28087, 35.7303, 36.17974,
    36.62918, 37.07862, 37.52806, 37.97749, 38.42693, 38.87637, 39.32581,
    39.77524, 40.22468, 40.67412, 41.12356, 41.57299, 42.02243, 42.47187,
    42.92131, 43.37074, 43.82018, 44.26962, 44.71906, 45.1685, 45.61793,
    46.06737, 46.51681, 46.96624, 47.41568, 47.86512, 48.31456, 48.76399,
    49.21343, 49.66287, 50.1123, 50.56174, 51.01118, 51.46062, 51.91006,
    52.35949, 52.80893, 53.25837, 53.70781, 54.15724, 54.60668, 55.05611,
    55.50555, 55.95499, 56.40443, 56.85386, 57.3033, 57.75274, 58.20218,
    58.65161, 59.10105, 59.55049, 59.99992, 60.44936, 60.8988, 61.34823,
    61.79767, 62.2471, 62.69654, 63.14598, 63.59542, 64.04485, 64.49429,
    64.94373, 65.39316, 65.8426, 66.29203, 66.74147, 67.1909, 67.64034,
    68.08978, 68.53922, 68.98865, 69.43808, 69.88752, 70.33695, 70.78639,
    71.23582, 71.68526, 72.13469, 72.58413, 73.03356, 73.48299, 73.93243,
    74.38187, 74.8313, 75.28073, 75.73016, 76.1796, 76.62903, 77.07846,
    77.52789, 77.97732, 78.42675, 78.87617, 79.32561, 79.77503, 80.22446,
    80.67389, 81.12331, 81.57273, 82.02216, 82.47157, 82.92099, 83.3704,
    83.81981, 84.26922, 84.71862, 85.16801, 85.61739, 86.06676, 86.51612,
    86.96545, 87.41474, 87.86398, 88.3131, 88.76199, 89.2103, 89.65597 ;

time = 243 ;

levsoi = 0.007100635, 0.027925, 0.06225858, 0.1188651, 0.2121934, 0.3660658,
    0.6197585, 1.038027, 1.727635, 2.864607 ;

levlak = 0.05, 0.6, 2.1, 4.6, 8.1, 12.6, 18.6, 25.6, 34.325, 44.775 ;
}



密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表