[solved!] Strange behaviour with GOTM-FABM-NPZD-OWS_PAPA

1) GOTM problem in temperature output

Different results with the same settings between GOTM 4.0.0 and 4.1.0.

Both models do NOT use T restoring (i.e. restoring time is set to some large value, e.g. 1e15), but the results are significantly different (Fig.1 4.0.0 on right 4.1.0 on left).

Figure 1: SST (bottom) and Time series of temperature in top 250 m

Nevertheless, I decided to restore T for 4.1.0 to fix physics. Now, T is restored alright (Fig.2 left), but the corresponding biology looks so strange (Fig.2 right).

Figure 2: Temperature (left) and Phytoplankton in top 250 m in 4.1.0.

• no swr.dat is present in 4.0.0.
• This is because heatflux.dat contains both swr (short wave radiation) and qh (surface heat flux) data.
• In 4.1.0, two variables are stored in different files swr.data and heatflux.dat.

I FOUND it... 'swr_method=0' was the problem. It should be 'swr_method=2'. Otherwise, no insolation, and no biological response.

The newly obtained results look promising! Temperature looks like the one in the Figure 2 without restoring, and the phytoplankton exists in the new run :D