jimmielin
commented
1 month ago Updated history field list with new CCPP field names:
***************** HISTORY FIELD LIST ******************
TROP_P Pa avg tropopause_air_pressure
TROP_T K avg tropopause_air_temperature
TROP_Z m avg tropopause_geopotential_height_wrt_surface
TROP_DZ m avg geopotential_height_difference_between_atmosphere_layer_and_tropopause
TROP_PD probability avg probability_distribution_of_tropopause_vertical_layer_index
TROP_FD probability avg tropopause_found
TROPP_P Pa avg tropopause_air_pressure_from_lapse_rate_method
TROPP_T K avg tropopause_air_temperature_from_lapse_rate_method
TROPP_Z m avg tropopause_geopotential_height_wrt_surface_from_lapse_rate_method
TROPP_DZ m avg geopotential_height_difference_between_atmosphere_layer_and_tropopause_from_lapse_rate_method
TROPP_PD probability avg probability_distribution_of_tropopause_vertical_layer_index_from_lapse_rate_method
TROPP_FD probability avg tropopause_found_from_lapse_rate_method
TROPF_P Pa avg tropopause_air_pressure_from_cold_point_method
TROPF_T K avg tropopause_air_temperature_from_cold_point_method
TROPF_Z m avg tropopause_geopotential_height_wrt_surface_from_cold_point_method
TROPF_DZ m avg geopotential_height_difference_between_atmosphere_layer_and_tropopause_from_cold_point_method
TROPF_PD probability avg probability_distribution_of_tropopause_vertical_layer_index_from_cold_point_method
TROPF_FD probability avg tropopause_found_from_cold_point_method
TROPC_P Pa avg tropopause_air_pressure_from_climatological_method
TROPC_T K avg tropopause_air_temperature_from_climatological_method
TROPC_Z m avg tropopause_geopotential_height_wrt_surface_from_climatological_method
TROPC_DZ m avg geopotential_height_difference_between_atmosphere_layer_and_tropopause_from_climatological_method
TROPC_PD probability avg probability_distribution_of_tropopause_vertical_layer_index_from_climatological_method
TROPC_FD probability avg tropopause_found_from_climatological_method
hstobie_trop fraction of ins vertical_layer_index_lower_bound_from_hybrid_stobie_linoz_with_climatological_backup_method_for_stratospheric_chemistry
hstobie_linoz fraction of ins vertical_layer_index_lower_bound_from_hybrid_stobie_linoz_with_climatological_backup_method_for_linearized_ozone_chemistry
hstobie_tropop fraction of ins tropopause_vertical_layer_index_from_hybrid_stobie_linoz_with_climatological_backup_method_for_chemistry
*************** END HISTORY FIELD LIST ****************
Implements CCPP-compliant tropopause_find.
Major changes to code
This is my first time doing this, so please feel free to suggest alternative approaches!
Standard CCPP-ization procedures, metadata, explicit use of fields instead of physics state, only passing active
:ncolfields, ... was made.tropopause_findaccepted "primary" and "backup" methods for finding the tropopause. All the methods known to be used in CAM are now calculated at the same time in thetropopause_find_runmain driver routine; it calls the underlying logic as appropriate with and populates the appropriate standard name physics fields. Because the same physical quantities for tropopause level, height, pressure, temperature, etc. using different methods are used throughout the CAM physics and simultaneously output in the history tape, the same standard names with a suffix are used to differentiate between these quantities computed by different methods. e.g.,The naming should probably be changed! I'm open to suggestions.
tropopause_read_file, which provides climatological data used as the default "backup" method to compute the tropopause (tropp_p_locandtropp_days), has been moved to a utility module in CAM-SIMA to read, regrid, and provide this climatology to the CCPP-ized physics; it also removeslchnkindexing since they're no longer used in CAM-SIMA.tropopause_find used to accept (
optional) arguments for returningtropT,tropZ,tropP(tropopause temperature, height, and pressure) which are no longer optional per the CAM-SIMA port guide.A "shim" in current CAM's
tropopause.F90was made to use the CCPP-ized routines, but preserve all existing behavior and functionality (b4b) as expected in current CAM. In CAM-SIMA all of this are not needed.Smaller changes
get_curr_calday()is now provided in the registry asfractional_calendar_days_on_end_of_current_timestepOutstanding issues/notes
Tropopause_find cannot be tested using snapshots within CAM-SIMA: In regular CAM,
tropopause_finditself does not update physics state/tendencies; its use in regular CAM is just providing the tropopause level, height, pressure, and temperature to the calling physics package; a final run oftropopause_outputwill then write out to history the computed tropopause properties using the primary method. This means if snapshots cannot verify whattropopause_findis doing within CAM-SIMA.How to make standard names for the same physical quantities computed using different methods that have to co-exist in one model run
Answer check
For existing CAM
QPC6(which uses CLUBB that usestropopause_find) answers match bit-for-bit for history (h0i) and snapshot files before and afterclubb_tend_cam. I can run more detailed regression tests for the CAM tag.For CAM-SIMA
Tricky way to verify answers
We can validate the answer fields of
TROP_P,TROP_T, andTROP_Zcan be against the output of CAM by running CAM-SIMA from a CAM snapshot and using history to write out these fields (instantaneously) at the first timestep. Note thatndenshas to be set to1in CAM (fordouble) output, and auser_setsnapshot should be taken for use beforetropopause_outputinphyspkg.F90for use inncdata.