anchor.diff:        spongeHTimeScale, spongeUTimeScale, spongeVTimeScale, &
         spongeH, spongeU, spongeV, &
         nx, ny, layers, RedGrav, DumpWind)
  -  print *, 'Execution ended normally'
     stop 0
   end program aronnax
   
  
positive.code:   allocate(wind_x(0:nx+1, 0:ny+1))
    allocate(wind_y(0:nx+1, 0:ny+1))
  
    nwrite = int(dumpFreq/dt)
    avwrite = int(avFreq/dt)
  
  
negative.code:   !   two older time steps.
    ! - The model then solves for the tendencies at the current step
    !   before solving for the fields at the next time step.
  
    print "(A)", "Initialized"
  
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    !!! MAIN LOOP OF THE MODEL STARTS HERE                                  !!!
  