Sort out tracer fields - Githubissues

hinnes97 / soc-rad-conv-fort

Fortran radiative-convective equilibrium code using met office SOCRATES radiation code

2 stars 0 forks source link

Sort out tracer fields #2

Closed hinnes97 closed 5 months ago

hinnes97 commented 1 year ago

Currently can handle 1 dry component and 1 condensing component (hard-coded as water). Should generalise so we can have multiple components that can vary. Need to keep an eye on how Socrates would handle the density and the gas constant

hinnes97 commented 1 year ago

Thoughts on how to do this:

Have a list of chem species with unique numbers (HITRAN index?)
Have a number, nq, of chemical species. q vector is length nq
At the beginning, create a vector length nq filled with the unique index of the chemical species
Can then loop over tracers doing and access quantity from phys file by doing

cp = 0. do i=1,nq call lookup_cp(tracer_index_list(i), cp_out) cp = cp + q(i) * cp_out enddo

hinnes97 commented 5 months ago

Now done with a limited number of tracers. Reads in dry composition from file and then adjusts the q of these species accordingly as the water vapour is changed due to saturation etc.