Changeset 14387


Ignore:
Timestamp:
2021-02-03T16:51:56+01:00 (6 months ago)
Author:
smasson
Message:

PPR: suppress real*8 ↔ real*16 compilation errors or implicit conversions, see #2617

Location:
vendors/PPR/src
Files:
15 edited

Legend:

Unmodified
Added
Removed
  • vendors/PPR/src/bfun1d.h90

    r14212 r14387  
    6060            select case (ndof) 
    6161            case (+1) 
    62                 bfun(1) = sval**1 / 1.d0 
    63                  
    64             case (+2) 
    65                 bfun(1) = sval**1 / 1.d0 
    66                 bfun(2) = sval**2 / 2.d0 
    67                  
    68             case (+3) 
    69                 bfun(1) = sval**1 / 1.d0 
    70                 bfun(2) = sval**2 / 2.d0 
    71                 bfun(3) = sval**3 / 3.d0 
    72                  
    73             case (+4) 
    74                 bfun(1) = sval**1 / 1.d0 
    75                 bfun(2) = sval**2 / 2.d0 
    76                 bfun(3) = sval**3 / 3.d0 
    77                 bfun(4) = sval**4 / 4.d0 
    78                  
    79             case (+5) 
    80                 bfun(1) = sval**1 / 1.d0 
    81                 bfun(2) = sval**2 / 2.d0 
    82                 bfun(3) = sval**3 / 3.d0 
    83                 bfun(4) = sval**4 / 4.d0 
    84                 bfun(5) = sval**5 / 5.d0 
    85                  
    86             case (+6) 
    87                 bfun(1) = sval**1 / 1.d0 
    88                 bfun(2) = sval**2 / 2.d0 
    89                 bfun(3) = sval**3 / 3.d0 
    90                 bfun(4) = sval**4 / 4.d0 
    91                 bfun(5) = sval**5 / 5.d0 
    92                 bfun(6) = sval**6 / 6.d0 
    93                  
    94             case (+7) 
    95                 bfun(1) = sval**1 / 1.d0 
    96                 bfun(2) = sval**2 / 2.d0 
    97                 bfun(3) = sval**3 / 3.d0 
    98                 bfun(4) = sval**4 / 4.d0 
    99                 bfun(5) = sval**5 / 5.d0 
    100                 bfun(6) = sval**6 / 6.d0 
    101                 bfun(7) = sval**7 / 7.d0 
     62                bfun(1) = sval**1 / 1.e0 
     63                 
     64            case (+2) 
     65                bfun(1) = sval**1 / 1.e0 
     66                bfun(2) = sval**2 / 2.e0 
     67                 
     68            case (+3) 
     69                bfun(1) = sval**1 / 1.e0 
     70                bfun(2) = sval**2 / 2.e0 
     71                bfun(3) = sval**3 / 3.e0 
     72                 
     73            case (+4) 
     74                bfun(1) = sval**1 / 1.e0 
     75                bfun(2) = sval**2 / 2.e0 
     76                bfun(3) = sval**3 / 3.e0 
     77                bfun(4) = sval**4 / 4.e0 
     78                 
     79            case (+5) 
     80                bfun(1) = sval**1 / 1.e0 
     81                bfun(2) = sval**2 / 2.e0 
     82                bfun(3) = sval**3 / 3.e0 
     83                bfun(4) = sval**4 / 4.e0 
     84                bfun(5) = sval**5 / 5.e0 
     85                 
     86            case (+6) 
     87                bfun(1) = sval**1 / 1.e0 
     88                bfun(2) = sval**2 / 2.e0 
     89                bfun(3) = sval**3 / 3.e0 
     90                bfun(4) = sval**4 / 4.e0 
     91                bfun(5) = sval**5 / 5.e0 
     92                bfun(6) = sval**6 / 6.e0 
     93                 
     94            case (+7) 
     95                bfun(1) = sval**1 / 1.e0 
     96                bfun(2) = sval**2 / 2.e0 
     97                bfun(3) = sval**3 / 3.e0 
     98                bfun(4) = sval**4 / 4.e0 
     99                bfun(5) = sval**5 / 5.e0 
     100                bfun(6) = sval**6 / 6.e0 
     101                bfun(7) = sval**7 / 7.e0 
    102102             
    103103            end select 
     
    107107            select case (ndof) 
    108108            case (+1) 
    109                 bfun(1) =           1.d0 
    110                  
    111             case (+2) 
    112                 bfun(1) =           1.d0 
    113                 bfun(2) = sval**1 * 1.d0 
    114                  
    115             case (+3) 
    116                 bfun(1) =           1.d0 
    117                 bfun(2) = sval**1 * 1.d0 
    118                 bfun(3) = sval**2 * 1.d0 
    119                  
    120             case (+4) 
    121                 bfun(1) =           1.d0 
    122                 bfun(2) = sval**1 * 1.d0 
    123                 bfun(3) = sval**2 * 1.d0 
    124                 bfun(4) = sval**3 * 1.d0 
    125                  
    126             case (+5) 
    127                 bfun(1) =           1.d0 
    128                 bfun(2) = sval**1 * 1.d0 
    129                 bfun(3) = sval**2 * 1.d0 
    130                 bfun(4) = sval**3 * 1.d0 
    131                 bfun(5) = sval**4 * 1.d0 
    132                  
    133             case (+6) 
    134                 bfun(1) =           1.d0 
    135                 bfun(2) = sval**1 * 1.d0 
    136                 bfun(3) = sval**2 * 1.d0 
    137                 bfun(4) = sval**3 * 1.d0 
    138                 bfun(5) = sval**4 * 1.d0 
    139                 bfun(6) = sval**5 * 1.d0 
    140                  
    141             case (+7) 
    142                 bfun(1) =           1.d0 
    143                 bfun(2) = sval**1 * 1.d0 
    144                 bfun(3) = sval**2 * 1.d0 
    145                 bfun(4) = sval**3 * 1.d0 
    146                 bfun(5) = sval**4 * 1.d0 
    147                 bfun(6) = sval**5 * 1.d0 
    148                 bfun(7) = sval**6 * 1.d0 
     109                bfun(1) =           1.e0 
     110                 
     111            case (+2) 
     112                bfun(1) =           1.e0 
     113                bfun(2) = sval**1 * 1.e0 
     114                 
     115            case (+3) 
     116                bfun(1) =           1.e0 
     117                bfun(2) = sval**1 * 1.e0 
     118                bfun(3) = sval**2 * 1.e0 
     119                 
     120            case (+4) 
     121                bfun(1) =           1.e0 
     122                bfun(2) = sval**1 * 1.e0 
     123                bfun(3) = sval**2 * 1.e0 
     124                bfun(4) = sval**3 * 1.e0 
     125                 
     126            case (+5) 
     127                bfun(1) =           1.e0 
     128                bfun(2) = sval**1 * 1.e0 
     129                bfun(3) = sval**2 * 1.e0 
     130                bfun(4) = sval**3 * 1.e0 
     131                bfun(5) = sval**4 * 1.e0 
     132                 
     133            case (+6) 
     134                bfun(1) =           1.e0 
     135                bfun(2) = sval**1 * 1.e0 
     136                bfun(3) = sval**2 * 1.e0 
     137                bfun(4) = sval**3 * 1.e0 
     138                bfun(5) = sval**4 * 1.e0 
     139                bfun(6) = sval**5 * 1.e0 
     140                 
     141            case (+7) 
     142                bfun(1) =           1.e0 
     143                bfun(2) = sval**1 * 1.e0 
     144                bfun(3) = sval**2 * 1.e0 
     145                bfun(4) = sval**3 * 1.e0 
     146                bfun(5) = sval**4 * 1.e0 
     147                bfun(6) = sval**5 * 1.e0 
     148                bfun(7) = sval**6 * 1.e0 
    149149             
    150150            end select 
     
    154154            select case (ndof) 
    155155            case (+1) 
    156                 bfun(1) =           0.d0 
    157                  
    158             case (+2) 
    159                 bfun(1) =           0.d0 
    160                 bfun(2) =           1.d0 
    161                  
    162             case (+3) 
    163                 bfun(1) =           0.d0 
    164                 bfun(2) =           1.d0 
    165                 bfun(3) = sval**1 * 2.d0 
    166                  
    167             case (+4) 
    168                 bfun(1) =           0.d0 
    169                 bfun(2) =           1.d0 
    170                 bfun(3) = sval**1 * 2.d0 
    171                 bfun(4) = sval**2 * 3.d0 
    172                  
    173             case (+5) 
    174                 bfun(1) =           0.d0 
    175                 bfun(2) =           1.d0 
    176                 bfun(3) = sval**1 * 2.d0 
    177                 bfun(4) = sval**2 * 3.d0 
    178                 bfun(5) = sval**3 * 4.d0 
    179                  
    180             case (+6) 
    181                 bfun(1) =           0.d0 
    182                 bfun(2) =           1.d0 
    183                 bfun(3) = sval**1 * 2.d0 
    184                 bfun(4) = sval**2 * 3.d0 
    185                 bfun(5) = sval**3 * 4.d0 
    186                 bfun(6) = sval**4 * 5.d0 
    187                  
    188             case (+7) 
    189                 bfun(1) =           0.d0 
    190                 bfun(2) =           1.d0 
    191                 bfun(3) = sval**1 * 2.d0 
    192                 bfun(4) = sval**2 * 3.d0 
    193                 bfun(5) = sval**3 * 4.d0 
    194                 bfun(6) = sval**4 * 5.d0 
    195                 bfun(7) = sval**5 * 6.d0 
     156                bfun(1) =           0.e0 
     157                 
     158            case (+2) 
     159                bfun(1) =           0.e0 
     160                bfun(2) =           1.e0 
     161                 
     162            case (+3) 
     163                bfun(1) =           0.e0 
     164                bfun(2) =           1.e0 
     165                bfun(3) = sval**1 * 2.e0 
     166                 
     167            case (+4) 
     168                bfun(1) =           0.e0 
     169                bfun(2) =           1.e0 
     170                bfun(3) = sval**1 * 2.e0 
     171                bfun(4) = sval**2 * 3.e0 
     172                 
     173            case (+5) 
     174                bfun(1) =           0.e0 
     175                bfun(2) =           1.e0 
     176                bfun(3) = sval**1 * 2.e0 
     177                bfun(4) = sval**2 * 3.e0 
     178                bfun(5) = sval**3 * 4.e0 
     179                 
     180            case (+6) 
     181                bfun(1) =           0.e0 
     182                bfun(2) =           1.e0 
     183                bfun(3) = sval**1 * 2.e0 
     184                bfun(4) = sval**2 * 3.e0 
     185                bfun(5) = sval**3 * 4.e0 
     186                bfun(6) = sval**4 * 5.e0 
     187                 
     188            case (+7) 
     189                bfun(1) =           0.e0 
     190                bfun(2) =           1.e0 
     191                bfun(3) = sval**1 * 2.e0 
     192                bfun(4) = sval**2 * 3.e0 
     193                bfun(5) = sval**3 * 4.e0 
     194                bfun(6) = sval**4 * 5.e0 
     195                bfun(7) = sval**5 * 6.e0 
    196196             
    197197            end select 
     
    201201            select case (ndof) 
    202202            case (+1) 
    203                 bfun(1) =           0.d0 
    204                  
    205             case (+2) 
    206                 bfun(1) =           0.d0 
    207                 bfun(2) =           0.d0 
    208                  
    209             case (+3) 
    210                 bfun(1) =           0.d0 
    211                 bfun(2) =           0.d0 
    212                 bfun(3) =           2.d0 
    213                  
    214             case (+4) 
    215                 bfun(1) =           0.d0 
    216                 bfun(2) =           0.d0 
    217                 bfun(3) =           2.d0 
    218                 bfun(4) = sval**1 * 6.d0 
    219                  
    220             case (+5) 
    221                 bfun(1) =           0.d0 
    222                 bfun(2) =           0.d0 
    223                 bfun(3) =           2.d0 
    224                 bfun(4) = sval**1 * 6.d0 
    225                 bfun(5) = sval**2 *12.d0 
    226                  
    227             case (+6) 
    228                 bfun(1) =           0.d0 
    229                 bfun(2) =           0.d0 
    230                 bfun(3) =           2.d0 
    231                 bfun(4) = sval**1 * 6.d0 
    232                 bfun(5) = sval**2 *12.d0 
    233                 bfun(6) = sval**3 *20.d0 
    234                  
    235             case (+7) 
    236                 bfun(1) =           0.d0 
    237                 bfun(2) =           0.d0 
    238                 bfun(3) =           2.d0 
    239                 bfun(4) = sval**1 * 6.d0 
    240                 bfun(5) = sval**2 *12.d0 
    241                 bfun(6) = sval**3 *20.d0 
    242                 bfun(7) = sval**4 *30.d0 
     203                bfun(1) =           0.e0 
     204                 
     205            case (+2) 
     206                bfun(1) =           0.e0 
     207                bfun(2) =           0.e0 
     208                 
     209            case (+3) 
     210                bfun(1) =           0.e0 
     211                bfun(2) =           0.e0 
     212                bfun(3) =           2.e0 
     213                 
     214            case (+4) 
     215                bfun(1) =           0.e0 
     216                bfun(2) =           0.e0 
     217                bfun(3) =           2.e0 
     218                bfun(4) = sval**1 * 6.e0 
     219                 
     220            case (+5) 
     221                bfun(1) =           0.e0 
     222                bfun(2) =           0.e0 
     223                bfun(3) =           2.e0 
     224                bfun(4) = sval**1 * 6.e0 
     225                bfun(5) = sval**2 *12.e0 
     226                 
     227            case (+6) 
     228                bfun(1) =           0.e0 
     229                bfun(2) =           0.e0 
     230                bfun(3) =           2.e0 
     231                bfun(4) = sval**1 * 6.e0 
     232                bfun(5) = sval**2 *12.e0 
     233                bfun(6) = sval**3 *20.e0 
     234                 
     235            case (+7) 
     236                bfun(1) =           0.e0 
     237                bfun(2) =           0.e0 
     238                bfun(3) =           2.e0 
     239                bfun(4) = sval**1 * 6.e0 
     240                bfun(5) = sval**2 *12.e0 
     241                bfun(6) = sval**3 *20.e0 
     242                bfun(7) = sval**4 *30.e0 
    243243             
    244244            end select 
  • vendors/PPR/src/ffsl1d.h90

    r14212 r14387  
    7878        integer :: head,tail,nprt 
    7979     
    80         head = +0 ; tail = +0 ; qedg = 0.d+0 
     80        head = +0 ; tail = +0 ; qedg = 0.e+0 
    8181     
    8282        do while (.true.) 
     
    251251        do  ipos = +2 , npos - 1 
    252252         
    253             if (uvel(ipos) .gt. +0.d0) then 
     253            if (uvel(ipos) .gt. +0.e0) then 
    254254             
    255255    !----------- integrate profile over upwind cell IPOS-1 ! 
     
    265265            end if 
    266266  
    267             ss11 = +1.d0 - 2.d0 * uCFL 
    268             ss22 = +1.d0 
     267            ss11 = +1.e0 - 2.e0 * uCFL 
     268            ss22 = +1.e0 
    269269 
    270270            call bfun1d(-1,mdof,ss11,vv11) 
     
    286286                       
    287287            else & 
    288         &   if (uvel(ipos) .lt. -0.d0) then 
     288        &   if (uvel(ipos) .lt. -0.e0) then 
    289289             
    290290    !----------- integrate profile over upwind cell IPOS+0 ! 
     
    300300            end if 
    301301 
    302             ss11 = -1.d0 - 2.d0 * uCFL 
    303             ss22 = -1.d0 
     302            ss11 = -1.e0 - 2.e0 * uCFL 
     303            ss22 = -1.e0 
    304304         
    305305            call bfun1d(-1,mdof,ss11,vv11) 
  • vendors/PPR/src/inv.h90

    r14212 r14387  
    140140    !-------------------------------- C = C + scal * A * B ! 
    141141 
    142         if (scal .eq. +1.d0) then 
     142        if (scal .eq. +1.e0) then 
    143143 
    144144        cmat(1,1) = cmat(1,1) &  
     
    157157 
    158158        else & 
    159         if (scal .eq. -1.d0) then 
     159        if (scal .eq. -1.e0) then 
    160160 
    161161        cmat(1,1) = cmat(1,1) &  
     
    211211    !-------------------------------- C = C + scal * A * B ! 
    212212 
    213         if (scal .eq. +1.d0) then 
     213        if (scal .eq. +1.e0) then 
    214214 
    215215        cmat(1,1) = cmat(1,1) & 
     
    253253 
    254254        else & 
    255         if (scal .eq. -1.d0) then 
     255        if (scal .eq. -1.e0) then 
    256256 
    257257        cmat(1,1) = cmat(1,1) & 
     
    500500    !---------------------------------------- L = C * A^-1 ! 
    501501 
    502         lmat(1,1) = +0.d0 
    503         lmat(1,2) = +0.d0 
    504         lmat(2,1) = +0.d0 
    505         lmat(2,2) = +0.d0 
     502        lmat(1,1) = +0.e0 
     503        lmat(1,2) = +0.e0 
     504        lmat(2,1) = +0.e0 
     505        lmat(2,2) = +0.e0 
    506506 
    507507        call mul_2x2(amat(3,1),adim,ainv,LDIM, & 
    508                     +1.d0,lmat,LDIM) 
     508                    +1.e0,lmat,LDIM) 
    509509         
    510510    !---------------------------------------- U = A^-1 * B !        
    511511 
    512         umat(1,1) = +0.d0 
    513         umat(1,2) = +0.d0 
    514         umat(2,1) = +0.d0 
    515         umat(2,2) = +0.d0 
     512        umat(1,1) = +0.e0 
     513        umat(1,2) = +0.e0 
     514        umat(2,1) = +0.e0 
     515        umat(2,2) = +0.e0 
    516516  
    517517        call mul_2x2(ainv,LDIM,amat(1,3),adim, & 
    518                     +1.d0,umat,LDIM) 
     518                    +1.e0,umat,LDIM) 
    519519         
    520520    !-------------------------------- S = D - C * A^-1 * B ! 
     
    526526 
    527527        call mul_2x2(lmat,LDIM,amat(1,3),adim, & 
    528                     -1.d0/adet,smat,LDIM) 
     528                    -1.e0/adet,smat,LDIM) 
    529529 
    530530        call inv_2x2(smat,LDIM,sinv,LDIM,sdet) 
     
    641641    !---------------------------------------- L = C * A^-1 ! 
    642642 
    643         lmat(1,1) = +0.d0 
    644         lmat(1,2) = +0.d0 
    645         lmat(1,3) = +0.d0 
    646         lmat(2,1) = +0.d0 
    647         lmat(2,2) = +0.d0 
    648         lmat(2,3) = +0.d0 
    649         lmat(3,1) = +0.d0 
    650         lmat(3,2) = +0.d0 
    651         lmat(3,3) = +0.d0 
     643        lmat(1,1) = +0.e0 
     644        lmat(1,2) = +0.e0 
     645        lmat(1,3) = +0.e0 
     646        lmat(2,1) = +0.e0 
     647        lmat(2,2) = +0.e0 
     648        lmat(2,3) = +0.e0 
     649        lmat(3,1) = +0.e0 
     650        lmat(3,2) = +0.e0 
     651        lmat(3,3) = +0.e0 
    652652 
    653653        call mul_3x3(amat(4,1),adim,ainv,LDIM, & 
    654                     +1.d0,lmat,LDIM) 
     654                    +1.e0,lmat,LDIM) 
    655655         
    656656    !---------------------------------------- U = A^-1 * B !        
    657657 
    658         umat(1,1) = +0.d0 
    659         umat(1,2) = +0.d0 
    660         umat(1,3) = +0.d0 
    661         umat(2,1) = +0.d0 
    662         umat(2,2) = +0.d0 
    663         umat(2,3) = +0.d0 
    664         umat(3,1) = +0.d0 
    665         umat(3,2) = +0.d0 
    666         umat(3,3) = +0.d0 
     658        umat(1,1) = +0.e0 
     659        umat(1,2) = +0.e0 
     660        umat(1,3) = +0.e0 
     661        umat(2,1) = +0.e0 
     662        umat(2,2) = +0.e0 
     663        umat(2,3) = +0.e0 
     664        umat(3,1) = +0.e0 
     665        umat(3,2) = +0.e0 
     666        umat(3,3) = +0.e0 
    667667  
    668668        call mul_3x3(ainv,LDIM,amat(1,4),adim, & 
    669                     +1.d0,umat,LDIM) 
     669                    +1.e0,umat,LDIM) 
    670670         
    671671    !-------------------------------- S = D - C * A^-1 * B ! 
     
    682682 
    683683        call mul_3x3(lmat,LDIM,amat(1,4),adim, & 
    684                     -1.d0/adet,smat,LDIM) 
     684                    -1.e0/adet,smat,LDIM) 
    685685 
    686686        call inv_3x3(smat,LDIM,sinv,LDIM,sdet) 
  • vendors/PPR/src/oscl1d.h90

    r14212 r14387  
    6767        do  ipos = +1, npos-1 
    6868        do  ivar = +1, nvar-0 
    69             oscl(1,ivar,ipos) = +0.d0 
    70             oscl(2,ivar,ipos) = +0.d0 
     69            oscl(1,ivar,ipos) = +0.e0 
     70            oscl(2,ivar,ipos) = +0.e0 
    7171        end do 
    7272        end do 
     
    144144        hhmm = hhll + hhcc + hhrr 
    145145 
    146         cmat(1,1) = -(hhcc+2.d0*hhrr)/(hhlc*hhmm) 
     146        cmat(1,1) = -(hhcc+2.e0*hhrr)/(hhlc*hhmm) 
    147147        cmat(1,2) = -(hhll-hhrr)* & 
    148         &          (3.d0*hhcc+2.d0*(hhll+hhrr))/& 
     148        &          (3.e0*hhcc+2.e0*(hhll+hhrr))/& 
    149149        &            (hhlc*hhrc*hhmm) 
    150         cmat(1,3) = +(hhcc+2.d0*hhll)/(hhrc*hhmm) 
    151  
    152         cmat(2,1) = +3.d0/(hhlc*hhmm) 
    153         cmat(2,2) = -3.d0*(2.d0*hhcc+hhll+hhrr)/& 
     150        cmat(1,3) = +(hhcc+2.e0*hhll)/(hhrc*hhmm) 
     151 
     152        cmat(2,1) = +3.e0/(hhlc*hhmm) 
     153        cmat(2,2) = -3.e0*(2.e0*hhcc+hhll+hhrr)/& 
    154154        &            (hhlc*hhrc*hhmm) 
    155         cmat(2,3) = +3.d0/(hhrc*hhmm) 
    156  
    157         do  ivar = 1, nvar 
    158  
    159             oscl(1,ivar,ipos) = +1.d0 * ( & 
     155        cmat(2,3) = +3.e0/(hhrc*hhmm) 
     156 
     157        do  ivar = 1, nvar 
     158 
     159            oscl(1,ivar,ipos) = +1.e0 * ( & 
    160160        & + cmat(1,1)*fdat(1,ivar,ipos-1) & 
    161161        & + cmat(1,2)*fdat(1,ivar,ipos+0) & 
    162162        & + cmat(1,3)*fdat(1,ivar,ipos+1) ) 
    163163 
    164             oscl(2,ivar,ipos) = +2.d0 * ( & 
     164            oscl(2,ivar,ipos) = +2.e0 * ( & 
    165165        & + cmat(2,1)*fdat(1,ivar,ipos-1) & 
    166166        & + cmat(2,2)*fdat(1,ivar,ipos+0) & 
     
    177177        hhrr = max(delx(head+2),dmin) 
    178178 
    179         cmat(1,1) = -2.d0 / (hhll+hhcc) 
    180         cmat(1,2) = +2.d0 / (hhll+hhcc) 
     179        cmat(1,1) = -2.e0 / (hhll+hhcc) 
     180        cmat(1,2) = +2.e0 / (hhll+hhcc) 
    181181 
    182182        do  ivar = 1, nvar 
     
    186186        & + cmat(1,2)*fdat(1,ivar,head+1) 
    187187              
    188             oscl(2,ivar,head) = +0.d0 
     188            oscl(2,ivar,head) = +0.e0 
    189189 
    190190        end do 
     
    196196        hhrr = max(delx(tail-0),dmin) 
    197197 
    198         cmat(1,2) = -2.d0 / (hhrr+hhcc) 
    199         cmat(1,3) = +2.d0 / (hhrr+hhcc) 
     198        cmat(1,2) = -2.e0 / (hhrr+hhcc) 
     199        cmat(1,3) = +2.e0 / (hhrr+hhcc) 
    200200 
    201201        do  ivar = 1, nvar 
     
    205205        & + cmat(1,3)*fdat(1,ivar,tail+0) 
    206206 
    207             oscl(2,ivar,tail) = +0.d0 
     207            oscl(2,ivar,tail) = +0.e0 
    208208 
    209209        end do 
     
    269269        & - .50d+0 * fdat(1,ivar,head+0) 
    270270         
    271             oscl(2,ivar,head) = +0.d0 
     271            oscl(2,ivar,head) = +0.e0 
    272272 
    273273        end do 
     
    281281        & - .50d+0 * fdat(1,ivar,tail-1) 
    282282             
    283             oscl(2,ivar,tail) = +0.d0 
     283            oscl(2,ivar,tail) = +0.e0 
    284284 
    285285        end do 
  • vendors/PPR/src/p1e.h90

    r14212 r14387  
    8484         
    8585            edge(ivar,1) = fdat(1,ivar,1) 
    86             dfdx(ivar,1) = 0.d0 
     86            dfdx(ivar,1) = 0.e0 
    8787             
    8888            edge(ivar,2) = fdat(1,ivar,1) 
    89             dfdx(ivar,2) = 0.d0 
     89            dfdx(ivar,2) = 0.e0 
    9090             
    9191        end do 
     
    104104    !--------------- reconstruction: constant grid-spacing ! 
    105105             
    106             dd10 = delx(+1) * 2.d0 
     106            dd10 = delx(+1) * 2.e0 
    107107             
    108108            do  ivar = +1, nvar 
     
    176176        &       fdat(+1,ivar,tail+0) 
    177177         
    178             dfdx(ivar,head-1) = 0.d0 
    179             dfdx(ivar,tail+1) = 0.d0 
     178            dfdx(ivar,head-1) = 0.e0 
     179            dfdx(ivar,tail+1) = 0.e0 
    180180 
    181181        end do 
  • vendors/PPR/src/p3e.h90

    r14212 r14387  
    8282         
    8383        integer, parameter :: NSIZ = +4 
    84         real*8 , parameter :: ZERO = 1.d-14   
     84        real*8 , parameter :: ZERO = 1.e-14   
    8585    
    8686        head = +3 ; tail = npos - 2 
     
    121121 
    122122                edge(ivar,ipos) = ( & 
    123         &     -      1.d0 * & 
     123        &     -      1.e0 * & 
    124124        &       fdat(1,ivar,ipos-2) & 
    125         &     +      7.d0 * & 
     125        &     +      7.e0 * & 
    126126        &       fdat(1,ivar,ipos-1) & 
    127         &     +      7.d0 * & 
     127        &     +      7.e0 * & 
    128128        &       fdat(1,ivar,ipos+0) & 
    129         &     -      1.d0 * & 
    130         &       fdat(1,ivar,ipos+1) ) / 12.d0 
     129        &     -      1.e0 * & 
     130        &       fdat(1,ivar,ipos+1) ) / 12.e0 
    131131             
    132132                dfdx(ivar,ipos) = ( & 
    133         &     +      1.d0 * & 
     133        &     +      1.e0 * & 
    134134        &       fdat(1,ivar,ipos-2) & 
    135         &     -     15.d0 * & 
     135        &     -     15.e0 * & 
    136136        &       fdat(1,ivar,ipos-1) & 
    137         &     +     15.d0 * & 
     137        &     +     15.e0 * & 
    138138        &       fdat(1,ivar,ipos+0) & 
    139         &     -      1.d0 * & 
    140         &       fdat(1,ivar,ipos+1) ) / 12.d0 
     139        &     -      1.e0 * & 
     140        &       fdat(1,ivar,ipos+1) ) / 12.e0 
    141141 
    142142                dfdx(ivar,ipos) = & 
     
    166166        &              +  delh(-1) ) / xhat 
    167167            xmap(-1) = -  delh(-1)   / xhat 
    168             xmap(+0) = +  0.d0 
     168            xmap(+0) = +  0.e0 
    169169            xmap(+1) = +  delh(+0)   / xhat 
    170170            xmap(+2) = +( delh(+0) & 
  • vendors/PPR/src/p5e.h90

    r14212 r14387  
    8282         
    8383        integer, parameter :: NSIZ = +6 
    84         real*8 , parameter :: ZERO = 1.d-14    
     84        real*8 , parameter :: ZERO = 1.e-14    
    8585 
    8686        head = +4 ; tail = npos - 3 
     
    121121 
    122122                edge(ivar,ipos) = & 
    123         &     + ( 1.d0 / 60.d0) * &  
     123        &     + ( 1.e0 / 60.e0) * &  
    124124        &       fdat(1,ivar,ipos-3) & 
    125         &     - ( 8.d0 / 60.d0) * & 
     125        &     - ( 8.e0 / 60.e0) * & 
    126126        &       fdat(1,ivar,ipos-2) & 
    127         &     + (37.d0 / 60.d0) * & 
     127        &     + (37.e0 / 60.e0) * & 
    128128        &       fdat(1,ivar,ipos-1) & 
    129         &     + (37.d0 / 60.d0) * & 
     129        &     + (37.e0 / 60.e0) * & 
    130130        &       fdat(1,ivar,ipos+0) & 
    131         &     - ( 8.d0 / 60.d0) * & 
     131        &     - ( 8.e0 / 60.e0) * & 
    132132        &       fdat(1,ivar,ipos+1) & 
    133         &     + ( 1.d0 / 60.d0) * & 
     133        &     + ( 1.e0 / 60.e0) * & 
    134134        &       fdat(1,ivar,ipos+2) 
    135135             
    136136                dfdx(ivar,ipos) = & 
    137         &     - ( 1.d0 / 90.d0) * &  
     137        &     - ( 1.e0 / 90.e0) * &  
    138138        &       fdat(1,ivar,ipos-3) & 
    139         &     + ( 5.d0 / 36.d0) * & 
     139        &     + ( 5.e0 / 36.e0) * & 
    140140        &       fdat(1,ivar,ipos-2) & 
    141         &     - (49.d0 / 36.d0) * & 
     141        &     - (49.e0 / 36.e0) * & 
    142142        &       fdat(1,ivar,ipos-1) & 
    143         &     + (49.d0 / 36.d0) * & 
     143        &     + (49.e0 / 36.e0) * & 
    144144        &       fdat(1,ivar,ipos+0) & 
    145         &     - ( 5.d0 / 36.d0) * & 
     145        &     - ( 5.e0 / 36.e0) * & 
    146146        &       fdat(1,ivar,ipos+1) & 
    147         &     + ( 1.d0 / 90.d0) * & 
     147        &     + ( 1.e0 / 90.e0) * & 
    148148        &       fdat(1,ivar,ipos+2) 
    149149 
     
    185185        &              +  delh(-1) ) / xhat 
    186186            xmap(-1) = -  delh(-1)   / xhat 
    187             xmap(+0) = +  0.d0 
     187            xmap(+0) = +  0.e0 
    188188            xmap(+1) = +  delh(+0)   / xhat 
    189189            xmap(+2) = +( delh(+0) & 
  • vendors/PPR/src/pbc.h90

    r14212 r14387  
    211211         
    212212        integer, parameter :: NSIZ = +3 
    213         real*8 , parameter :: ZERO = +1.d-14   
     213        real*8 , parameter :: ZERO = +1.e-14   
    214214 
    215215        head = +2; tail = npos - 2 
     
    245245        xmap(-1) =-(delh(-1) + & 
    246246        &           delh(+0)*0.5d0)/xhat 
    247         xmap(+0) = -1.d0 
    248         xmap(+1) = +1.d0 
     247        xmap(+0) = -1.e0 
     248        xmap(+1) = +1.e0 
    249249        xmap(+2) = (delh(+1) + & 
    250250        &           delh(+0)*0.5d0)/xhat 
     
    454454 
    455455            eval(-1) = & 
    456         &   fdat(1,ivar,head-1) * 1.d0 
     456        &   fdat(1,ivar,head-1) * 1.e0 
    457457            eval(+0) = & 
    458458        &   fdat(1,ivar,head-1) * .5d0 + & 
     
    542542 
    543543        integer, parameter :: NSIZ = +3 
    544         real*8 , parameter :: ZERO = +1.d-14 
     544        real*8 , parameter :: ZERO = +1.e-14 
    545545 
    546546        head = +2; tail = npos - 2 
     
    576576        xmap(-1) =-(delh(-1) + & 
    577577        &           delh(+0)*0.5d0)/xhat 
    578         xmap(+0) = -1.d0 
    579         xmap(+1) = +1.d0 
     578        xmap(+0) = -1.e0 
     579        xmap(+1) = +1.e0 
    580580        xmap(+2) = (delh(+1) + & 
    581581        &           delh(+0)*0.5d0)/xhat 
     
    791791        &   fdat(1,ivar,tail+1) * .5d0 
    792792            eval(+2) = & 
    793         &   fdat(1,ivar,tail+1) * 1.d0 
     793        &   fdat(1,ivar,tail+1) * 1.e0 
    794794    
    795795            gval(+0) = & 
  • vendors/PPR/src/plm.h90

    r14212 r14387  
    127127            fhat(1,ivar,1) = & 
    128128        &   fdat(1,ivar,1) 
    129             fhat(2,ivar,1) = 0.d+0 
     129            fhat(2,ivar,1) = 0.e+0 
    130130        end do 
    131131        end if 
     
    235235            fhat(1,ivar,1) = & 
    236236        &   fdat(1,ivar,1) 
    237             fhat(2,ivar,1) = 0.d+0 
     237            fhat(2,ivar,1) = 0.e+0 
    238238        end do 
    239239        end if 
     
    325325        real*8  :: fell,ferr,scal 
    326326 
    327         real*8 , parameter :: ZERO = 1.d-14 
     327        real*8 , parameter :: ZERO = 1.e-14 
    328328 
    329329    !---------------------------- 2nd-order approximations ! 
     
    401401        real*8  :: fell,ferr,scal 
    402402 
    403         real*8 , parameter :: ZERO = 1.d-14 
     403        real*8 , parameter :: ZERO = 1.e-14 
    404404 
    405405    !---------------------------- 2nd-order approximations ! 
  • vendors/PPR/src/ppm.h90

    r14212 r14387  
    9595            fhat(1,ivar,+1) = & 
    9696        &   fdat(1,ivar,+1) 
    97             fhat(2,ivar,+1) = 0.d0 
    98             fhat(3,ivar,+1) = 0.d0 
     97            fhat(2,ivar,+1) = 0.e0 
     98            fhat(3,ivar,+1) = 0.e0 
    9999        end do 
    100100        end if 
     
    104104    !------------------- reconstruct function on each cell ! 
    105105 
    106         uhat = +0.d+0 
    107         lhat = +0.d+0 
     106        uhat = +0.e+0 
     107        lhat = +0.e+0 
    108108 
    109109        do  ipos = +1 , npos-1 
     
    154154    !----------------------------- pref. unlimited profile ! 
    155155 
    156             wval(1) = +1.d+0 
    157             wval(2) = +0.d+0 
     156            wval(1) = +1.e+0 
     157            wval(2) = +0.e+0 
    158158             
    159159            case (mono_limit) 
     
    167167    !----------------------------- pref. monotonic profile ! 
    168168 
    169             wval(1) = +0.d+0 
    170             wval(2) = +1.d+0 
     169            wval(1) = +0.e+0 
     170            wval(2) = +1.e+0 
    171171             
    172172            case (weno_limit) 
     
    190190    !----------------------------- pref. unlimited profile ! 
    191191 
    192             wval(1) = +1.d+0 
    193             wval(2) = +0.d+0 
     192            wval(1) = +1.e+0 
     193            wval(2) = +0.e+0 
    194194             
    195195            end if 
     
    269269         
    270270        if((ffrr - ff00) * &  
    271     &      (ff00 - ffll) .lt. 0.d+0) then 
     271    &      (ff00 - ffll) .lt. 0.e+0) then 
    272272 
    273273    !----------------------------------- "flatten" extrema ! 
     
    276276 
    277277            lhat(1) = ff00 
    278             lhat(2) = 0.d0 
    279             lhat(3) = 0.d0 
     278            lhat(2) = 0.e0 
     279            lhat(3) = 0.e0 
    280280              
    281281            return 
     
    286286     
    287287        if((ffll - fell) * & 
    288     &      (fell - ff00) .le. 0.d+0) then 
     288    &      (fell - ff00) .le. 0.e+0) then 
    289289 
    290290            mono = +1 
     
    295295 
    296296        if((ffrr - ferr) * & 
    297     &      (ferr - ff00) .le. 0.d+0) then 
     297    &      (ferr - ff00) .le. 0.e+0) then 
    298298 
    299299            mono = +1 
     
    322322    &                  / lhat(3) 
    323323 
    324         if ((turn .ge. -1.d+0)& 
    325     &  .and.(turn .le. +0.d+0)) then 
     324        if ((turn .ge. -1.e+0)& 
     325    &  .and.(turn .le. +0.e+0)) then 
    326326 
    327327        mono =   +2 
     
    342342 
    343343        else & 
    344     &   if ((turn .gt. +0.d+0)& 
    345     &  .and.(turn .le. +1.d+0)) then 
     344    &   if ((turn .gt. +0.e+0)& 
     345    &  .and.(turn .le. +1.e+0)) then 
    346346 
    347347        mono =   +2 
  • vendors/PPR/src/pqm.h90

    r14212 r14387  
    100100        do  ivar = +1, nvar 
    101101            fhat(1,ivar,+1) = fdat(1,ivar,+1) 
    102             fhat(2,ivar,+1) = 0.d0 
    103             fhat(3,ivar,+1) = 0.d0 
    104             fhat(4,ivar,+1) = 0.d0 
    105             fhat(5,ivar,+1) = 0.d0 
     102            fhat(2,ivar,+1) = 0.e0 
     103            fhat(3,ivar,+1) = 0.e0 
     104            fhat(4,ivar,+1) = 0.e0 
     105            fhat(5,ivar,+1) = 0.e0 
    106106        end do 
    107107        end if 
     
    169169    !----------------------------- pref. unlimited profile ! 
    170170             
    171             wval(1) = +1.d+0 
    172             wval(2) = +0.d+0 
     171            wval(1) = +1.e+0 
     172            wval(2) = +0.e+0 
    173173             
    174174            case (mono_limit) 
     
    183183    !----------------------------- pref. monotonic profile ! 
    184184             
    185             wval(1) = +0.d+0 
    186             wval(2) = +1.d+0 
     185            wval(1) = +0.e+0 
     186            wval(2) = +1.e+0 
    187187             
    188188            case (weno_limit) 
     
    207207    !----------------------------- pref. unlimited profile ! 
    208208             
    209             wval(1) = +1.d+0 
    210             wval(2) = +0.d+0 
     209            wval(1) = +1.e+0 
     210            wval(2) = +0.e+0 
    211211             
    212212            end if 
     
    286286         
    287287        uhat(1) = & 
    288     & + (30.d+0 / 16.d+0) * ff00 & 
    289     & - ( 7.d+0 / 16.d+0) *(ferr+fell) & 
    290     & + ( 1.d+0 / 16.d+0) *(derr-dell) 
     288    & + (30.e+0 / 16.e+0) * ff00 & 
     289    & - ( 7.e+0 / 16.e+0) *(ferr+fell) & 
     290    & + ( 1.e+0 / 16.e+0) *(derr-dell) 
    291291        uhat(2) = & 
    292     & + ( 3.d+0 /  4.d+0) *(ferr-fell) & 
    293     & - ( 1.d+0 /  4.d+0) *(derr+dell) 
     292    & + ( 3.e+0 /  4.e+0) *(ferr-fell) & 
     293    & - ( 1.e+0 /  4.e+0) *(derr+dell) 
    294294        uhat(3) = & 
    295     & - (30.d+0 /  8.d+0) * ff00 & 
    296     & + (15.d+0 /  8.d+0) *(ferr+fell) & 
    297     & - ( 3.d+0 /  8.d+0) *(derr-dell) 
     295    & - (30.e+0 /  8.e+0) * ff00 & 
     296    & + (15.e+0 /  8.e+0) *(ferr+fell) & 
     297    & - ( 3.e+0 /  8.e+0) *(derr-dell) 
    298298        uhat(4) = & 
    299     & - ( 1.d+0 /  4.d+0) *(ferr-fell  & 
     299    & - ( 1.e+0 /  4.e+0) *(ferr-fell  & 
    300300    &                      -derr-dell) 
    301301        uhat(5) = & 
    302     & + (30.d+0 / 16.d+0) * ff00 & 
    303     & - (15.d+0 / 16.d+0) *(ferr+fell) & 
    304     & + ( 5.d+0 / 16.d+0) *(derr-dell) 
     302    & + (30.e+0 / 16.e+0) * ff00 & 
     303    & - (15.e+0 / 16.e+0) *(ferr+fell) & 
     304    & + ( 5.e+0 / 16.e+0) *(derr-dell) 
    305305   
    306306    !-------------------------------- "mono" slope-limiter ! 
    307307     
    308308        if((ffrr - ff00) * &  
    309     &      (ff00 - ffll) .le. 0.d+0) then 
     309    &      (ff00 - ffll) .le. 0.e+0) then 
    310310 
    311311    !----------------------------------- "flatten" extrema ! 
     
    314314 
    315315            lhat(1) = ff00 
    316             lhat(2) = 0.d0 
    317             lhat(3) = 0.d0 
    318             lhat(4) = 0.d0 
    319             lhat(5) = 0.d0 
     316            lhat(2) = 0.e0 
     317            lhat(3) = 0.e0 
     318            lhat(4) = 0.e0 
     319            lhat(5) = 0.e0 
    320320              
    321321            return 
     
    326326 
    327327        if((ffll - fell) * & 
    328     &      (fell - ff00) .le. 0.d+0) then 
     328    &      (fell - ff00) .le. 0.e+0) then 
    329329 
    330330            mono = +1 
     
    334334        end if 
    335335 
    336         if (dell * dfds(0) .lt. 0.d+0) then 
     336        if (dell * dfds(0) .lt. 0.e+0) then 
    337337 
    338338            mono = +1 
     
    343343 
    344344        if((ffrr - ferr) * & 
    345     &      (ferr - ff00) .le. 0.d+0) then 
     345    &      (ferr - ff00) .le. 0.e+0) then 
    346346 
    347347            mono = +1 
     
    351351        end if 
    352352 
    353         if (derr * dfds(0) .lt. 0.d+0) then 
     353        if (derr * dfds(0) .lt. 0.e+0) then 
    354354 
    355355            mono = +1 
     
    362362     
    363363        lhat(1) = & 
    364     & + (30.d+0 / 16.d+0) * ff00 & 
    365     & - ( 7.d+0 / 16.d+0) *(ferr+fell) & 
    366     & + ( 1.d+0 / 16.d+0) *(derr-dell) 
     364    & + (30.e+0 / 16.e+0) * ff00 & 
     365    & - ( 7.e+0 / 16.e+0) *(ferr+fell) & 
     366    & + ( 1.e+0 / 16.e+0) *(derr-dell) 
    367367        lhat(2) = & 
    368     & + ( 3.d+0 /  4.d+0) *(ferr-fell) & 
    369     & - ( 1.d+0 /  4.d+0) *(derr+dell) 
     368    & + ( 3.e+0 /  4.e+0) *(ferr-fell) & 
     369    & - ( 1.e+0 /  4.e+0) *(derr+dell) 
    370370        lhat(3) = & 
    371     & - (30.d+0 /  8.d+0) * ff00 & 
    372     & + (15.d+0 /  8.d+0) *(ferr+fell) & 
    373     & - ( 3.d+0 /  8.d+0) *(derr-dell) 
     371    & - (30.e+0 /  8.e+0) * ff00 & 
     372    & + (15.e+0 /  8.e+0) *(ferr+fell) & 
     373    & - ( 3.e+0 /  8.e+0) *(derr-dell) 
    374374        lhat(4) = & 
    375     & - ( 1.d+0 /  4.d+0) *(ferr-fell  & 
     375    & - ( 1.e+0 /  4.e+0) *(ferr-fell  & 
    376376    &                      -derr-dell) 
    377377        lhat(5) = & 
    378     & + (30.d+0 / 16.d+0) * ff00 & 
    379     & - (15.d+0 / 16.d+0) *(ferr+fell) & 
    380     & + ( 5.d+0 / 16.d+0) *(derr-dell) 
     378    & + (30.e+0 / 16.e+0) * ff00 & 
     379    & - (15.e+0 / 16.e+0) *(ferr+fell) & 
     380    & + ( 5.e+0 / 16.e+0) *(derr-dell) 
    381381 
    382382    !------------------ calc. inflexion via 2nd-derivative ! 
    383383         
    384         call roots_2(12.d+0 * lhat(5), & 
    385     &                 6.d+0 * lhat(4), & 
    386     &                 2.d+0 * lhat(3), & 
     384        call roots_2(12.e+0 * lhat(5), & 
     385    &                 6.e+0 * lhat(4), & 
     386    &                 2.e+0 * lhat(3), & 
    387387    &                 iflx  , haveroot ) 
    388388 
     
    391391        turn = +0 
    392392 
    393         if ( ( iflx(1) .gt. -1.d+0 ) & 
    394     &  .and. ( iflx(1) .lt. +1.d+0 ) ) then 
     393        if ( ( iflx(1) .gt. -1.e+0 ) & 
     394    &  .and. ( iflx(1) .lt. +1.e+0 ) ) then 
    395395 
    396396    !------------------ check for non-monotonic inflection ! 
    397397 
    398398        grad = lhat(2) & 
    399     &+ (iflx(1)**1) * 2.d+0* lhat(3) & 
    400     &+ (iflx(1)**2) * 3.d+0* lhat(4) & 
    401     &+ (iflx(1)**3) * 4.d+0* lhat(5) 
    402  
    403         if (grad * dfds(0) .lt. 0.d+0) then 
     399    &+ (iflx(1)**1) * 2.e+0* lhat(3) & 
     400    &+ (iflx(1)**2) * 3.e+0* lhat(4) & 
     401    &+ (iflx(1)**3) * 4.e+0* lhat(5) 
     402 
     403        if (grad * dfds(0) .lt. 0.e+0) then 
    404404 
    405405            if (abs(dfds(-1)) & 
     
    418418        end if 
    419419             
    420         if ( ( iflx(2) .gt. -1.d+0 ) & 
    421     &  .and. ( iflx(2) .lt. +1.d+0 ) ) then 
     420        if ( ( iflx(2) .gt. -1.e+0 ) & 
     421    &  .and. ( iflx(2) .lt. +1.e+0 ) ) then 
    422422 
    423423    !------------------ check for non-monotonic inflection ! 
    424424                 
    425425        grad = lhat(2) & 
    426     &+ (iflx(2)**1) * 2.d+0* lhat(3) & 
    427     &+ (iflx(2)**2) * 3.d+0* lhat(4) & 
    428     &+ (iflx(2)**3) * 4.d+0* lhat(5) 
    429  
    430         if (grad * dfds(0) .lt. 0.d+0) then 
     426    &+ (iflx(2)**1) * 2.e+0* lhat(3) & 
     427    &+ (iflx(2)**2) * 3.e+0* lhat(4) & 
     428    &+ (iflx(2)**3) * 4.e+0* lhat(5) 
     429 
     430        if (grad * dfds(0) .lt. 0.e+0) then 
    431431 
    432432            if (abs(dfds(-1)) & 
     
    454454 
    455455        derr = & 
    456     &- ( 5.d+0 /  1.d+0) * ff00 & 
    457     &+ ( 3.d+0 /  1.d+0) * ferr & 
    458     &+ ( 2.d+0 /  1.d+0) * fell 
     456    &- ( 5.e+0 /  1.e+0) * ff00 & 
     457    &+ ( 3.e+0 /  1.e+0) * ferr & 
     458    &+ ( 2.e+0 /  1.e+0) * fell 
    459459        dell = & 
    460     &+ ( 5.d+0 /  3.d+0) * ff00 & 
    461     &- ( 1.d+0 /  3.d+0) * ferr & 
    462     &- ( 4.d+0 /  3.d+0) * fell 
    463  
    464         if (dell*dfds(+0) .lt. 0.d+0) then 
    465  
    466         dell =   0.d+0 
     460    &+ ( 5.e+0 /  3.e+0) * ff00 & 
     461    &- ( 1.e+0 /  3.e+0) * ferr & 
     462    &- ( 4.e+0 /  3.e+0) * fell 
     463 
     464        if (dell*dfds(+0) .lt. 0.e+0) then 
     465 
     466        dell =   0.e+0 
    467467 
    468468        ferr = & 
    469     &+ ( 5.d+0 /  1.d+0) * ff00 & 
    470     &- ( 4.d+0 /  1.d+0) * fell 
     469    &+ ( 5.e+0 /  1.e+0) * ff00 & 
     470    &- ( 4.e+0 /  1.e+0) * fell 
    471471        derr = & 
    472     &+ (10.d+0 /  1.d+0) * ff00 & 
    473     &- (10.d+0 /  1.d+0) * fell 
     472    &+ (10.e+0 /  1.e+0) * ff00 & 
     473    &- (10.e+0 /  1.e+0) * fell 
    474474 
    475475        else & 
    476     &   if (derr*dfds(+0) .lt. 0.d+0) then 
    477  
    478         derr =   0.d+0 
     476    &   if (derr*dfds(+0) .lt. 0.e+0) then 
     477 
     478        derr =   0.e+0 
    479479 
    480480        fell = & 
    481     &+ ( 5.d+0 /  2.d+0) * ff00 & 
    482     &- ( 3.d+0 /  2.d+0) * ferr 
     481    &+ ( 5.e+0 /  2.e+0) * ff00 & 
     482    &- ( 3.e+0 /  2.e+0) * ferr 
    483483        dell = & 
    484     &- ( 5.d+0 /  3.d+0) * ff00 & 
    485     &+ ( 5.d+0 /  3.d+0) * ferr 
     484    &- ( 5.e+0 /  3.e+0) * ff00 & 
     485    &+ ( 5.e+0 /  3.e+0) * ferr 
    486486 
    487487        end if 
    488488 
    489489        lhat(1) = & 
    490     &+ (30.d+0 / 16.d+0) * ff00 & 
    491     &- ( 7.d+0 / 16.d+0) *(ferr+fell) & 
    492     &+ ( 1.d+0 / 16.d+0) *(derr-dell) 
     490    &+ (30.e+0 / 16.e+0) * ff00 & 
     491    &- ( 7.e+0 / 16.e+0) *(ferr+fell) & 
     492    &+ ( 1.e+0 / 16.e+0) *(derr-dell) 
    493493        lhat(2) = & 
    494     &+ ( 3.d+0 /  4.d+0) *(ferr-fell) & 
    495     &- ( 1.d+0 /  4.d+0) *(derr+dell) 
     494    &+ ( 3.e+0 /  4.e+0) *(ferr-fell) & 
     495    &- ( 1.e+0 /  4.e+0) *(derr+dell) 
    496496        lhat(3) = & 
    497     &- (30.d+0 /  8.d+0) * ff00 & 
    498     &+ (15.d+0 /  8.d+0) *(ferr+fell) & 
    499     &- ( 3.d+0 /  8.d+0) *(derr-dell) 
     497    &- (30.e+0 /  8.e+0) * ff00 & 
     498    &+ (15.e+0 /  8.e+0) *(ferr+fell) & 
     499    &- ( 3.e+0 /  8.e+0) *(derr-dell) 
    500500        lhat(4) = & 
    501     &- ( 1.d+0 /  4.d+0) *(ferr-fell  & 
     501    &- ( 1.e+0 /  4.e+0) *(ferr-fell  & 
    502502    &                     -derr-dell) 
    503503        lhat(5) = & 
    504     &+ (30.d+0 / 16.d+0) * ff00 & 
    505     &- (15.d+0 / 16.d+0) *(ferr+fell) & 
    506     &+ ( 5.d+0 / 16.d+0) *(derr-dell) 
     504    &+ (30.e+0 / 16.e+0) * ff00 & 
     505    &- (15.e+0 / 16.e+0) *(ferr+fell) & 
     506    &+ ( 5.e+0 / 16.e+0) *(derr-dell) 
    507507 
    508508        end if 
     
    515515 
    516516        derr = & 
    517     &- ( 5.d+0 /  3.d+0) * ff00 & 
    518     &+ ( 4.d+0 /  3.d+0) * ferr & 
    519     &+ ( 1.d+0 /  3.d+0) * fell 
     517    &- ( 5.e+0 /  3.e+0) * ff00 & 
     518    &+ ( 4.e+0 /  3.e+0) * ferr & 
     519    &+ ( 1.e+0 /  3.e+0) * fell 
    520520        dell = & 
    521     &+ ( 5.d+0 /  1.d+0) * ff00 & 
    522     &- ( 2.d+0 /  1.d+0) * ferr & 
    523     &- ( 3.d+0 /  1.d+0) * fell 
    524  
    525         if (dell*dfds(+0) .lt. 0.d+0) then 
    526  
    527         dell =   0.d+0 
     521    &+ ( 5.e+0 /  1.e+0) * ff00 & 
     522    &- ( 2.e+0 /  1.e+0) * ferr & 
     523    &- ( 3.e+0 /  1.e+0) * fell 
     524 
     525        if (dell*dfds(+0) .lt. 0.e+0) then 
     526 
     527        dell =   0.e+0 
    528528 
    529529        ferr = & 
    530     &+ ( 5.d+0 /  2.d+0) * ff00 & 
    531     &- ( 3.d+0 /  2.d+0) * fell 
     530    &+ ( 5.e+0 /  2.e+0) * ff00 & 
     531    &- ( 3.e+0 /  2.e+0) * fell 
    532532        derr = & 
    533     &+ ( 5.d+0 /  3.d+0) * ff00 & 
    534     &- ( 5.d+0 /  3.d+0) * fell 
     533    &+ ( 5.e+0 /  3.e+0) * ff00 & 
     534    &- ( 5.e+0 /  3.e+0) * fell 
    535535 
    536536        else &  
    537     &   if (derr*dfds(+0) .lt. 0.d+0) then 
    538  
    539         derr =   0.d+0 
     537    &   if (derr*dfds(+0) .lt. 0.e+0) then 
     538 
     539        derr =   0.e+0 
    540540 
    541541        fell = & 
    542     &+ ( 5.d+0 /  1.d+0) * ff00 & 
    543     &- ( 4.d+0 /  1.d+0) * ferr 
     542    &+ ( 5.e+0 /  1.e+0) * ff00 & 
     543    &- ( 4.e+0 /  1.e+0) * ferr 
    544544        dell = & 
    545     &- (10.d+0 /  1.d+0) * ff00 & 
    546     &+ (10.d+0 /  1.d+0) * ferr 
     545    &- (10.e+0 /  1.e+0) * ff00 & 
     546    &+ (10.e+0 /  1.e+0) * ferr 
    547547 
    548548        end if 
    549549 
    550550        lhat(1) = & 
    551     &+ (30.d+0 / 16.d+0) * ff00 & 
    552     &- ( 7.d+0 / 16.d+0) *(ferr+fell) & 
    553     &+ ( 1.d+0 / 16.d+0) *(derr-dell) 
     551    &+ (30.e+0 / 16.e+0) * ff00 & 
     552    &- ( 7.e+0 / 16.e+0) *(ferr+fell) & 
     553    &+ ( 1.e+0 / 16.e+0) *(derr-dell) 
    554554        lhat(2) = & 
    555     &+ ( 3.d+0 /  4.d+0) *(ferr-fell) & 
    556     &- ( 1.d+0 /  4.d+0) *(derr+dell) 
     555    &+ ( 3.e+0 /  4.e+0) *(ferr-fell) & 
     556    &- ( 1.e+0 /  4.e+0) *(derr+dell) 
    557557        lhat(3) = & 
    558     &- (30.d+0 /  8.d+0) * ff00 & 
    559     &+ (15.d+0 /  8.d+0) *(ferr+fell) & 
    560     &- ( 3.d+0 /  8.d+0) *(derr-dell) 
     558    &- (30.e+0 /  8.e+0) * ff00 & 
     559    &+ (15.e+0 /  8.e+0) *(ferr+fell) & 
     560    &- ( 3.e+0 /  8.e+0) *(derr-dell) 
    561561        lhat(4) = & 
    562     &- ( 1.d+0 /  4.d+0) *(ferr-fell  & 
     562    &- ( 1.e+0 /  4.e+0) *(ferr-fell  & 
    563563    &                     -derr-dell) 
    564564        lhat(5) = & 
    565     &+ (30.d+0 / 16.d+0) * ff00 & 
    566     &- (15.d+0 / 16.d+0) *(ferr+fell) & 
    567     &+ ( 5.d+0 / 16.d+0) *(derr-dell) 
     565    &+ (30.e+0 / 16.e+0) * ff00 & 
     566    &- (15.e+0 / 16.e+0) *(ferr+fell) & 
     567    &+ ( 5.e+0 / 16.e+0) *(derr-dell) 
    568568 
    569569        end if 
  • vendors/PPR/src/rmap1d.h90

    r14212 r14387  
    7676        &   intent(inout) , optional :: tCPU 
    7777 
    78         real*8 , parameter :: RTOL = +1.d-14 
     78        real*8 , parameter :: RTOL = +1.e-14 
    7979 
    8080    !------------------------------------------- variables ! 
     
    245245    !------------------------------------- initializations ! 
    246246 
    247         vvlo(+1:+5) = 0.d0 
    248         vvhi(+1:+5) = 0.d0 
     247        vvlo(+1:+5) = 0.e0 
     248        vvhi(+1:+5) = 0.e0 
    249249 
    250250    !------------- remap FDAT from XPOS to XNEW using FHAT ! 
     
    286286            do  ivar = +1,nvar 
    287287                 
    288                 fnew(idof,ivar,kpos) = 0.d0 
     288                fnew(idof,ivar,kpos) = 0.e0 
    289289             
    290290            end do 
     
    375375    !   account for FP roundoff. 
    376376 
    377         sdat = 0.d0; serr = 0.d0 
     377        sdat = 0.e0; serr = 0.e0 
    378378        do  ipos = +1, npos-1 
    379379        do  ivar = +1, nvar-0 
     
    408408        sdat =  sdat + serr 
    409409 
    410         snew = 0.d0; serr = 0.d0 
     410        snew = 0.e0; serr = 0.e0 
    411411        do  ipos = +1, nnew-1 
    412412        do  ivar = +1, nvar-0 
     
    451451        do  ivar = +1, nvar-0 
    452452 
    453             if (serr(ivar) .gt. 0.d0) then 
     453            if (serr(ivar) .gt. 0.e0) then 
    454454 
    455455                vmin = kmin(ivar) 
     
    460460 
    461461            else & 
    462         &   if (serr(ivar) .lt. 0.d0) then 
     462        &   if (serr(ivar) .lt. 0.e0) then 
    463463 
    464464                vmax = kmax(ivar) 
  • vendors/PPR/src/root1d.h90

    r14212 r14387  
    5252        real*8 :: sq,ia,a0,b0,c0,x0 
    5353 
    54         real*8, parameter :: rt = +1.d-14 
     54        real*8, parameter :: rt = +1.e-14 
    5555 
    5656        a0 = abs(aa) 
  • vendors/PPR/src/util1d.h90

    r14212 r14387  
    111111            call random_number (rand) 
    112112 
    113             rand = 2.d0 * (rand-.5d0) 
     113            rand = 2.e0 * (rand-.5d0) 
    114114             
    115115            move = rand *  move 
  • vendors/PPR/src/weno1d.h90

    r14212 r14387  
    7171        real*8 :: omin,omax,wsum 
    7272         
    73         real*8 , parameter :: ZERO = +1.d-16 
     73        real*8 , parameter :: ZERO = +1.e-16 
    7474        
    7575        if (size(delx).gt.+1) then 
     
    112112        wval(1) = wval(1) / wsum 
    113113    !   wval(2) = wval(2) / wsum 
    114         wval(2) =-wval(1) + 1.d0 ! wval(2)/wsum but robust ! 
     114        wval(2) =-wval(1) + 1.e0 ! wval(2)/wsum but robust ! 
    115115 
    116116        return 
     
    169169        if (ipos-halo.lt.head) then 
    170170         
    171             omax = 1.d0  
    172             omin = 0.d0 ; return 
     171            omax = 1.e0  
     172            omin = 0.e0 ; return 
    173173         
    174174        end if 
     
    176176        if (ipos+halo.gt.tail) then 
    177177             
    178             omax = 1.d0  
    179             omin = 0.d0 ; return 
     178            omax = 1.e0  
     179            omin = 0.e0 ; return 
    180180         
    181181        end if 
     
    204204    !---------------------------------------- "lower" part ! 
    205205       
    206         delh = 0.d0 
     206        delh = 0.e0 
    207207 
    208208        do  hpos = ipos-1, imin, -1 
     
    236236    !---------------------------------------- "upper" part !      
    237237       
    238         delh = 0.d0 
     238        delh = 0.e0 
    239239         
    240240        do  hpos = ipos+1, imax, +1 
     
    320320        if (ipos-halo.lt.head) then 
    321321         
    322             omax = 1.d0  
    323             omin = 0.d0 ; return 
     322            omax = 1.e0  
     323            omin = 0.e0 ; return 
    324324         
    325325        end if 
     
    327327        if (ipos+halo.gt.tail) then 
    328328             
    329             omax = 1.d0  
    330             omin = 0.d0 ; return 
     329            omax = 1.e0  
     330            omin = 0.e0 ; return 
    331331         
    332332        end if 
     
    344344        dfx2 = oscl(2,ivar,ipos) 
    345345       
    346         oval = (2.d0**1*dfx1)**2 & 
    347         &    + (2.d0**2*dfx2)**2 
     346        oval = (2.e0**1*dfx1)**2 & 
     347        &    + (2.e0**2*dfx2)**2 
    348348       
    349349        omin = oval 
     
    352352    !---------------------------------------- "lower" part ! 
    353353       
    354         delh = 0.d0 
     354        delh = 0.e0 
    355355 
    356356        do  hpos = ipos-1, imin, -1 
     
    358358    !------------------ calc. derivatives centred on IPOS. !      
    359359 
    360         delh = delh + 2.d0 
     360        delh = delh + 2.e0 
    361361         
    362362        dfx1 = oscl(1,ivar,hpos) 
     
    367367    !------------------ indicator: NORM(H^N * D^N/DX^N(F)) ! 
    368368 
    369         oval = (2.d0**1*dfx1)**2 & 
    370     &        + (2.d0**2*dfx2)**2 
     369        oval = (2.e0**1*dfx1)**2 & 
     370    &        + (2.e0**2*dfx2)**2 
    371371 
    372372        if (oval .lt. omin) then 
     
    381381    !---------------------------------------- "upper" part !      
    382382       
    383         delh = 0.d0 
     383        delh = 0.e0 
    384384         
    385385        do  hpos = ipos+1, imax, +1 
     
    387387    !------------------ calc. derivatives centred on IPOS. !      
    388388 
    389         delh = delh - 2.d0 
     389        delh = delh - 2.e0 
    390390         
    391391        dfx1 = oscl(1,ivar,hpos) 
     
    396396    !------------------ indicator: NORM(H^N * D^N/DX^N(F)) ! 
    397397 
    398         oval = (2.d0**1*dfx1)**2 & 
    399     &        + (2.d0**2*dfx2)**2 
     398        oval = (2.e0**1*dfx1)**2 & 
     399    &        + (2.e0**2*dfx2)**2 
    400400   
    401401        if (oval .lt. omin) then 
Note: See TracChangeset for help on using the changeset viewer.