loadModel( 'LTD.eml' )

PKC = createEntityStub( 'Variable:/CELL/CYTOPLASM:PKC' )

aDict = {}
aList = ( "AMPAR", "AMPAR_P", "PKC", "PKCactive" )

for aVariableName in aList:
    aDict[ aVariableName ]  = createLoggerStub( "Variable:/CELL/CYTOPLASM:" + aVariableName + ":Value" )
    aDict[ aVariableName ].create()

message( "start simulation......" )
duration = 1000
run( duration )

PKC[ 'Value' ] = PKC[ 'Value' ] * 100
message( "multiply PKC value by 100" )

run ( 500 )

from ecell.ECDDataFile import *

for aVariableName in aList:
    aDataFile = ECDDataFile( aDict[ aVariableName ].getData(500,1900,0.1) )
    aDataFile.save( aVariableName + ".ecd" )
    del aDataFile