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

xuyufufzu

#IF Tap
* U! J9 O6 q# @, ^ #STORAGE STOR:26
0 ]* q% V5 h- e, r+ w #IF SAT==TAP
  #LOCAL REAL RVD1_1
% W' `" z6 A; u6 @. D- i      RVD1_1 = $TapI
7 ^! }+ D: U8 i: E+ J$ _, @      IF (RVD1_1 .LT. 1.0E-8) RVD1_1 = 1.0
#ENDIF
#ELSE
+ I& W) H9 T5 Z/ k) O5 T) R #STORAGE STOR:21
#ENDIF
!  TRANSFORMER SATURATION SUBROUTINE
#SUBROUTINE TSAT21 Transformer Saturation Subroutine
8 L6 y; E: c0 n5 i" J#LOCAL REAL X_KNEE
7 S) L/ X  W8 N+ x7 G4 F# d! |#IF Txk>0
# f  l: R- f& I8 e0 g$ {5 N; W      IF (TIME.GT.$Txk) THEN
) c& E1 u5 x3 J         X_KNEE = $Xknee
      ELSE
- T& b* K% ?' j( r4 x* }& g  ?; R, y         X_KNEE = -$Xknee
      ENDIF
#ELSE
      X_KNEE = $Xknee
#ENDIF
#LOCAL REAL TURNS_V
* K3 _! I- _* c+ a#IF SAT==1
- Z$ m8 e: {" L6 R) J  F      TURNS_V = $V1
#ELSE
( l7 k* [! L% R      TURNS_V = $V2
#ENDIF
#IF SAT==TAP
     + * RVD1_1
7 }) {3 H9 P/ |% x2 u2 f( O#ENDIF
9 i/ F# x; Q/ R9 O1 }1 B/ Y#IF SAT==1
#IF CuL>0.0
5 V6 ]/ U, c1 Q0 Y' y8 u  c& Y! a& I      CALL TSAT21($X1,$B1,
#ELSE
8 c& b6 u1 O5 n5 @/ b" C0 S# Q+ }      CALL TSAT21($A1,$B1,
#ENDIF
#ELSE
" b/ F3 ~) R& z* ~ #IF CuL>0.0
9 d4 S( b; w9 w: k# p4 A0 l; @' |      CALL TSAT21($X2,$B2,
#ELSE
      CALL TSAT21($A2,$B2,
#ENDIF
" h: g0 k! Q' y" u2 ~1 @" k#ENDIF
! N, E0 q# A5 H; V     + 0,0, 0,0, $SS,$Tmva,TURNS_V,$Xair,X_KNEE,$f,$Tdc,$I_m1, 0.0, $Enab)
; a0 s, ^5 d* k% d; X' R#OUTPUT REAL IMAG {-STOR(NEXC-21+17)}
3 p: d* p# ~& C5 h8 j#OUTPUT REAL FLUX {STOR(NEXC-21+11)}
#IF Tap
8 c% E% |1 O0 }: D. C" {: e1 i. B      CALL TRTAP5($#[1], $Tap, $SS, $TapI)
, _6 M9 P) Z* V$ P$ l+ v#ENDIF
!

风残冷月

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

drzyx

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

suner1986

什么意思啊，看不懂fortran

sgkevin

完全看不懂。。

小小鸟儿

这是什么程序里面的?

huafeinaojin

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

buddam

要学习学习了

gpbaby199

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