关于变压器代码，高手必看？

xuyufufzu

#IF Tap
#STORAGE STOR:26
#IF SAT==TAP
0 J1 M$ {$ h( \1 G  #LOCAL REAL RVD1_1
      RVD1_1 = $TapI
' Z5 _, D6 z' i( A      IF (RVD1_1 .LT. 1.0E-8) RVD1_1 = 1.0
. _+ }$ k9 ]. Z  e! b #ENDIF
#ELSE
4 T9 a* T5 m$ v! ~! O1 C" ^$ b #STORAGE STOR:21
& f2 S* F  o- H' n- Q& K, a6 x) W2 A1 h#ENDIF
!  TRANSFORMER SATURATION SUBROUTINE
; j1 Z) i# R) y4 d#SUBROUTINE TSAT21 Transformer Saturation Subroutine
2 d5 q& E' S1 B- A+ G$ r#LOCAL REAL X_KNEE
! `  L  F. |6 T0 A#IF Txk>0
+ B* d* q- b2 N$ j5 H      IF (TIME.GT.$Txk) THEN
% x6 e: o  k* r- f$ {- }         X_KNEE = $Xknee
0 R8 i# C3 b, F) f: D8 s! \, l      ELSE
         X_KNEE = -$Xknee
      ENDIF
#ELSE
9 `- Z1 @2 t& x- j/ x      X_KNEE = $Xknee
#ENDIF
, c% d. a8 x  _1 \- O; e0 b& e. G#LOCAL REAL TURNS_V
( ~) O  V8 q. l7 O$ B#IF SAT==1
      TURNS_V = $V1
0 t, O" ?9 K# z$ n3 o#ELSE
      TURNS_V = $V2
' a( M  W: O2 K* j" j6 g#ENDIF
/ }, Z3 r: _4 m  I4 {) X. T( b8 e#IF SAT==TAP
1 w5 L4 w2 Y0 D0 p     + * RVD1_1
/ b( P7 f3 a; v" Y#ENDIF
/ L$ D. t% _% z# B2 m2 I#IF SAT==1
#IF CuL>0.0
% i; \( F" i- x" R. G      CALL TSAT21($X1,$B1,
#ELSE
      CALL TSAT21($A1,$B1,
2 N2 @. N* |$ i/ b; I0 B8 D #ENDIF
: q9 I# x& n  f% q& k  Q) |% `#ELSE
4 y7 |. M$ Q, } #IF CuL>0.0
      CALL TSAT21($X2,$B2,
#ELSE
2 H; p4 c9 ^& K0 e& |      CALL TSAT21($A2,$B2,
6 S% o, s& A: H6 Z+ \ #ENDIF
; Q; w# C$ [" l6 V#ENDIF
' Y  K) p& w2 G     + 0,0, 0,0, $SS,$Tmva,TURNS_V,$Xair,X_KNEE,$f,$Tdc,$I_m1, 0.0, $Enab)
#OUTPUT REAL IMAG {-STOR(NEXC-21+17)}
3 k' n1 m8 n8 |7 f#OUTPUT REAL FLUX {STOR(NEXC-21+11)}
#IF Tap
      CALL TRTAP5($#[1], $Tap, $SS, $TapI)
#ENDIF
!

风残冷月

俺也是低手，过来学习一下

drzyx

看不懂，不过还是顶一下。

suner1986

什么意思啊，看不懂fortran

sgkevin

完全看不懂。。

小小鸟儿

这是什么程序里面的?

huafeinaojin

看英文help里的内容，可能会有帮助

buddam

要学习学习了

gpbaby199

这个是变压器里面script里面DSDYN的饱和子程序吧？