[mesa-users] Atmosphere structure missing in profiles

[Bill Paxton]

Hi Ehsan,

The profile is for information about the internal structure of the mesa model on a zone-by-zone basis.

The added atmosphere is NOT part of the mesa model and so info about it doesn't go in the profile.

If you'd like to output information about the added atmosphere, you can use this (from star_data).

      ! atmosphere structure information (only created on demand)
      integer :: atm_structure_num_pts
      real(dp), pointer :: atm_structure(:,:) ! (num_results_for_create_atm,atm_structure_num_pts)
         ! defined at points in atmosphere
         ! atm_structure(:,1) is base of atmosphere
         ! atm_structure(:,atm_structure_num_pts) is top of atmosphere

You can see this being used in star/private/pulsation_info.   E.g. to get grada at a point k in the atmosphere, 
	s% atm_structure(atm_grada,k)
For the radius, you need to add a delta to the model radius:
               s% r(1) + s% atm_structure(atm_delta_r,k)

The value "num_results_for_create_atm" is defined in atm/public/atm_def along with def's for what is saved:
      
      ! info about structure of atmosphere
      integer, parameter :: atm_xm = 1 ! mass of atm exterior to this point (g)
      integer, parameter :: atm_delta_r = atm_xm+1 ! radial distance above base of envelope (cm)
      integer, parameter :: atm_lnP = atm_delta_r+1
      integer, parameter :: atm_lnd = atm_lnP+1
      integer, parameter :: atm_lnT = atm_lnd+1
      integer, parameter :: atm_gradT = atm_lnT+1
      integer, parameter :: atm_kap = atm_gradT+1
      integer, parameter :: atm_gamma1 = atm_kap+1
      integer, parameter :: atm_grada = atm_gamma1+1
      integer, parameter :: atm_chiT = atm_grada+1
      integer, parameter :: atm_chiRho = atm_chiT+1
      integer, parameter :: atm_cv = atm_chiRho+1
      integer, parameter :: atm_cp = atm_cv+1
      integer, parameter :: atm_lnfree_e = atm_cp+1
      integer, parameter :: atm_dlnkap_dlnT = atm_lnfree_e+1
      integer, parameter :: atm_dlnkap_dlnd = atm_dlnkap_dlnT+1
      integer, parameter :: atm_lnPgas = atm_dlnkap_dlnd+1
      integer, parameter :: atm_tau = atm_lnPgas+1
      integer, parameter :: atm_gradr = atm_tau+1

      integer, parameter :: num_results_for_create_atm = atm_gradr 


If you have already created the atm_structure info for the pulsation analysis, then it should still be okay to use it without needing to call do_create_atm.

Good luck!

Bill




On Jan 22, 2014, at 8:38 AM, Ehsan Moravveji wrote:

> Hello all,
> I would like to store the atmosphere structure in the profile files as they are stored in pulsation files. Indeed, many quantities will be kept non-varying in the atmosphere, and some will. 
> I currently use version 5548, and would like to know if "there could be a patch" to take care of this. I have no idea if such a possible patch/update need be tailored to specific version, or can be made available for the most recent version (that I do not want to use!).
> 
> Kind regards.
> Ehsan.
> ------------------------------------------------------------------------------
> CenturyLink Cloud: The Leader in Enterprise Cloud Services.
> Learn Why More Businesses Are Choosing CenturyLink Cloud For
> Critical Workloads, Development Environments & Everything In Between.
> Get a Quote or Start a Free Trial Today. 
> http://pubads.g.doubleclick.net/gampad/clk?id=119420431&iu=/4140/ostg.clktrk
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users