Main menu

Navigation

Slab ocean gx3v7 grid issue

6 posts / 0 new
Last post
shineng.hu@...
Slab ocean gx3v7 grid issue

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@...

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.nc

domain1 = /glade/p/cesm/cseg/inputdata/ocn/docn7/SOM/pop_frc.gx3v7.110128.nc

file1 = /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@...

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 =      1

48:(shr_stream_set) filename = /glade/p/cesm/cseg/inputdata/lnd/clm2/ndepdata/fndep_clm_hist_simyr1849-2006_1.9x2.5_c100428.nc

59:(shr_stream_set) size of filename =      1

59:(shr_stream_set) filename = /glade/p/cesm/cseg/inputdata/lnd/clm2/ndepdata/fndep_clm_hist_simyr1849-2006_1.9x2.5_c100428.nc

0: Reading setup_nml

0: Reading grid_nml

0: Reading ice_nml

0: Reading tracer_nml

0:CalcWorkPerBlock: Total blocks:    50 Ice blocks:    50 IceFree blocks:     0 Land blocks:     0

0:MCT::m_Router::initp_: GSMap indices not increasing...Will correct

0:MCT::m_Router::initp_: RGSMap indices not increasing...Will correct

0:MCT::m_Router::initp_: RGSMap indices not increasing...Will correct

0:MCT::m_Router::initp_: GSMap indices not increasing...Will correct

14: (domain_check_grid_mct) ERROR: incompatible domain grid coordinates

15: (domain_check_grid_mct) ERROR: incompatible domain grid coordinates

36: (domain_check_grid_mct) ERROR: incompatible domain grid coordinates

16: (domain_check_grid_mct) ERROR: incompatible domain grid coordinates

.....

70:(domain_check_grid_mct) : incompatible domain grid coordinates

55:(domain_check_grid_mct) : incompatible domain grid coordinates

40: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 coordinates

17: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 coordinates

32: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.exe

31:     MPT Version: SGI MPT 2.15  09/03/16 04:15:54

31:

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.exe

40:     MPT Version: SGI MPT 2.15  09/03/16 04:15:54

40:

40:MPT: --------stack traceback-------

45:(domain_check_grid_mct) : incompatible domain grid coordinates

45:02D.MCT(MPEU)::die.: from (domain_check_grid_mct)()

27:(domain_check_grid_mct) : incompatible domain grid coordinates

27:01B.MCT(MPEU)::die.: from (domain_check_grid_mct)()

33:(domain_check_grid_mct) : incompatible domain grid coordinates

33:021.MCT(MPEU)::die.: from (domain_check_grid_mct)()

=========

Shineng

 

Log in or register to post comments

Who's new

  • 20171204300@...
  • poornadurga.g@...
  • lina.boljka@...
  • nuistwangjing@...
  • vineetm@...