Different Climatology in NCAR's AMIP and My Own

[FuchengY]

What version of the code are you using?
CESM 2.1.5


Have you made any changes to files in the source tree?
No


Describe every step you took leading up to the problem:
Create new case with f09_f09_mg17 & FHIST_BGC, and run for 10 years on Derecho.

In CaseStatus:

2024-10-04 09:22:30: xmlchange success <command> ./xmlchange STOP_N=2 </command>
---------------------------------------------------
2024-10-04 09:22:48: xmlchange success <command> ./xmlchange RESUBMIT=4 </command>
---------------------------------------------------
2024-10-04 09:25:49: case.build starting
---------------------------------------------------
CESM version is release-cesm2.1.5
Processing externals description file : Externals.cfg (/glade/u/home/fuchengy/my_cesm_code)
Processing externals description file : Externals_CAM.cfg (/glade/u/home/fuchengy/my_cesm_code/components/cam)
Processing externals description file : Externals_CISM.cfg (/glade/u/home/fuchengy/my_cesm_code/components/cism)
Processing externals description file : Externals_CLM.cfg (/glade/u/home/fuchengy/my_cesm_code/components/clm)
Processing externals description file : Externals_POP.cfg (/glade/u/home/fuchengy/my_cesm_code/components/pop)
Checking local status of required & optional components: cam, chem_proc, carma, clubb, cosp2, cice, cime, cism, source_cism, clm, fates, mosart, pop, cvmix, marbl, rtm, ww3,
./cime
clean sandbox, on cime5.6.49
./components/cam
clean sandbox, on cam_cesm2_1_rel_60
./components/cam/chem_proc
clean sandbox, on tools/proc_atm/chem_proc/release_tags/chem_proc5_0_03_rel
./components/cam/src/physics/carma/base
clean sandbox, on carma/release_tags/carma3_49_rel
./components/cam/src/physics/clubb
clean sandbox, on vendor_clubb_r8099_n03
./components/cam/src/physics/cosp2/src
clean sandbox, on v2.1.4cesm
./components/cice
clean sandbox, on cice5_cesm2_1_1_20231220
./components/cism
clean sandbox, on cism-release-cesm2.1.2_04
./components/cism/source_cism
clean sandbox, on release-cism2.1.04
./components/clm
clean sandbox, on release-clm5.0.37
./components/clm/src/fates
clean sandbox, on sci.1.30.0_api.8.0.0
./components/mosart
clean sandbox, on release-cesm2.0.04
./components/pop
clean sandbox, on pop2_cesm2_1_rel_n15
./components/pop/externals/CVMix
clean sandbox, on v0.93-beta
./components/pop/externals/MARBL
clean sandbox, on cesm2.1-n00
./components/rtm
clean sandbox, on release-cesm2.0.04
./components/ww3
clean sandbox, on ww3_181001
2024-10-04 09:37:17: case.build success
---------------------------------------------------
2024-10-04 09:37:17: case.submit starting
---------------------------------------------------
2024-10-04 09:37:30: case.submit success case.run:6269912.desched1, case.st_archive:6269913.desched1


Describe your problem or question:

I created an AMIP experiment and run for 10 years. The output shows that the precipitation (PRECT) is quite low compared with AMIP simulation on NCAR's server (/glade/campaign/collections/cmip/CMIP6/timeseries-cmip6/f.e21.FHIST_BGC.f09_f09_mg17.CMIP6-AMIP.001/atm/proc/tseries/month_1). The figure shows the precipitation difference (My Own - NCAR) during 1979-1985 (unit: mm/day). Is it a problem of my case? Or can I expect a similar climatology?

 

[fischer]

Hi Fucheng,

Can you point me to your case directory? Looking at your CaseStatus, it looks like you only ran for 10 days instead of 10 years.

Chris

 

[FuchengY]

Hi Fucheng,

Can you point me to your case directory? Looking at your CaseStatus, it looks like you only ran for 10 days instead of 10 years.

Chris

Hi Chris,

I did set STOP_OPTION='nyears'. Please see the attached file. If it is a result for 10 days, we will see a lot of weather patterns.

The NCAR's AMIP simulations starts from 1950 but I run from 1979, so is there a sense of spin-up for AMIP? I don't think atmopshere needs a long time to adjust.

Best

Fucheng

 

[FuchengY]

Hi Fucheng,

Can you point me to your case directory? Looking at your CaseStatus, it looks like you only ran for 10 days instead of 10 years.

Chris

My namelist is like:

avgflag_pertape = 'A','I'
nhtfrq = 0, -3
mfilt = 12, 1
empty_htapes=.True.
fincl1='PRECT'
fincl2='U'

Although mfilt was set to 12, the files were still output every month.

 

[fischer]

Hi Fucheng,

I'm going to move this to the CAM forum, were someone more familiar with the atmosphere model can help you.

Chris

 

[islas]

It may be helpful to also point us toward the script you're using to make the plot. This difference in precipitation is larger than anything I would expect. Is it possible there's a plotting issue? Anyway, pointing to the script you're using for plotting if it's on casper/derecho might help.

 

[FuchengY]

It may be helpful to also point us toward the script you're using to make the plot. This difference in precipitation is larger than anything I would expect. Is it possible there's a plotting issue? Anyway, pointing to the script you're using for plotting if it's on casper/derecho might help.

I did not plot on server, but I just moved all files needed to my directory: /glade/derecho/scratch/fuchengy/test. If you are interested in, please feel free to have a look at.

 

[FuchengY]

It may be helpful to also point us toward the script you're using to make the plot. This difference in precipitation is larger than anything I would expect. Is it possible there's a plotting issue? Anyway, pointing to the script you're using for plotting if it's on casper/derecho might help.

I just figured out. I did made a mistake on calculation. The climatological precipitation is similar in NCAR's AMIP and my own's.