[MITgcm-support] Qnet and Qsw heat fluxes after ceckpointc54

[Martin Losch]

Hi Qentin,

the model variable Qnet is really the net heat flux. If you have a closer look at the code that you copied from external_forcing_surf.F, you'll realize that Qsw is actually subtracted (and not added) from Qnet (because Qsw can penetrate water, it will be added later and distributed over more than just the surface layer).

If you want shortwave heat to penetrate into the ocean, you need to specify a qswfile, independtly of using the surfQnetFile or the obsolete surfQfile. 

Martin

On Jun 18, 2013, at 12:29 PM, quentin jamet <quentinjamett at gmail.com> wrote:

> Hi all,
> 
> I am running a forced ocean model on the CS-grid shape at low resolution.
> I have questioning about heat fluxes.
> When looking at the definition of Qnet (surface net heat flux) in FFIELD.h, I can read :
> Qnet = latent + sensible + net longwave + net shortwave
> 
> So I concluded that shortwave heating has to be included into the QnetFile. 
> This is coherent with external_fields_load.F where Qsw data are added to Q data to obtain Qnet data.
> So if I directly prescribe a QnetFile, it is the same that prescribing a Qfile and a QswFile.
> But in external_forcing_surf.F, Qsw is also added to Qnet : 
> 
> C     Net heat flux Qnet:
>           surfaceForcingT(i,j,bi,bj) = surfaceForcingT(i,j,bi,bj)
>      &       - ( Qnet(i,j,bi,bj)
> #ifdef SHORTWAVE_HEATING
>      &          -Qsw(i,j,bi,bj)
> #endif
>      &         ) *recip_Cp*mass2rUnit
> 
> I so wonder which file do I have to prescibe to obtain a total heat flux containing shortwave heating?
> 
> Thanks in advance
> 
> ---------------------
> Quentin Jamet
> 'Le vent nous pousse toujours, il s'agit de savoir dans quel sens!'
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support