シミュレーションの演習のコード

GS9.BGSのコードを以下に示します。このコードは BGSP用です。第1行の

BOND SYMBOL EXPR OUT DEFIN PARM INIT

はコードの一部ですから省略してはいけません。また、BGSPではカラム指定が厳格で1カラムでもずれるとエラーになります。また、許される文字は英数字、コロン(:)、スラッシュ(/)、加減乗算記号(+-*)、半角空白のみであり、制御文字が入るとエラーになります。倍角の空白はエラーを生じますが、ソースをプリントしても気づきにくいので注意が必要です。

GS9のコード

BOND SYMBOL EXPR OUT DEFIN PARM INIT
SF1    11
ZJ     -11   13    14
ZJ     -14   19    12
CI1    13
TF1    19    20
OJ     -20   22    23    24
II1    24
CI2    23
RE1    22
RE2    12
/
BC: P13   EF13
BC: Q12   FL12
/
SF1
PA: PE1    400.0           inflow
      F=PE1
EOD
CI1
PA: PC1    1.0             capacitive load
      C=Z/PC1
EOD
TF1
PA: PT1    2.0
      TF=PT1
EOD
RE2
PA: PR2    600.0
FB: XXR    DP23
RV: XXF1 
      IF(XXR.LT.0.0) THEN
       XXF1=0.0
      ELSE
       XXF1=XXR
      ENDIF
      R=Z*PR2*XXF1
EOD
II1
PA: PM1    3.0D-06         piston mass
      L=Z/PM1
EOD
CI2
PA: PC2    1.0
      C=Z/PC2
       IF(C.LT.0.0) THEN
        C=200.0*C+5.0
       ELSE
        IF(C.LT.0.4) THEN
         C=C+5.0
        ELSE
          C=200.0*C-74.6
        ENDIF
       ENDIF
EOD
RE1
PA: PR1    0.01
       R=Z*PR1
EOD
/
0.00              0.10              0.0001             1000
/