[Mitgcm-support] Re: advection code

[mitgcm-support at dev.mitgcm.org]

Alistair,

that's good news!
Could you please send the code as well? 

Arne

Alistair Adcroft wrote:
> 
> Arne,
> 
> Surprise! Here's the first attempt at third order upwind.
> 
> I haven't compiled the code (hence not debugged it)!
> 
> I changed the CPP macro  first_order  to  UPWIND  and
> introduced another called  THIRD_ORDER.
> 
> Truth table:
>                           UPWIND    THIRD_ORDER
>   Centered second order   #undef    #undef
>   First-order upwind      #define   #undef
>   Third-order upwind      #define   #define
> 
> Can you try the following code, calc_gt.F, and if it works
> do the same for calc_gs.F  . calc_mom_rhs.F is much more
> work so stay with the current version for now.
> 
> Let me know what happens,
> 
> Alistair.

-- 

Dr. Arne Biastoch
Scripps Institution of Oceanography        phone: +1-858-822-3787
Physical Oceanography Research Division    fax  : +1-858-534-9820
MS: 0230                                   email: abiastoch at ucsd.edu
8605 La Jolla Shores Dr.            
La Jolla, CA 92093-0230, U.S.A.   http://www.ecco.ucsd.edu/~biastoch