[MITgcm-devel] missing barrier causing problems for packing

Matthew Mazloff mmazloff at ucsd.edu
Mon Jan 30 17:13:19 EST 2012


Hello

We have come across the problem where the packing initiates before all  
processors finish writing adxx_theta, and thus sometimes ecco_cost is  
incorrect (zeros were read/written instead of proper gradients.)

We need to have a barrier put in either

around line 261 at the beginning of the_main_loop.F before       CALL  
INITIALISE_VARIA( mythid )

or in the_model_main.F  around line 600 between
       CALL THE_MAIN_LOOP( myCurrentTime, myCurrentIter, myThid )
and
          CALL CTRL_PACK( .FALSE. , mythid )


Let me know what is preferable -- and if someone wants to go ahead and  
check this in that would be great!

Thanks
Matt





More information about the MITgcm-devel mailing list