LCOV - code coverage report
Current view: top level - core - nwrfst.f (source / functions) Hit Total Coverage
Test: combined.info Lines: 0 34 0.0 %
Date: 2019-09-08 04:53:50 Functions: 0 1 0.0 %

          Line data    Source code
       1           0 :       SUBROUTINE nwrfst(mrad,nsol,is,it,nmatch,nzero,ferro,ec,rc,
       2           0 :      +                  pow,piw,gc,err,var,dv,varnew,errnew)
       3             : 
       4             : c..........................................................nwrfst
       5             : c starting values for E, A, A_out, A_inw
       6             : c
       7             :       IMPLICIT NONE
       8             : C     ..
       9             : C     .. Scalar Arguments ..
      10             :       INTEGER, INTENT (IN) :: mrad
      11             :       REAL ec
      12             :       INTEGER is,it,nmatch,nsol,nzero
      13             :       LOGICAL ferro
      14             : C     ..
      15             : C     .. Array Arguments ..
      16             :       REAL dv(4),err(4),errnew(4),gc(2,2,mrad),piw(2,2),pow(2,2),
      17             :      +     rc(mrad),var(4),varnew(4)
      18             : C     ..
      19             : C     .. Local Scalars ..
      20             :       REAL ratt,rr,trymix
      21             :       INTEGER iv,j,n
      22             : C     ..
      23             : C     .. Local Arrays ..
      24             :       REAL niw(2),now(2)
      25             : C     ..
      26             : C     .. Data statements ..
      27             : C
      28             :       DATA trymix/0.010/
      29             : C     ..
      30             : C                                      --------------------
      31             : C                                       START VALUES FOR
      32             : C                                           PARAMETERS
      33             : C                                      --------------------
      34           0 :       var(1) = ec
      35           0 :       var(2) = pow(is,is)/piw(is,is)
      36             : c
      37           0 :       IF ((nsol.EQ.2) .AND. ferro) THEN
      38           0 :          DO 10 j = 1,nsol
      39           0 :             now(j) = 0.00
      40           0 :    10    CONTINUE
      41           0 :          DO 30 n = 1,nmatch - 1
      42           0 :             rr = rc(n)**3
      43           0 :             DO 20 j = 1,nsol
      44           0 :                now(j) = now(j) + gc(j,j,n)**2*rr
      45           0 :    20       CONTINUE
      46           0 :    30    CONTINUE
      47           0 :          DO 40 j = 1,nsol
      48           0 :             niw(j) = 0.00
      49           0 :    40    CONTINUE
      50           0 :          DO 60 n = nmatch,nzero - 1
      51           0 :             rr = rc(n)**3
      52           0 :             DO 50 j = 1,nsol
      53           0 :                niw(j) = niw(j) + gc(j,j,n)**2*rr
      54           0 :    50       CONTINUE
      55           0 :    60    CONTINUE
      56           0 :          ratt = pow(it,it)/piw(it,it)
      57             :          var(3) = trymix* (now(is)+niw(is)*var(2))/
      58           0 :      +            (now(it)+niw(it)*ratt)
      59           0 :          var(4) = ratt*var(3)/var(2)
      60             :       ELSE
      61           0 :          DO 70 iv = 3,4
      62           0 :             err(iv) = 0.00
      63           0 :             errnew(iv) = 0.00
      64           0 :             var(iv) = 0.00
      65           0 :             varnew(iv) = 0.00
      66           0 :             dv(iv) = 0.00
      67           0 :    70    CONTINUE
      68             :       END IF
      69             : c
      70           0 :       RETURN
      71             :       END

Generated by: LCOV version 1.13