[MITgcm-devel] ocean_emissivity consistency

[Jean-Michel Campin]

Hi Martin,

I don't much how thing goes regading long-wave in pkg/seaice,
but in pkg/exf, when ocean_emissivity < 1, we should also 
apply the ocean_emissivity factor to the downward long-wave
flux when computing the absorbed lw flux:
LW_net(+=up) = ocean_emissivity*( sigma*SST^4 - lw_down )

In Large & Yeager 2004, eq.12, they recommand to use 
ocean_emissivity=1 for this reason. And over sea-ice,
their equation 23 is clear, with ice_emissivity on both
side. Do we have the same issue in pkg/seaice ?

Cheers,
Jean-Michel

On Mon, Dec 19, 2011 at 06:49:01AM -0500, Martin Losch wrote:
> Update of /u/gcmpack/MITgcm/pkg/seaice
> In directory forge:/tmp/cvs-serv17544/pkg/seaice
> 
> Modified Files:
> 	seaice_readparms.F 
> Log Message:
> make ocean_emissivity*stefanBoltzmann the default for
> SEAICE_emissivity, if useEXF; requires setting a wrong parameter value
> in global_ocean.cs32x15/input.seaice/data.seaice for backward
> compatiblity
> 
> 
> _______________________________________________
> MITgcm-cvs mailing list
> MITgcm-cvs at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-cvs