Branch to provide coupling with WaveWatchIII

[jcastill]

It will allow sending the total ice fraction and the surface current to the wave program, and receiving the significant wave height and the wave to ocean energy flux. The key key_ww3 will be needed to activate this coupling.

[jcastill]

The name of the branch will be dev_r5518_ww3_coupling

[jcastill]

The revision used should be the @5847. Revision @5851 merges the mslp branch dev_r5518_ukv_mslp

[jcastill]

There was a problem when sending the sea surface height to the wave model when ln_apr_dyn is true (pressure gradient added to the equations): the program crashes in the routine sbc_cpl_snd, in the lines for ssnd(jps_ssh)%laction. When ln_apr_dyn is true, the program updates ztmp1 with the values of ssh_ib and ssh_ibb, but at the first time step these variables are not allocated. The previous variables are allocated in the routine sbc_apr, which is called only after sbc_cpl_snd is called. I think this happens because at some point the call to sbc_cpl_snd was moved from the end of routine stp to the beginning.

The problem was solved by modifying a bit the problematic lines in the sbc_cpl_snd routine:

IF( ssnd(jps_ssh )%laction ) THEN

! ! removed inverse barometer ssh when Patm
! forcing is used (for sea-ice dynamics)
IF( ln_apr_dyn ) THEN

IF( kt /= nit000 ) THEN

ztmp1(:,:) = sshb(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) )

ELSE

ztmp1(:,:) = sshb(:,:)

ENDIF

ELSE

ztmp1(:,:) = sshn(:,:)

ENDIF

[jcastill]

Ticket reopened because further science changes are needed for ocean/wave coupling

[jcastill]

The latest developments related to this ticket have already been merged into the trunk. Further modifications will be taken into account on new tickets.

[nemo]

Remove 'Development' type