Attempting parallelization in ocean, need help in science

hi I was attempting to optimize the ocean code and ended doing the following for a particular reason do k=1,km //main k loop ....   if(k==1)then   do kk=1,km           call hdifft(kk, WORKN_TEMP(:,:,:,kk), TMIX, UMIX, VMIX, this_block)   enddo   endif WORKN = WORKN_TEMP(:,:,:,k)...enddo That is call all of the hdifft routine over all the verticals at once(Assuming mix_coeffs_kpp is completed for the first vertical and vmix_type_kpp is true).

 Can the final answer/NC file give wrong results, if else how?
 I have found that my NC file do differ and not able to find out how or why? there are no dependencies. Can anyone help me on why I see these errors?
 I am running a pure MPI based run and OpenMP is off.
 
Back
Top