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

xuyufufzu

#IF Tap
4 ~* O$ P* s) x( N #STORAGE STOR:26
4 g' l% X+ `) e #IF SAT==TAP
/ ^" L! a8 C/ a" G2 N  #LOCAL REAL RVD1_1
# B* S& W5 i4 m, e      RVD1_1 = $TapI
      IF (RVD1_1 .LT. 1.0E-8) RVD1_1 = 1.0
. N6 T- Y0 a2 G3 n( M: q #ENDIF
#ELSE
9 d2 B9 [* i  U& ?$ n5 w% s1 R/ v #STORAGE STOR:21
7 [3 L& R4 y  k  L2 R" e, X#ENDIF
: P3 ?' W9 L% e6 T!  TRANSFORMER SATURATION SUBROUTINE
#SUBROUTINE TSAT21 Transformer Saturation Subroutine
" A1 E* p" }5 E) S' k6 M5 n#LOCAL REAL X_KNEE
/ H* d* d! U* ~" |- }4 K: ]#IF Txk>0
+ ^9 A1 v  D( @3 O6 F      IF (TIME.GT.$Txk) THEN
         X_KNEE = $Xknee
      ELSE
         X_KNEE = -$Xknee
7 B0 }# d% R5 n/ k4 P8 E! C      ENDIF
2 w5 I! Q" L5 w6 M#ELSE
! `: S" j7 V7 z9 v8 i3 S3 R8 u" ?" ^      X_KNEE = $Xknee
0 \" O0 S, q  d. @#ENDIF
7 {9 ~% }* D6 h#LOCAL REAL TURNS_V
#IF SAT==1
. t7 j2 P  O/ s" d, B& [      TURNS_V = $V1
#ELSE
      TURNS_V = $V2
#ENDIF
# l6 k% t; w  L2 b7 M' l/ N#IF SAT==TAP
     + * RVD1_1
6 D7 j& w" z3 v* ]3 k3 k) o8 e8 U#ENDIF
- K2 j; D9 X8 M8 ~#IF SAT==1
" T+ I0 ^- {1 V, } #IF CuL>0.0
* K+ E- O, L* l- g- ~      CALL TSAT21($X1,$B1,
#ELSE
      CALL TSAT21($A1,$B1,
$ t  w* E( }; `6 K8 u) P6 T #ENDIF
+ p1 t/ Z, |+ b. y. r- E# [4 Z/ Q#ELSE
" F: F1 P/ A# @/ u3 e( U #IF CuL>0.0
0 C, w6 B' a% A- c. C! o      CALL TSAT21($X2,$B2,
#ELSE
      CALL TSAT21($A2,$B2,
#ENDIF
2 `) h( \; d' Y& }  V$ Q' T- I#ENDIF
     + 0,0, 0,0, $SS,$Tmva,TURNS_V,$Xair,X_KNEE,$f,$Tdc,$I_m1, 0.0, $Enab)
; ~7 }, l; Q0 E! \, K7 @! `#OUTPUT REAL IMAG {-STOR(NEXC-21+17)}
#OUTPUT REAL FLUX {STOR(NEXC-21+11)}
#IF Tap
3 N# e( H% I) ^( x      CALL TRTAP5($#[1], $Tap, $SS, $TapI)
7 B$ |2 W" H+ s, P5 b% G. Q# x#ENDIF
!

风残冷月

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

drzyx

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

suner1986

什么意思啊，看不懂fortran

sgkevin

完全看不懂。。

小小鸟儿

这是什么程序里面的?

huafeinaojin

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

buddam

要学习学习了

gpbaby199

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