Hi Kelly!

I had totally forgot that I had put FLX 1,2 and 3 in my "myoutpufields.txt". It was before I put the variables in the Registry and recompiled.

Apparently the file said: +:h:10:FLX1,FLX2,FLX3, so that is where output stream number 10 came from. I don't know why I chose outputstream 10...

Anyway, I made the file empty and ran the model again. The variables are still 0 in the historyfile, but not when I look at the ouput from the print statements. So, the problem still exists....

What I am trying to figure out in a larger perspective is my energy balance at the surface. According to this I use following lines in a Python script to calculate it myself(TIDSLANGD means number of timesteps):

[GLW = Filobjekt.variables['GLW'][0:TIDSLANGD, 1, 1]

SWDOWN = Filobjekt.variables['SWDOWN'][0:TIDSLANGD, 1, 1]

GRDFLX = Filobjekt.variables['GRDFLX'][0:TIDSLANGD, 1, 1]

HFX = Filobjekt.variables['HFX'][0:TIDSLANGD, 1, 1]

LH = Filobjekt.variables['LH'][0:TIDSLANGD, 1, 1]

EMISS = Filobjekt.variables['EMISS'][0:TIDSLANGD, 1, 1]

ALBEDO = Filobjekt.variables['ALBEDO'][0:TIDSLANGD, 1, 1]

FLX1 = Filobjekt.variables['FLX1'][0:TIDSLANGD, 1, 1]

FLX2 = Filobjekt.variables['FLX2'][0:TIDSLANGD, 1, 1]

FLX3 = Filobjekt.variables['FLX3'][0:TIDSLANGD, 1, 1]

SWUP = ALBEDO*SWDOWN

GUPLW = EMISS*5.67e-8*(TSK+273.15)**4

RESIDUAL = GLW + SWDOWN - SWUP - GUPLW - HFX - LH + GRDFLX - FLX1 - FLX2 - FLX3

I want to compare my RESIDUAL to the NOAHRES variable from the Noah Land surface scheme:

NOAHRES = Filobjekt.variables['NOAHRES'][0:TIDSLANGD, 1, 1]

NOAHRES is defined like this in Noah lsm driver(with my write statements below):

! Residual of surface energy balance equation terms

!

noahres(i,j) = ( solnet + lwdn ) - sheat + ssoil - eta &

- ( emissi * STBOLT * (t1**4) ) - flx1 - flx2 - flx3

WRITE(*,*)"MH FLX1=",FLX1

WRITE(*,*)"MH FLX2=",FLX2

WRITE(*,*)"MH FLX3=",FLX3

RESIDUAL and NOAHRES are not the same, which is strange. It is exactly the same signs in my equation as in the model, where I have seen that LH is the same as ETA, SSOIL is the same as GRDFLX and HFX is the same as SHEAT. I suppose that is right, or?

So the first thing I want to see is if the FLX-variables can be responsible for this... I don't think so, because they seem too small. Anyway, at the moment they are ZERO at every timestep in my equation, so it is difficult to really verify this hypotheses.

Another reflection :

In the description of ETA and SHEAT in module_sf_noahlsm.F it says that SHEAT is positive UP FROM SURFACE and that ETA is NEGATIVE up from surface, but the variables has THE SAME sign in the equation in the model:

noahres(i,j) = ( solnet + lwdn ) - sheat + ssoil - eta &

- ( emissi * STBOLT * (t1**4) ) - flx1 - flx2 - flx3

SHEAT SENSIBLE HEAT FLUX (W M-2: POSITIVE, IF UPWARD FROM

! SURFACE)'

ETA ACTUAL LATENT HEAT FLUX (W m-2: NEGATIVE, IF UP FROM

! SURFACE)

Thanks for your time Kelly!

//Martin