Something wrong with the boundary input datasets

[dfeijat@126_com]

Dear CAM3.0 developers,
When I try to plot the file :sst_HadOIBl_bc_64x128_clim_c020411.nc, some problems encounter me like the followings:

ga-> sdfopen d:sst_HadOIBl_bc_64x128_clim_c020411.nc
Scanning self-describing file: d:sst_HadOIBl_bc_64x128_clim_c020411.nc
SDF file has no discernable time coordinate.
Found 2 data variables in SDF file.
Found displayable variable sst_cpl with 0 levels in SDF file.
Found displayable variable ice_cov with 0 levels in SDF file.
SDF file d:sst_HadOIBl_bc_64x128_clim_c020411.nc is open as file 1
LON set to 0 360
LAT set to -87.8638 87.8638
LEV set to 0 0
Time values set: 1:1:1:0 1:1:1:0
ga-> q file
File 1 :
Descriptor: d:sst_HadOIBl_bc_64x128_clim_c020411.nc
Binary: d:sst_HadOIBl_bc_64x128_clim_c020411.nc
Type = Gridded
Xsize = 128 Ysize = 64 Zsize = 1 Tsize = 1
Number of Variables = 2
sst_cpl 0 -999 BCS Pseudo SST
ice_cov 0 -999 BCS Pseudo Sea-ice concentration
ga-> q ctlinfo
dset d:sst_HadOIBl_bc_64x128_clim_c020411.nc
title
undef 1e+36
xdef 128 linear 0 2.8125
ydef 64 levels -87.8638 -85.0965 -82.3129 -79.5256 -76.7369 -73.9475 -71.1578 -6 8.3678
-65.5776 -62.7874 -59.997 -57.2066 -54.4162 -51.6257 -48.8352 -46.0447 -43.2542 -40.4636
-37.6731 -34.8825 -32.0919 -29.3014 -26.5108 -23.7202 -20.9296 -18.139 -15.3484 -12.5578
-9.76715 -6.97653 -4.18592 -1.39531 1.39531 4.18592 6.97653 9.76715 12.5578 15. 3484
18.139 20.9296 23.7202 26.5108 29.3014 32.0919 34.8825 37.6731 40.4636 43.2542
46.0447 48.8352 51.6257 54.4162 57.2066 59.997 62.7874 65.5776 68.3678 71.1578
73.9475 76.7369 79.5256 82.3129 85.0965 87.8638
zdef 1 linear 0 0
tdef 1 linear 00Z01JAN1 1mn
vars 2
sst_cpl 0 -999 BCS Pseudo SST
ice_cov 0 -999 BCS Pseudo Sea-ice concentration
ga-> q dim
Default file number is: 1
X is varying Lon = 0 to 360 X = 1 to 129
Y is varying Lat = -87.8638 to 87.8638 Y = 1 to 64
Z is fixed Lev = 0 Z = 1
T is fixed Time = 00Z01JAN1 T = 1

I use the software GrADS1.8sL11. Trying to sdfopen the file, I fould it`s hard to reach my destination.
Why the Time is only set to 1 month which should be 12 months for the climatological SST.
The same problem encounters with the file :sst_HadOIBl_bc_64x128_1949-2001_c020411.nc which the Tsize should be 636 months.

How can I correct it?

Thanks and Regards.

-----------------------

Yours!
defeijat@126.com

cams.cma

 

[eaton]

I'm not familiar with the GrADS software, but from scanning the output from the sdfopen command it appears that the only problem is not being able to recognize a time coordinate. This is a problem with the file metadata. The reason this file works in CAM is that only the data stored in the date and datesec variables are used for time interpolation. The file contains a time variable, but its metadata does not follow the COARDS or CF conventions and I believe that's why GrADS is unable to recognize it. We are aware of this problem and will fix it in future datasets.

 

[dfeijat@126_com]

Dear eaton,
Thanks for your reply.
It is appreciative that you said that you would fix the problem in future datasets. But, could you tell me how can I deal with the present datasets in advance? Cauze I want to analyze it to reveal the information it contains which I need in my essay before I run it on the CAM model.
Maybe it is too fulsome for you to answer my question. But any of your hints would help me greatly.

Thanks in advance!

----------------------
Yours!

dfeijat@126.com

cams.cma

 

[dfeijat@126_com]

Or could you tell me which software could recognize the time coordinate contains in the files?
By the way, when would the new-meliorative datasets appear?

 

[dfeijat@126_com]

I'm not familiar with the GrADS software, but from scanning the output from the sdfopen command it appears that the only problem is not being able to recognize a time coordinate. This is a problem with the file metadata. The reason this file works in CAM is that only the data stored in the date and datesec variables are used for time interpolation. The file contains a time variable, but its metadata does not follow the COARDS or CF conventions and I believe that's why GrADS is unable to recognize it. We are aware of this problem and will fix it in future datasets.

------------------------------------------------------------------------
========================================================

Dear CAM3.0 developers,
I am trying to run the cam3.0 with the monthly observed SSTs, sst_HadOIBl_bc_128x256_1949_2001_c020812.nc, but when I submit the job, the model stop with the followings errors:
+16382 0: Successfully initialized time variant variables
+16383 0:
+16384 0: Successfully initialized the land model
+16385 0: begin initial run at:
+16386 0: nstep= 0 year= 0 month= 9 day= 1 seconds= 0
+16387 0:
+16388 0:******************************************************
+16389 0:
+16390 0: SSTINI: Failed to find dates bracketing ncdate, ncsec= 831 85800
+16391 0: ENDRUN: called without a message string


My question is, is there something wrong with the dataset, sst_HadOIBl_bc_128x256_1949_2001_c020812.nc, or somewhere in my namelist?

To show the errors that encountered me clearly, I give the namelist with which I cam verbosely as following:

+1 &camexp
+2 absems_data = '/u/zuozhy/du/t85/atm/cam2/rad/abs_ems_factors_fastvx.c030508.nc'
+3 aeroptics = '/u/zuozhy/du/t85/atm/cam2/rad/AerosolOptics_c040105.nc'
+4 archive_dir = '/archive/00/zuozhy/case.t85.dom.sst.52'
+5 bndtvaer = '/u/zuozhy/du/t85/atm/cam2/rad/AerosolMass_V_128x256_clim_c031022.nc'
+6 bndtvdms = '/u/zuozhy/du/t85/atm/cam2/scyc/DMS_emissions_128x256_clim_c040122.nc'
+7 bndtvghg = '/u/zuozhy/du/t85/atm/cam2/ggas/ghg_1870_2100_c040122.nc'
+8 bndtvo = '/u/zuozhy/du/t85/atm/cam2/ozone/pcmdio3.r8.64x1_L60_clim_c970515.nc'
+9 bndtvoxid = '/u/zuozhy/du/t85/atm/cam2/scyc/oxid_128x256_L26_clim_c040112.nc'
+10 bndtvs = '/u/zuozhy/du/t85/atm/cam2/sst/sst_HadOIBl_bc_128x256_1949_2001_c020812.nc'
+11 bndtvscon = '/u/zuozhy/du/t85/atm/cam2/rad/scon_1870_2100_c040122.nc'
+12 bndtvsox = '/u/zuozhy/du/t85/atm/cam2/scyc/SOx_emissions_128x256_L2_1870-1871_c040520.nc'
+13 bndtvvolc = '/u/zuozhy/du/t85/atm/cam2/rad/VolcanicMass_1870-1999_64x1_L18_c040115.nc'
+14 caseid = 'case.t85.dom.sst.52'
+15 dtime = 600
+16 empty_htapes = .TRUE.
+17 fexcl1 = 'CLDHGH', 'CLDICE', 'CLDLIQ', 'CLDLOW', 'CLDMED', 'CLDTOT', 'CLOUD', 'CONCLD',
+18 'FICE', 'GCLDLWP', 'ICLDIWP', 'ICLDLWP', 'LWCF', 'PBLH', 'SWCF', 'TGCLDIWP',
+19 'TGCLDLWP'
+20 fincl1 = 'CMFDQ', 'CMFDQR', 'CMFDT', 'DCQ', 'DTCOND', 'DTH', 'DTV', 'ETADOT', 'FLN200',
+21 'FLN200C', 'FLNS', 'FLNSC', 'FLNSOI', 'FLNT', 'FLNTC', 'FLUT', 'FLUTC', 'FSDS',
+22 'FSDSC', 'FSN200', 'FSN200C', 'FSNS', 'FSNSC', 'FSNSOI', 'FSNT', 'FSNTC', 'FSNTOA',
+23 'FSNTOAC', 'FU', 'FV', 'LANDFRAC', 'LHFLX', 'LHFLXOI', 'LPSTEN', 'NSTEP', 'OCNFRAC',
+24 'OMEGA', 'OMEGA500', 'OMEGA850', 'OMEGAT', 'OMEGAU', 'PBOT', 'PDELDRY', 'PHIS', 'PRECC',
+25 'PRECCFRQ', 'PRECCINT', 'PRECCav', 'PRECL', 'PRECLFRQ', 'PRECLINT', 'PRECLav', 'PRECSC', 'PRECSH',
+26 'PRECSL', 'PRECT', 'PRECTMX', 'PS', 'PSDRY', 'PSL', 'Q', 'Q200', 'Q850',
+27 'QC', 'QFLX', 'QRL', 'QRS', 'RELHUM', 'SFCLDICE', 'SFCLDLIQ', 'SFQ', 'SGH',
+28 'SHFLX', 'SHFLXOI', 'SNOWHICE', 'SNOWHLND', 'SOLIN', 'SRFRAD', 'SST', 'T', 'T300',
+29 'T850', 'TAUGWX', 'TAUGWY', 'TAUTMSX', 'TAUTMSY', 'TAUX', 'TAUY', 'TBOT', 'TMQ',
+30 'TREFHT', 'TS', 'TSMN', 'TSMX', 'TTEND', 'U', 'U200', 'U850', 'UBOT',
+31 'US', 'USTAR', 'UTEND', 'UU', 'V', 'V200', 'V850', 'VBOT', 'VD01',
+32 'VQ', 'VT', 'VU', 'VV', 'WSPEED', 'Z050', 'Z3', 'Z300', 'Z500',
+33 'Z700', 'ZBOT'
+34 iyear_ad = 1950
+35 ncdata = '/u/zuozhy/du/t85/atm/cam2/inic/gaus/cami_0000-09-01_128x256_L26_c040422.nc'
+36 nelapse = -18980
+37 nhtfrq = 0
+38 nrefrq = 0
+39 nsrest = 0
+40 sstcyc = .false.
+41 /
+42 &clmexp
+43 finidat = '/u/zuozhy/du/t85/lnd/clm2/inidata_2.1/cam/clmi_0000-09-01_128x256_c040422.nc'
+44 fpftcon = '/u/zuozhy/du/t85/lnd/clm2/pftdata/pft-physiology'
+45 fsurdat = '/u/zuozhy/du/t85/lnd/clm2/srfdata/cam/clms_128x256_c031031.nc'
+46 /
"

Shall you tell me how can I correct the errors?

Thanks in advance!
---------------------------------------
Yours!
dfeijat@126.com
cams.cma
Friday, 10, June, 2005


BNDTVS='/u/zuozhy/du/t85/atm/cam2/sst/sst_HadOIBl_bc_128x256_1949_2001_c020812.nc'

=======================================================

Dear eaton,

Is this the same problem with the file sst_HadOIBl_bc_128x256_1949_2001_c020812.nc that I couldn`t plot out by GrADS, while in this time, it doesn`t work in CAM3.0.

Could you give me some clues to accomplish my job.

Thank you in advance!

-----------------------
Yours!
dfeijat@126.com
cams.cma
Friday,10/06/2005

 

[eaton]

The problem is in the namelist. By default the initial time of the run is read from the initial file which in this case is 0000-09-01. This time is not contained in the SST dataset you specified. Set an appropriate initial date using the namelist variable start_ymd.

 

[dfeijat@126_com]

The problem is in the namelist. By default the initial time of the run is read from the initial file which in this case is 0000-09-01. This time is not contained in the SST dataset you specified. Set an appropriate initial date using the namelist variable start_ymd.

=======================================================

Dear eaton,

Thank you for your help first.
When I set the namelist variable start_ymd=194909 & et al. correspondingly, it really works.

But I am confused about why should I set the start_ymd et al., could the model initial time of the run be set to be 0000-01-01 ?

Or in other words, why the model set to start run from 0000-09-01 but not 0000-01-01 ? Could you tell me the reason ? Or could you give me some conceptual hints about ?

Thank you in advance !

------------------------------------------------------

Yours!
dfeijat@126.com
cams.cma
Saturday,11/06/2005

 

[dfeijat@126_com]

Dear eaton,

I got another error with the inputdataset like the following:

+16278 0:-----------------------------------------------------------
+16279 0:
+16280 0: Initializing variables for time accumulation .....
+16281 0:------------------------------------------------------------
+16282 0: Accumulated fields
+16283 0: No Name Units Type Period Inival Description
+16284 0:_______________________________________________________________________
+16285 0: 1 TREFAV K timeavg 6 0. average over an hour of 2-m temperature
+16286 0:_______________________________________________________________________
+16287 0:
+16288 0:------------------------------------------------------------
+16289 0: Successfully initialized variables for accumulation
+16290 0:
+16291 0: Attempting to initialize time variant variables
+16292 0: Reading initial data from initial dataset
+16293 0: (GETFIL): attempting to find local file clmi_0000-09-01_128x256_c040422.nc
+16294 0: (GETFIL): using /u/zuozhy/du/t85/lnd/clm2/inidata_2.1/cam/clmi_0000-09-01_128x256_c040422.nc
+16295 0: Attempting to read monthly vegetation data .....
+16296 0: nstep = 0 month = 49 day = 9
+16297 0: (GETFIL): attempting to find local file clms_128x256_c031031.nc
+16298 0: (GETFIL): using /u/zuozhy/du/t85/lnd/clm2/srfdata/cam/clms_128x256_c031031.nc
+16299 0: netcdf error from readMonthlyVegetation
+16300 0: ENDRUN: called without a message string

-------------------------------------------------------------

Could you help me corret the error ?

The namelist I used included as following:

+1 &camexp
+2 absems_data = '/u/zuozhy/du/t85/atm/cam2/rad/abs_ems_factors_fastvx.c030508.nc'
+3 aeroptics = '/u/zuozhy/du/t85/atm/cam2/rad/AerosolOptics_c040105.nc'
+4 bndtvaer = '/u/zuozhy/du/t85/atm/cam2/rad/AerosolMass_V_128x256_clim_c031022.nc'
+5 bndtvdms = '/u/zuozhy/du/t85/atm/cam2/scyc/DMS_emissions_128x256_clim_c040122.nc'
+6 bndtvghg = '/u/zuozhy/du/t85/atm/cam2/ggas/ghg_1870_2100_c040122.nc'
+7 bndtvo = '/u/zuozhy/du/t85/atm/cam2/ozone/pcmdio3.r8.64x1_L60_clim_c970515.nc'
+8 bndtvoxid = '/u/zuozhy/du/t85/atm/cam2/scyc/oxid_128x256_L26_clim_c040112.nc'
+9 bndtvs = '/u/zuozhy/du/t85/atm/cam2/sst/sst_HadOIBl_bc_128x256_1949_2001_c020812.nc'
+10 bndtvscon = '/u/zuozhy/du/t85/atm/cam2/rad/scon_1870_2100_c040122.nc'
+11 bndtvsox = '/u/zuozhy/du/t85/atm/cam2/scyc/SOx_emissions_128x256_L2_1870-1871_c040520.nc'
+12 bndtvvolc = '/u/zuozhy/du/t85/atm/cam2/rad/VolcanicMass_1870-1999_64x1_L18_c040115.nc'
+13 caseid = 'case.t85.dom.sst.52'
+14 ctitle = 'multi-year194909-200109 SST dataset'
+15 dtime = 600
+16 empty_htapes = .TRUE.
+17 fexcl1 = 'CLDHGH', 'CLDICE', 'CLDLIQ', 'CLDLOW', 'CLDMED', 'CLDTOT', 'CLOUD', 'CONCLD',
+18 'FICE', 'GCLDLWP', 'ICLDIWP', 'ICLDLWP', 'LWCF', 'PBLH', 'SWCF', 'TGCLDIWP',
+19 'TGCLDLWP'
+20 fincl1 = 'CMFDQ', 'CMFDQR', 'CMFDT', 'DCQ', 'DTCOND', 'DTH', 'DTV', 'ETADOT', 'FLN200',
+21 'FLN200C', 'FLNS', 'FLNSC', 'FLNSOI', 'FLNT', 'FLNTC', 'FLUT', 'FLUTC', 'FSDS',
+22 'FSDSC', 'FSN200', 'FSN200C', 'FSNS', 'FSNSC', 'FSNSOI', 'FSNT', 'FSNTC', 'FSNTOA',
+23 'FSNTOAC', 'FU', 'FV', 'LANDFRAC', 'LHFLX', 'LHFLXOI', 'LPSTEN', 'NSTEP', 'OCNFRAC',
+24 'OMEGA', 'OMEGA500', 'OMEGA850', 'OMEGAT', 'OMEGAU', 'PBOT', 'PDELDRY', 'PHIS', 'PRECC',
+25 'PRECCFRQ', 'PRECCINT', 'PRECCav', 'PRECL', 'PRECLFRQ', 'PRECLINT', 'PRECLav', 'PRECSC', 'PRECSH',
+26 'PRECSL', 'PRECT', 'PRECTMX', 'PS', 'PSDRY', 'PSL', 'Q', 'Q200', 'Q850',
+27 'QC', 'QFLX', 'QRL', 'QRS', 'RELHUM', 'SFCLDICE', 'SFCLDLIQ', 'SFQ', 'SGH',
+28 'SHFLX', 'SHFLXOI', 'SNOWHICE', 'SNOWHLND', 'SOLIN', 'SRFRAD', 'SST', 'T', 'T300',
+29 'T850', 'TAUGWX', 'TAUGWY', 'TAUTMSX', 'TAUTMSY', 'TAUX', 'TAUY', 'TBOT', 'TMQ',
+30 'TREFHT', 'TS', 'TSMN', 'TSMX', 'TTEND', 'U', 'U200', 'U850', 'UBOT',
+31 'US', 'USTAR', 'UTEND', 'UU', 'V', 'V200', 'V850', 'VBOT', 'VD01',
+32 'VQ', 'VT', 'VU', 'VV', 'WSPEED', 'Z050', 'Z3', 'Z300', 'Z500',
+33 'Z700', 'ZBOT'
+34 inithist = 'none'
+35 iyear_ad = 1950
+36 mss_irt = 0
+37 ncdata = '/u/zuozhy/du/t85/atm/cam2/inic/gaus/cami_0000-09-01_128x256_L26_c040422.nc'
+38 nelapse = -18980
+39 nhtfrq = 0
+40 nrefrq = 0
+41 nsrest = 0
+42 rest_pfile = '/dgpfs/fs1/zuozhy/t85.dom.sst.52/case.t85.dom.sst.52/cam3.rpointer'
+43 sstcyc = .false.
+44 start_tod = 0
+45 start_ymd = 194909
+46 stop_tod = 0
+47 stop_ymd = 20020101
+48 /
+49 &clmexp
+50 finidat = '/u/zuozhy/du/t85/lnd/clm2/inidata_2.1/cam/clmi_0000-09-01_128x256_c040422.nc'
+51 fpftcon = '/u/zuozhy/du/t85/lnd/clm2/pftdata/pft-physiology'
+52 fsurdat = '/u/zuozhy/du/t85/lnd/clm2/srfdata/cam/clms_128x256_c031031.nc'
+53 /


Thanks in advance !
-------------------------------------
Yours!
dfeijat@126.com
cams.cma
Saturday,11/06/2005

 

[dfeijat@126_com]

It is my incaution that made the so-called error encountered me.

I check it more and more times and finally found the mistake I had made which I give as following:


+45 start_ymd = 194909

should be set to
start_ymd = 19490901

It is ok now !



But why the initial run time should be set to 0000-09-01 still confuse me.

I am looking forward to hearing from you as soon as possible.

Thanks and regards !

-----------------------------------------------------------
Sincerely !
dfeijat@126.com
cams.cma
Sunday,12/06/2005

 

[louisgarcia723@yahoo_com]

The BOUNDARY= data set provides the boundary information for the sequential test. At stage

, the data set is usually created with an ODS OUTPUT statement from the "Boundary Information" table created by the SEQDESIGN procedure. At each subsequent stage, the data set is usually created with an ODS OUTPUT statement from the "Test Information" table that was created by the SEQTEST procedure at the previous stage. See the term paper writing service section: SEQTEST Procedure for an illustration of the BOUNDARY= data set option.The BOUNDARY= data set contains the following variables:[*]_Scale_, the boundary scale, with the value MLE for the maximum likelihood estimate, STDZ for the standardized Z, SCORE for the score statistic, or PVALUE for the nominal

-value. Note that for a two-sided design, the nominal

-value is the one-sided fixed-sample

-value under the null hypothesis with a lower alternative hypothesis.[*]_Stop_, the stopping criterion, with the value REJECT for rejecting the null hypothesis

, ACCEPT for accepting

, or BOTH for both rejecting and accepting

[*]_ALT_, the type of alternative hypothesis, with the value UPPER for an upper alternative, LOWER for a lower alternative, or TWOSIDED for a two-sided alternative[*]_Stage_, the stage number[*]the boundary variables, a subset of Bound_LA for lower

boundary, Bound_LB for lower

boundary, Bound_UB for upper

boundary, and Bound_UA for upper

boundary[*]AltRef_L, the lower alternative reference, if ALT=LOWER or ALT=TWOSIDED[*]AltRef_U, the upper alternative reference, if ALT=UPPER or ALT=TWOSIDED[*]_InfoProp_, the information proportion at each stage[/list]