atmospheric pressure in POP2

[jjakacki@iopan_gda_pl]

 We found, that atmospheric pressure in the active ocean model (POP2) is ignored. We work on small shelf sea - Baltic the model provides insufficient dense water transport during Baltic Sea inflows. So we looked at the data transfer from atmosphere (datm). 
The model has in barotropic.F90 lines (because I modified the code, in original version line numbers could be different): 449 !-----------------------------------------------------------------------
450 !
451 !     calculate negative gradient of surface atmospheric pressure
452 !     and add it to r.h.s. forcing
453 !
454 !-----------------------------------------------------------------------
455
456      if (ap_data_type /= 'none') then
457
458          call grad(1, WORK1, WORK2, ATM_PRESS(:,:,iblock), this_block)
459
460          WORK3 = WORK3 - c2dtp*WORK1
461          WORK4 = WORK4 - c2dtp*WORK2
462      endif I commented out lines 456 and 462 and the differences between integration with and without commented lines is attached. The image presents sea level differences between integration with atmospheric pressure and without it. I also added image with atmospheric pressure for the same day and contours on sea level differences are plotted. In my point of view it looks like a bug, not like a conscious. Regards,Jaromir