How to run F compset using 2023 intial data

[fybi]

Dear All,

I am trying to run CESM "F_2000_SPCAM_m2005_" compset using 2023 meterology field as the initial data and 2023 SST inputdata, but I find that the MERRA2 data that has been interpolated to the CESM requirement format is only up to 2022-12-31. I have tried to interpolate the MERRA2 forcing data horizontally and vertically, but I encountered an error:
pio_support::pio_die:: myrank= -1 : ERROR: nf_mod.F90: 1224 :
NetCDF: Invalid dimension ID or name
forrtl: severe (174): SIGSEGV, segmentation fault occurred
The cesm.log is attached and below is the intial data info：
ncdump -h MERRA2_0.23x0.31_20230726.nc
netcdf MERRA2_0.23x0.31_20230726 {
dimensions:
time = 8 ;
lev = 30 ;
lat = 768 ;
lon = 1152 ;
variables:
float U(time, lev, lat, lon) ;
U:_FillValue = 9.96921e+36f ;
U:long_name = "eastward_wind" ;
U:units = "m s-1" ;
U:info = "NCL function hyi2hyo used to interpolate" ;
float time(time) ;
time:standard_name = "time" ;
time:long_name = "time" ;
time:units = "days since 1970-01-01 00:00:00" ;
time:calendar = "gregorian" ;
time:axis = "T" ;
double lev(lev) ;
lev:positive = "down" ;
lev:standard_name = "atmosphere_hybrid_sigma_pressure_coordinate" ;
lev:formula_terms = "a: hyam b: hybm p0: P0 ps: PS" ;
lev:units = "level" ;
lev:long_name = "hybrid level at midpoints (1000*(A+B))" ;
double lat(lat) ;
lat:standard_name = "latitude" ;
lat:long_name = "latitude" ;
lat:units = "degrees_north" ;
lat:axis = "Y" ;
double lon(lon) ;
lon:standard_name = "longitude" ;
lon:long_name = "longitude" ;
lon:units = "degrees_east" ;
lon:axis = "X" ;
float V(time, lev, lat, lon) ;
V:_FillValue = 9.96921e+36f ;
V:long_name = "northward_wind" ;
V:units = "m s-1" ;
V:info = "NCL function hyi2hyo used to interpolate" ;
float T(time, lev, lat, lon) ;
T:_FillValue = 9.96921e+36f ;
T:long_name = "air_temperature" ;
T:units = "K" ;
T:info = "NCL function hyi2hyo used to interpolate" ;
float Q(time, lev, lat, lon) ;
Q:_FillValue = 9.96921e+36f ;
Q:long_name = "specific_humidity" ;
Q:units = "kg kg-1" ;
Q:info = "NCL function hyi2hyo used to interpolate" ;
int date(time) ;
date:long_name = "current date (YYYYMMDD)" ;
int datesec(time) ;
datesec:long_name = "current seconds of current date" ;
float PS(time, lat, lon) ;
PS:long_name = "surface_pressure" ;
PS:units = "Pa" ;
float PHIS(time, lat, lon) ;
PHIS:long_name = "surface geopotential height" ;
PHIS:units = "m+2 s-2" ;
}

Can anyone please suggest me how can I solve this problem? Or is there any tutorial available explaining the CESM inputdata format.

Thanks!

 

[peverley]

Hi,

Can you clarify how you are interpreting the MERRA2 forcing data?

Thanks,
Courtney

 

[fybi]

Hi,

Can you clarify how you are interpreting the MERRA2 forcing data?

Thanks,
Courtney

Thank you very much for your reply!
Do you mean how I interpolated the MERRA2 forcing data？I used cdo to interpolate the MERRA2 forcing data (from https://rda.ucar.edu/datasets/ds313.3/) horizontally, and the vertical interpolation used the function hyi2hyo_Wrap.

 

[peverley]

Yes, thank you for translating my typo!

Can you reattach your log files (cesm and atm) so I can see more of what's happening?

Thanks!
Courtney

 

[fybi]

Of course, the two log files are attached here. Thanks so much and hope to get your guidance!

Best

 

[minimax]

Of course, the two log files are attached here. Thanks so much and hope to get your guidance!

Best

it seems that you have a problem with the following file, check it
/data/linguangxing/CESM/inputdata/atm/cam/topo/USGS_gtopo30_0.23x0.31_remap_c06
1107.nc
( NetCDF: Invalid dimension ID or name)

 

[minimax]

try to interpolate using cdo:
cdo griddes <original merra file> > grid_merra.txt
cdo remapbil,grid_merra.txt <merra original file> <new merra file>
may be it helps you.

 

[fybi]

try to interpolate using cdo:
cdo griddes <original merra file> > grid_merra.txt
cdo remapbil,grid_merra.txt <merra original file> <new merra file>
may be it helps you.

Dear minimax，
Thank you for your helpful advise. I have reinterpolated the original MERRA2 data, and it seems to have worked. But I met a new error:
" ENDRUN: called without a message string
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 1"

Do you have any suggestions for this question? and I modified the inputdata by changing the ncdata and bnd_topo in user_nl_cam, is it right or do I need to change something else? The interpolated data information and files are attached here.
Thank you very much!