Open gelianxing opened 4 months ago
This is in CER_Growth.for, line 737. TVR1 uses values SENLFG and SENLFGRS which are calculated later in the routine. Maybe this TVR1 calculation could be moved later in the routine to use today's senescence values. I don't think it should go to integrate because it is still a rate calculation. @palderman @GerritHoogenboom @fabiooliveira72 - if you have time, please have a look.
Linking directly for convenience:
Just did a quick search and we'd need to be careful about where we move it, because TVR1
is used to calculate RTWTGRS
just below:
and RTWTGRS
is then used subsequently:
The dependency graph seems to balloon rather quickly after that. Would it be easier to move senescence calculations earlier?
@palderman Thanks for looking. Yes, maybe. But it might be circular logic and we just need to live with the fact that this is a daily model.
I agree move senescence calculations earlier
IF (((LFWT+GROLF-SENLFG-SENLFGRS) ! Prevent divide by zero & +(STWT+GROST)+(RSWT+GRORS)) .GT. 0.0)THEN TVR1 = ! Conc & (RSWT+GRORS-SENRS)/ & ((LFWT+GROLF-SENLFG-SENLFGRS) & +(STWT+GROST)+(RSWT+GRORS)) ELSE TVR1 = 0.0 END IF
SENLFG、SENLFGRS why these variables use yesterday’s variables,maybe can move to integr