[MITgcm-devel] [MITgcm-cvs] MITgcm/verification/seaice_itd/input.melt CVS Commit

Wed Apr 3 20:10:07 EDT 2013

Torge, not sure I am answering correct question as I have not followed discussion,
but if you are asking how to do the exchange between tiles,
look at following lines at bottom of seaice_read_pickup.F for an example:

C--    Update overlap regions
       CALL EXCH_UV_XY_RL( uIce, vIce,.TRUE.,myThid)
       _EXCH_XY_RL( HEFF, myThid )
       _EXCH_XY_RL( AREA, myThid )
       CALL EXCH_3D_RL   ( TICES, MULTDIM, myThid )
       _EXCH_XY_RL(TICE , myThid )
c      _EXCH_XY_RL(YNEG , myThid )
       _EXCH_XY_RL(HSNOW, myThid )
#if (defined(SEAICE_CGRID) && defined(SEAICE_ALLOW_EVP))
       IF ( SEAICEuseEVP ) THEN
          _EXCH_XY_RL(seaice_sigma1 , myThid )
          _EXCH_XY_RL(seaice_sigma2 , myThid )
          _EXCH_XY_RL(seaice_sigma12, myThid )
       ENDIF
#endif /* SEAICE_CGRID SEAICE_ALLOW_EVP */

Dimitris Menemenlis

On Apr 2, 2013, at 11:40 PM, Torge Martin wrote:

I found that the overlap extent in AREAITD and HEFFITD is filled with 0 but is filled with the same value of the one center coordinate in AREA and HEFF when running without SEAICE_ITD (using "print*,AREA" etc. during runtime). Would that make a difference and how do I fill the overlap extent with the correct numbers?

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-devel/attachments/20130404/4ae49f59/attachment.htm>