TUV-x photolysis #1471
TUV-x photolysis #1471
Conversation
modified: .gitmodules modified: cime_config/buildlib modified: libraries/tuv-x (new commits)
| call mee_ion_final() | ||
| call rate_diags_final() | ||
| call species_sums_final() | ||
| call tuvx_finalize() |
There was a problem hiding this comment.
Could maybe wrap this within if (tuvx_active) to be consistent with the other calls
| !----------------------------------------------------------------- | ||
| ! ... get calculated photolysis rates from TUV-x | ||
| !----------------------------------------------------------------- | ||
| call tuvx_get_photo_rates( state, pbuf, ncol, lchnk, zmid, zint, & |
There was a problem hiding this comment.
It looks like the only place tuvx_get_photo_rates uses state is to get state%ncol in line ~679, but ncol is already passed in here. Maybe state does not have to be passed in to call tuvx_get_photo_rates
| subroutine gas_phase_chemdr(lchnk, ncol, imozart, q, & | ||
| phis, zm, zi, calday, & | ||
| tfld, pmid, pdel, pint, rpdel, rpdeldry, & | ||
| cldw, troplev, troplevchem, & | ||
| ncldwtr, ufld, vfld, & | ||
| delt, ps, & | ||
| fsds, ts, asdir, ocnfrac, icefrac, & | ||
| precc, precl, snowhland, ghg_chem, latmapback, & | ||
| drydepflx, wetdepflx, cflx, fire_sflx, fire_ztop, nhx_nitrogen_flx, noy_nitrogen_flx, & | ||
| use_hemco, qtend, pbuf) | ||
| drydepflx, wetdepflx, cflx, fire_sflx, fire_ztop, & | ||
| nhx_nitrogen_flx, noy_nitrogen_flx, & | ||
| use_hemco, qtend, pbuf, state) |
There was a problem hiding this comment.
These changes appear to be unnecessary as state was already previously passed in as the first argument, and was moved to the last one for some reason.
| use mo_chm_diags, only : chm_diags, het_diags | ||
| use perf_mod, only : t_startf, t_stopf | ||
| use gas_wetdep_opts, only : gas_wetdep_method | ||
| use physics_types, only : physics_state |
There was a problem hiding this comment.
This line is a duplicate of line 318 and does not need to be added again
| call gas_phase_chemdr(lchnk, ncol, imozart, state%q, & | ||
| state%phis, state%zm, state%zi, calday, & | ||
| state%t, state%pmid, state%pdel, state%pint, state%rpdel, state%rpdeldry, & | ||
| cldw, tropLev, tropLevChem, ncldwtr, state%u, state%v, chem_dt, state%ps, & | ||
| fsds, cam_in%ts, cam_in%asdir, cam_in%ocnfrac, cam_in%icefrac, & | ||
| cam_out%precc, cam_out%precl, cam_in%snowhland, ghg_chem, state%latmapback, & | ||
| drydepflx, wetdepflx, cam_in%cflx, cam_in%fireflx, cam_in%fireztop, & | ||
| nhx_nitrogen_flx, noy_nitrogen_flx, use_hemco, ptend%q, pbuf ) | ||
| nhx_nitrogen_flx, noy_nitrogen_flx, use_hemco, ptend%q, pbuf, state ) |
There was a problem hiding this comment.
These changes appear to be unnecessary as state was already previously passed in as the first argument, and was moved to the last one for some reason.
| integer, parameter :: PROFILE_INDEX_SO2 = 7 ! Sulfur dioxide profile index | ||
| integer, parameter :: PROFILE_INDEX_NO2 = 8 ! Nitrogen dioxide profile index |
There was a problem hiding this comment.
These appear to be non-public and unused, can I ask if these are intended to be not implemented?
| ! Definition of the MS93 wavelength grid TODO add description of this | ||
| integer, parameter :: NUM_BINS_MS93 = 4 | ||
| real(kind=r8), parameter :: WAVELENGTH_EDGES_MS93(NUM_BINS_MS93+1) = & | ||
| (/ 181.6_r8, 183.1_r8, 184.6_r8, 190.2_r8, 192.5_r8 /) |
There was a problem hiding this comment.
Missed TODO here about the description here
| deallocate( cam_grids ) | ||
| deallocate( cam_profiles ) | ||
| deallocate( cam_radiators ) | ||
| deallocate( wavelength ) |
There was a problem hiding this comment.
Since wavelength => cam_grids%get_grid(...) I would suggest switching the deallocate order for safety
| deallocate( cam_grids ) | |
| deallocate( cam_profiles ) | |
| deallocate( cam_radiators ) | |
| deallocate( wavelength ) | |
| deallocate( wavelength ) | |
| deallocate( cam_grids ) | |
| deallocate( cam_profiles ) | |
| deallocate( cam_radiators ) |
|
|
||
| host_profile => profiles%get_profile( "temperature", "K" ) | ||
| this%profiles_( PROFILE_INDEX_TEMPERATURE ) = this%core_%get_updater( host_profile, found ) | ||
| call assert( 418735162, found ) |
There was a problem hiding this comment.
Here and throughout the file, there are many calls to musica_assert::assert with seemingly random numbers for unique error codes. It seems to be consistent with what's done elsewhere in MUSICA, but I would request to at least have a comment somewhere that these numbers are just unique error codes and don't have any particular meaning other than being unique, and this is how the MUSICA library handles errors. Thanks!
| ! || || | ||
| ! || | ||
| ! ----------------- i_int = pver | ||
| ! ================= i_imd = pver ------------------ i_int = 2 |
There was a problem hiding this comment.
typo here:
| ! ================= i_imd = pver ------------------ i_int = 2 | |
| ! ================= i_mid = pver ------------------ i_int = 2 |
2 participants
This provides the option to use the TUV-x library to compute photo-chemical rates in-line as an alternative to the traditional table look up method.
Closes #757