get_whitelight_image Subroutine

public subroutine get_whitelight_image(qunit, fl)

Uses

  • proc~~get_whitelight_image~~UsesGraph proc~get_whitelight_image get_whitelight_image module~mod_global_parameters mod_global_parameters proc~get_whitelight_image->module~mod_global_parameters module~mod_basic_types mod_basic_types module~mod_global_parameters->module~mod_basic_types module~mod_connectivity mod_connectivity module~mod_global_parameters->module~mod_connectivity module~mod_constants mod_constants module~mod_global_parameters->module~mod_constants module~mod_physicaldata mod_physicaldata module~mod_global_parameters->module~mod_physicaldata module~mod_variables mod_variables module~mod_global_parameters->module~mod_variables mpi mpi module~mod_global_parameters->mpi module~mod_variables->module~mod_basic_types

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: qunit
type(te_fluid), intent(in) :: fl

Calls

proc~~get_whitelight_image~~CallsGraph proc~get_whitelight_image get_whitelight_image proc~get_image get_image proc~get_whitelight_image->proc~get_image proc~mpistop mpistop proc~get_whitelight_image->proc~mpistop proc~get_image->proc~mpistop mpi_allreduce mpi_allreduce proc~get_image->mpi_allreduce proc~get_cor_image get_cor_image proc~get_image->proc~get_cor_image proc~get_line_info get_line_info proc~get_image->proc~get_line_info proc~init_vectors_cartesian init_vectors_cartesian proc~get_image->proc~init_vectors_cartesian proc~init_vectors_spherical init_vectors_spherical proc~get_image->proc~init_vectors_spherical proc~integrate_emission_cartesian integrate_emission_cartesian proc~get_image->proc~integrate_emission_cartesian proc~integrate_emission_spherical integrate_emission_spherical proc~get_image->proc~integrate_emission_spherical proc~integrate_whitelight_spherical integrate_whitelight_spherical proc~get_image->proc~integrate_whitelight_spherical proc~output_data output_data proc~get_image->proc~output_data mpi_abort mpi_abort proc~mpistop->mpi_abort proc~dot_product_loc dot_product_loc proc~get_cor_image->proc~dot_product_loc proc~get_line_info->proc~mpistop proc~cross_product_loc cross_product_loc proc~init_vectors_cartesian->proc~cross_product_loc proc~cartesian_to_spherical cartesian_to_spherical proc~init_vectors_spherical->proc~cartesian_to_spherical proc~init_vectors_spherical->proc~cross_product_loc proc~spherical_to_cartesian spherical_to_cartesian proc~init_vectors_spherical->proc~spherical_to_cartesian proc~integrate_emission_cartesian->proc~get_cor_image proc~integrate_emission_cartesian->proc~get_line_info proc~get_euv get_EUV proc~integrate_emission_cartesian->proc~get_euv proc~get_sxr get_SXR proc~integrate_emission_cartesian->proc~get_sxr proc~integrate_emission_spherical->proc~get_line_info dsin dsin proc~integrate_emission_spherical->dsin proc~integrate_emission_spherical->proc~dot_product_loc proc~get_cor_image_spherical get_cor_image_spherical proc~integrate_emission_spherical->proc~get_cor_image_spherical proc~integrate_emission_spherical->proc~get_euv proc~integrate_emission_spherical->proc~get_sxr proc~get_unit_vector_spherical get_unit_vector_spherical proc~integrate_emission_spherical->proc~get_unit_vector_spherical proc~integrate_emission_spherical->proc~spherical_to_cartesian proc~integrate_whitelight_spherical->dsin proc~integrate_whitelight_spherical->proc~dot_product_loc proc~integrate_whitelight_spherical->proc~get_cor_image_spherical proc~get_thomson_parameters get_Thomson_parameters proc~integrate_whitelight_spherical->proc~get_thomson_parameters proc~integrate_whitelight_spherical->proc~get_unit_vector_spherical proc~get_whitelight_thomson get_whitelight_Thomson proc~integrate_whitelight_spherical->proc~get_whitelight_thomson proc~output_data->proc~mpistop proc~write_image_vticc write_image_vtiCC proc~output_data->proc~write_image_vticc proc~write_image_vtucc write_image_vtuCC proc~output_data->proc~write_image_vtucc dacos dacos proc~cartesian_to_spherical->dacos proc~get_cor_image_spherical->proc~dot_product_loc proc~get_cor_image_spherical->proc~spherical_to_cartesian proc~get_euv->proc~mpistop dexp dexp proc~get_sxr->dexp proc~get_unit_vector_spherical->dsin dcos dcos proc~get_unit_vector_spherical->dcos proc~spherical_to_cartesian->dsin proc~spherical_to_cartesian->dcos proc~write_image_vticc->proc~get_line_info proc~write_image_vtucc->proc~get_line_info

Source Code

    subroutine get_whitelight_image(qunit,fl)
      use mod_global_parameters

      integer, intent(in) :: qunit
      type(te_fluid), intent(in) :: fl
      character(20) :: datatype
      double precision :: LASCO_rsl

      if (mype==0) print *,&
          '###################################################'

      if (whitelight_instrument=='LASCO/C1') then
        LASCO_rsl=5.6d0/instrument_resolution_factor
        if (mype==0) print *,&
            'Systhesizing white light image (observed by LASCO/C1).'
      else if (whitelight_instrument=='LASCO/C2') then
        LASCO_rsl=11.4d0/instrument_resolution_factor
        if (mype==0) print *,&
            'Systhesizing white light image (observed by LASCO/C2).'
      else if (whitelight_instrument=='LASCO/C3') then
        LASCO_rsl=56.d0/instrument_resolution_factor
        if (mype==0) print *,&
            'Systhesizing white light image (observed by LASCO/C3).'
      else
        call MPISTOP('Whitelight synthesis: instrument is not supported!')
      endif

      if (mype==0) write(*,'(a,f5.1,a,f5.1,a,f5.1,a,f5.1,a)') ' Pixel: ',&
         LASCO_rsl*0.725,' Mm x ',LASCO_rsl*0.725, ' Mm  (', LASCO_rsl,&
          ' arcsec x ', LASCO_rsl, ' arcsec) '
      if (mype==0) print *, 'Unit of white light flux: average Sun brightness'

      datatype='image_whitelight'

      if (mype==0) then
        if (activate_unit_arcsec) then
          print *, 'Unit of length: arcsec (~725 km)'
        else
          if (SI_unit) then
            if (mype==0) write(*,'(a,f8.1,a)') ' Unit of length: ',&
               unit_length/1.d6,' Mm'
          else
            if (mype==0) write(*,'(a,f8.1,a)') ' Unit of length: ',&
               unit_length/1.d8,' Mm'
          endif
        endif
      endif

      if (coordinate==spherical) then
        if (mype==0) write(*,'(a,f6.3,f8.3,f8.3,a)') &
' Mapping: R=0 of the simulation box is located at [X=0,Y=0] of the image'
        call get_image(qunit,datatype,fl)
      else
        call MPISTOP("Whitelight synthesis: this coordinate is not supported!")

      endif

      if (mype==0) print *,&
          '###################################################'

    end subroutine get_whitelight_image