JSOC/proj/vfisv/apps/lim_init.f90 - view

Return to lim_init.f90 CVS log

Up to [Development] / JSOC / proj / vfisv / apps

File: [Development] / JSOC / proj / vfisv / apps / lim_init.f90 (download)
Revision: 1.1, Tue Apr 10 21:17:08 2012 UTC (11 years, 1 month ago) by keiji
Branch: MAIN
CVS Tags: Ver_LATEST, Ver_9-5, Ver_9-41, Ver_9-4, Ver_9-3, Ver_9-2, Ver_9-1, Ver_9-0, Ver_8-8, Ver_8-7, Ver_8-6, Ver_8-5, Ver_8-4, Ver_8-3, Ver_8-2, Ver_8-12, Ver_8-11, Ver_8-10, Ver_8-1, Ver_8-0, Ver_7-1, Ver_7-0, Ver_6-4, Ver_6-3, Ver_6-2, HEAD

*** empty log message ***

SUBROUTINE LIM_INIT (CONT)
! Initialize limits and normalization factors

  USE INV_PARAM
  IMPLICIT NONE
  REAL(DP), INTENT(IN)           :: CONT

! Used to be in inv_utils.f90
! Limits for the model parameters. Used in FINE_TUNE_MODEL
! Order: eta0, inclination (deg), azimuth (deg), damping, Doppler width, 
! field strength (gauss), line-of-sight velocity (cm/s), source function, 
! source function gradient, filling factor

ICONT=CONT

  LOWER_LIMIT = (/1D0, 0D0 , 0D0 , 1D-4,1D0,5D0,-7D5, 1.5E-1*ICONT, 1.5E-1*ICONT, 0D0/)
  UPPER_LIMIT  = (/1D3, 180D0,180D0,5D0, 5D2,5D3, 7D5, 1.2D0*ICONT,  1.2D0*ICONT,  1D0/)

  !NORM(:)=(/25D0,90D0,90D0,1D0,50D0,1500D0,1E5_DP,0.5D0*ICONT,0.5D0*ICONT,0.5D0/)
  NORM(:)=(/0.2D0,7D0,50D0,1D0,0.4D0,50D0,5000D0,0.003D0*ICONT,0.004D0*ICONT,0.5D0/)

! Limits on DMODEL
! I have tried with more generous limits at this is the best combination I could find
  DLIMIT(:)=(/10D0,25D0,25D0,0.1D0,30D0,500D0,1E5_DP,0.5D0*ICONT,0.5D0*ICONT,0.25D0/)
! Relative limits. 0 means don't use
  RLIMIT(:)=(/0.0D0,0D0,0D0,0D0,0D0,0D0,0D0,0D0,0D0,0D0/)
 
END SUBROUTINE LIM_INIT
!CVSVERSIONINFO "$Id: lim_init.f90,v 1.1 2012/04/10 22:17:08 keiji Exp $"

Karen Tian