13/09/14 21:43:33.37 4CwOwE5K
>>11
>そういう解釈でいいのかな?
それでいいよ。ところで実際のプログラムをちょっと見てみた。
FUNCTION STKSIM (RK, R, POP, CATCH)
C Calculate the stock trajectory with parameter RK and R: return
C the current stock size
C RK = notional carrying capacity; R = productivity parameter * 1.4184
COMMON /MANDAT/ ISTART, IYEAR, NS, NZ, RKLO, RKHI
INTEGER ISTART, IYEAR, NS, NZ
REAL RKLO, RKHI
REAL CATCH(0:*), POP(0:*)
POP(ISTART) = RK
DO 10 I = ISTART + 1, IYEAR
D = POP(I - 1) / RK
POP(I) = POP (I - 1) * (1. + R * (1. - D * D)) - CATCH(I-1)
IF (POP(I) .LE. 0) GOTO 20
10 CONTINUE
STKSIM = POP(IYEAR)
RETURN
20 STKSIM = 0.
END
これはRPMプログラムのRK(環境収容力)とR(μ * 1.4184)を元に、最初の年からIYEAR年後の資源量を計算するサブルーチン。
配列POPに各年の資源量を格納する。CATCHは毎年の捕獲量。