PERFECT_RESTART oddity

Message

kate · #1 Unread post by **kate** » Fri Nov 02, 2012 4:39 pm

In get_state.F, there's code to read in nstp and friends on a PERFECT_RESTART. I watched in the debugger and this code was never executed. It doesn't matter because nstp is being set in main3d based on iic.

The reason the code isn't executed is here:

Code: Select all

IF ((model.eq.0).and.(nrrec(ng).ne.0)) THEN

For my domain, model was set to one. To find out why, go to initial.F and see:

Code: Select all

#ifdef INI_FILE
      CALL get_state (ng, iNLM, 1, INI(ng)%name, IniRec, Tindex)
#else
      IF (nrrec(ng).ne.0) THEN
        CALL get_state (ng, 0, 1, INI(ng)%name, IniRec, Tindex)
      END IF
#endif

In other words, model will only get set to zero on restart in the case of ANA_INITIAL.