--- trunk/libf/dyn3d/leapfrog.f90 2012/04/20 14:58:43 61 +++ trunk/libf/dyn3d/leapfrog.f90 2012/07/26 14:37:37 62 @@ -106,7 +106,7 @@ call assert(shape(ucov) == (/iim + 1, jjm + 1, llm/), "leapfrog") itaufin = nday * day_step - ! "day_step" is a multiple of "iperiod", therefore "itaufin" is one too + ! "day_step" is a multiple of "iperiod", therefore so is "itaufin". dq = 0.