Scheduled Downtime
On Tuesday 24 October 2023 @ 5pm MT the forums will be in read only mode in preparation for the downtime. On Wednesday 25 October 2023 @ 5am MT, this website will be down for maintenance and expected to return online later in the morning.
Normal Operations
The forums are back online with normal operations. If you notice any issues or errors related to the forums, please reach out to help@ucar.edu

Energy Conservation Within CAM

I am using CAM3.0 coupled with CLM and the Slab Ocean Model. I have added a surface flux parameterization and want to check to make sure the energy within the model is still being conserved. Is there a diagnostic already output that I can use? I have noticed that ENGYCORR is a possibility.... :?
 

pjr

Member
scapps said:
I am using CAM3.0 coupled with CLM and the Slab Ocean Model. I have added a surface flux parameterization and want to check to make sure the energy within the model is still being conserved. Is there a diagnostic already output that I can use? I have noticed that ENGYCORR is a possibility.... :?

ENGYCORR is a diagnostic of the tendencies required by some dynamical cores (the eulerian and SLD cores) to conserve energy. The diagnostic is only relevant to the dynamics, not to a physically based flux parameterization.

Which fluxes are you modifying. If you are changing the heat, moisture and momentum fluxes you should:
1) look at those fluxes (called SHFLX, LHFLX, TAUX, TAUY) on the history file
2) you can also run the standard CAM diagnostic package and look at the set 1 diagnostics. Once the model has reached an equilibrium you should see the numbers for the first two rows (labelled RESTOM, and RESTOA) agree to a few percent for a five year run or longer. You cannot expect better agreement than that because the atmosphere and slab ocean are always slightly out of equilibrium with each other.

Note that you should first do this in a standalong CAM run (no SOM). Without knowing exactly what you are doing and why I am hesitant to provide more advice.

Phil
 
Thank you for your response. I am modifying sensible, latent and momentum fluxes. Using a probability density function to represent subgrid-scale windspeeds, the fluxes are weighted according to the probability of that particular bin wind speed. I have already plotted a time series of TOM net solar minus TOM net longwave and noticed no trend throughout the run period....
 
Top