~ubuntu-branches/debian/jessie/eso-midas/jessie

« back to all changes in this revision

Viewing changes to applic/fit/libsrc/ftdfpr.for

  • Committer: Package Import Robot
  • Author(s): Ole Streicher
  • Date: 2014-04-22 14:44:58 UTC
  • Revision ID: package-import@ubuntu.com-20140422144458-okiwi1assxkkiz39
Tags: upstream-13.09pl1.2+dfsg
ImportĀ upstreamĀ versionĀ 13.09pl1.2+dfsg

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
C @(#)ftdfpr.for        19.1 (ES0-DMD) 02/25/03 13:17:33
 
2
C===========================================================================
 
3
C Copyright (C) 1995 European Southern Observatory (ESO)
 
4
C
 
5
C This program is free software; you can redistribute it and/or 
 
6
C modify it under the terms of the GNU General Public License as 
 
7
C published by the Free Software Foundation; either version 2 of 
 
8
C the License, or (at your option) any later version.
 
9
C
 
10
C This program is distributed in the hope that it will be useful,
 
11
C but WITHOUT ANY WARRANTY; without even the implied warranty of
 
12
C MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 
13
C GNU General Public License for more details.
 
14
C
 
15
C You should have received a copy of the GNU General Public 
 
16
C License along with this program; if not, write to the Free 
 
17
C Software Foundation, Inc., 675 Massachusetss Ave, Cambridge, 
 
18
C MA 02139, USA.
 
19
C
 
20
C Corresponding concerning ESO-MIDAS should be addressed as follows:
 
21
C       Internet e-mail: midas@eso.org
 
22
C       Postal address: European Southern Observatory
 
23
C                       Data Management Division 
 
24
C                       Karl-Schwarzschild-Strasse 2
 
25
C                       D 85748 Garching bei Muenchen 
 
26
C                       GERMANY
 
27
C===========================================================================
 
28
C
 
29
      SUBROUTINE FTDFPR(PARAM,LEN,VAL,LENV,ISTAT)                               
 
30
C++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++                   
 
31
C.COPYRIGHT: Copyright (c) 1987 European Southern Observatory,                  
 
32
C                                         all rights reserved                   
 
33
C                                                                               
 
34
C.VERSION: 1.0  ESO-FORTRAN Conversion, AA  17:25 - 13 JAN 1988                 
 
35
C                                                                               
 
36
C.LANGUAGE: F77+ESOext                                                          
 
37
C                                                                               
 
38
C.AUTHOR: J.D.PONZ                                                              
 
39
C                                                                               
 
40
C.IDENTIFICATION                                                                
 
41
C                                                                               
 
42
C  FITLIB.FOR   VERSION 1.0  27 MAR 1984                                        
 
43
C                                                                               
 
44
C.PURPOSE                                                                       
 
45
C                                                                               
 
46
C  INTERFACE ROUTINES FOR THE FITTING STRUCTURES                                
 
47
C DEFINES NEW PARAMETER GUESS                                                   
 
48
C                                                                               
 
49
C.ALGORITHM                                                                     
 
50
C                                                                               
 
51
C  USE MIDAS I/O INTERFACES TO FRAMES AND TABLES                                
 
52
C                                                                               
 
53
C.KEYWORDS                                                                      
 
54
C                                                                               
 
55
C  NON LINEAR FITTING                                                           
 
56
C                                                                               
 
57
C                                                                               
 
58
C----------------------------------------------------------------               
 
59
C                                                                               
 
60
C      IMPLICIT NONE                                                            
 
61
C                                                                               
 
62
C INPUT ARGUMENTS                                                               
 
63
C PARAM CHAR PARAMETER NAME                                                     
 
64
C LEN INTG PARAMETER LENGTH                                                     
 
65
C VAL CHAR PARAMETER GUESS                                                      
 
66
C LENV INTG VALUE LENGTH                                                        
 
67
C                                                                               
 
68
C OUTPUT ARGUMENTS                                                              
 
69
C ISTAT INTG STATUS RETURN                                                      
 
70
C                                                                               
 
71
      INTEGER LEN,LENV,ISTAT,I,J,II,II1,L,IFUN,IOFF
 
72
      CHARACTER*(*) PARAM,VAL                                                   
 
73
      CHARACTER NAME*12,FSPEC*80,CVAL*30                                        
 
74
       INCLUDE 'MID_INCLUDE:FITI.INC/NOLIST'                                     
 
75
       INCLUDE 'MID_INCLUDE:FITC.INC/NOLIST'                                     
 
76
C                                                                               
 
77
C SEARCH FOR PARAMETER (I IS THE PARAMETER NUMBER)                              
 
78
C                                                                               
 
79
      DO 10 I = 1,FZNPTOT                                                     
 
80
          NAME   = FZPTOKEN(I)                                                
 
81
          IF (NAME(1:FZPLEN(I)).EQ.PARAM(1:LEN)) GO TO 20                     
 
82
   10 CONTINUE                                                                  
 
83
      ISTAT  = FERPFU                                                       
 
84
      RETURN                                                                    
 
85
C                                                                               
 
86
C FIND FUNCTION  (J IS THE FUNCTION NUMBER)                                     
 
87
C                                                                               
 
88
   20 CONTINUE                                                                  
 
89
      IF (FZFIXED(I).GT.0) THEN                                               
 
90
          ISTAT  = FERPFU                                                   
 
91
          RETURN                                                                
 
92
                                                                                
 
93
      END IF                                                                    
 
94
                                                                                
 
95
      IOFF   = 0                                                                
 
96
      DO 30 J = 1,FZNFUN                                                      
 
97
          IF (I.LE.IOFF+FZACTPAR(J)) THEN                                     
 
98
              IFUN   = FZFCODE(J)                                             
 
99
              GO TO 40                                                          
 
100
                                                                                
 
101
          END IF                                                                
 
102
                                                                                
 
103
          IOFF   = IOFF + FZACTPAR(J)                                         
 
104
   30 CONTINUE                                                                  
 
105
      ISTAT  = FERFUN                                                       
 
106
      RETURN                                                                    
 
107
C                                                                               
 
108
C SEARCH FOR INITIAL VALUE IN FSPEC                                             
 
109
C                                                                               
 
110
   40 FSPEC  = FZSPEC(J)                                                      
 
111
      NAME   = PARAM(1:LEN)//'='                                                
 
112
      II     = INDEX(FSPEC,NAME(1:LEN+1))                                       
 
113
      IF (II.EQ.0) THEN                                                         
 
114
          II     = INDEX(FSPEC,')') + 2                                         
 
115
                                                                                
 
116
      ELSE                                                                      
 
117
          II1    = INDEX(FSPEC(II:),' ')                                        
 
118
          FSPEC(II:) = FZSPEC(J) (II1:)                                       
 
119
      END IF                                                                    
 
120
C                                                                               
 
121
C INSERT VAL IN POSITION II                                                     
 
122
C                                                                               
 
123
      IF (FZFIXED(I).EQ.0) THEN                                               
 
124
          CVAL   = VAL(1:LENV)//'@'                                             
 
125
          L      = LENV + 1                                                     
 
126
                                                                                
 
127
      ELSE                                                                      
 
128
          CVAL   = VAL                                                          
 
129
          L      = LENV                                                         
 
130
      END IF                                                                    
 
131
                                                                                
 
132
      FZSPEC(J) = FSPEC(1:II-1)//CVAL(1:L)//FSPEC(II:)                        
 
133
      RETURN                                                                    
 
134
                                                                                
 
135
      END