Hi, I am currently using CESM2.1.5 to set up a HS94-Povalni-Kushner case. The final aim is to place a planetary-scale wave forcing at the bottom to generate sudden warmings in the stratosphere. However, after successfully setting up the original HS94 case and trying to modify the stratospheric settings following the guidance on the CESM website (Running with a different analytical relaxation temperature profile and damping parameters e.g., the Polvani and Kushner (2002) setup | Community Earth System Model), I encountered a problem possibly associated with the IC file. The commands I used are: The commands I used to set up the case are:
./create_newcase --case /data/jiachengye/CESM_exps/ssw1 --compset FTJ16 --res T42z30_T42_mg17 --compiler intel --run-unsupported
./xmlchange STOP_OPTION=ndays,STOP_N=20
./xmlchange --append --file env_build.xml --id CAM_CONFIG_OPTS --val="-nlev=40" (The default setting for PK02)
./case.setup
./case.build
./case.submit.
The content in the user_nl_cam is simply:"
&pkstrat_nl
pkstrat=.True.
vgamma=2.
/
ncdata='/data/jiachengye/CESM_new/dabiic/HS1994.ne16.L40.nc'
"
I placed all modified .f90 files (i.e., runtime_opts.F90, held_suarez.F90, held_suarez_cam.F90) in the directory /data/jiachengye/CESM_exps/ssw1/SourceMods/src.cam/. The namelist file namelist_definition.xml under the directory $CESM/components/cam/bld/namelist_files/ is also correspondingly updated to add the additional pkstrat statement before setting up the case. According to the error report ($OUTPUTDIR/dry_held_suarez/ssw1/run/cesm.log.250311-014419):
"
read_inidat: ERROR: model parameters do not match initial dataset parameters
Model Parameters: plon = 128 plat = 64
Dataset Parameters: dlon = 123435928 dlat = 0
ERROR:
read_inidat: ERROR: model parameters do not match initial dataset parameters
"
I suspect the problem is associated with the IC data file because the resolution (T42z30_T42_mg17) I chose to set up the case is not compatible with the default resolution of the IC data. In the ncl file (/dabiic/makeic.ncl) used to generate initial data, I set
res = "ne16",
which seems to be the error source. However, so far I don't know how to generate the corresponding initial data, as changing from 'ne16' to '64x128' appears to be prohibited. Any guidance would be much appreciated! I have attached the relevant files for reference. I apologize for this seemingly easy question.
./create_newcase --case /data/jiachengye/CESM_exps/ssw1 --compset FTJ16 --res T42z30_T42_mg17 --compiler intel --run-unsupported
./xmlchange STOP_OPTION=ndays,STOP_N=20
./xmlchange --append --file env_build.xml --id CAM_CONFIG_OPTS --val="-nlev=40" (The default setting for PK02)
./case.setup
./case.build
./case.submit.
The content in the user_nl_cam is simply:"
&pkstrat_nl
pkstrat=.True.
vgamma=2.
/
ncdata='/data/jiachengye/CESM_new/dabiic/HS1994.ne16.L40.nc'
"
I placed all modified .f90 files (i.e., runtime_opts.F90, held_suarez.F90, held_suarez_cam.F90) in the directory /data/jiachengye/CESM_exps/ssw1/SourceMods/src.cam/. The namelist file namelist_definition.xml under the directory $CESM/components/cam/bld/namelist_files/ is also correspondingly updated to add the additional pkstrat statement before setting up the case. According to the error report ($OUTPUTDIR/dry_held_suarez/ssw1/run/cesm.log.250311-014419):
"
read_inidat: ERROR: model parameters do not match initial dataset parameters
Model Parameters: plon = 128 plat = 64
Dataset Parameters: dlon = 123435928 dlat = 0
ERROR:
read_inidat: ERROR: model parameters do not match initial dataset parameters
"
I suspect the problem is associated with the IC data file because the resolution (T42z30_T42_mg17) I chose to set up the case is not compatible with the default resolution of the IC data. In the ncl file (/dabiic/makeic.ncl) used to generate initial data, I set
res = "ne16",
which seems to be the error source. However, so far I don't know how to generate the corresponding initial data, as changing from 'ne16' to '64x128' appears to be prohibited. Any guidance would be much appreciated! I have attached the relevant files for reference. I apologize for this seemingly easy question.