[MITgcm-devel] ALLOW_SEAICE_COST_AREASST

[Matthew Mazloff]

Hi Patrick,

You may want to remove the
   IF ( useSeaice) THEN
in ecco/the_main_loop.F
just so it isn't forgotten about when one defines  
ALLOW_SEAICE_COST_AREASST

I still turn set useSeaice = FALSE for the adjoint run but I do want  
the adjoint to see the cost term (correct?).  so I commented this out:

#ifdef ALLOW_SEAICE
#ifndef DISABLE_DEBUGMODE
       IF ( debugLevel .GE. debLevB )
      &    CALL DEBUG_CALL('seaice_cost_driver',myThid)
#endif
CMM      IF ( useSeaice) THEN
       call timer_start('seaice_cost_driver [ECCO SPIN-DOWN]', mythid)
       call seaice_cost_driver( myiter, mytime, mythid )
       call timer_stop ('seaice_cost_driver [ECCO SPIN-DOWN]', mythid)
CMM      ENDIF
#endif


On another note, could we set  SEAICEmultidimadvection in data.seaice  
instead of in seaice_advdiff.F as it is currently done::
       SEAICEmultidimadvection = .TRUE.
       IF ( SEAICEadvScheme.EQ.ENUM_CENTERED_2ND
      & .OR.SEAICEadvScheme.EQ.ENUM_UPWIND_3RD
      & .OR.SEAICEadvScheme.EQ.ENUM_CENTERED_4TH ) THEN
        SEAICEmultiDimAdvection = .FALSE.
       ENDIF
If so, this would allow us to use various advection schemes with the  
adjoint (e.g. 33), instead of just the ones above (especially since 3  
and 4 aren't even implemented yet...meaning we are really stuck with 2)

-Matt

ps> SOSE_2005-2007 is now running!  just in time for 2009 :)