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

« back to all changes in this revision

Viewing changes to applic/fit/libsrc/ftimag.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  @(#)ftimag.for       19.1 (ESO-DMD) 02/25/03 13:17:34 
 
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 Massachusetts Ave, Cambridge, 
 
18
C MA 02139, USA.
 
19
C
 
20
C Correspondence 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 FTIMAG(NAME,MASK,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 INITIALIZE DATA STRUCTURE TO FIT IMAGE DATA                                   
 
48
C DEPENDENT VARIABLE IS THE IMAGE DATA                                          
 
49
C INDEPENDENT VARIABLES ARE IMAGE AXES                                          
 
50
C                                                                               
 
51
C.ALGORITHM                                                                     
 
52
C                                                                               
 
53
C  USE MIDAS I/O INTERFACES TO FRAMES AND TABLES                                
 
54
C                                                                               
 
55
C.KEYWORDS                                                                      
 
56
C                                                                               
 
57
C  NON LINEAR FITTING                                                           
 
58
C                                                                               
 
59
C.VERSION
 
60
C 021122        last modif
 
61
C                                                                               
 
62
C----------------------------------------------------------------               
 
63
C                                                                               
 
64
C                                                                               
 
65
C INPUT PARAMETERS                                                              
 
66
C NAME CHAR IMAGE NAME                                                          
 
67
C MASK CHAR OPTIONAL WEIGHTING MASK                                             
 
68
C   BLANKS IF NOT USED                                                          
 
69
C                                                                               
 
70
C OUTPUT PARAMETERS                                                             
 
71
C ISTAT INTG STATUS RETURN                                                      
 
72
C                                                                               
 
73
      IMPLICIT NONE
 
74
 
 
75
      CHARACTER*(*) NAME,MASK                                                   
 
76
      CHARACTER UNIT*72,IDENT*72                                                
 
77
 
 
78
      INTEGER NPIX(8), I,ISTAT,NAXIS,IMASK
 
79
      INTEGER*8  IPTR
 
80
 
 
81
      DOUBLE PRECISION STR(3), STP(3)
 
82
 
83
      INCLUDE 'MID_INCLUDE:ST_DEF.INC'
 
84
      INCLUDE 'MID_INCLUDE:FITI.INC/NOLIST'  
 
85
      INCLUDE 'MID_INCLUDE:FITC.INC/NOLIST'  
 
86
      INCLUDE 'MID_INCLUDE:ST_DAT.INC'
 
87
C                                                                               
 
88
C READ OPTIONAL MASK                                                            
 
89
C                                                                               
 
90
      IF (MASK(1:1).EQ.' ') THEN                                                
 
91
          FZPTRM = 0                                                          
 
92
          FZWEIGHT = 0                                                        
 
93
      ELSE                                                                      
 
94
          CALL STIGET(MASK,D_R4_FORMAT,F_I_MODE,F_IMA_TYPE,
 
95
     +    3,NAXIS,NPIX,STR,STP,IDENT,UNIT,IPTR,IMASK,
 
96
     +                ISTAT)                                                    
 
97
          FZPTRM = IPTR                                                       
 
98
          FZWEIGHT = 1                                                        
 
99
      END IF                                                                    
 
100
C                                                                               
 
101
C READ IMAGE                                                                    
 
102
C                                                                               
 
103
      CALL STIGET(NAME,D_R4_FORMAT,F_I_MODE,F_IMA_TYPE,3,NAXIS,
 
104
     +            NPIX,FZSTART,FZSTEP,IDENT,
 
105
     +            UNIT,IPTR,FZIDEN,ISTAT)          
 
106
      FZPTRI = IPTR                                                           
 
107
      FZNAXIS = NAXIS                                                         
 
108
      FZDVAR = 0                                                              
 
109
      DO 10, I = 1,3 
 
110
          FZIVAR(I) = I                                                       
 
111
          FZNPIX(I) = MAX(1,NPIX(I))                                          
 
112
   10 CONTINUE                                                                  
 
113
C                                                                               
 
114
C INITIALIZE VARIABLES IN THE COMMON AREA                                       
 
115
C                                                                               
 
116
      FZNAME = NAME                                                           
 
117
      FZTYPE = 'BDF '                                                         
 
118
      FZNFUN = 0                                                              
 
119
      FZNDAT = 0                                                              
 
120
      FZNITER = 0                                                             
 
121
      FZNPTOT = 0                                                             
 
122
      FZRELAX = 0.                                                            
 
123
      FZCHISQ = 0.                                                            
 
124
      FZCCHIS = 0.                                                            
 
125
      DO 20, I = 1,FZFUNMAX                                                    
 
126
          FZFCODE(I) = 0                                                      
 
127
          FZACTPAR(I) = 0                                                     
 
128
          FZSPEC(I) = ' '                                                     
 
129
   20 CONTINUE                                                                  
 
130
      DO 30, I = 1,FZPARMAX                                                    
 
131
          FZERROR(I) = 0.D0                                                   
 
132
          FZVALUE(I) = 0.D0                                                   
 
133
          FZGUESS(I) = 0.D0                                                   
 
134
          FZUNCER(I) = 0.                                                     
 
135
          FZFIXED(I) = -1                                                     
 
136
   30 CONTINUE                                                                  
 
137
      FZMAPPED = 1                                                            
 
138
      RETURN                                                                    
 
139
      END