Hello,
I was looking at the code to understand what is the difference between climatology and nudging cppkeys.
In step2d.F, there is something strange in the treatment of SSH vs velocities.
More precisely, line 419 and following
#ifdef ZNUDGING
...
zeta_new(i,j)=zeta_new(i,j) + dtfast*Znudgcof(i,j)
& *(ssh(i,j)-zeta_new(i,j))
#endif /* ZNUDGING */
While, line 1441, there is a second condition in the if statement
#if defined M2NUDGING && defined M2CLIMATOLOGY
...
DUnew = dtfast*M2nudgcof(i,j)*(ubclm(i,j)-ubar(i,j,knew))
...
ubar(i,j,knew)=ubar(i,j,knew) + DUnew
#endif
In step3d_t, temperature is treated as the velocities (with both TCLIMATOLOGY and TNUDGING, while SSH has only ZNUDGING).
- So why SSH is treated differently?