File: N:\mfix\model\init_resid.f

1     !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
2     !                                                                      C
3     !  Module name: INIT_RESID(IER)                                        C
4     !  Purpose: Initialize residuals                                       C
5     !                                                                      C
6     !                                                                      C
7     !  Author: M. Syamlal                                 Date: 8-JUL-96   C
8     !  Reviewer:                                          Date:            C
9     !                                                                      C
10     !                                                                      C
11     !  Literature/Document References:                                     C
12     !                                                                      C
13     !  Variables referenced:                                               C
14     !  Variables modified:                                                 C
15     !                                                                      C
16     !  Local variables:                                                    C
17     !                                                                      C
18     !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
19     !
20           SUBROUTINE INIT_RESID()
21     !...Translated by Pacific-Sierra Research VAST-90 2.06G5  12:17:31  12/09/98
22     !...Switches: -xf
23     !
24     !  Include param.inc file to specify parameter values
25     !
26     !-----------------------------------------------
27     !   M o d u l e s
28     !-----------------------------------------------
29           USE param
30           USE param1
31           USE physprop
32           USE residual
33           IMPLICIT NONE
34     !-----------------------------------------------
35     !   G l o b a l   P a r a m e t e r s
36     !-----------------------------------------------
37     !-----------------------------------------------
38     !   D u m m y   A r g u m e n t s
39     !-----------------------------------------------
40     !
41     !                      residual number
42           INTEGER          L
43     !
44     !                      Phase index
45           INTEGER          M
46     !
47     !-----------------------------------------------
48     !
49           L = 1
50           IF (NRESID > 0) THEN
51              M = 0
52              IF (MMAX + 1 > 0) THEN
53                 RESID(:NRESID,:MMAX) = ZERO
54                 M = MMAX + 1
55              ENDIF
56              L = NRESID + 1
57           ENDIF
58           RETURN
59           END SUBROUTINE INIT_RESID
60