I managed to trace the error to ray_effect.f
! INITIALISATION
WRITE(*,*) "Test2b"
CALL OV( 'X=C ' ,KE,KE,KE,0.D0,NPOIN )
WRITE(*,*) "TEst2c"
CALL OV( 'X=C ' ,IH,IH,IH,0.D0,NPOIN )
WRITE(*,*) "Test2d"
CALL OS( 'X=0 ' ,X=EFF)
WRITE(*,*) "Test2"
!
The error happens after "Test2d". I tried writing EFF to screen to see the output, but the script doesn't compile, i.e.,
WRITE(*,*) "Test2d", EFF
EFF is listed as "SUNSHINE EFFECT ON ALGAE GROWTH". It feels like I'm missing a keyword in the cas file for this, but I'm not clear on which one needs to be activated. I have attached my draft cas files in case.
Any help appreciated.