Slab ocean gx3v7 grid issue

[shineng_hu@gmail_com]

I was trying to run CESM 1.0.6 (compset E1850CN; resolution T31_g37), but the run failed due to a model grid issue. The default 'pop_frc.1x1d.090130.nc' worked fine, but it failed when I changed to use 'pop_frc.gx3v7.110128.nc'. I have pasted the error message at the end of message. I tried a similar case but for the f19_g16 resolution, and it worked fine with either 'pop_frc.1x1d.090130.nc' or 'pop_frc.gx1v6.180416.nc'. Does anyone know how to fix the problem?Eventually I would like to diganose the q-flux from a fully coupled run (compset B1850CN; resolution T31_g37), and use it for slab-ocean simulation (compset E1850CN; resolution T31_g37). If the grid issue is not solved, I wonder would it be scientifically valid to interpolate the g37 q-flux to 1x1 degree so that the slab ocean can be smoothly running? Any information regarding those two questions will be appreciated!Shineng===============(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.stream.txt(shr_file_freeUnit) unit   95 was not in use(docn_comp_init)  ocn mode = SOM(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 increasing       11  2.900000000000      359.3000000000    (shr_map_checkGrids_global)  ERROR Xsrc not regular lat,lon        1       6  323.3000000000      323.3000000000(shr_map_checkGrids_global)  ERROR Ysrc not regular lat,lon        2      50 0.2997398353064     0.2997382175421(shr_sys_abort) ERROR: ('shr_map_checkGrids_global')  ERROR(shr_sys_abort) WARNING: calling shr_mpi_abort() and stopping================

 

[dbailey]

What is the resolution combination you are using for your E compset? It is trying to interpolate from the curvilinear gx3 grid, which the data ocean cannot do.
Dave

 

[dbailey]

This is puzzling. I wonder if the domain file is inconsistent with the SOM forcing file. What is your domainfile in docn_ocn_in?

 

[shineng_hu@gmail_com]

Hi Dave, here are the domain files I used. The resolution is T31_g37. DOMAINFILE = /glade/p/cesm/cseg/inputdata/ocn/docn7/domain.ocn.gx3v7.090903.ncdomain1 = /glade/p/cesm/cseg/inputdata/ocn/docn7/SOM/pop_frc.gx3v7.110128.ncfile1 = /glade/p/cesm/cseg/inputdata/ocn/docn7/SOM/pop_frc.gx3v7.110128.nc Shineng 

 

[dbailey]

Try using the pop_frc.gx3v7.110128.nc as your domain file as well. I believe the model is somehow thinking it needs to interpolate and it shouldn't be doing that. Dave

 

[shineng_hu@gmail_com]

Hi Dave,I tried to point the domain file to pop_frc.gx3v7.110128.nc as you suggested, but it failed since it doesn't have the variable 'frac'. Next I tried to add 'frac' to the forcing file, but it still didn't work. Below are the reported errors.=========48:(shr_stream_set) size of filename =      148:(shr_stream_set) filename = /glade/p/cesm/cseg/inputdata/lnd/clm2/ndepdata/fndep_clm_hist_simyr1849-2006_1.9x2.5_c100428.nc59:(shr_stream_set) size of filename =      159:(shr_stream_set) filename = /glade/p/cesm/cseg/inputdata/lnd/clm2/ndepdata/fndep_clm_hist_simyr1849-2006_1.9x2.5_c100428.nc0: Reading setup_nml0: Reading grid_nml0: Reading ice_nml0: Reading tracer_nml0:CalcWorkPerBlock: Total blocks:    50 Ice blocks:    50 IceFree blocks:     0 Land blocks:     00:MCT::m_Router::initp_: GSMap indices not increasing...Will correct0:MCT::m_Router::initp_: RGSMap indices not increasing...Will correct0:MCT::m_Router::initp_: RGSMap indices not increasing...Will correct0:MCT::m_Router::initp_: GSMap indices not increasing...Will correct14: (domain_check_grid_mct) ERROR: incompatible domain grid coordinates15: (domain_check_grid_mct) ERROR: incompatible domain grid coordinates36: (domain_check_grid_mct) ERROR: incompatible domain grid coordinates16: (domain_check_grid_mct) ERROR: incompatible domain grid coordinates.....70:(domain_check_grid_mct) : incompatible domain grid coordinates55:(domain_check_grid_mct) : incompatible domain grid coordinates40:028.MCT(MPEU)::die.: from (domain_check_grid_mct)()42:02A.MCT(MPEU)::die.: from (domain_check_grid_mct)()17:(domain_check_grid_mct) : incompatible domain grid coordinates17:011.MCT(MPEU)::die.: from (domain_check_grid_mct)()39:027.MCT(MPEU)::die.: from (domain_check_grid_mct)()55:037.MCT(MPEU)::die.: from (domain_check_grid_mct)()32:(domain_check_grid_mct) : incompatible domain grid coordinates32:020.MCT(MPEU)::die.: from (domain_check_grid_mct)()56:038.MCT(MPEU)::die.: from (domain_check_grid_mct)()31:MPT ERROR: Rank 31(g:31) is aborting with error code 2.31:     Process ID: 11561, Host: r1i2n13, Program: /glade2/scratch2/shinengh/testCTL2_so_T31_g37/bld/testCTL2_so_T31_g37.ccsm.exe31:     MPT Version: SGI MPT 2.15  09/03/16 04:15:5431:31:MPT: --------stack traceback-------59:03B.MCT(MPEU)::die.: from (domain_check_grid_mct)()60:03C.MCT(MPEU)::die.: from (domain_check_grid_mct)()61:03D.MCT(MPEU)::die.: from (domain_check_grid_mct)()62:03E.MCT(MPEU)::die.: from (domain_check_grid_mct)()69:045.MCT(MPEU)::die.: from (domain_check_grid_mct)()43:02B.MCT(MPEU)::die.: from (domain_check_grid_mct)()40:MPT ERROR: Rank 40(g:40) is aborting with error code 2.40:     Process ID: 62129, Host: r1i2n18, Program: /glade2/scratch2/shinengh/testCTL2_so_T31_g37/bld/testCTL2_so_T31_g37.ccsm.exe40:     MPT Version: SGI MPT 2.15  09/03/16 04:15:5440:40:MPT: --------stack traceback-------45:(domain_check_grid_mct) : incompatible domain grid coordinates45:02D.MCT(MPEU)::die.: from (domain_check_grid_mct)()27:(domain_check_grid_mct) : incompatible domain grid coordinates27:01B.MCT(MPEU)::die.: from (domain_check_grid_mct)()33:(domain_check_grid_mct) : incompatible domain grid coordinates33:021.MCT(MPEU)::die.: from (domain_check_grid_mct)()=========Shineng