Maybe one can give me a hint to find the right place. I used the module_sf_noahmpdrv.F to modify SWE.
JLOOP : DO J=jts,jte
ILOOP : DO I = its, ite
......
SWE = MAX(0.0, SWE - 0.01) !FERO(I,J) * DT ) !+ grid%fdep(I,J) * DT)
SNOWERACC(I,J) = SNOWERACC(I,J) - 0.01 !FERO(I,J) *DT* (1/rhosn) *1000 ![m]
.....
ENDDO ILOOP ! of I loop
ENDDO JLOOP ! of J loop
!------------------------------------------------------
END SUBROUTINE noahmplsm
That was a test: to reduce SWE about 0.01 every time step.. but after one hour there was just a difference of 0.01 between the normal and the test run?? But SNOWERACC was much higher.
Is there somewhere in the code a line that overwrites SWE ?
I really tried my best to find it out by myself, but now I am running out of time for my PhD..

If anybody is familiar with the code.. please help me

Christina