Ocean Modeling Discussion

The idea of the new layer thickness is that is a way to preserve the stratigraphy of the bed after deposition has occurred and new layers need to be created. You might need to fine tune the choice of new_layerthick so that it makes sense to your application.

Well, I think it might be because you have newlayer_thick set to zero, so it will create new layers of zero thickness. Could you change that?

Could you please modify a line in the sed_bed_cohesive.F code? Could you change line 705 from: DO k=2,Nbed to DO k=2,Nbed-1 That specific loop should now look like: DO k=2,Nbed IF (thck_avail.lt.thck_to_add) THEN thck_avail=thck_avail+bed(i,j,k,ithck) Ksed=k END IF END DO That should deal with the ...

Looks like an interesting problem. In the cohesive bed routine, computations happen in both the Ksed layer and in the Ksed+1, but that should not exceed the number of bed layers. “! Ksed is used for grain props, and Ksed+1 for bed_tau_crit at bottom of layer Ksed” The error message seems to be hap ...

MUD_TAUCR_TIME == 28800.0d0 is in seconds (as the rest of the model). That value implies a consolidation nudging time scale (tendency toward equilibrium critical stress profile) of 8 hours. As the Swelling time scale is 100 times slower, then swelling is order 800 hours in your simulation. bed_tau ...