sean_m_davis@noaa_gov
Member
I am attempting to force an f-compset (F_2000_CAM5) run with the SST/sea ice output from a b-compset (B_2000) run using the procedure described here
http://www.cesm.ucar.edu/models/cesm1.2/cesm/doc/usersguide/x2288.html
The problem I'm running into is related to the fact that my B-compset output is on an irregular 384x320 grid --
If I follow the recipe above (i.e., I leave alone the OCN_DOMAIN_* and ICE_DOMAIN_* variables in env_run.xml), then my run crashes on the OCN initialization:
[Sean.M.Davis@fe4 run]$ cat ocn.log.130911-181900
(docn_comp_init) ocn_in = docn_ocn_in
(docn_comp_init) decomp = 1d
(docn_comp_init) restfilm = undefined
(docn_comp_init) restfils = undefined
(shr_strdata_readnml) reading input namelist file: docn_ocn_in
(shr_stream_init) Reading file docn.streams.txt.prescribed
(shr_file_freeUnit) unit 94 was not in use
(docn_comp_init) ocn mode = SSTDATA
(shr_strdata_init) calling shr_dmodel_mapSet for fill
(shr_strdata_init) calling shr_dmodel_mapSet for remap
(shr_map_checkGrids_global) ERROR Xsrc not monotonic 1 321.6875089518 320.5625089205
(shr_map_checkGrids_global) ERROR Xsrc not regular lat,lon 1 188 320.5625114233 320.5625114233
(shr_map_checkGrids_global) ERROR Ysrc not regular lat,lon 36 1 -79.22052259842 -79.22052260746
(shr_sys_abort) ERROR: ('shr_map_checkGrids_global') ERROR
(shr_sys_abort) WARNING: calling shr_mpi_abort() and stopping
[Sean.M.Davis@fe4 run]$
If I then modify the OCN_DOMAIN_* to point to the same domain file as SSTICE, the OCN initialization completes, but the model fails on the ICE initialization:
--- FYI, the OCN_DOMAIN and SSTICE_GRID files are as follows (in env_run)
[Sean.M.Davis@fe4 run]$ tail ocn.log.130911-174453
(shr_strdata_print)
(shr_strdata_print) nvectors = 0
(shr_strdata_print) ----------------------------------------------------------
(docn_comp_init) initialize gsmaps
(docn_comp_init) copy domains
(docn_comp_init) allocate AVs
(shr_dmodel_readLBUB) reading file: /scratch2/portfolios/BMC/chimera/sdavis/cesm-1.1.0/scripts/F_2000_CAM5cesm1.1_dom_bdo1979_1981/B_2000_CAM5cesm1.1bdo3_1979_1981_avg.pop.nc 1200
(shr_dmodel_readLBUB) reading file: /scratch2/portfolios/BMC/chimera/sdavis/cesm-1.1.0/scripts/F_2000_CAM5cesm1.1_dom_bdo1979_1981/B_2000_CAM5cesm1.1bdo3_1979_1981_avg.pop.nc 1
(docn_comp_run) ocn: model date 1510101 0s
(docn_comp_init) docn_comp_init done
(see attached ice.log.130911-174453 for details on what crashed during ice initialization)
I thought that the maybe the cause of the crash was that my ICE_DOMAIN_* variables needed to be set similarly to the OCN_DOMAIN_* variables, so i tried that -- i.e.,
But when I did that, I got a different error in the ice initialization:
[Sean.M.Davis@fe4 run]$ tail ice.log.130911-165544
Distribution type: blkrobin
Distribution weight: latitude
{min,max}Blocks = 1 5
Number of ghost cells: 1
Active processors: 192 192
Block size: nx_block = 6
ny_block = 6
[Sean.M.Davis@fe4 run]$ tail cpl.log.130911-165544
(seq_mct_drv) : Initialize each component: atm, lnd, rof, ocn, ice, and glc
(seq_mct_drv) : Initialize atm component ATM
(seq_mct_drv) : Initialize lnd component LND
(seq_mct_drv) : Initialize rof component ROF
(seq_mct_drv) : Initialize ocn component OCN
(seq_mct_drv) : Initialize ice component ICE
(shr_sys_abort) ERROR: latlongrid: ni,nj not equal to nx_global,ny_global
(shr_sys_abort) WARNING: calling shr_mpi_abort() and stopping
(see also attached cpl/ice logs with ending 165544)
I notice in the CICE namelist documentation (http://www.cesm.ucar.edu/models/cesm1.2/cesm/doc/modelnl/nl_cice.html) that there is a group of variables called CICE grid settings (grid_nml). Currently, the grid_type is 'lonlatgrid'. I'm wondering if some of the "grid_nml" settings in my ice_in need to be changed? The documentation online is not really helpful for figuring out what I need to change though.
I've attached my ice_in for reference.
http://www.cesm.ucar.edu/models/cesm1.2/cesm/doc/usersguide/x2288.html
The problem I'm running into is related to the fact that my B-compset output is on an irregular 384x320 grid --
If I follow the recipe above (i.e., I leave alone the OCN_DOMAIN_* and ICE_DOMAIN_* variables in env_run.xml), then my run crashes on the OCN initialization:
[Sean.M.Davis@fe4 run]$ cat ocn.log.130911-181900
(docn_comp_init) ocn_in = docn_ocn_in
(docn_comp_init) decomp = 1d
(docn_comp_init) restfilm = undefined
(docn_comp_init) restfils = undefined
(shr_strdata_readnml) reading input namelist file: docn_ocn_in
(shr_stream_init) Reading file docn.streams.txt.prescribed
(shr_file_freeUnit) unit 94 was not in use
(docn_comp_init) ocn mode = SSTDATA
(shr_strdata_init) calling shr_dmodel_mapSet for fill
(shr_strdata_init) calling shr_dmodel_mapSet for remap
(shr_map_checkGrids_global) ERROR Xsrc not monotonic 1 321.6875089518 320.5625089205
(shr_map_checkGrids_global) ERROR Xsrc not regular lat,lon 1 188 320.5625114233 320.5625114233
(shr_map_checkGrids_global) ERROR Ysrc not regular lat,lon 36 1 -79.22052259842 -79.22052260746
(shr_sys_abort) ERROR: ('shr_map_checkGrids_global') ERROR
(shr_sys_abort) WARNING: calling shr_mpi_abort() and stopping
[Sean.M.Davis@fe4 run]$
If I then modify the OCN_DOMAIN_* to point to the same domain file as SSTICE, the OCN initialization completes, but the model fails on the ICE initialization:
--- FYI, the OCN_DOMAIN and SSTICE_GRID files are as follows (in env_run)
[Sean.M.Davis@fe4 run]$ tail ocn.log.130911-174453
(shr_strdata_print)
(shr_strdata_print) nvectors = 0
(shr_strdata_print) ----------------------------------------------------------
(docn_comp_init) initialize gsmaps
(docn_comp_init) copy domains
(docn_comp_init) allocate AVs
(shr_dmodel_readLBUB) reading file: /scratch2/portfolios/BMC/chimera/sdavis/cesm-1.1.0/scripts/F_2000_CAM5cesm1.1_dom_bdo1979_1981/B_2000_CAM5cesm1.1bdo3_1979_1981_avg.pop.nc 1200
(shr_dmodel_readLBUB) reading file: /scratch2/portfolios/BMC/chimera/sdavis/cesm-1.1.0/scripts/F_2000_CAM5cesm1.1_dom_bdo1979_1981/B_2000_CAM5cesm1.1bdo3_1979_1981_avg.pop.nc 1
(docn_comp_run) ocn: model date 1510101 0s
(docn_comp_init) docn_comp_init done
(see attached ice.log.130911-174453 for details on what crashed during ice initialization)
I thought that the maybe the cause of the crash was that my ICE_DOMAIN_* variables needed to be set similarly to the OCN_DOMAIN_* variables, so i tried that -- i.e.,
But when I did that, I got a different error in the ice initialization:
[Sean.M.Davis@fe4 run]$ tail ice.log.130911-165544
Distribution type: blkrobin
Distribution weight: latitude
{min,max}Blocks = 1 5
Number of ghost cells: 1
Active processors: 192 192
Block size: nx_block = 6
ny_block = 6
[Sean.M.Davis@fe4 run]$ tail cpl.log.130911-165544
(seq_mct_drv) : Initialize each component: atm, lnd, rof, ocn, ice, and glc
(seq_mct_drv) : Initialize atm component ATM
(seq_mct_drv) : Initialize lnd component LND
(seq_mct_drv) : Initialize rof component ROF
(seq_mct_drv) : Initialize ocn component OCN
(seq_mct_drv) : Initialize ice component ICE
(shr_sys_abort) ERROR: latlongrid: ni,nj not equal to nx_global,ny_global
(shr_sys_abort) WARNING: calling shr_mpi_abort() and stopping
(see also attached cpl/ice logs with ending 165544)
I notice in the CICE namelist documentation (http://www.cesm.ucar.edu/models/cesm1.2/cesm/doc/modelnl/nl_cice.html) that there is a group of variables called CICE grid settings (grid_nml). Currently, the grid_type is 'lonlatgrid'. I'm wondering if some of the "grid_nml" settings in my ice_in need to be changed? The documentation online is not really helpful for figuring out what I need to change though.
I've attached my ice_in for reference.