Hi,
I am running CESM2.1.2 with slab ocean aquaplanet simulation (Qflx is zero) with the seasonal cycle. I am encountering an asymmetry in the north-south SOLIN with a value of 0.74w/m^2 increase in the southern hemisphere. This is very similar to this queries North ITCZ. I have used different compiler gcc/7.5m,gcc/9 and gcc/12 to compile the model all of them gives exactly same asymmetry. I also tried with cesm 2.1.0. Does anyone have a comment on this please let me know?
The changes I done are :
./create_newcase --case Testing_SOLIN_intel --compset QSC6 --res f19_f19_mg17 --compiler intel --mpilib impi
I introduced seasonality by including orbital parameters in user_nl_cpl:
orb_eccen = 0.
orb_obliq = 23.5
orb_mvelp = 0.
orb_mode = "fixed_parameters"
And commented these two lines in the following block in seq_infodata_mod.F90 file and put it into SourceMods/src.drv/:
if ( infodata%aqua_planet ) then
infodata%aqua_planet_sst = 1
!infodata%perpetual = .true.
!infodata%perpetual_ymd = aqua_perpetual_ymd
endif
For sanity check. I made the obliquity to -23.5 which make the asymmetry 0.74w/m^2 increase to northern hemisphere.
If anyone has encountered this issue, any suggestions or resolutions would be greatly appreciated.
I am moving this to the simpler model forum. I really have not idea on this one.