Examples/calcium1/code1.f

   1 C Copyright (C) 2009-2012  EDF R&D
   2 C
   3 C This library is free software; you can redistribute it and/or
   4 C modify it under the terms of the GNU Lesser General Public
   5 C License as published by the Free Software Foundation; either
   6 C version 2.1 of the License.
   7 C
   8 C This library is distributed in the hope that it will be useful,
   9 C but WITHOUT ANY WARRANTY; without even the implied warranty of
  10 C MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11 C Lesser General Public License for more details.
  12 C
  13 C You should have received a copy of the GNU Lesser General Public
  14 C License along with this library; if not, write to the Free Software
  15 C Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  16 C
  17 C See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  18 C
  19
  20        SUBROUTINE SERV1(compo,a,b,c)
  21        include 'calcium.hf'
  22        integer compo(2)
  23        integer i, nval, info, z(10), l
  24        integer*8 lz(10)
  25        integer*4 z4(10),lo(10)
  26
  27        real*8 dd(10),a,b,c,ti,tf,t
  28        real*4 u(20)
  29        real*4 tti,ttf,tt
  30        character*10 s(3)
  31        character*20 rs(3)
  32        character*64 instance
  33
  34        write(6,*)a,b
  35        call cpcd(compo,instance,info)
  36        write(6,*)"instance name=",instance
  37
  38 C  write
  39        tt=0.
  40        t=0.
  41        dd(1)=125.45
  42        dd(2)=8.8
  43        i=1
  44        l=10
  45        CALL cpedb(compo,CP_ITERATION,t,i,'ba',l,dd,info)
  46        CALL cpedb(compo,CP_ITERATION,t,2,'ba',l,dd,info)
  47        CALL cpedb(compo,CP_ITERATION,t,3,'ba',l,dd,info)
  48        write(6,*)'info=',info
  49        call flush(6)
  50
  51        s(1)="titi"
  52        s(2)="tututu"
  53        s(3)="tatatata"
  54        write(6,*)'s=',s
  55        l=3
  56        CALL cpech(compo,CP_ITERATION,tt,i,'bb',l,s,info)
  57        write(6,*)'info=',info
  58        call flush(6)
  59
  60        z(1)=1
  61        z(2)=8
  62        z(3)=0
  63        write(6,*)'z=',z(1)
  64        write(6,*)'z=',z(2)
  65        write(6,*)'z=',z(3)
  66        l=10
  67        CALL cpeen(compo,CP_ITERATION,tt,i,'bc',l,z,info)
  68        write(6,*)'info=',info
  69        call flush(6)
  70
  71        u(1)=1
  72        u(2)=8
  73        u(3)=4
  74        u(4)=4
  75        u(5)=5
  76        u(6)=5
  77        write(6,*)'u=',u(1)
  78        write(6,*)'u=',u(2)
  79        write(6,*)'u=',u(3)
  80        write(6,*)'u=',u(4)
  81        write(6,*)'u=',u(5)
  82        write(6,*)'u=',u(6)
  83        CALL cpecp(compo,CP_ITERATION,tt,i,'bd',l,u,info)
  84        write(6,*)'info=',info
  85        call flush(6)
  86
  87        u(1)=1.1
  88        u(2)=8.8
  89        u(3)=4.4
  90        write(6,*)'u=',u(1)
  91        write(6,*)'u=',u(2)
  92        write(6,*)'u=',u(3)
  93        CALL cpere(compo,CP_ITERATION,tt,i,'be',l,u,info)
  94        write(6,*)'info=',info
  95        call flush(6)
  96
  97        lo(1)=1
  98        lo(2)=0
  99        lo(3)=1
 100        write(6,*)'lo=',lo(1)
 101        write(6,*)'lo=',lo(2)
 102        write(6,*)'lo=',lo(3)
 103        CALL cpelo(compo,CP_ITERATION,tt,i,'bf',l,lo,info)
 104        write(6,*)'info=',info
 105        call flush(6)
 106
 107        lz(1)=11
 108        lz(2)=22
 109        lz(3)=33
 110        write(6,*)'lz=',lz(1)
 111        write(6,*)'lz=',lz(2)
 112        write(6,*)'lz=',lz(3)
 113        CALL cpeln(compo,CP_ITERATION,tt,i,'bg',l,lz,info)
 114        write(6,*)'info=',info
 115        call flush(6)
 116
 117        z4(1)=1
 118        z4(2)=8
 119        z4(3)=0
 120        write(6,*)'z4=',z4(1)
 121        write(6,*)'z4=',z4(2)
 122        write(6,*)'z4=',z4(3)
 123        CALL cpein(compo,CP_ITERATION,tt,i,'bh',l,z4,info)
 124        write(6,*)'info=',info
 125        call flush(6)
 126
 127        lz(1)=11
 128        lz(2)=22
 129        lz(3)=2**30
 130        lz(3)=2**20*lz(3)
 131        write(6,*)'lz=',lz(1)
 132        write(6,*)'lz=',lz(2)
 133        write(6,*)'lz=',lz(3)
 134        CALL cpelg(compo,CP_ITERATION,tt,i,'bi',l,lz,info)
 135        write(6,*)'info=',info
 136        call flush(6)
 137
 138 C  read
 139        ti=0.
 140        tf=1.
 141        i=1
 142        dd(1)=0.
 143        dd(2)=0.
 144        dd(3)=0.
 145        l=3
 146        CALL cpldb(compo,CP_ITERATION,ti,tf,i,'aa',l,nval,dd,info)
 147        write(6,*)'info=',info
 148        write(6,*)'dd=',dd(1)
 149        write(6,*)'dd=',dd(2)
 150        write(6,*)'dd=',dd(3)
 151        write(6,*)'nval=',nval
 152        call flush(6)
 153
 154        tti=0.
 155        ttf=1.
 156        i=1
 157        CALL cplch(compo,CP_ITERATION,tti,ttf,i,'ab',l,nval,rs,info)
 158        write(6,*)'info=',info
 159        write(6,*)'rs=',rs
 160        write(6,*)'nval=',nval
 161        call flush(6)
 162
 163        z(1)=0
 164        z(2)=0
 165        z(3)=0
 166        CALL cplen(compo,CP_ITERATION,tti,ttf,i,'ac',l,nval,z,info)
 167        write(6,*)'info=',info
 168        write(6,*)'nval=',nval
 169        write(6,*)'z=',z(1)
 170        write(6,*)'z=',z(2)
 171        write(6,*)'z=',z(3)
 172        call flush(6)
 173
 174        u(1)=0
 175        u(2)=0
 176        u(3)=0
 177        u(4)=0
 178        u(5)=0
 179        u(6)=0
 180        CALL cplcp(compo,CP_ITERATION,tti,ttf,i,'ad',l,nval,u,info)
 181        write(6,*)'info=',info
 182        write(6,*)'nval=',nval
 183        write(6,*)'u=',u(1)
 184        write(6,*)'u=',u(2)
 185        write(6,*)'u=',u(3)
 186        write(6,*)'u=',u(4)
 187        write(6,*)'u=',u(5)
 188        write(6,*)'u=',u(6)
 189        call flush(6)
 190
 191        u(1)=0
 192        u(2)=0
 193        u(3)=0
 194        CALL cplre(compo,CP_ITERATION,tti,ttf,i,'ae',l,nval,u,info)
 195        write(6,*)'info=',info
 196        write(6,*)'nval=',nval
 197        write(6,*)'u=',u(1)
 198        write(6,*)'u=',u(2)
 199        write(6,*)'u=',u(3)
 200        call flush(6)
 201
 202        lo(1)=0
 203        lo(2)=0
 204        lo(3)=0
 205        CALL cpllo(compo,CP_ITERATION,tti,ttf,i,'af',l,nval,lo,info)
 206        write(6,*)'info=',info
 207        write(6,*)'nval=',nval
 208        write(6,*)'lo=',lo(1)
 209        write(6,*)'lo=',lo(2)
 210        write(6,*)'lo=',lo(3)
 211        call flush(6)
 212
 213        lz(1)=0
 214        lz(2)=0
 215        lz(3)=0
 216        CALL cplln(compo,CP_ITERATION,tti,ttf,i,'ag',l,nval,lz,info)
 217        write(6,*)'info=',info
 218        write(6,*)'nval=',nval
 219        write(6,*)'lz=',lz(1)
 220        write(6,*)'lz=',lz(2)
 221        write(6,*)'lz=',lz(3)
 222        call flush(6)
 223
 224        z4(1)=0
 225        z4(2)=0
 226        z4(3)=0
 227        CALL cplin(compo,CP_ITERATION,tti,ttf,i,'ah',l,nval,z4,info)
 228        write(6,*)'info=',info
 229        write(6,*)'nval=',nval
 230        write(6,*)'z4=',z4(1)
 231        write(6,*)'z4=',z4(2)
 232        write(6,*)'z4=',z4(3)
 233        call flush(6)
 234
 235        lz(1)=0
 236        lz(2)=0
 237        lz(3)=0
 238        CALL cpllg(compo,CP_ITERATION,tti,ttf,i,'ai',l,nval,lz,info)
 239        write(6,*)'info=',info
 240        write(6,*)'nval=',nval
 241        write(6,*)'lz=',lz(1)
 242        write(6,*)'lz=',lz(2)
 243        write(6,*)'lz=',lz(3)
 244        call flush(6)
 245
 246        call cpfini(compo,'aa',1,info)
 247        write(6,*)'info=',info
 248        call flush(6)
 249
 250        i=1
 251        l=3
 252        CALL cpldb(compo,CP_ITERATION,ti,tf,i,'aa',l,nval,dd,info)
 253        write(6,*)'info=',info
 254        call flush(6)
 255
 256        call cpeffi(compo,'aa',3,info)
 257        write(6,*)'info=',info
 258        call flush(6)
 259
 260        i=3
 261        l=3
 262        CALL cpldb(compo,CP_ITERATION,ti,tf,i,'aa',l,nval,dd,info)
 263        write(6,*)'info=',info
 264        call flush(6)
 265
 266        CALL cpfin(compo,CP_ARRET,info)
 267
 268        c=a+b
 269        return
 270        end