39 |
|
|
40 |
c == end of interface == |
c == end of interface == |
41 |
|
|
42 |
|
jtlo = mybylo(mythid) |
43 |
|
jthi = mybyhi(mythid) |
44 |
|
itlo = mybxlo(mythid) |
45 |
|
ithi = mybxhi(mythid) |
46 |
|
jmin = 1 |
47 |
|
jmax = sny |
48 |
|
imin = 1 |
49 |
|
imax = snx |
50 |
|
|
51 |
do bj = jtlo,jthi |
do bj = jtlo,jthi |
52 |
do bi = itlo,ithi |
do bi = itlo,ithi |
53 |
do j = jmin,jmax |
do j = jmin,jmax |
54 |
do i = imin,imax |
do i = imin,imax |
55 |
c #1: |
c #1: |
56 |
gencost_modfld(i,j,bi,bi,1) = qnet(i,j,bi,bi) |
c gencost_modfld(i,j,bi,bi,1) = qnet(i,j,bi,bi) |
57 |
c #2: |
c #2: |
58 |
gencost_modfld(i,j,bi,bi,2) = empmr(i,j,bi,bi) |
gencost_modfld(i,j,bi,bi,1) = empmr(i,j,bi,bi) |
59 |
c #3: |
c #3: |
60 |
gencost_modfld(i,j,bi,bi,3) = fu(i,j,bi,bi) |
c gencost_modfld(i,j,bi,bi,3) = fu(i,j,bi,bi) |
61 |
c #4: |
c #4: |
62 |
gencost_modfld(i,j,bi,bi,4) = fv(i,j,bi,bi) |
c gencost_modfld(i,j,bi,bi,4) = fv(i,j,bi,bi) |
63 |
c |
c |
64 |
enddo |
enddo |
65 |
enddo |
enddo |