integer levoff
levoff = mod(modelstartdate(1)/100,100)
levmon = (irectmp-1) + levoff
levmon = (irec-1) + levoff